Commit cc860c04 authored by brunner's avatar brunner
Browse files

No commit message

No commit message
parent f50ec4c5
......@@ -148,8 +148,8 @@ class ObservationOperator(object):
repeat(self.forecast_nmembers))
]
with Pool(3) as pool:
pool.starmap(self.extract_model_data, args)
# with Pool(3) as pool:
# pool.starmap(self.extract_model_data, args)
for i in range(0,self.forecast_nmembers):
idx = str(i).zfill(3)
......@@ -176,7 +176,6 @@ class ObservationOperator(object):
time_stamp = str((dacycle['time.start']+timedelta(hours=hstart)).strftime('%Y%m%d%H'))+'_'+str((dacycle['time.start']+timedelta(hours=hstop)).strftime('%Y%m%d%H'))
sites = ("lhw","brm","jfj","ssl")
self.dacycle = dacycle
cosmo_start = dacycle['time.start'].strftime('%Y%m%d%H') #+timedelta(hours=168)
cosmo_out = "/scratch/snx3000/parsenov/ctdas/"+cosmo_start+"_"+str(hstart)+"_"+str(hstop+1)+"/cosmo/output/"
......@@ -209,6 +208,7 @@ class ObservationOperator(object):
cdo.remapnn("lon=7.92_lat=47.92,", input = cosmo_out+"CO2_"+ens+"_"+time_stamp+".nc", output = cosmo_out+"CO2_ssl_"+ens+"_"+time_stamp+".nc")
cdo.remapnn("lon=7.92_lat=47.92,", input = cosmo_out+"hhl.nc", output = cosmo_out+"hhl_ssl.nc")
sites = ("lhw","brm","jfj","ssl")
for s,ss in enumerate(sites):
site_height.main(cosmo_out, str(ens), ss, time_stamp)
......
......@@ -347,8 +347,8 @@ class ObsPackObservations(Observations):
for obs in self.datalist: # first loop over all available data points to set flags correctly
obs.mdm = 0.1 # default is very high model-data-mismatch, until explicitly set by script
obs.flag = 1
obs.mdm = 1 # default is very high model-data-mismatch, until explicitly set by script
obs.flag = 0
identifier = obs.code
# species, site, method, lab, datasetnr = identifier.split('_')
# nr_obs_per_day = len([c.code for c in self.datalist if c.code == obs.code and c.xdate.day == obs.xdate.day and c.flag == 0])
......
......@@ -10,7 +10,7 @@ def main(dir, ens, station, times):
rlat = ifile_co2.variables['lat'][:]
rlon = ifile_co2.variables['lon'][:]
time = ifile_co2.variables['time'][:]
co2 = ifile_co2.variables['CO2'][:]*29./44.
co2 = ifile_co2.variables['CO2'][:] #*29./44. # not needed, it's in the observationoperator class
hhl = ifile_hhl.variables['HHL'][:]
ifile_co2.close()
ifile_hhl.close()
......
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