diff --git a/PE_CBL.py b/PE_CBL.py
index 5db83c7d35027e694825952991c09acd6cf38f10..96a95492a52d36f7dec7d6d99cfbf2c53d1edf6f 100644
--- a/PE_CBL.py
+++ b/PE_CBL.py
@@ -49,7 +49,7 @@ 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 = True  # CHECKED OK # Comparison EAKF vs LETKF
+    fig09 = False  # CHECKED OK # Comparison EAKF vs LETKF
     fig10 = True  # CHECKED OK # DA diagnostics (profiles) for perfect-obs experiments
     
     # Other figures
@@ -605,13 +605,13 @@ if __name__ == '__main__':
         plot_diagnostics(exp1,show='errors',axes=[ax1,ax2,ax3,ax4],label='B$_4$ with PE',linecolor='#1f77b4',zmax=1500)
         plot_diagnostics(exp2,show='errors',axes=[ax1,ax2,ax3,ax4],label='B$_4$ without PE',linecolor='#9EC9E9',zmax=1500)
         ax1.set_title('a) Analysis error')
-        ax1.set_xlabel(r'RMSE$_{y_a}$')
+        ax1.set_xlabel(r'RMSE$_{\overline{y}_a}$')
         ax2.set_title('b) First-guess error')
-        ax2.set_xlabel(r'RMSE$_{y_b}$')
+        ax2.set_xlabel(r'RMSE$_{\overline{y}_b}$')
         ax3.set_title('c) Error reduction')
-        ax3.set_xlabel(r'RMSE$_{y_b}$-RMSE$_{y_a}$')
+        ax3.set_xlabel(r'RMSE$_{\overline{y}_b}$-RMSE$_{\overline{y}_a}$')
         ax4.set_title('d) Consistency of innovations')
-        ax4.set_xlabel(r'$(\sigma_{y_b}+\sigma_{y_o})$/RMSE$_{y_b}$')
+        ax4.set_xlabel(r'$(\sigma_{y_b}+\sigma_{y_o})$/RMSE$_{\overline{y}_b}$')
         ax1.set_ylabel('height (m)')
         ax3.set_ylabel('height (m)')
         ax1.legend(frameon=False,fontsize=8,loc=7)