Global Index (short | long) | Local contents | Local Index (short | long)
Look at heat content regressions
| This script calls | |
|---|---|
clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP9_hctend_ALL.mat
load LP9_hctend_89.mat
[latu, lonu] = getll('u', lims);
[latw, lonw] = getll('wl', lims);
nfrm = length(lags);
clear tem47 tem89
for i = 1:nfrm;
tem47(i,:,:) = interp2(lonw, latw, squeeze(wbtp47(i,:,:)), lonu, latu');
tem89(i,:,:) = interp2(lonw, latw, squeeze(wbtp89(i,:,:)), lonu, latu');
end
totbar =ubtp89 + vbtp89 + tem89;% ubtp47 + vbtp47 + tem47 + ...
ubtp89 + vbtp89 + tem89;
figure(1); figure_orient;
get_global;
XAX = lonw; YAX = latw; FRAME = [105 300 -60 60];
XAX = lonu; YAX = latu; FRAME = [105 300 -60 60];
tem = totbar;
cint = 0.2;
lagind = [-3:3];
for i = 1:6;
lind = find(lags == lagind(i));
subplot(3,2,i);
gcont(tem(lind,:,:), cint);
dc2(squeeze(tem(lind,:,:)));
ylabel(['Lag = ' num2str(lags(lind))]);
end
clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP9_detrend_L1-7_EOF_yr101-1000.mat
lims = [105 300 -60 60];
hc = getheat(lims, 8:10, 101:1000);
[lat, lon] = getll('temp', lims);
hc = detrend(hc);
ct = getnc('temp', [180 270 -6 6], 1, 101:1000);
ct = mean2(mean2(shiftdim(ct, 1)));
ct = squeeze(detrend(ct));
[b, a] = butter(9, 2/9);
lphc = filtfilt(b, a, hc);
lpct = filtfilt(b, a, ct);
lags = [-30:30];
hcreg = regress_eof(hc, pcs, lags);
lims2 = [180 210 -60 60];
[xk, yk] = keep_var(lims2, lon, lat);
tem = 1e-6*hcreg(:, yk, xk);
tem = squeeze(mean2(shiftdim(tem, 2)))';
timind = 200:2:400;
tem = squeeze(mean2(shiftdim(1e-6*lphc(timind, yk, xk), 2)))';
figure(2); figure_orient(1);
sd(1)
cint = 25; clev = sort([-cint:-cint:min(min(tem)) cint:cint:max(max(tem)) 0]);
contourf(timind, lat, tem, clev); shading flat
axis([min(timind) max(timind) -60 60]);
caxis([-500 500])
hold on;
vline([200:10:400], '--k');
hline(0, '-k')
hold off
sd(2);
plot(timind, lpct(timind));
hold on;
vline([200:10:400], '--k');
hold off
axis([min(timind) max(timind) -0.5 0.5]);