Consolidate

Created by: lkugler

config/cfg.py

 from dartwrf import utils
 
 exp = utils.Experiment()
-exp.expname = "obs-cold_localization-narrow"
+exp.expname = "test_WV73_many_inf3"
 exp.model_dx = 2000
-exp.n_ens = 10
+exp.n_ens = 20
 
 exp.filter_kind = 1
-exp.prior_inflation = 0
-exp.post_inflation = 4
+exp.prior_inflation = 3
+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 = False #(2, 10)
 exp.superob_km = False  # False or int (spatial averaging of observations)
 exp.adjust_obs_impact = False
 
 exp.use_existing_obsseq = False  # False or pathname (use precomputed obs_seq.out files)
-exp.use_existing_obsseq = '/users/students/lehre/advDA_s2023/dartwrf_tutorial/very_cold_observation.out'
+#exp.use_existing_obsseq = '/users/students/lehre/advDA_s2023/dartwrf_tutorial/very_cold_observation.out'
 
 # path to the nature run, where we take observations from
 exp.nature = '/users/students/lehre/advDA_s2023/data/sample_nature/'
@@ -39,9 +40,9 @@ 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., 
-            cov_loc_radius_km=20)
+            n_obs=256, obs_locations='square_array_evenly_on_grid',
+            error_generate=2., error_assimilate=2., 
+            cov_loc_radius_km=10)
 
 ir108 = dict(plotname='Brightness temperature IR 10.8µm', plotunits='[K]',
              kind='MSG_4_SEVIRI_TB', sat_channel=9, 
@@ -58,29 +59,30 @@ radar = dict(plotname='Radar reflectivity', plotunits='[dBz]',
 
 t = dict(plotname='Temperature', plotunits='[K]',
          kind='RADIOSONDE_TEMPERATURE', 
-         #n_obs=22500, obs_locations='square_array_evenly_on_grid',
+         #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=50)
+         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]',
-           kind='SYNOP_TEMPERATURE', n_obs=1, 
-           error_generate=0.1, error_assimilate=0.1,
-           cov_loc_radius_km=40)
+           kind='SYNOP_TEMPERATURE', 
+           #n_obs=1, obs_locations=[(45., 0.)],
+           n_obs=256, obs_locations='square_array_evenly_on_grid',
+           error_generate=1, error_assimilate=1,
+           cov_loc_radius_km=10)
 
 psfc = dict(plotname='SYNOP Pressure', plotunits='[Pa]',
             kind='SYNOP_SURFACE_PRESSURE', n_obs=1, 
             error_generate=50., error_assimilate=100.,
             cov_loc_radius_km=32)
 
-exp.observations = [t2m]
+exp.observations = [wv73]
 exp.update_vars = ['U', 'V', 'W', 'THM', 'PH', 'MU', 'QVAPOR', 'QCLOUD', 'QICE', 'PSFC']
-#exp.update_vars = ['U', 'V', 'W', 'T', 'PH', 'MU', 'QVAPOR', 'PSFC']