% plot time-pressure series for the whole cruise
clf
b2rcolormap(24);
subplot(4,1,1,'align')
[c,h]=contourf(yday(i1:i2),pint,thetaint(:,i1:i2),286:2:320);
set(h,'edgecolor','none')
set(gca,'ylim',[500 1020])
set(gca,'ydir','reverse')
colorbar
hold on
title({'STRATUS 2007 time-pressure section' 'potential temperature (K)'})
subplot(4,1,2,'align')
[c,h]=contourf(yday(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
title('specific humidity (g/kg)');
subplot(4,1,3,'align')
[c,h]=contourf(yday(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
title('zonal velocity (m s ^{-1})')
ylabel('pressure (hPa)')
subplot(4,1,4,'align')
[c,h]=contourf(yday(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
title('meridional velocity (m s^{-1})')
xlabel('year day (UTC, Oct 18-Nov 3)')
orient tall
%print('-depsc2',[way_raw_images_balloon 'timepres_thetaquv.eps'])
%print('-djpeg',[way_raw_images_balloon 'timepres_thetaquv.jpg'])
