Global Index (short | long) | Local contents | Local Index (short | long)
sst = sst(360:1800,:,:); oro = oro(360:1800,:,:);
| This script calls | This script is called by |
|---|---|
clear
cd /home/disk/hayes2/dvimont/csm/300yr/srf
load sst1-150.mat
load oro1-150.mat
load lims1-150.mat
whos
lims = [0 360 -90 90];
sst = subset(sst, nxa, nya, lims);
oro = subset(oro, nxa, nya, lims);
nxkeep = nxa(find(nxa >= lims(1) & nxa <= lims(2)));
nykeep = nya(find(nya >= lims(3) & nya <= lims(4)));
[ntmstp, nypts, nxpts] = size(sst);
whos
sst = reshape(sst, ntmstp, nypts*nxpts);
oro = reshape(oro, ntmstp, nypts*nxpts);
oro = mean(oro);
dokeep = find(oro == 0);
sst = sst(:,dokeep);
[sst,clim] = annave(sst);
temp = zeros(12, nypts*nxpts);
temp(1:12,dokeep) = clim;
temp = reshape(temp, 12, nypts, nxpts);
limpac = [115 285 -60 65];
temp = subset(temp, nxa, nya, limpac);
figure(1)
orient landscape
subplot(2,2,1)
gcont(squeeze(temp(1,:,:)),limpac,[273:2:305])
dc
title('SST: January Climatology')
ylabel('Contour Interval: 3 deg C')
subplot(2,2,2)
gcont(squeeze(temp(4,:,:)),limpac,[273:2:305])
dc
title('SST: April Climatology')
subplot(2,2,3)
gcont(squeeze(temp(7,:,:)),limpac,[273:2:305])
dc
xlabel('SST: July Climatology')
ylabel('Contour Interval: 3 deg C')
subplot(2,2,4)
gcont(squeeze(temp(10,:,:)),limpac,[273:2:305])
dc
xlabel('SST: September Climatology')
ylabel('Contour Interval: 3 deg C')