From 9604f73ca98ad275f35207586464c0fc1b75ee2d Mon Sep 17 00:00:00 2001
From: lkugler <lukas.kugler@gmail.com>
Date: Wed, 7 Jun 2023 15:01:15 +0200
Subject: [PATCH] .
---
dartwrf/exp_config.py | 25 ++++++++++++++-----------
1 file changed, 14 insertions(+), 11 deletions(-)
diff --git a/dartwrf/exp_config.py b/dartwrf/exp_config.py
index 71650ab..6d90a43 100644
--- a/dartwrf/exp_config.py
+++ b/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']
-
--
GitLab