def convert_tz_xr(ds, input_tz='UTC', output_tz='US/Eastern', time_coord='time'):
    """
    Converts an xarray dataset from one timezone to another.

    Note that this will not work if you have problematic date features within your 
    simulation (e.g., time changes). This is why xarray hasn't integrated this feature yet. 

    Parameters
    ----------
    :param ds: `xarray.DataSet`
        DataSet that you want converted from one time zone to another.
    :param input_tx: string
        String identifying the current time zone of the DataSet.
    :param output_tx: string
        String identifying the desired time zone of the DataSet.
    :param time_coord: string
        Name of the coordinate holding date/time information within `xarray.DataSet`
        that you passed in with `ds`.
    :return: `xarray.DataSet`
        DataSet with updated time zone reflected in the time coordinate. 
    """
    tidx_in = ds.time.to_index().tz_localize(tz=input_tz)
    ds.coords[time_coord] = tidx_in.tz_convert(output_tz).tz_localize(None)
    return ds