diff --git a/PE_CBL.py b/PE_CBL.py
index 750dc0a7783f5b75b54fdfeeae046cceec8fa278..5e346ac8fe5789cd44e7feaa12a751ad9d4b5b91 100644
--- a/PE_CBL.py
+++ b/PE_CBL.py
@@ -95,20 +95,6 @@ if __name__ == '__main__':
     assert np.mod(default_da_settings['assimilation_interval'],integration_dt)==0,\
         "Assimilation interval must be an integer multiplier of model dt"
 
-    # Decide what figures to plot
-    fig01 = False
-    fig02 = False
-    fig03 = False
-    fig04 = False
-    fig05 = False
-    fig06 = False
-    fig07 = False
-    fig08 = True
-    
-    # Other possible, optional, plots
-    opt01 = False # assimilation of a single observation
-    opt02 = False # assimilation of profiles at two times
-
     # Whether or not to run experiments without parameter estimation
     # Applies only to sets B&C; no-PE run is always computed for experiments A&D
     noPE_runs = False
@@ -423,6 +409,20 @@ if __name__ == '__main__':
 
     ########################################################################
 
+    # Decide what figures to plot
+    fig01 = False
+    fig02 = False
+    fig03 = False
+    fig04 = False
+    fig05 = False
+    fig06 = False
+    fig07 = False
+    fig08 = True
+    
+    # Other possible, optional, plots
+    opt01 = False # assimilation of a single observation
+    opt02 = False # assimilation of profiles at two times
+
     if fig01:
 
         # Create a copy of the default settings