% plot 20S lontitude-height section
x=repmat(lo,[3 1]); y=repmat([4.8;5.0;NaN],[1 length(lo)]);
y(1,hour==23 | hour==0)=4.6;
j=(hour==11 | hour==12); % midday for oct day labels
iz=hz<=5e3;
icbst=find(time_cb>301.75,1,'first');

clf
b2rcolormap(24);
subplot(4,1,1,'align')
[c,h]=contourf(lo(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(lon_cb(icbst:end-1),cb2_85(icbst:end)/1e3,'y.','markersize',2.5)
plot(lo(i1:i2),hinvbase(i1:i2)/1e3,'m.')
plot(lo(i1:i2),hinvtop(i1:i2)/1e3,'m.')
plot(x(:,ii),y(:,ii),'w')
text(x(1,ii&j),4.8+zeros(1,sum(ii&j)),num2str(floor(octday(ii&j))),'color','w','horizontalalignment','right');
title({[cruise ' ' yyyy ' eastward 20\circS section'] 'potential temperature (K)'})
subplot(4,1,2,'align')
[c,h]=contourf(lo(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(:,ii),y(:,ii),'w')
text(x(1,ii&j),4.8+zeros(1,sum(ii&j)),num2str(floor(octday(ii&j))),'color','w','horizontalalignment','right');
title('specific humidity (g/kg)');
subplot(4,1,3,'align')
[c,h]=contourf(lo(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(:,ii),y(:,ii),'k')
text(x(1,ii&j),4.8+zeros(1,sum(ii&j)),num2str(floor(octday(ii&j))),'color','w','horizontalalignment','right');
title('zonal velocity (m s ^{-1})')
ylabel('height (km)')
subplot(4,1,4,'align')
[c,h]=contourf(lo(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(:,ii),y(:,ii),'k')
text(x(1,ii&j),4.8+zeros(1,sum(ii&j)),num2str(floor(octday(ii&j))),'color','w','horizontalalignment','right');
title('meridional velocity (m s^{-1})')
xlabel('longitude (\circ)')
orient tall
