% plot 85W latitude-height section
x=repmat(la(i1:i2),[3 1]); y=repmat([4.8;5.0;NaN],[1 i2-i1+1]);
y(1,hour(i1:i2)==23 | hour(i1:i2)==0)=4.6;
j=(hour(i1:i2)==11 | hour(i1:i2)==12); % midday for oct day labels
iz=hz<=5e3;

clf
b2rcolormap(24);
subplot(4,1,1,'align')
[c,h]=contourf(la(i1:i2),hz(iz)*1e-3,thetaz(iz,i1:i2),286:2:320);
set(h,'edgecolor','none')
set(gca,'ylim',[0 5])
colorbar
hold on
plot(lat_cb(1:end-1),cb2_85/1e3,'y.','markersize',2.5)
plot(la(i1:i2),hinvbase(i1:i2)/1e3,'m.')
plot(la(i1:i2),hinvtop(i1:i2)/1e3,'m.')
plot(x,y,'w')
text(x(1,j),4.8+zeros(1,sum(j)),num2str(floor(octday(j))),'color','w','horizontalalignment','right');
title({[cruise ' ' yyyy ' southward section'] 'potential temperature (K)'})
subplot(4,1,2,'align')
[c,h]=contourf(la(i1:i2),hz(iz)*1e-3,qz(iz,i1:i2),0:1:15);
set(h,'edgecolor','none')
set(gca,'ylim',[0 5])
colorbar
hold on
plot(x,y,'w')
text(x(1,j),4.8+zeros(1,sum(j)),num2str(floor(octday(j))),'color','w','horizontalalignment','right');
title('specific humidity (g/kg)');
subplot(4,1,3,'align')
[c,h]=contourf(la(i1:i2),hz(iz)*1e-3,uz(iz,i1:i2),-12:12);
set(h,'edgecolor','none')
colorbar
set(gca,'ylim',[0 5]);
colorbar
hold on
plot(x,y,'k')
text(x(1,j),4.8+zeros(1,sum(j)),num2str(floor(octday(j))),'color','k','horizontalalignment','right');
title('zonal velocity (m s ^{-1})')
ylabel('height (km)')
subplot(4,1,4,'align')
[c,h]=contourf(la(i1:i2),hz(iz)*1e-3,vz(iz,i1:i2),-12:12);
set(h,'edgecolor','none')
colorbar
set(gca,'ylim',[0 5])
colorbar
hold on
plot(x,y,'k')
text(x(1,j),4.8+zeros(1,sum(j)),num2str(floor(octday(j))),'color','k','horizontalalignment','right');
title('meridional velocity (m s^{-1})')
xlabel('latitude (\circ)')
orient tall
