Global Index (short | long) | Local contents | Local Index (short | long)
[temp, lat, lon, depth, middepth] = getheat(lims, lev, tim, varn);
| This script calls | |
|---|---|
clear
varn = 'v';
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_detrend_L1-7_yr101-1000.mat
lev = 1:10;
tim = 101:1000;
lims = [100 180 25 40];
lags = -10:5; nfrm = lags;
temp = getnc(varn, lims, lev, tim);
[lat, lon, depth] = getll(varn, lims);
% Get western boundary
[ntim, nlev, ny, nx] = size(temp);
wb = NaN * ones(ny, 1);
for i = 1:ny;
land = find(isnan(squeeze(temp(1,1,i,:))));
wb(i) = max(land)+1;
end
temp2 = NaN * ones(ntim, nlev, ny, 5);
for i = 1:ny;
temp2(:,:,i,:) = temp(:,:,i,[wb(i):(wb(i)+4)]);
end
temp2 = shiftdim(squeeze(mean2(shiftdim(temp2, 2))), 1);
treg = regress_eof(temp2, pcs, nfrm);
lon = lon(min(wb):(min(wb)+4));
cint = 0.01; clev = [-.5:cint:-cint cint:cint:.5];
figure(3); figure_landscape;
for i = 1:4;
subplot(4,4,4*i-3);
pncont(lon, -1*depth, squeeze(treg(i,:,:)), clev, 0, 'k');
axis([120 150 -550 0]);
ylabel(['Lag = ' num2str(lags(i))]);
set(gca, 'YTick', -500:100:0);
subplot(4,4,4*i-2);
pncont(lon, -1*depth, squeeze(treg(i+4,:,:)), clev, 0, 'k');
axis([120 150 -550 0]);
ylabel(['Lag = ' num2str(lags(i+4))]);
set(gca, 'YTick', -500:100:0);
subplot(4,4,4*i-1);
pncont(lon, -1*depth, squeeze(treg(i+8,:,:)), clev, 0, 'k');
axis([120 150 -550 0]);
ylabel(['Lag = ' num2str(lags(i+8))]);
set(gca, 'YTick', -500:100:0);
subplot(4,4,4*i);
pncont(lon, -1*depth, squeeze(treg(i+12,:,:)), clev, 0, 'k');
axis([120 150 -550 0]);
ylabel(['Lag = ' num2str(lags(i+12))]);
set(gca, 'YTick', -500:100:0);
end
% Look at tx index as defined by the taux anomalies in the NPac
lims1 = [160 180 30 37.5];
lims2 = [170 190 45 52.5];
tx1 = getnc('taux', lims1, 1, tim);
tx2 = getnc('taux', lims2, 1, tim);
tx = squeeze(mean(mean((tx1 - tx2), 2), 3));
[b, a] = butter(6, 2/10);
tx = filtfilt(b, a, detrend(tx));
pcs = tx;