diff --git a/config/cfg.py b/config/cfg.py
index b643e1bf2670580b399293f6d8f8f2e2286ed4f2..3ae1873ac4db8c449e5abb28ce16fb79f5903c62 100755
--- a/config/cfg.py
+++ b/config/cfg.py
@@ -9,10 +9,15 @@ class ExperimentConfiguration(object):
exp = ExperimentConfiguration()
-exp.expname = "exp_v1.17_P1-1_WV_10z"
+exp.expname = "exp_v1.18_P1_nature"
exp.model_dx = 2000
-exp.n_ens = 40
-exp.n_nodes = 10
+exp.n_ens = 4
+exp.n_nodes = 1
+
+#exp.nature_wrfout = '/home/fs71386/lkugler/data/sim_archive/exp_v1.17_P1_nature/2008-07-30_06: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.<iens>.wrfprof'
+#exp.input_profile = '/home/fs71386/lkugler/wrf_profiles/data/wrf/ens/2021-05-04/raso.fc.<iens>.wrfprof'
+
# localize vertically, if it has a vertical position
# needs a horizontal scale too, to calculate the vertical normalization
@@ -53,7 +58,7 @@ psfc = dict(plotname='SYNOP Pressure', plotunits='[dBz]',
cov_loc_radius_km=32)
-exp.observations = [wv73] #radar] # 108, wv73, vis]
+exp.observations = [] #wv73] #radar] # 108, wv73, vis]
#exp.update_vars = ['T', 'QVAPOR', 'QCLOUD', 'QICE','CLDFRA']
exp.update_vars = ['U', 'V', 'T', 'PH', 'MU', 'QVAPOR', 'QCLOUD', 'QICE', 'TSK', 'CLDFRA']
diff --git a/config/clusters.py b/config/clusters.py
index 8d2709972303374e701dcfcd5238407532ce882f..b8679cec3c1166dd04a57c52b385ccc924c54f77 100755
--- a/config/clusters.py
+++ b/config/clusters.py
@@ -35,10 +35,6 @@ vsc.dart_srcdir = '/home/fs71386/lkugler/DART/DART-9.11.9/models/wrf/work'
vsc.dartrundir = '/gpfs/data/fs71386/lkugler/run_DART'
vsc.scriptsdir = '/home/fs71386/lkugler/DART-WRF/scripts/'
-vsc.nature_wrfout = '/home/fs71386/lkugler/data/sim_archive/exp_v1.17_P1_nature/2008-07-30_06:00/1/wrfout_d01_%Y-%m-%d_%H:%M:%S'
-#vsc.input_profile = '/home/fs71386/lkugler/wrf_profiles/data/wrf/ens/from_LMU/raso.nat.<iens>.wrfprof'
-vsc.input_profile = '/home/fs71386/lkugler/wrf_profiles/data/wrf/ens/2021-05-04/raso.fc.<iens>.wrfprof'
-
vsc.ideal = vsc.userdir+'/compile/bin/ideal-v4.2.2_v1.16.exe'
vsc.wrfexe = vsc.userdir+'/compile/bin/wrf-v4.3_v1.16.exe'
vsc.namelist = vsc.scriptsdir+'/../templates/namelist.input'
@@ -58,7 +54,6 @@ jet.srcdir = '/jetfs/home/lkugler/compile/WRF/WRF-4.1.5/run'
jet.scriptsdir = ''
jet.archive_base = '/jetfs/home/lkugler/data_jetfs/sim_archive/'
jet.dartrundir = '/jetfs/home/lkugler/DART-WRF/rundir'
-jet.nature_wrfout = '/raid61/scratch/lkugler/VSC/sim_archive/OSSE_v1.10_LMU+shear/2/single/wrfout_d01_%Y-%m-%d_%H:%M:%S'
jet.ideal = jet.userdir+'/compile/bin/ideal.exe'
jet.wrfexe = jet.userdir+'/compile/bin/wrf-v4.2_v1.10.dmpar.exe'
diff --git a/scripts/assim_synth_obs.py b/scripts/assim_synth_obs.py
index 794abfd93d840ef52179a752eff6a88d4f558211..fdbdc44ac113d30aeb94564d61775c61f6650920 100755
--- a/scripts/assim_synth_obs.py
+++ b/scripts/assim_synth_obs.py
@@ -237,11 +237,11 @@ def obs_operator_nature(time):
def link_nature_to_dart_truth(time):
# get wrfout_d01 from nature run
- shutil.copy(time.strftime(cluster.nature_wrfout),
+ shutil.copy(time.strftime(exp.nature_wrfout),
cluster.dartrundir+'/wrfout_d01')
# DART may need a wrfinput file as well, which serves as a template for dimension sizes
symlink(cluster.dartrundir+'/wrfout_d01', cluster.dartrundir+'/wrfinput_d01')
- print('linked', time.strftime(cluster.nature_wrfout), 'to', cluster.dartrundir+'/wrfout_d01')
+ print('linked', time.strftime(exp.nature_wrfout), 'to', cluster.dartrundir+'/wrfout_d01')
def prepare_nature_dart(time):
diff --git a/scripts/assim_synth_obs_separate_obstypes.py b/scripts/assim_synth_obs_separate_obstypes.py
index 6689af0979ad3db854c57dd8874b0d3fa7a47dd7..9ae4a740210bd71f9c14d85a878fe1e0baaf41c7 100755
--- a/scripts/assim_synth_obs_separate_obstypes.py
+++ b/scripts/assim_synth_obs_separate_obstypes.py
@@ -189,7 +189,7 @@ def obs_operator_nature(time):
def link_nature_to_dart_truth(time):
# get wrfout_d01 from nature run
- shutil.copy(time.strftime(cluster.nature_wrfout),
+ shutil.copy(time.strftime(exp.nature_wrfout),
cluster.dartrundir+'/wrfout_d01')
# DART may need a wrfinput file as well, which serves as a template for dimension sizes
symlink(cluster.dartrundir+'/wrfout_d01', cluster.dartrundir+'/wrfinput_d01')
diff --git a/scripts/prepare_nature.py b/scripts/prepare_nature.py
index 549bd4b4d53a49460563fde1fddc001291ae79d9..b0ed8ccac3e8ebe572b629ba17f36edad7ac035b 100755
--- a/scripts/prepare_nature.py
+++ b/scripts/prepare_nature.py
@@ -6,7 +6,7 @@ from utils import symlink, copy
time = dt.datetime.strptime(sys.argv[1], '%Y-%m-%d_%H:%M')
# get wrfout_d01 from nature run
-shutil.copy(time.strftime(cluster.nature_wrfout),
+shutil.copy(time.strftime(exp.nature_wrfout),
cluster.dartrundir+'/wrfout_d01')
import wrfout_add_geo
diff --git a/scripts/prepare_wrfinput.py b/scripts/prepare_wrfinput.py
index 1a8aca58f97b752f03728f186652a17c9fe7447e..981c31fc4da5062fbed52676d83a36e99007ae5e 100755
--- a/scripts/prepare_wrfinput.py
+++ b/scripts/prepare_wrfinput.py
@@ -9,7 +9,7 @@ init_time = dt.datetime.strptime(sys.argv[1], '%Y-%m-%d_%H:%M')
for iens in range(1, exp.n_ens+1):
print('preparing ens', iens)
- input_prof = (cluster.input_profile).replace('<iens>', str(iens).zfill(3))
+ input_prof = (exp.input_profile).replace('<iens>', str(iens).zfill(3))
rundir = cluster.wrf_rundir(iens)
os.makedirs(rundir, exist_ok=True)