From 56d5232263c046a2ac4f973d190083e046f1eeb1 Mon Sep 17 00:00:00 2001
From: lkugler <lukas.kugler@gmail.com>
Date: Tue, 5 Jul 2022 10:44:15 +0200
Subject: [PATCH] .
---
config/cfg.py | 18 ++++++++++--------
templates/input.nml | 4 ++--
2 files changed, 12 insertions(+), 10 deletions(-)
diff --git a/config/cfg.py b/config/cfg.py
index f48bf9f..3046f68 100755
--- a/config/cfg.py
+++ b/config/cfg.py
@@ -7,15 +7,17 @@ class ExperimentConfiguration(object):
pass
exp = ExperimentConfiguration()
-exp.expname = "exp_v1.21_P3_wbub7_WV+VIS_obs10_loc20"
+exp.expname = "exp_v1.21_P2_rr_REFL_obs2-10_loc20_oe10"
exp.model_dx = 2000
exp.n_ens = 40
exp.n_nodes = 10
+exp.inflation = True
+exp.sec = True
#exp.nature_wrfout = '/home/fs71386/lkugler/data/sim_archive/exp_v1.19_P5+su_nat2/2008-07-30_07:00/1/wrfout_d01_%Y-%m-%d_%H:%M:%S'
-exp.nature_wrfout = '/home/fs71386/lkugler/data/sim_archive/exp_v1.19_P3_wbub7_nat/2008-07-30_12:00/1/wrfout_d01_%Y-%m-%d_%H:%M:%S'
+#exp.nature_wrfout = '/home/fs71386/lkugler/data/sim_archive/exp_v1.19_P3_wbub7_nat/2008-07-30_12:00/1/wrfout_d01_%Y-%m-%d_%H:%M:%S'
#exp.nature_wrfout = '/home/fs71386/lkugler/data/sim_archive/exp_v1.19_Pwbub5_nat/2008-07-30_12:00/1/wrfout_d01_%Y-%m-%d_%H:%M:%S'
-#exp.nature_wrfout = '/home/fs71386/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_wrfout = '/home/fs71386/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_wrfout = '/home/fs71386/lkugler/data/sim_archive/exp_v1.19_P4_nat/2008-07-30_07:00/1/wrfout_d01_%Y-%m-%d_%H:%M:%S'
#exp.input_profile = '/home/fs71386/lkugler/wrf_profiles/data/wrf/ens/2021-05-04/raso.nat.001.wrfprof'
@@ -32,9 +34,9 @@ exp.input_profile = '/gpfs/data/fs71386/lkugler/initial_profiles/wrf/ens/2022-03
# needs a horizontal scale too, to calculate the vertical normalization
# since you can not specify different vertical localizations for diff. variables
exp.cov_loc_vert_km_horiz_km = (2, 20)
-#exp.superob_km = 10
+exp.superob_km = 10
-n_obs = 961 # 22500: 2km, 5776: 4km, 121: 30km, 256:16x16 (20km); 961: 10km resoltn # radar: n_obs for each observation height level
+n_obs = 22500 # 22500: 2km, 5776: 4km, 121: 30km, 256:16x16 (20km); 961: 10km resoltn # radar: n_obs for each observation height level
vis = dict(plotname='VIS 0.6µm', plotunits='[1]',
kind='MSG_4_SEVIRI_BDRF', sat_channel=1, n_obs=n_obs,
@@ -53,7 +55,7 @@ ir108 = dict(plotname='Brightness temperature IR 10.8µm', plotunits='[K]',
radar = dict(plotname='Radar reflectivity', plotunits='[dBz]',
kind='RADAR_REFLECTIVITY', n_obs=n_obs,
- error_generate=2.5, error_assimilate=5,
+ error_generate=2.5, error_assimilate=10,
heights=np.arange(2000, 14001, 2000),
cov_loc_radius_km=20)
@@ -74,8 +76,8 @@ psfc = dict(plotname='SYNOP Pressure', plotunits='[dBz]',
cov_loc_radius_km=32)
-exp.observations = [vis, wv73]
-exp.update_vars = ['U', 'V', 'W', 'T', 'THM', 'PH', 'MU', 'QVAPOR', 'QCLOUD', 'QICE', 'PSFC']
+exp.observations = [radar]
+exp.update_vars = ['U', 'V', 'W', 'THM', 'PH', 'MU', 'QVAPOR', 'QCLOUD', 'QICE', 'PSFC']
#exp.update_vars = ['U', 'V', 'W', 'T', 'PH', 'MU', 'QVAPOR', 'PSFC']
# directory paths depend on the name of the experiment
diff --git a/templates/input.nml b/templates/input.nml
index 9f7de60..97968eb 100644
--- a/templates/input.nml
+++ b/templates/input.nml
@@ -57,7 +57,7 @@
output_sd = .true.
write_all_stages_at_end = .false.
- inf_flavor = 0, 4,
+ inf_flavor = 0, <post_inflation>,
inf_initial_from_restart = .true., .false.,
inf_sd_initial_from_restart = .true., .false.,
inf_initial = 1.00, 0.90,
@@ -107,7 +107,7 @@
cutoff = <cov_loc_radian>,
sort_obs_inc = .false.,
spread_restoration = .false.,
- sampling_error_correction = .true.,
+ sampling_error_correction = <sampling_error_correction>,
adaptive_localization_threshold = -1,
output_localization_diagnostics = .false.,
localization_diagnostics_file = 'localization_diagnostics',
--
GitLab