new function to plot vertical profiles of obs and background ensemble

graphics.py

@@ -690,8 +690,8 @@ def plot_overview(exp,label,filename,show_spread=False,show_transformed_paramete
         ax3.set_title(r'd) Exp. %s, evolution of $p\prime\prime$'%label)
     else:
         ax3.set_title(r'd) Exp. %s, evolution of $p$'%label)
+    ax3.set_yticks(np.arange(5+1))
     ax3.set_xlabel('Time (h)')
-    ax3.set_yticks(np.arange(int(maxtime/3600)))
     ax0.set_xlim([0,maxtime/3600])
     ax1.set_xlim([0,maxtime/3600])
     ax2.set_xlim([0,maxtime/3600])
@@ -706,3 +706,29 @@ def plot_overview(exp,label,filename,show_spread=False,show_transformed_paramete
     #
     fig.savefig(filename,format='png',dpi=300)
     p.close(fig)
+
+def plot_profiles(exp,settings,filename):
+
+    # Make plots
+    fig, [ax1, ax2] = p.subplots(1,2,constrained_layout=True)
+    fig.set_size_inches(6,3)
+
+    which_cycle = 0
+    tt1 = (which_cycle+1)*settings['assimilation_interval']
+    ax1 = plot_CBL_assimilation(exp,None,which_cycle=which_cycle,ax=ax1)
+    ax1.set_ylim([0,2000])
+    ax1.set_xlim([288,298])
+    ax1.set_ylabel(r'height (m)')
+    ax1.set_title(r'a) Analysis at $t=$%u s'%tt1)
+    ax1.legend(frameon=False)
+
+    which_cycle = 31
+    tt2 = (which_cycle+1)*settings['assimilation_interval']
+    ax2 = plot_CBL_assimilation(exp,None,which_cycle=which_cycle,ax=ax2)
+    ax2.set_ylim([0,2000])
+    ax2.set_xlim([288,298])
+    ax2.set_ylabel(r'height (m)')
+    ax2.set_title(r'b) Analysis at $t=$%u s'%tt2)
+
+    fig.savefig(filename,format='png',dpi=300)
+    p.close(fig)
\ No newline at end of file