alles guet

da/cosmo/observationoperator_octe.py

@@ -108,16 +108,17 @@ class ObservationOperator(object):
 
         self.lambda_file = os.path.join(self.outputdir, 'lambda.%s.nc' % self.dacycle['time.sample.stamp'])
         ofile = Dataset(self.lambda_file, mode='w')
-        opar = ofile.createDimension('param', self.nparams)
-        omem = ofile.createDimension('member', self.forecast_nmembers)#len(members.nmembers))
+        opar = ofile.createDimension('nparam', self.nparams)
+        omem = ofile.createDimension('nensembles', self.forecast_nmembers)#len(members.nmembers))
 
-        l = ofile.createVariable('lambda', np.float32, ('member','param'),fill_value=-999.99)
+        l = ofile.createVariable('lambda', np.float32, ('nensembles','nparam'),fill_value=-999.99)
         co2 = np.empty(shape=(self.forecast_nmembers,self.nparams))
 
         for m in range(0,20):
             co2[m,:] = members[m].param_values
         l[:] = co2
         ofile.close()
+        os.system('cp '+self.lambda_file+' '+dacycle['da.vprm']+'/lambdas.nc')
 
         os.chdir(dacycle['da.obsoperator.home'])
 
@@ -127,7 +128,9 @@ class ObservationOperator(object):
             else:
                 os.rename(dacycle['dir.da_run']+"/"+absolute_start_time+"_"+str(starth+lag*168)+"_"+str(endh+lag*168), dacycle['dir.da_run']+"/non_opt_"+dacycle['time.start'].strftime('%Y%m%d%H')+"_"+str(starth+lag*168)+"_"+str(endh+lag*168))
 
-        os.system('python run_chain.py '+self.dacycle['run.name']+' '+absolute_start_time_ch+' '+str(starth+lag*168)+' '+str(endh+lag*168)+' -j meteo icbc emissions int2lm post_int2lm octe online_vprm cosmo')
+        os.system('python run_chain.py '+self.dacycle['run.name']+' '+absolute_start_time_ch+' '+str(starth+lag*168)+' '+str(endh+lag*168)+' -j octe online_vprm cosmo -f')
+        sys.exit()
+        #os.system('python run_chain.py '+self.dacycle['run.name']+' '+absolute_start_time_ch+' '+str(starth+lag*168)+' '+str(endh+lag*168)+' -j meteo icbc emissions int2lm post_int2lm octe online_vprm cosmo')
         logging.info('COSMO done!')
         os.chdir(dacycle['dir.da_run'])
 

da/cosmo/statevector.py

@@ -291,7 +291,7 @@ class StateVector(object):
         self.ensemble_members[lag].append(newmember)
 
         # Create members 1:nmembers and add to ensemble_members list
-        for member in range(1, self.nmembers):
+        for member in range(1, self.nmembers+1):
 #            rands = np.random.uniform(low=-1., high=1., size=self.nparams-1)
  #           rands_bg = np.random.uniform(low=-0.05, high=0.05, size=1)
 #            rands = np.random.randn(self.nparams-1)

octe.rc

@@ -99,10 +99,8 @@ da.obsoperator         : cosmo
 ! be used as observation operator in this experiment.
 !
 
-!da.obsoperator.home    : /store/empa/em05/parsenov/cosmo_my_prc_chain
 da.obsoperator.home    : /store/empa/em05/parsenov/cosmo_processing_chain
-da.bio.input           : /store/empa/em05/parsenov/cosmo_input/vprm/processed
-da.bg.input           : /store/empa/em05/parsenov/cosmo_input/icbc/processed
+da.vprm                : /store/empa/em05/parsenov/cosmo_input/online_vprm
 da.obsoperator.rc      : ${da.obsoperator.home}/tm5-ctdas-ei-zoom.rc
 
 !forward.savestate.exceptsam : TRUE