diff --git a/config/cfg.py b/config/cfg.py index c0340bd94405b672c4c9761ef9c6a4aba8c18f64..0d0f93647640b616b18a5511bdbd5ff694d7831d 100755 --- a/config/cfg.py +++ b/config/cfg.py @@ -2,13 +2,14 @@ from dartwrf.utils_workflow import ExperimentConfiguration from config import clusters # from . = problem in archivedir cluster = clusters.srvx1 # change cluster configuration here -exp = ExperimentConfiguration() -exp.expname = "test_srvx1" #"exp_v1.22_P3_wbub7_WV62_obs10_loc20_oe1" +exp = utils.ExperimentConfiguration() +exp.expname = "exp_v1.22_P2_rr_VIS+WV73_obs10_loc20_inf2" exp.model_dx = 2000 exp.n_ens = 40 exp.filter_kind = 1 -exp.inflation = True +exp.prior_inflation = 2 +exp.post_inflation = 0 exp.sec = True exp.reject_smallFGD = False exp.cov_loc_vert_km_horiz_km = (3, 20) @@ -95,7 +96,7 @@ psfc = dict(plotname='SYNOP Pressure', plotunits='[Pa]', cov_loc_radius_km=32) -exp.observations = [t] +exp.observations = [vis,wv73] exp.update_vars = ['U', 'V', 'W', 'THM', 'PH', 'MU', 'QVAPOR', 'QCLOUD', 'QICE', 'PSFC'] #exp.update_vars = ['U', 'V', 'W', 'T', 'PH', 'MU', 'QVAPOR', 'PSFC'] diff --git a/dartwrf/assim_synth_obs.py b/dartwrf/assim_synth_obs.py index ce777c3c1f8db76b74f81358a4b7c63b9b94f074..3b66917bb3e150b071ad3ea1c78cb32634614f7b 100755 --- a/dartwrf/assim_synth_obs.py +++ b/dartwrf/assim_synth_obs.py @@ -42,7 +42,8 @@ def set_DART_nml(just_prior_values=False): "<adjust_obs_impact>": '.true.' if exp.adjust_obs_impact else '.false.', "<filter_kind>": str(int(exp.filter_kind)), "<sampling_error_correction>": '.true.' if exp.sec else '.false.', - "<post_inflation>": '4' if exp.inflation else '0', + "<prior_inflation>": str(exp.prior_inflation), + "<post_inflation>": str(exp.post_inflation), "<n_ens>": str(int(exp.n_ens)), "<cov_loc_radian>": "0.00000001", # dummy value, used for types not mentioned below "<list_obstypes>": "'" + "','".join(list_obstypes) + "'", diff --git a/templates/input.nml b/templates/input.nml index ee5c2a72ccdee5cb98f989789b8b2e1fabd6ee7c..3cc9b63416f67b10e1660c02c3ec9aa3ce1c083e 100644 --- a/templates/input.nml +++ b/templates/input.nml @@ -57,7 +57,7 @@ output_sd = .true. write_all_stages_at_end = .false. - inf_flavor = 0, <post_inflation>, + inf_flavor = <prior_inflation>, <post_inflation>, inf_initial_from_restart = .true., .false., inf_sd_initial_from_restart = .true., .false., inf_initial = 1.00, 0.90,