From 14d2581fbe5c18f1f1774fd5aaff88062d76527a Mon Sep 17 00:00:00 2001
From: Stefano Serafin <serafin@jet01.img.univie.ac.at>
Date: Thu, 27 Feb 2025 10:36:47 +0100
Subject: [PATCH] Fixes in Fig 1
---
PE_CBL.py | 20 +++++++++++---------
1 file changed, 11 insertions(+), 9 deletions(-)
diff --git a/PE_CBL.py b/PE_CBL.py
index 96a9549..7fd1011 100644
--- a/PE_CBL.py
+++ b/PE_CBL.py
@@ -41,7 +41,7 @@ if __name__ == '__main__':
########################################################################
# Manuscript figures
- fig01 = False # CHECKED OK # Sensitivity to p
+ fig01 = True # CHECKED OK # Sensitivity to p
fig02 = False # CHECKED OK # Successful PE experiment
fig03 = False # CHECKED OK # Consistency checks
fig04 = False # CHECKED OK # Experiments with varying weights
@@ -49,14 +49,14 @@ if __name__ == '__main__':
fig06 = False # CHECKED OK # Analysis increments in model error experiments
fig07 = False # CHECKED OK # Experiment with error growth
fig08 = False # CHECKED OK # Comparison successful vs unsuccessful PE
- fig09 = False # CHECKED OK # Comparison EAKF vs LETKF
- fig10 = True # CHECKED OK # DA diagnostics (profiles) for perfect-obs experiments
+ fig09 = False # CHECKED OK # Comparison EAKF vs LETKF
+ fig10 = False # CHECKED OK # DA diagnostics (profiles) for perfect-obs experiments
# Other figures
- opt01 = False # CHECKED OK # assimilation of a single observation
- opt02 = False # CHECKED OK # assimilation of profiles at two times
- opt03 = False # CHECKED OK # assimilation of real observations
- opt04 = False # CHECKED OK # assimilation of real observations, different spinup times
+ opt01 = False # CHECKED OK # Assimilation of a single observation
+ opt02 = False # CHECKED OK # Assimilation of profiles at two times
+ opt03 = False # CHECKED OK # Assimilation of real observations
+ opt04 = False # CHECKED OK # Assimilation of real observations, different spinup times
opt05 = False # CHECKED OK # Check impact of localization cutoff
opt06 = False # CHECKED OK # Check impact of RTPS variance inflation factor for parameters
opt07 = False # CHECKED OK # Check impact of RTPS variance inflation factor for state vector
@@ -78,7 +78,7 @@ if __name__ == '__main__':
maxtime = 21600
cbl_settings['maxtime'] = maxtime
- # Panel a) deterministic run using default settings
+ # Panel c) deterministic run using default settings
cbl_det = CBL(cbl_settings)
cbl_det.initialize(1)
cbl_det.run(output_full_history=True)
@@ -93,18 +93,20 @@ if __name__ == '__main__':
theta_profiles = []
for pfac in p_factors:
cbl_settings['pfac'] = pfac
+ cbl_settings['do_parameter_estimation'] = False
cbl_pf = CBL(cbl_settings)
cbl_pf.initialize(1)
cbl_pf.run()
theta_profiles.append(cbl_pf.x[:cbl_pf.nz])
- # Panel c) spread induced by p
+ # Panel d) spread induced by p
# Do a free ensemble run (ensemble size set expliclity)
# Set size to 200 members
with open(workdir+'/default_cbl.json', 'r') as fp:
cbl_settings_free = json.load(fp)
cbl_settings_free['maxtime'] = maxtime
cbl_settings_free['perturb_ensemble_state'] = False
+ cbl_settings_free['do_parameter_estimation'] = False
cbl_free = CBL(cbl_settings_free)
cbl_free.initialize(200)
cbl_free.run(output_full_history=True)
--
GitLab