Documentation of hc_sst_lagged_regression

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

Help text

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_detrend_L1-7_yr550-1000.mat; pcs = -1*pcs;

Cross-Reference Information

Listing of script hc_sst_lagged_regression

clear

lims = [100 300 -65 65];
tim = 101:1000;
lev = 4:7;

ct = getnc('temp', [180 270 -6 6], 1, tim);
ct = squeeze(mean(mean(pcs, 2), 3));
ct = detrend(pcs);
[b1, a1] = butter(9, 2/15);
[b2, a2] = butter(9, 2/30);
pcs = filtfilt(b1, a1, ct) - filtfilt(b2, a2, ct);

heat = getheat(lims, lev, tim);
temp = getnc('temp', lims, 1, tim);
[lat, lon] = getll('temp', lims);

lags = [-12:0];
hreg = regress_eof(heat, pcs, lags);
treg = regress_eof(temp, pcs, lags);

default_global; FRAME = [105 290 -60 60];
figure(1); figure_orient;

ind = 2;
for i = 1:4;
  subplot(4,2,2*i-1);
  gcont(1e-8*hreg(i+ind,:,:), [-5:.2:-.2 .2:.2:5]);
  dc2(hreg(i+ind,:,:))
  ylabel(['Lag = ' num2str(lags(i+ind))]);
end
ind = 6;
for i = 1:4;
  subplot(4,2,2*i-0);
  gcont(1e-8*hreg(i+ind,:,:), [-5:.2:-.2 .2:.2:5]);
  dc2(hreg(i+ind,:,:))
  ylabel(['Lag = ' num2str(lags(i+ind))]);
end


subplot(4,2,1);
  title('LP10:  80-270m Heat Content');
subplot(4,2,7)
  xlabel('Contour Interval:  1x10^7 J m^-^2 std^-^1');

for i = 1:4;
  subplot(4,2,2*i);
  gcont(treg(i,:,:), [-1:.025:-.025 .025:.025:1]);
  dc2(treg(i,:,:))
  ylabel(['Lag = ' num2str(lags(i))]);
end
subplot(4,2,2);
  title('LP10:  12.5m Temperature');
subplot(4,2,8)
  xlabel('Contour Interval:  0.025 K std^-^1');

cd /home/disk/tao/dvimont/matlab/CSIRO/IUGG_talk/Plots


clear

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_detrend_L1-7_yr101-1000.mat; pcs = -1*pcs;

lims = [100 300 0 65];
tim = 101:1000;
lev = 1:7;

heat = getheat(lims, lev, tim);
[lath, lonh] = getll('temp', lims);
v = getheat(lims, lev, tim, 'v');
[latv, lonv] = getll('v', lims);

[strf, lat_out, lon_out] = vel_to_strfnc(squeeze(mean(v)), latv, lonv);

lags = [-13:-6];
hreg = regress_eof(heat, pcs, lags);

figure(1); figure_orient;

default_global; FRAME = [105 290 2.5 62.5];

rlim = 0.4;
for i = 1:4;
subplot(4,2,2*i-1);
  XAX = lon_out; YAX = lat_out;
  gcont(1e-6*strf, [-20:20]);
  hh = dc2(strf);
  hold on;
  XAX = lonh; YAX = lath;
  color_shade(squeeze(1e-8*hreg(i,:,:)), rlim, 'r');
  color_shade(squeeze(-1e-8*hreg(i,:,:)), rlim, 'b');
  hold off
  ylabel(['Lag = ' num2str(lags(i))]);
end


for i = 1:4;
subplot(4,2,2*i);
  XAX = lon_out; YAX = lat_out;
  gcont(1e-6*strf, [-20:20]);
  hh = dc2(strf);
  hold on;
  XAX = lonh; YAX = lath;
  color_shade(squeeze(1e-8*hreg(i+4,:,:)), rlim, 'r');
  color_shade(squeeze(-1e-8*hreg(i+4,:,:)), rlim, 'b');
  hold off
  ylabel(['Lag = ' num2str(lags(i+4))]);
end


h = patch([150 200 200 150], [30 30 40 40], [-100 -100 -100 -100], 'r');