Documentation of lagged_1000yr_regression

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

Cross-Reference Information

Listing of script lagged_1000yr_regression

clear
lims = [106 312 -63 63];
lev = 1:7;
tim = 101:1000;
lags = [-8:6];
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
load LP10_detrend_L1-7_yr101-1000.mat
temp = getnc('temp', lims, 1, tim);
[latt, lont, deptht] = getll('temp', lims);
[heat, lat, lon, depth, middepth] = getheat(lev, tim, lims);
[treg, ccoef] = regress_eof(temp, pcs, lags);
[hreg, hcoef] = regress_eof(heat, pcs, lags);
cint = 0.05; clev = [-5:cint:-cint cint:cint:5];
default_global;  XAX = lont;  YAX = latt;  FRAME = [110 300 -60 60];
figure(3); figure_orient; offset = 9;
cint = 0.05; clev = [-5:cint:-cint cint:cint:5];
for i = 1:4;
  rind = (i-1)+offset;
  subplot(4,2,2*i-1);
    gcont(-1*treg(rind,:,:), clev);
    dc2(treg(rind,:,:));
    ylabel(['Lag ' num2str(lags(rind))]);
    if i == 4;
      xlabel(['Contour Interval:  ' num2str(cint) ' K std^-^1']);
    end
end
cint = 0.25; clev = [-5:cint:-cint cint:cint:5];
for i = 1:4;
  rind = (i-1)+offset;
  subplot(4,2,2*i);
    gcont(-1e-8*hreg(rind,:,:), clev);
    dc2(treg(rind,:,:));
    ylabel(['Lag ' num2str(lags(rind))]);
    if i == 4;
      xlabel(['Contour Interval:  ' num2str(cint) ' x 10^8 J m^-^2 std^-^1']);
    end
end