move paths in config

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']
 

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'

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):

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')

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

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)