fix: vertloc, .000

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
 

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

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

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)