% plot latitude-pressure (vertical) section
clf
b2rcolormap(24);
subplot(4,1,1,'align')
[c,h]=contourf(la(i1:i2),pint,thetaint(:,i1:i2),286:2:320);
set(h,'edgecolor','none')
set(gca,'ylim',[500 1020])
set(gca,'ydir','reverse')
colorbar
hold onx

plot(x,y,'w')
for j=2:4:length(octday);
    text(x(1,j),530,num2str(round(octday(j))),'color','w','horizontalalignment','right');
end
title({[cruise ' ' yyyy ' southward section'] 'potential temperature (K)'})
subplot(4,1,2,'align')
[c,h]=contourf(la(i1:i2),pint,qint(:,i1:i2),0:1:13);
set(h,'edgecolor','none')
set(gca,'ylim',[500 1020])
set(gca,'ydir','reverse')
colorbar
hold on
plot(x,y,'w')
for j=2:4:length(yday); text(x(1,j),530,num2str(round(octday(j))),'color','w','horizontalalignment','right'); end
title('specific humidity (g/kg)');
subplot(4,1,3,'align')
[c,h]=contourf(la(i1:i2),pint,uint(:,i1:i2),-10:10);
set(h,'edgecolor','none')
set(gca,'ydir','reverse')
colorbar
set(gca,'ylim',[500 1020]);
set(gca,'ydir','reverse')
colorbar
hold on
plot(x,y,'k')
for j=2:4:length(yday); text(x(1,j),530,num2str(round(octday(j))),'color','k','horizontalalignment','right'); end
title('zonal velocity (m s ^{-1})')
ylabel('pressure (hPa)')
subplot(4,1,4,'align')
[c,h]=contourf(la(i1:i2),pint,vint(:,i1:i2),-10:10);
set(h,'edgecolor','none')
set(gca,'ydir','reverse')
colorbar
set(gca,'ylim',[500 1020])
set(gca,'ydir','reverse')
colorbar
hold on
plot(x,y,'k')
for j=2:4:length(octday); text(x(1,j),530,num2str(round(octday(j))),'color','k','horizontalalignment','right'); end
title('meridional velocity (m s^{-1})')
xlabel('latitude (\circ)')
orient tall
