h1=plot((qr_diff(1,:)+qs_diff(1,:)+qg_diff(1,:)+qh_diff(1,:)).*1000,[1:50],'k-','LineWidth',4);
Unrecognized function or variable 'qr_diff'.
text(-0.06,42,'d) total q','FontSize',42,'FontWeight','bold');
xticklabels({'-0.8';'-0.4';'0';'0.4';'0.8'});
xlabel('g kg^{-1}','FontSize',22,'FontWeight','bold');
aa = get(gca,'XTickLabel'); set(gca,'XTickLabel',aa,'FontSize',22,'FontWeight','bold');
yticklabels({'0';'10';'20';'30';'40';'50'});
ylabel('Model Level','FontSize',22,'FontWeight','bold');
bb = get(gca,'YTickLabel'); set(gca,'YTickLabel',bb,'FontSize',22,'FontWeight','bold');
set(gca,'TickDir','out');
set(gcf,'units','pixels','outerposition',[1 1 900 900],'windowstyle','normal')
set(gca,'units','normalized','position',[0.1300 0.1100 0.7750 0.8150],'plotboxaspectratio',[1 1 1])
saveas(gcf,[times(t,:),'_stincre_totq_vprof_',exps{nex},'.png'])