Documentation of lagged_regressions

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

Cross-Reference Information

Listing of script lagged_regressions

clear
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
load LP10_detrend_L1-7_yr101-1000.mat
lims = [100 300 -45 65];
tim = 101:1000;
lags = [-10:0];
[temp, taux] = getnc('temp', 'taux', lims, 1, tim);
heat = getheat(lims, 4:7, tim);
[lat, lon, depth] = getll('temp', lims);
[latx, lonx] = getll('taux', lims);
[treg, tcoef] = regress_eof(temp, pcs, lags);
[hreg, hcoef] = regress_eof(heat, pcs, lags);
[xreg, xcoef] = regress_eof(taux, pcs, lags);
figure(1); figure_landscape;
default_global; FRAME = [110 300 -40 60];
cint1 = 0.15;  clev1 = [-1.5:cint1:-cint1 cint1:cint1:1.5];
cint2 = 0.15;  clev2 = [-1.5:cint2:-cint2 cint2:cint2:1.5];
cint3 = 0.15;  clev3 = [-1.5:cint3:-cint3 cint3:cint3:1.5];
global GRDX_SPACING; GRDX_SPACING = 45;
for i = 1:3;
  subplot(3,3,3*i-2);
    tem = zeros(size(hcoef(4*i-1,:,:)));
    tem(find(hcoef(i+3,:,:).^2 > 0.1)) = NaN;
    eval(['clev = clev' num2str(i) ';']);
    gcont(-1e-8*hreg(i+3,:,:), clev);
    dc2(hreg(i+3,:,:));
    shade_ccoef(tem);
    ylabel(['Lag = ' num2str(lags(i+3)) ' Years']);
end
xlabel(['Contour:  ' num2str(100*cint1) 'x10^6 J m^-^2 std^-^1'])
cint1 = 0.2;  clev1 = [-10:cint1:-cint1 cint1:cint1:10];
cint2 = 0.2;  clev2 = [-10:cint2:-cint2 cint2:cint2:10];
cint3 = 0.2;  clev3 = [-10:cint3:-cint3 cint3:cint3:10];
for i = 1:3;
  subplot(3,3,3*i-1);
    tem = zeros(size(tcoef(i+3-1,:,:)));
    tem(find(tcoef(i+3,:,:).^2 > 0.1)) = NaN;
    eval(['clev = clev' num2str(i) ';']);
    gcont(-1e1*treg(i+3,:,:), clev);
    dc2(treg(i+3,:,:));
    shade_ccoef(tem);
end
xlabel(['Contour:  ' num2str(0.1*cint1) ' K std^-^1'])
cint1 = 5;  clev1 = [-100:cint1:-cint1 cint1:cint1:100];
cint2 = 5;  clev2 = [-100:cint2:-cint2 cint2:cint2:100];
cint3 = 5;  clev3 = [-100:cint3:-cint3 cint3:cint3:100];
default_global; XAX = lonx; YAX = latx; FRAME = [110 300 -40 60];
for i = 1:3;
  subplot(3,3,3*i);
    default_global; XAX = lonx; YAX = latx; FRAME = [110 300 -40 60];
    tem = zeros(size(xcoef(i+3,:,:)));
    tem(find(xcoef(i+3,:,:).^2 > 0.1)) = NaN;
    eval(['clev = clev' num2str(i) ';']);
    gcont(-1e4*(xreg(i+3,:,:)/10), clev);
    default_global; FRAME = [110 300 -40 60];
    dc2(treg(i+3,:,:));
    default_global; XAX = lonx; YAX = latx; FRAME = [110 300 -40 60];
    shade_ccoef(tem);
end
xlabel(['Contour:  ' num2str(cint1) 'x10^-^4 N m^-^2 std^-^1'])
subplot(3,3,1);
  title('80-270m Heat Content');
subplot(3,3,2);
  title('12.5m Temperature');
subplot(3,3,3);
  title('Zonal Wind Stress');
cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/GRL_Figs
print -dps2 figure3.ps