def mooring_1dvar(file_in): """Standard EcoFOCI Mooring .nc files with one dimensional parameters as a function of time such as temperature at a particular level Timestep of data is assumed to be in fractions of a day""" #dir_path = os.path.dirname(os.path.abspath(__file__)) ncfile = file_in ###nc readin/out nchandle = ncutil.ncopen(ncfile) params = constants.nc_vars_moor() ncdata = ncutil.ncreadfile(nchandle, params) ncutil.ncclose(nchandle) ###data massaging #time_all = ncdata[:,0] + ncdata[:,1] xx = ncutil.nc_missing(ncdata[:, 9], flag=1e35, setfill='Zero') pytime = util.EPICdate2udunits(ncdata[:, 0], ncdata[:, 1]) dt = 1. / pytime['interval_min'] #data is 4 times daily time_base = 'days' time = pytime['timeint'] #time = util.subsample(time, int(pytime.get('interval_min')) / 4) variance = np.var(xx) #normalize print 'Variance = %s ' % (variance) x = (xx - np.mean(xx)) / np.sqrt(variance) variance = np.var(x) return (xx, x, dt, np.array(time), variance, time_base)
def mooring_1dvar(file_in): """ TODO: make file choice and var choice more transparent""" #dir_path = os.path.dirname(os.path.abspath(__file__)) ncfile = file_in ###nc readin/out nchandle = ncutil.ncopen(ncfile) params = constants.nc_vars_moor() ncdata = ncutil.ncreadfile(nchandle, params) ncutil.ncclose(nchandle) ###data massaging #time_all = ncdata[:,0] + ncdata[:,1] xx = ncutil.nc_missing(ncdata[:,9], flag=1e35, setfill='Zero') pytime = util.EPICdate2udunits(ncdata[:,0], ncdata[:,1]) dt = 1. / pytime['interval_min'] #data is 4 times daily time_base = 'days' time = pytime['timeint'] #time = util.subsample(time, int(pytime.get('interval_min')) / 4) variance = np.var(xx) #normalize print 'Variance = %s ' % (variance) x = (xx - np.mean(xx)) / np.sqrt(variance) variance = np.var(x) return (xx, x,dt,np.array(time), variance, time_base)
def mooring_1dvar(file_in): """Standard EcoFOCI Mooring .nc files with one dimensional parameters as a function of time such as temperature at a particular level Timestep of data is assumed to be in fractions of a day""" #dir_path = os.path.dirname(os.path.abspath(__file__)) ncfile = file_in ###nc readin/out nchandle = ncutil.ncopen(ncfile) params = constants.nc_vars_moor() ncdata = ncutil.ncreadfile(nchandle, params) ncutil.ncclose(nchandle) ###data massaging #time_all = ncdata[:,0] + ncdata[:,1] xx = ncutil.nc_missing(ncdata[:,9], flag=1e35, setfill='Zero') pytime = util.EPICdate2udunits(ncdata[:,0], ncdata[:,1]) dt = 1. / pytime['interval_min'] #data is 4 times daily time_base = 'days' time = pytime['timeint'] #time = util.subsample(time, int(pytime.get('interval_min')) / 4) variance = np.var(xx) #normalize print 'Variance = %s ' % (variance) x = (xx - np.mean(xx)) / np.sqrt(variance) variance = np.var(x) return (xx, x,dt,np.array(time), variance, time_base)
def mooring_1dvar(file_in): """ TODO: make file choice and var choice more transparent""" #dir_path = os.path.dirname(os.path.abspath(__file__)) ncfile = file_in ###nc readin/out nchandle = ncutil.ncopen(ncfile) params = constants.nc_vars_moor() ncdata = ncutil.ncreadfile(nchandle, params) ncutil.ncclose(nchandle) ###data massaging #time_all = ncdata[:,0] + ncdata[:,1] xx = ncutil.nc_missing(ncdata[:, 9], flag=1e35, setfill='Zero') pytime = util.EPICdate2udunits(ncdata[:, 0], ncdata[:, 1]) dt = 1. / pytime['interval_min'] #data is 4 times daily time_base = 'days' time = pytime['timeint'] #time = util.subsample(time, int(pytime.get('interval_min')) / 4) variance = np.var(xx) #normalize print 'Variance = %s ' % (variance) x = (xx - np.mean(xx)) / np.sqrt(variance) variance = np.var(x) return (xx, x, dt, np.array(time), variance, time_base)