diff --git a/scripts/assim_synth_obs.py b/scripts/assim_synth_obs.py index 37dde003a196e46bf24bc456fa548e3c80e6b02a..305c553886fa53d8f498619c1ccdced9df689958 100755 --- a/scripts/assim_synth_obs.py +++ b/scripts/assim_synth_obs.py @@ -266,7 +266,7 @@ def prepare_prior_ensemble(assim_time, prior_init_time, prior_valid_time, prior_ os.system('rm -rf '+cluster.dartrundir+'/perfect_output_*') os.system('rm -rf '+cluster.dartrundir+'/obs_seq.fina*') - os.system(cluster.python+' '+cluster.scriptsdir+'/link_rttov.py') + os.system(cluster.python+' '+cluster.scriptsdir+'/link_dart_rttov.py') def calc_obserr_WV73(Hx_nature, Hx_prior): diff --git a/scripts/link_rttov.py b/scripts/link_dart_rttov.py similarity index 100% rename from scripts/link_rttov.py rename to scripts/link_dart_rttov.py diff --git a/scripts/pre_assim.py b/scripts/pre_assim.py deleted file mode 100755 index bb2e41180af321eb185d988ace265c736c24e666..0000000000000000000000000000000000000000 --- a/scripts/pre_assim.py +++ /dev/null @@ -1,76 +0,0 @@ -import os, sys, shutil -import datetime as dt -from config.cfg import exp, cluster -from utils import symlink, copy_scp_srvx8, copy, sed_inplace, try_remove -import wrfout_add_geo - -def run(assim_time, prior_init_time, prior_valid_time, prior_path_exp): - """Prepares DART files for running filter - i.e. - - links first guess state to DART first guess filenames - - creates wrfinput_d01 files - - adds geo-reference (xlat,xlon) coords so that DART can deal with the files - - writes txt files so DART knows what input and output is - - removes probably pre-existing files which could lead to problems - """ - os.makedirs(cluster.dartrundir, exist_ok=True) - - print('prepare prior state estimate') - for iens in range(1, exp.n_ens+1): - print('link wrfout file to DART background file') - wrfout_run = prior_path_exp \ - +prior_init_time.strftime('/%Y-%m-%d_%H:%M/') \ - +str(iens) \ - +prior_valid_time.strftime('/wrfout_d01_%Y-%m-%d_%H:%M:%S') - dart_ensdir = cluster.dartrundir+'/advance_temp'+str(iens) - wrfout_dart = dart_ensdir+'/wrfout_d01' - - os.makedirs(dart_ensdir, exist_ok=True) - print('copy', wrfout_run, 'to', wrfout_dart) - copy(wrfout_run, wrfout_dart) - symlink(wrfout_dart, dart_ensdir+'/wrfinput_d01') - - # ensure prior time matches assim time (can be off intentionally) - if assim_time != prior_valid_time: - print('overwriting time in prior from nature wrfout') - os.system(cluster.ncks+' -A -v XTIME ' - +cluster.dartrundir+'/wrfout_d01 '+wrfout_dart) - - # this seems to be necessary (else wrong level selection) - wrfout_add_geo.run(cluster.dartrundir+'/../geo_em.d01.nc', wrfout_dart) - - fpath = cluster.dartrundir+'/input_list.txt' - print('writing', fpath) - try_remove(fpath) - with open(fpath, 'w') as f: - for iens in range(1, exp.n_ens+1): - f.write('./advance_temp'+str(iens)+'/wrfout_d01') - f.write('\n') - - fpath = cluster.dartrundir+'/output_list.txt' - print('writing', fpath) - try_remove(fpath) - with open(fpath, 'w') as f: - for iens in range(1, exp.n_ens+1): - f.write('./filter_restart_d01.'+str(iens).zfill(4)) - f.write('\n') - - - print('removing preassim and filter_restart') - os.system('rm -rf '+cluster.dartrundir+'/preassim_*') - os.system('rm -rf '+cluster.dartrundir+'/filter_restart*') - os.system('rm -rf '+cluster.dartrundir+'/output_mean*') - os.system('rm -rf '+cluster.dartrundir+'/output_sd*') - os.system('rm -rf '+cluster.dartrundir+'/perfect_output_*') - os.system('rm -rf '+cluster.dartrundir+'/obs_seq.fina*') - - os.system(cluster.python+' '+cluster.scriptsdir+'/link_rttov.py') - - -if __name__ == '__main__': - assim_time = dt.datetime.strptime(sys.argv[1], '%Y-%m-%d_%H:%M') - prior_init_time = dt.datetime.strptime(sys.argv[2], '%Y-%m-%d_%H:%M') - prior_valid_time = dt.datetime.strptime(sys.argv[3], '%Y-%m-%d_%H:%M') - prior_path_exp = str(sys.argv[4]) - - run(assim_time, prior_init_time, prior_valid_time, prior_path_exp) diff --git a/scripts/run_ens.jet.sh b/scripts/run_ens.jet.sh deleted file mode 100755 index 57cf7b5df386c19451433fba389a399850a87dfa..0000000000000000000000000000000000000000 --- a/scripts/run_ens.jet.sh +++ /dev/null @@ -1,27 +0,0 @@ -export SLURM_STEP_GRES=none - -## $SLURM_ARRAY_TASK_ID -echo "SLURM_ARRAY_TASK_ID:"$SLURM_ARRAY_TASK_ID - -EXPNAME="OSSE_v1.10_test" -NAMELIST="namelist.input" - -USERDIR=/jetfs/home/lkugler/ -DATADIR=$USERDIR -SRC_DIR=/jetfs/home/lkugler/compile/WRF/WRF-v4.2/run/ -IDEAL_EXE=/jetfs/home/lkugler/compile/bin/ideal.exe -WRF_EXE=/jetfs/home/lkugler/compile/bin/wrf-v4.2_v1.10.dmpar.exe - - -# VSC support: be careful with correct pinning ! -pinning=(0-11 12-23 24-35 36-47) -for ((n=1; n<=4; n++)) -do - RUNDIR=$USERDIR/run_WRF/$EXPNAME/$IENS - cd $RUNDIR - rm -r wrfout_d01_* - echo 'mpirun -genv I_MPI_PIN_PROCESSOR_LIST=${pinning[$n]} -np 12 ./wrf.exe >/dev/null 2>&1' - mpirun -genv I_MPI_PIN_PROCESSOR_LIST=${pinning[$n]} -np 12 ./wrf.exe >/dev/null 2>&1 - cd ../ -done -wait