diff --git a/PE_CBL.py b/PE_CBL.py
index 8f89ef585e53134641a3339523fd67cf4b8a6b5f..4acb54260673465a0f832e3b20958e3d2af65396 100644
--- a/PE_CBL.py
+++ b/PE_CBL.py
@@ -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)