T=0:1:30;
dt=1.5;
r=.8;
Le=2500;
cp=1000;
R=287;
b=cp/Le*dt./qsea(T)./(1-r+.62*Le*1000/R./(T+273.16).^2*dt);
dq=qsea(T).*(1-r+.61*Le*1000/R*dt./(T+273.16).^2);
rat=(1+1.61*r*qsea(T)/1000).*(1+Le/cp./(T+273.16).*b*r.*qsea(T));
rat1=(1+1.61*qsea(T)/1000).*(1+dt./(T+273.16)*r.*qsea(T)./dq);

figure;plot(T,b);
figure;plot(T,rat);
