From feff66d872faca7c54f44de18a5b3f8355d64338 Mon Sep 17 00:00:00 2001 From: lkugler <lukas.kugler@gmail.com> Date: Mon, 22 May 2023 18:20:50 +0200 Subject: [PATCH] rename dart_rundir --- dartwrf/assim_synth_obs.py | 8 ++++---- tests/test_dart-rttov.py | 8 ++++---- 2 files changed, 8 insertions(+), 8 deletions(-) diff --git a/dartwrf/assim_synth_obs.py b/dartwrf/assim_synth_obs.py index 14be283..39ac427 100755 --- a/dartwrf/assim_synth_obs.py +++ b/dartwrf/assim_synth_obs.py @@ -170,7 +170,7 @@ def archive_filteroutput(time): for ftype in ftypes: for iens in range(1, exp.n_ens + 1): fname = "/"+ftype+"_member_" + str(iens).zfill(4) + ".nc" - copy(cluster.dartrundir + fname, archive_assim + fname) + copy(cluster.dart_rundir + fname, archive_assim + fname) for f in ["output_mean.nc", "output_sd.nc"]: # copy mean and sd to archive copy(cluster.dart_rundir + "/" + f, archive_assim + "/" + f) @@ -226,7 +226,7 @@ def set_obserr_assimilate_in_obsseqout(oso, outfile="./obs_seq.out"): where_oso_iskind = oso.df.kind == kind if obscfg["error_assimilate"] == False: - osf_prior = obsseq.ObsSeq(cluster.dartrundir + "/obs_seq.final") # this file will be generated by `evaluate()` + osf_prior = obsseq.ObsSeq(cluster.dart_rundir + "/obs_seq.final") # this file will be generated by `evaluate()` where_osf_iskind = osf_prior.df.kind == kind @@ -240,7 +240,7 @@ def set_obserr_assimilate_in_obsseqout(oso, outfile="./obs_seq.out"): oso.to_dart(outfile) def qc_obs(time, oso): - osf_prior = obsseq.ObsSeq(cluster.dartrundir + "/obs_seq.final") + osf_prior = obsseq.ObsSeq(cluster.dart_rundir + "/obs_seq.final") # obs should be superobbed already! for i, obscfg in enumerate(exp.observations): @@ -555,7 +555,7 @@ def main(time, prior_init_time, prior_valid_time, prior_path_exp): evaluate(time, output_format="%Y-%m-%d_%H:%M_obs_seq.final-eval_prior_allobs") print(" assign observation-errors for assimilation ") - set_obserr_assimilate_in_obsseqout(oso, outfile=cluster.dartrundir + "/obs_seq.out") + set_obserr_assimilate_in_obsseqout(oso, outfile=cluster.dart_rundir + "/obs_seq.out") if do_QC: print(" 2.3) reject observations? ") diff --git a/tests/test_dart-rttov.py b/tests/test_dart-rttov.py index 4f617c7..4d502b3 100644 --- a/tests/test_dart-rttov.py +++ b/tests/test_dart-rttov.py @@ -35,7 +35,7 @@ else: datadir = '/jetfs/home/lkugler/data/sim_archive/' img_dir = '/jetfs/home/lkugler/analysis_plots/' -dartrundir = datadir+'/../run_DART/test_jet' +dart_rundir = datadir+'/../run_DART/test_jet' def test_rttov(): @@ -59,9 +59,9 @@ def test_rttov(): shutil.copy(datadir+'/'+case+'/2008-07-30_12:00/1/'+time.strftime('wrfout_d01_2008-07-30_%H:%M:00'), - dartrundir + "/wrfout_d01") + dart_rundir + "/wrfout_d01") - wrfout_add_geo.run(datadir+ "/geo_em.d01.nc", dartrundir + "/wrfout_d01") + wrfout_add_geo.run(datadir+ "/geo_em.d01.nc", dart_rundir + "/wrfout_d01") osq.create_obsseqin_alltypes(time, [obscfg]) @@ -69,7 +69,7 @@ def test_rttov(): aso.run_perfect_model_obs(nproc=1) - obs = obsseq.ObsSeq(dartrundir+'/obs_seq.out') + obs = obsseq.ObsSeq(dart_rundir+'/obs_seq.out') obs_map = obs.df.truth.values.reshape(150,150) -- GitLab