Skip to content
Snippets Groups Projects
Commit 392b2e27 authored by Stefano Serafin's avatar Stefano Serafin
Browse files

cleared useless bits of code

parent 06f0df1c
Branches
Tags
No related merge requests found
......@@ -25,24 +25,7 @@ def load_or_run(settings):
if __name__ == '__main__':
# Global settings
nens = 20
sigma_b_init = 0.1 # increased by a factor 10 in experiment B
sigma_o_gn = 0.1 # set equal to sigma_o_as in experiments C and F
sigma_o_as = 0.1 # increased by a factor 10 in experiments C and F
# Define seed for random number generation, to ensure reproducibility
rnseed = 181612
# Safety check on time steps
#integration_dt = 0.25*default_cbl_settings["dz"]**2/default_cbl_settings["Kmax"]
#assert np.mod(default_da_settings['assimilation_interval'],integration_dt)==0,\
# "Assimilation interval must be an integer multiplier of model dt"
# Whether or not to run experiments without parameter estimation
# Applies only to sets B-C-D; no-PE run is always computed for experiments A
noPE_runs = False
##################################################################################
# Check if working directory exists, if not, create it.
workdir = "./runs/seed=181612_enssize=20_EAKF_6hrs/"
if not os.path.exists(workdir):
......@@ -91,6 +74,7 @@ if __name__ == '__main__':
da_settings_D2 = json.load(fp)
exp_D2 = load_or_run(da_settings_D2)
noPE_runs = False
if noPE_runs:
with open('./config_files/da_exp_B1_noPE.json', 'r') as fp:
da_settings_B1_noPE = json.load(fp)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment