Documentation of hc_sst_lagged_regression_color

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

Cross-Reference Information

Listing of script hc_sst_lagged_regression_color

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_detrend_L1-7_yr550-1000.mat; pcs = -1*pcs;
lims = [100 300 -65 65];
tim = 551:1000;
lev = 4:7;
heat = getheat(lims, lev, tim);
temp = getnc('temp', lims, 1, tim);
[lat, lon] = getll('temp', lims);
lags = [-12 -10 -8 -7];
hreg = regress_eof(heat, pcs, lags);
treg = regress_eof(temp, pcs, lags);
default_global; FRAME = [105 290 -60 60];
figure(1); figure_orient;
dx = mean(diff(XAX));
dy = mean(diff(YAX));
XAX2 = [XAX(1)-dx; XAX; XAX(nx)+dx];
YAX2 = [YAX(1)-dy; YAX; YAX(ny)+dy];
for i = 1:4;
  subplot(4,2,2*i);
    hreg3 = squeeze(hreg(i,:,:));
    [ny, nx] = size(hreg3);
    hreg3 = [ones(1, nx+2); [ones(ny, 1) hreg3 ones(ny, 1)]; ones(1, nx+2)];
    contourf2(XAX2, YAX2, 1e-7*hreg3, [-17:12]);
    hold on;
      pncont(XAX, YAX, squeeze(1e-7*hreg(i,:,:)), [-18:2:12], 0, 'k');
    hold off;
    ylabel(['Lag = ' num2str(lags(i))]);
    caxis([-17 12])
    colorbar2
    dc3(0);
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');