clc;clear all;close all;
load exp1.mat;
rho_c_top1 = dim.m_show(:,end);
rho_r_LCL1 = 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;
x_line = linspace(0,1.5,20);
y_line = ones(20,1) * 1;
h2=plot(x_line,y_line,'--');
set(h2,'linewidth',2,'color','black');

xlim([0 1.5]);
ylim([0 4]);
xlabel('t(h)');
ylabel('g m^{-3}');
grid on;
legend('\rho_c','\rho_th');
title('bubble top \rho_c');
set(gca,'Fontweight','bold','Fontsize',15);
set(gca,'Linewidth',2);
filename_png=['cloud rain time LCL single','.png'];
saveas(gcf,filename_png,'png');