From d04078822276996ca97e5f74556059aa8a140b4e Mon Sep 17 00:00:00 2001
From: lkugler <lukas.kugler@gmail.com>
Date: Wed, 29 Sep 2021 16:38:31 +0200
Subject: [PATCH] move paths in config
---
config/cfg.py | 13 +++++++++----
config/clusters.py | 5 -----
scripts/assim_synth_obs.py | 4 ++--
scripts/assim_synth_obs_separate_obstypes.py | 2 +-
scripts/prepare_nature.py | 2 +-
scripts/prepare_wrfinput.py | 2 +-
6 files changed, 14 insertions(+), 14 deletions(-)
diff --git a/config/cfg.py b/config/cfg.py
index b643e1b..3ae1873 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 8d27099..b8679ce 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 794abfd..fdbdc44 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 6689af0..9ae4a74 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 549bd4b..b0ed8cc 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 1a8aca5..981c31f 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)
--
GitLab