From 823eaee9cc6a1bfa645d75cd9001b557846b9c81 Mon Sep 17 00:00:00 2001
From: lkugler <lukas.kugler@gmail.com>
Date: Tue, 31 Jan 2023 16:07:20 +0100
Subject: [PATCH] change in config
---
config/cfg.py | 63 +++++++++++++++++++++++++++------------------------
1 file changed, 34 insertions(+), 29 deletions(-)
diff --git a/config/cfg.py b/config/cfg.py
index d2c4c8e..c7272f4 100755
--- a/config/cfg.py
+++ b/config/cfg.py
@@ -1,4 +1,4 @@
-from config import clusters # from . = problem in archivedir
+from config import clusters
cluster = clusters.jet # change cluster configuration here
class ExperimentConfiguration(object):
@@ -6,34 +6,31 @@ class ExperimentConfiguration(object):
pass
exp = ExperimentConfiguration()
-exp.expname = "test_jet" #"exp_v1.22_P3_wbub7_WV62_obs10_loc20_oe1"
+exp.expname = "exp_v1.22_P2_rr_WV62_obs10_loc20_oe2"
exp.model_dx = 2000
exp.n_ens = 40
-exp.n_nodes = 40
+exp.size_jobarray = 40
exp.filter_kind = 1
exp.inflation = True
exp.sec = True
exp.reject_smallFGD = False
-exp.cov_loc_vert_km_horiz_km = False #(3, 20)
+exp.cov_loc_vert_km_horiz_km = (3, 20)
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.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 = '/jetfs/home/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 = '/jetfs/home/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 = '/jetfs/home/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'
-#exp.input_profile = '/home/fs71386/lkugler/wrf_profiles/data/wrf/ens/2021-05-04/raso.nat.<iens>.wrfprof'
-#exp.input_profile = '/home/fs71386/lkugler/wrf_profiles/data/wrf/ens/2021-05-04/raso.fc.<iens>.wrfprof'
-#exp.input_profile = '/home/fs71386/lkugler/data/initial_profiles/wrf/ens/large_mean_error/raso.nat.<iens>.wrfprof'
exp.input_profile = '/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'
@@ -43,53 +40,61 @@ exp.input_profile = '/jetfs/home/lkugler/data/initial_profiles/wrf/ens/2022-03-3
# needs a horizontal scale too, to calculate the vertical normalization
# since you can not specify different vertical localizations for diff. variables
-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: 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,
+ kind='MSG_4_SEVIRI_BDRF', sat_channel=1,
+ n_obs=961, obs_locations='square_array_evenly_on_grid',
+ # n_obs=1, obs_locations=[(44.141, -0.99)],
error_generate=0.03, error_assimilate=0.03,
cov_loc_radius_km=20)
wv62 = dict(plotname='Brightness temperature WV 6.2µm', plotunits='[K]',
- kind='MSG_4_SEVIRI_TB', sat_channel=5, n_obs=n_obs,
- error_generate=1., error_assimilate=1.,
+ kind='MSG_4_SEVIRI_TB', sat_channel=5,
+ n_obs=961, obs_locations='square_array_evenly_on_grid',
+ error_generate=1., error_assimilate=2.,
cov_loc_radius_km=20)
wv73 = dict(plotname='Brightness temperature WV 7.3µm', plotunits='[K]',
- kind='MSG_4_SEVIRI_TB', sat_channel=6, n_obs=n_obs,
- error_generate=1., error_assimilate=1.,
- cov_loc_radius_km=10)
+ 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)
ir108 = dict(plotname='Brightness temperature IR 10.8µm', plotunits='[K]',
- kind='MSG_4_SEVIRI_TB', sat_channel=9, n_obs=n_obs,
+ kind='MSG_4_SEVIRI_TB', sat_channel=9,
+ n_obs=1, obs_locations='square_array_evenly_on_grid',
error_generate=5., error_assimilate=10.,
cov_loc_radius_km=32)
radar = dict(plotname='Radar reflectivity', plotunits='[dBz]',
- kind='RADAR_REFLECTIVITY', n_obs=n_obs,
- error_generate=2.5, error_assimilate=5,
+ kind='RADAR_REFLECTIVITY',
+ n_obs=5776, obs_locations='square_array_evenly_on_grid',
+ error_generate=2.5, error_assimilate=2.5,
heights=range(2000, 14001, 2000),
- cov_loc_radius_km=20)
+ cov_loc_radius_km=1)
t = dict(plotname='Temperature', plotunits='[K]',
- kind='RADIOSONDE_TEMPERATURE', n_obs=n_obs,
- error_generate=0.2, error_assimilate=0.4,
- heights=[5000,], #range(1000, 17001, 2000),
- cov_loc_radius_km=4)
+ kind='RADIOSONDE_TEMPERATURE',
+ n_obs=22500, 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=1.5)
q = dict(plotname='Specific humidity', plotunits='[kg/kg]',
- kind='RADIOSONDE_SPECIFIC_HUMIDITY', n_obs=n_obs,
+ kind='RADIOSONDE_SPECIFIC_HUMIDITY', n_obs=1,
error_generate=0., error_assimilate=5*1e-5,
heights=[1000], #range(1000, 17001, 2000),
cov_loc_radius_km=0.1)
t2m = dict(plotname='SYNOP Temperature', plotunits='[K]',
- kind='SYNOP_TEMPERATURE', n_obs=n_obs,
+ kind='SYNOP_TEMPERATURE', n_obs=1,
error_generate=0.1, error_assimilate=1.,
cov_loc_radius_km=20)
psfc = dict(plotname='SYNOP Pressure', plotunits='[Pa]',
- kind='SYNOP_SURFACE_PRESSURE', n_obs=n_obs,
+ kind='SYNOP_SURFACE_PRESSURE', n_obs=1,
error_generate=50., error_assimilate=100.,
cov_loc_radius_km=32)
--
GitLab