diff --git a/PE_CBL.py b/PE_CBL.py
index 788a9a0436bd15fca7431c693dd7485775f27ee8..f87cf84c88bad6083e598e1cfed4fc1a6b6c4135 100644
--- a/PE_CBL.py
+++ b/PE_CBL.py
@@ -14,7 +14,7 @@ if __name__ == '__main__':
 
     # Global settings
     nobs = 16
-    nens = 20
+    nens = 200
     sigma_b_init = 0.1 # increased by a factor 10 in experiment B
     sigma_o_gn = 0.0 # 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
@@ -227,6 +227,7 @@ if __name__ == '__main__':
         ax3 = plot_CBL_PE(exp_A,None,ax=ax3)
         ax3.set_title(r'd) Exp. A, evolution of $p$')
         ax3.set_xlabel('Time (h)')
+        ax3.set_yticks([0,1,2,3,4,5])
         p.colorbar(c0,orientation='horizontal')
         p.colorbar(c1,orientation='horizontal')
         p.colorbar(c2,orientation='horizontal')
@@ -532,6 +533,7 @@ if __name__ == '__main__':
         ax3 = plot_CBL_PE(exp_D,None,ax=ax3)
         ax3.set_title(r'd) Exp. D, evolution of $p$')
         ax3.set_xlabel('Time (h)')
+        ax3.set_yticks([0,1,2,3,4,5])
         p.colorbar(c0,orientation='horizontal')
         p.colorbar(c1,orientation='horizontal')
         p.colorbar(c2,orientation='horizontal')
@@ -590,6 +592,7 @@ if __name__ == '__main__':
         ax3 = plot_CBL_PE(exp_E,None,ax=ax3)
         ax3.set_title(r'd) Exp. E, evolution of $p$')
         ax3.set_xlabel('Time (h)')
+        ax3.set_yticks([0,1,2,3,4,5])
         p.colorbar(c0,orientation='horizontal')
         p.colorbar(c1,orientation='horizontal')
         p.colorbar(c2,orientation='horizontal')