tutorial

analysis_only.py

@@ -6,9 +6,9 @@ import datetime as dt
 from dartwrf.workflows import WorkFlows
 
 
-prior_path_exp = '/mnt/jetfs/scratch/lkugler/data/sim_archive/exp_v1.19_P3_wbub7_noDA'
+prior_path_exp = '/users/students/lehre/advDA_s2023/data/sample_ensemble/'
 prior_init_time = dt.datetime(2008,7,30,12)
-prior_valid_time = dt.datetime(2008,7,30,12,30)
+prior_valid_time = dt.datetime(2008,7,30,13)
 assim_time = prior_valid_time
 
 w = WorkFlows(exp_config='cfg.py', server_config='srvx1.py')

config/cfg.py

 from dartwrf import utils
 
 exp = utils.Experiment()
-exp.expname = "test_newcode"
+exp.expname = "obs-cold_localization-narrow"
 exp.model_dx = 2000
 exp.n_ens = 10
 
@@ -9,30 +9,20 @@ exp.filter_kind = 1
 exp.prior_inflation = 0
 exp.post_inflation = 4
 exp.sec = True
-exp.reject_smallFGD = False
-exp.cov_loc_vert_km_horiz_km = (3, 20)
+exp.cov_loc_vert_km_horiz_km = (4, 40)
 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 = '/jetfs/home/lkugler/data/sim_archive/NoImpactFactors/obs_seq_out/2008-07-30_%H:%M_obs_seq.out'
-#exp.use_existing_obsseq = '/jetfs/home/lkugler/data/sim_archive/exp_v1.21_P3_wbub7_VIS_obs10_loc20/obs_seq_out/2008-07-30_%H:%M_obs_seq.out'  
-#exp.use_existing_obsseq = '/gpfs/data/fs71386/lkugler/sim_archive/exp_v1.21_P3_wbub7_REFL2D_obs10_loc20_oe5/obs_seq_out/2008-07-30_%H:%M_obs_seq.out'
-#exp.use_existing_obsseq = '/gpfs/data/fs71386/lkugler/sim_archive/exp_v1.21_P2_rr_VIS_obs20_loc4/obs_seq_out/2008-07-30_%H:%M_obs_seq.out'
+exp.use_existing_obsseq = '/users/students/lehre/advDA_s2023/dartwrf_tutorial/very_cold_observation.out'
 
-#exp.nature = '/mnt/jetfs/scratch/lkugler/data/sim_archive/exp_v1.19_P3_wbub7_nat/2008-07-30_12:00/1'
-exp.nature = '/mnt/jetfs/scratch/lkugler/data/sim_archive/exp_v1.18_P1_nature/2008-07-30_06:00/1'
+# path to the nature run, where we take observations from
+exp.nature = '/users/students/lehre/advDA_s2023/data/sample_nature/'
 
 exp.input_profile = '/mnt/jetfs/home/lkugler/data/initial_profiles/wrf/ens/2022-03-31/raso.fc.<iens>.wrfprof'
-#exp.input_profile = '/gpfs/data/fs71386/lkugler/initial_profiles/wrf/ens/2022-03-31/raso.nat.<iens>.wrfprof'
-#exp.input_profile = '/gpfs/data/fs71386/lkugler/initial_profiles/wrf/ens/2022-05-18/raso.fc.<iens>.wrfprof'
 
-
-# localize vertically, if it has a vertical position
-# needs a horizontal scale too, to calculate the vertical normalization
-# since you can not specify different vertical localizations for diff. variables
-
-# n_obs= 22500: 2km, 5776: 4km, 121: 30km, 256:16x16 (20km); 961: 10km resoltn # radar: n_obs for each observation height level
+# n_obs can be 22500: 2km, 5776: 4km, 121: 30km, 256:16x16 (20km); 961: 10km resoltn 
+# if radar: then n_obs is for each observation height level
 
 vis = dict(plotname='VIS 0.6µm', plotunits='[1]',
            kind='MSG_4_SEVIRI_BDRF', sat_channel=1, 
@@ -82,15 +72,15 @@ q = dict(plotname='Specific humidity', plotunits='[kg/kg]',
 
 t2m = dict(plotname='SYNOP Temperature', plotunits='[K]',
            kind='SYNOP_TEMPERATURE', n_obs=1, 
-           error_generate=0.1, error_assimilate=1.,
-           cov_loc_radius_km=20)
+           error_generate=0.1, error_assimilate=0.1,
+           cov_loc_radius_km=40)
 
 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 = [t]
+exp.observations = [t2m]
 exp.update_vars = ['U', 'V', 'W', 'THM', 'PH', 'MU', 'QVAPOR', 'QCLOUD', 'QICE', 'PSFC']
 #exp.update_vars = ['U', 'V', 'W', 'T', 'PH', 'MU', 'QVAPOR', 'PSFC']
 

config/srvx1.py

@@ -63,7 +63,7 @@ cluster.archive_base = utils.userhome+'/data/sim_archive/'
 cluster.srcdir = '/users/staff/lkugler/AdvDA23/DART/WRF-4.3/run'
 cluster.dart_srcdir = '/users/students/lehre/advDA_s2023/DART/models/wrf/work'
 cluster.rttov_srcdir = '/users/students/lehre/advDA_s2023/RTTOV13/rtcoef_rttov13/'
-cluster.scriptsdir = utils.userhome+'/DART-WRF/dartwrf/'
+cluster.scriptsdir = utils.userhome+'/AdvDA23/DART-WRF/dartwrf/'
 cluster.geo_em = '/users/students/lehre/advDA_s2023/data/geo_em.d01.nc'
 
 # templates/run scripts