clear
close all
Znight=300;
Zday=2500;
time=[0:0.1:24];
PBL=Znight;
PBL=Znight+(Zday-Znight)*exp(-(time-14).^2/20);
%PBL=Zday*sin(time/6-pi/2);
subplot(3,1,1)
plot(time,PBL)
ylabel('PBLH')
axis([0 24 -50 5000])
hold on
T=[0,0, time, 24, 24,0];
Z=[0 , Znight,PBL, Znight,0,0];
fill(T,Z,'g')
subplot(3,1,2)
AOD=0.1;
EXT=AOD./PBL*1000;
plot(time,EXT)
set(gca,'xlim',[0 24])
ylabel('EXT')

subplot(3,1,3)
SSC=3;
PM10=EXT/3*1000;
plot(time,PM10)
set(gca,'xlim',[0 24])
xlabel('time')
ylabel('PM10')
subplotfix
SetFontSize(13)