def __init__(self, name=None, units='kg/m^3', temperature=None, salinity=None): if (temperature is None or salinity is None or not isinstance(temperature, TemperatureTS) or not isinstance(salinity, SalinityTS)): raise ValueError('Must provide temperature and salinity ' 'time series Environment objects') if len(temperature.time.time) > len(salinity.time.time): density_times = temperature.time else: density_times = salinity.time dummy_pt = np.array([ [0, 0], ]) import gsw from gnome import constants data = [ gsw.rho(salinity.at(dummy_pt, t), temperature.at(dummy_pt, t, units='C'), constants.atmos_pressure * 0.0001) for t in density_times.time ] TimeseriesData.__init__(self, name, units, time=density_times, data=data)
def __init__(self, name=None, units='kg/m^3', temperature=None, salinity=None): if (temperature is None or salinity is None or not isinstance(temperature, TemperatureTS) or not isinstance(salinity, SalinityTS)): raise ValueError('Must provide temperature and salinity ' 'time series Environment objects') if len(temperature.time.time) > len(salinity.time.time): density_times = temperature.time else: density_times = salinity.time dummy_pt = np.array([[0, 0], ]) import gsw from gnome import constants data = [gsw.rho(salinity.at(dummy_pt, t), temperature.at(dummy_pt, t, units='C'), constants.atmos_pressure * 0.0001) for t in density_times.time] TimeseriesData.__init__(self, name, units, time=density_times, data=data)
def __init__(self, name=None, units='K', time=None, data=None, **kwargs): if 'timeseries' in kwargs: ts = kwargs['timeseries'] time = map(lambda e: e[0], ts) data = np.array(map(lambda e: e[1], ts)) TimeseriesData.__init__(self, name, units, time, data=data)