.

dartwrf/exp_config.py

@@ -4,13 +4,13 @@ exp = utils.Experiment()
 exp.expname = "template_experiment"
 exp.model_dx = 2000
 exp.n_ens = 40
-exp.superob_km = False  # False or int (spatial averaging of observations, unit: km)
+exp.superob_km = False  # False or int (spatial averaging of observations)
 
 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 = '/jetfs/home/lkugler/data/sim_archive//exp_v1.18_P1_nature+1//obs_seq_out/%Y-%m-%d_%H:%M_obs_seq.out'
 
 # path to the nature run, where we take observations from
-exp.nature_wrfout = '/mnt/jetfs/scratch/lkugler/data/sim_archive/exp_v1.18_P1_nature/2008-07-30_06:00/1/wrfout_d01_%Y-%m-%d_%H:%M:%S'
+exp.nature_expname = 'exp_v1.18_P1_nature+1'
 
 exp.input_profile = '/mnt/jetfs/home/lkugler/data/initial_profiles/wrf/ens/2022-03-31/raso.fc.<iens>.wrfprof'
 
@@ -30,6 +30,9 @@ exp.dart_nml = {'&assim_tools_nml':
                             output_sd='.true.',
                             stages_to_write='output',
                         ),
+                '&quality_control_nml':
+                    dict(outlier_threshold='-1',
+                        ),
                 '&location_nml':
                     dict(horiz_dist_only='.true.',
                         ),
@@ -70,9 +73,9 @@ wv62 = dict(var_name='Brightness temperature WV 6.2µm', unit='[K]',
 
 wv73 = dict(var_name='Brightness temperature WV 7.3µm', unit='[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., 
-            loc_horiz_km=20)
+            n_obs=256, obs_locations='square_array_evenly_on_grid',
+            error_generate=1., error_assimilate=1., 
+            loc_horiz_km=10)
 
 ir108 = dict(var_name='Brightness temperature IR 10.8µm', unit='[K]',
              kind='MSG_4_SEVIRI_TB', sat_channel=9, 
@@ -102,16 +105,16 @@ q = dict(var_name='Specific humidity', unit='[kg/kg]',
          loc_horiz_km=0.1, loc_vert_km=2.5)
 
 t2m = dict(var_name='SYNOP Temperature', unit='[K]',
-           kind='SYNOP_TEMPERATURE', n_obs=1, 
-           error_generate=0.1, error_assimilate=0.1,
-           loc_horiz_km=40, loc_vert_km=2.5)
+           kind='SYNOP_TEMPERATURE', 
+           n_obs=256, obs_locations='square_array_evenly_on_grid',
+           error_generate=0.3, error_assimilate=0.3,
+           loc_horiz_km=10, loc_vert_km=2)
 
 psfc = dict(var_name='SYNOP Pressure', unit='[Pa]',
             kind='SYNOP_SURFACE_PRESSURE', n_obs=1, 
             error_generate=50., error_assimilate=100.,
             loc_horiz_km=32, loc_vert_km=5)
 
-exp.observations = [vis]
+exp.observations = [wv73]
 exp.update_vars = ['U', 'V', 'W', 'THM', 'PH', 'MU', 'QVAPOR', 'QCLOUD', 'QICE', 'PSFC']
 
-