diff --git a/config/cfg.py b/config/cfg.py index 0b1d7379b12ff97336fdf865e471714c20773666..1dabc00b39d65a9ed4f467406bdd813a801cfb92 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 b292d40e43c1d23f1c610d0db2a29167c1d3419b..d574098705f2cfc9ea6187e1643b7c3e349d3a4f 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 c40a7c73b69bbc7f9c8ffb0c6f7ef198496888b0..fb5945e564d26e680f24cf399bce357d93df39c3 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 41a8c260174425e80b1c727871ce71b1c66ae077..e9df5d5c9478f5f55a00ff38db8a5cb9b6050d40 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)