
clf;

subplot(221)
q=find(lat>19.85 & isfinite(temp(5,:))');
contourf(lon(q),p,temp(:,q),5:.5:30)
set(gca,'ydir','reverse')
axis([min(lon(q)) max(lon(q)) 0 500])
shading flat
xlabel('longitude')
ylabel('depth (m)')

subplot(223)
contourf(lon(q),p,sal(:,q),35:.1:37.5);
set(gca,'ydir','reverse')
axis([min(lon(q)) max(lon(q)) 0 500])
shading flat
xlabel('longitude')
ylabel('depth (m)')

q=find(lon>-23.5 & isfinite(temp(5,:))');
subplot(222)
contourf(lat(q),p,temp(:,q),5:.5:30)
set(gca,'ydir','reverse')
set(gca,'xdir','reverse')
axis([min(lat(q)) max(lat(q)) 0 500])
shading flat
xlabel('latitude (23W)')
ylabel('depth (m)')
colorbarvert;

subplot(224)
contourf(lat(q),p,sal(:,q),35:.1:37.5);
set(gca,'ydir','reverse')
set(gca,'xdir','reverse')
axis([min(lat(q)) max(lat(q)) 0 500])
shading flat
xlabel('latitude (23W)')
ylabel('depth (m)')
colorbarvert;
