No commit message

No commit message

da/cosmo/observationoperator_octe.py

@@ -141,7 +141,7 @@ class ObservationOperator(object):
             pool.starmap(self.extract_model_data, args)
 
         for i in range(0,self.forecast_nmembers):
-            idx = str(i).zfill(3)
+            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.01)*1E6   # in ppm