clc;clear all;close all;
load exp1.mat;
rho_c_top1 = dim.m_show(:,end);
rho_r_LCL1 = dim.rho_r_LCL;
load exp2.mat;
rho_c_top2 = dim.m_show(:,end);
rho_r_LCL2 = dim.rho_r_LCL;
load exp3.mat;
rho_c_top3 = dim.m_show(:,end);
rho_r_LCL3 = dim.rho_r_LCL;
tpos_hit = dim.tpos_hit';   % falling rain time coordinate
tpos = dim.tpos;    % cloud water time coordinate

tpos_hit = tpos_hit - 1e3/dim.VT;

figure(1);
h1=plot(tpos/3600,rho_c_top1*1e3);
set(h1,'linewidth',2,'color','black');
hold on;
h1=plot(tpos/3600,rho_c_top2*1e3);
set(h1,'linewidth',2,'color','red');
hold on;
h1=plot(tpos/3600,rho_c_top3*1e3);
set(h1,'linewidth',2,'color','green');
hold on;

h2=plot(tpos_hit(1:end-1)/3600,rho_r_LCL1(1:end-1)*1e3,':');
set(h2,'linewidth',2,'color','black');
hold on;
h2=plot(tpos_hit(1:end-1)/3600,rho_r_LCL2(1:end-1)*1e3,':');
set(h2,'linewidth',2,'color','red');
hold on;
h2=plot(tpos_hit(1:end-1)/3600,rho_r_LCL3(1:end-1)*1e3,':');
set(h2,'linewidth',2,'color','green');

xlim([0 1.5]);
ylim([0 4]);
xlabel('t(h)');
ylabel('g m^{-3}');
grid on;
legend('\rho_{c} ref','\rho_{c} \rho_{th}=2gm^{-3}','\rho_{c} \epsilon=0.5km^{-1}','\rho_r ref','\rho_r \rho_{th}=2gm^{-3}','\rho_r \epsilon=0.5km^{-1}');
set(gca,'Fontweight','bold','Fontsize',15);
set(gca,'Linewidth',2);
filename_png=['Figure3_cloud rain time LCL','.png'];
filename_eps=['Figure3_cloud rain time LCL','.eps'];
saveas(gcf,filename_png,'png');
saveas(gcf,filename_eps,'epsc');