From ec2d7ad3c303df8930c5fa8cecb5eeed0d0b74de Mon Sep 17 00:00:00 2001
From: Kugler Lukas <lukas.kugler@univie.ac.at>
Date: Wed, 21 Oct 2020 15:30:23 +0200
Subject: [PATCH] fix: vertloc, .000
---
config/cfg.py | 6 +++---
scripts/create_obsseq.py | 2 +-
scripts/gen_synth_obs.py | 4 ++--
scripts/run_obs_diag.py | 2 +-
4 files changed, 7 insertions(+), 7 deletions(-)
diff --git a/config/cfg.py b/config/cfg.py
index 0b1d737..1dabc00 100755
--- a/config/cfg.py
+++ b/config/cfg.py
@@ -9,11 +9,11 @@ class ExperimentConfiguration(object):
exp = ExperimentConfiguration()
-exp.expname = "exp_v1.12_LMU_radar_vertloc1km"
+exp.expname = "exp_v1.12_LMU_so_radar_vertloc1km"
exp.model_dx = 2000
exp.timestep = 10
-exp.n_ens = 40
-exp.n_nodes = 10
+exp.n_ens = 20
+exp.n_nodes = 5
n_obs = 64 # radar: n_obs for each observation height level
diff --git a/scripts/create_obsseq.py b/scripts/create_obsseq.py
index b292d40..d574098 100755
--- a/scripts/create_obsseq.py
+++ b/scripts/create_obsseq.py
@@ -140,7 +140,7 @@ obs_kind_definitions
lon = coords[i_obs-1][1]
lat = coords[i_obs-1][0]
- hgt_m = str(coords[i_obs-1][2])+'.000'
+ hgt_m = str(coords[i_obs-1][2])
lon_rad = str(degr_to_rad(lon))
lat_rad = str(degr_to_rad(lat))
diff --git a/scripts/gen_synth_obs.py b/scripts/gen_synth_obs.py
index c40a7c7..fb5945e 100755
--- a/scripts/gen_synth_obs.py
+++ b/scripts/gen_synth_obs.py
@@ -88,7 +88,7 @@ def set_input_nml(sat_channel=False, just_prior_values=False,
'<cov_loc_radian>': str(cov_loc_radian)}
if cov_loc_vert_km:
- cov_loc_vert_rad = cov_loc_vert_km/cov_loc_radian
+ cov_loc_vert_rad = cov_loc_vert_km*1000/cov_loc_radian
options['<horiz_dist_only>'] = '.false.'
options['<vert_norm_hgt>'] = str(cov_loc_vert_rad)
else:
@@ -113,7 +113,7 @@ if __name__ == "__main__":
fpath_obs_coords = cluster.archivedir()+time.strftime('/%Y-%m-%d_%H:%M/obs_coords.pkl')
# remove any existing observation files
- os.system('rm -f '+cluster.dartrundir+'/obs_seq_*.out')
+ os.chdir(cluster.dartrundir); os.system('rm -f obs_seq_*.out obs_seq.in obs_seq.final')
# loop over observation types
for i_obs, obscfg in enumerate(exp.observations):
diff --git a/scripts/run_obs_diag.py b/scripts/run_obs_diag.py
index 41a8c26..e9df5d5 100644
--- a/scripts/run_obs_diag.py
+++ b/scripts/run_obs_diag.py
@@ -19,7 +19,7 @@ def run(folder_obs_seq_final):
for obserr_iszero in ['.true.', '.false.']:
print('ensure correct input.nml')
- copy(cluster.scriptsdir+'/../templates/input.nml',
+ copy(cluster.scriptsdir+'/../templates/input.prioronly.nml',
rundir_program+'/input.nml')
sed_inplace(rundir_program+'/input.nml', '<n_ens>', str(int(exp.n_ens)))
sed_inplace(rundir_program+'/input.nml', '<zero_error_obs>', obserr_iszero)
--
GitLab