.

config/cfg.py

 from dartwrf import utils
 
 exp = utils.Experiment()
-exp.expname = "obs-T_inflation-1.1"
+exp.expname = "advDA_Raso-T_inf2-4"
 exp.model_dx = 2000
-exp.n_ens = 20
+exp.n_ens = 40
 
 exp.filter_kind = 1
-exp.prior_inflation = 3
-exp.inf_initial = 1.1
+exp.prior_inflation = 2
+exp.inf_initial = 4
 exp.post_inflation = 0
 exp.sec = True
-exp.cov_loc_vert_km_horiz_km = (4, 40)
+exp.cov_loc_vert_km_horiz_km = (2, 100)
 exp.superob_km = False  # False or int (spatial averaging of observations)
 exp.adjust_obs_impact = False
 
@@ -41,7 +41,7 @@ wv62 = dict(plotname='Brightness temperature WV 6.2µm', plotunits='[K]',
 wv73 = dict(plotname='Brightness temperature WV 7.3µm', plotunits='[K]',
             kind='MSG_4_SEVIRI_TB', sat_channel=6, 
             n_obs=961, obs_locations='square_array_evenly_on_grid',
-            error_generate=1., error_assimilate=3., 
+            error_generate=1., error_assimilate=1., 
             cov_loc_radius_km=20)
 
 ir108 = dict(plotname='Brightness temperature IR 10.8µm', plotunits='[K]',
@@ -59,16 +59,16 @@ radar = dict(plotname='Radar reflectivity', plotunits='[dBz]',
 
 t = dict(plotname='Temperature', plotunits='[K]',
          kind='RADIOSONDE_TEMPERATURE', 
-         n_obs=961, obs_locations='square_array_evenly_on_grid',
-         #n_obs=1, obs_locations=[(45., 0.)],
+         #n_obs=961, obs_locations='square_array_evenly_on_grid',
+         n_obs=1, obs_locations=[(45., 0.)],
          error_generate=0.2, error_assimilate=0.2,
-         heights=[1000,], #range(1000, 17001, 2000),
-         cov_loc_radius_km=40)
+         heights=range(700, 17701, 500),
+         cov_loc_radius_km=100)
 
 q = dict(plotname='Specific humidity', plotunits='[kg/kg]',
          kind='RADIOSONDE_SPECIFIC_HUMIDITY', n_obs=1,
          error_generate=0., error_assimilate=5*1e-5,
-         heights=[1000], #range(1000, 17001, 2000),
+         heights=[1500], #range(1000, 17001, 2000),
          cov_loc_radius_km=0.1)
 
 t2m = dict(plotname='SYNOP Temperature', plotunits='[K]',

dartwrf/assim_synth_obs.py

@@ -476,7 +476,7 @@ def prepare_inflation_2(time, prior_init_time):
     if os.path.isfile(f_prior):
         copy(f_prior, f_new)
         print(f_prior, 'copied to', f_new)
-    else:
+    else:  # no prior inflation file at the first assimilation
         warnings.warn(f_prior + ' does not exist. Using default file instead.')
         copy(f_default, f_new)
 

templates/input.nml

@@ -63,7 +63,7 @@
    inf_initial                 = <inf_initial>,                0.90,
    inf_lower_bound             = 1.0,                     1.0,
    inf_upper_bound             = 1000000.0,               1000000.0,
-   inf_damping                 = 0.9,                     1.0,
+   inf_damping                 = 1.0,                     1.0,
    inf_sd_initial              = 0.6,                     0.0,
    inf_sd_lower_bound          = 0.6,                     0.0,
    inf_sd_max_change           = 1.05,                    1.05,