Commit e2099517 authored by brunner's avatar brunner
Browse files

No commit message

No commit message
parent ac0f0d98
......@@ -136,7 +136,7 @@ class ObservationOperator(object):
idx = str(i+1).zfill(3)
cosmo_file = os.path.join('/store/empa/em05/parsenov/cosmo_data/model_'+idx+'_%s.nc' % dacycle['time.sample.stamp'])
ifile = Dataset(cosmo_file, mode='r')
model_data[i,:] = np.squeeze(ifile.variables['CO2'][:])*29./44.)*1E6 # in ppm
model_data[i,:] = (np.squeeze(ifile.variables['CO2'][:])*29./44.)*1E6 # in ppm
ifile.close()
for j,data in enumerate(zip(ids,obs,mdm)):
......@@ -149,14 +149,14 @@ class ObservationOperator(object):
logging.info('ObservationOperator finished successfully, output file written (%s)' % self.simulated_file)
def run_forecast_model(self,dacycle)
def run_forecast_model(self,dacycle):
self.prepare_run()
self.run(dacycle)
def extract_model_data(self,dacycle,hts,ensnum):
time_stamp = dacycle['time.sample.stamp']
sites = ("lhw","brm","jfj","ssl")
self.dacycle = dacycle
cosmo_time_stamp = dacycle['time.start'].strftime('%Y%m%d%H') #+timedelta(hours=168)
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment