%cbl_plot_gm.m
%plotting program for cbl_growth model
%x = fetch (km); h = cbl height; thetam= mixed layer potential temp.; 
%hs = heat flux

 xmin = 0; xmax = 50; 			%km
 hmin = 0; hmax = 800; 			%m
 thetamin = 250; thetamax = 260; 	%K
 hsmin = 0; hsmax = 250; 		%W/m2 

 subplot(3,1,1), plot(x,h ,'b-');
hold on; 
 set(gca, 'fontsize', 12);
 xlabel('fetch (km)');
 ylabel('CIBL height (m)');
 axis([xmin xmax hmin hmax]);

 subplot(3,1,2), plot(x,thetam,'b-');
hold on;
 set(gca, 'fontsize', 12);
 xlabel('fetch (km)');
 ylabel('potential temperature (K)');
 axis([xmin xmax thetamin thetamax]);

 subplot(3,1,3), plot(x,hs,'b-');
hold on;
% subplot(3,1,3), plot(x,lhf,'b--');
 set(gca, 'fontsize', 12);
 xlabel('fetch (km)');
 ylabel('heat fluxes (W m^{-2})');
 axis([xmin xmax hsmin hsmax]);
 orient tall