Documentation of tcline_regression

Global Index (short | long) | Local contents | Local Index (short | long)

Cross-Reference Information

Listing of script tcline_regression

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load HP10_detrend_L1-7_EOF_yr101-1000.mat
lims = [100 300 -3 3];
tim = 101:1000;
temp = getnc('temp', lims, 1:10, tim);
w = getnc('wl', lims, 1:9, tim);
temp = squeeze(shiftdim(mean2(shiftdim(temp, 2)), 2));
w = squeeze(shiftdim(mean2(shiftdim(w, 2)), 2));
nfrm = -2:2;
treg = regress_eof(temp, pcs, nfrm);
wreg = regress_eof(w, pcs, nfrm);
[latt, lont, deptht] = getll('temp', lims);
[latw, lonw, depthw] = getll('wl', lims);
cintt = 0.025; clevt = sort([0:-cintt:min(min(min(treg))) ...
                            cintt:cintt:max(max(max(treg)))]);
cintw = 0.05; clevw = sort([0:-cintw:min(min(min(wreg))) ...
                            cintw:cintw:max(max(max(wreg)))]);
xt = 120:30:270;
xtlab = ['120E'; '150E'; '180 '; '150W'; '120W'; '90W '];
for i = 1:5;
  subplot(7,2,2*i-1);
    pncont(lont, -1*deptht, squeeze(treg(i,:,:)), clevt, 'k');
    axis([110 280 -500 0]);
    dcv(lont, -1*deptht, squeeze(treg(i,:,:)));
    set(gca, 'XTick', xt, 'XTickLabel', xtlab);
    set(gca, 'YTick', [-500:100:0], 'YTickLabel', [500:-100:0]);
    ylabel(['Lag = ' num2str(nfrm(i))]);
  subplot(7,2,2*i);
    pncont(lonw, -1*depthw, squeeze(wreg(i,:,:)), clevw, 'k');
    axis([110 280 -500 0]);
    dcv(lont, -1*deptht, squeeze(treg(i,:,:)));
    set(gca, 'XTick', xt, 'XTickLabel', xtlab);
    set(gca, 'YTick', [-500:100:0], 'YTickLabel', [500:-100:0]);
end
subplot(7,2,1);
  title('< HPPC1, Temperature >');
subplot(7,2,2);
  title('< HPPC1, Vertical Velocity >');
cd ~/Thesis/Chap3