Skip to content
Snippets Groups Projects
Commit 5a611a1e authored by brunner's avatar brunner
Browse files

No commit message 

No commit message
parent 6da331f4
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
da/cosmo/observationoperator.py
+ 1
− 1
View file @ 5a611a1e
...@@ -121,7 +121,7 @@ class ObservationOperator(object): ...@@ -121,7 +121,7 @@ class ObservationOperator(object):
# UNCOMMENT FROM HERE # UNCOMMENT FROM HERE
for dt in rrule.rrule(rrule.HOURLY, dtstart=dacycle['time.start']+timedelta(hours=24*lag*int(dacycle['time.cycle'])), until=dacycle['time.start']+timedelta(hours=(lag+1)*24*int(dacycle['time.cycle']))): for dt in rrule.rrule(rrule.HOURLY, dtstart=dacycle['time.start']+timedelta(hours=24*lag*int(dacycle['time.cycle'])), until=dacycle['time.start']+timedelta(hours=(lag+1)*24*int(dacycle['time.cycle']))):
print('time time',lag,dt.strftime('%Y%m%d%H')) logging.info('Multiplying emissions with parameters for lag %d, date %s' % (lag, dt.strftime('%Y%m%d%H')))
for ens in range(0,self.forecast_nmembers): for ens in range(0,self.forecast_nmembers):
ens = str(ens).zfill(3) ens = str(ens).zfill(3)
cdo.setunit("'kg m-2 s-1' -expr,GPP_"+ens+"_F=CO2_GPP_F*parametermap -merge "+os.path.join(dacycle['da.bio.input'], 'gpp_%s.nc' % dt.strftime('%Y%m%d%H')), input = os.path.join(dacycle['restartmap.dir'],"parameters_lag"+str(lag)+"."+ens+".nc"), output = os.path.join(dacycle['da.bio.input'], 'ensemble', "gpp_"+ens+"_%s.nc" % dt.strftime('%Y%m%d%H'))) cdo.setunit("'kg m-2 s-1' -expr,GPP_"+ens+"_F=CO2_GPP_F*parametermap -merge "+os.path.join(dacycle['da.bio.input'], 'gpp_%s.nc' % dt.strftime('%Y%m%d%H')), input = os.path.join(dacycle['restartmap.dir'],"parameters_lag"+str(lag)+"."+ens+".nc"), output = os.path.join(dacycle['da.bio.input'], 'ensemble', "gpp_"+ens+"_%s.nc" % dt.strftime('%Y%m%d%H')))
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment