Documentation of hc_eofs

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

Cross-Reference Information

Listing of script hc_eofs

clear
loadpcs
lims = [-0.1 360 -90 90];
tim = 101:1000;
hc1 = getheat(lims, 1:7, tim);
hc1 = detrend(hc1);
[b1, a1] = butter(9, 2/9);
[b2, a2] = butter(9, 2/10);
[reg1, c1] = regress_eof(hc1-filtfilt(b2, a2, hc1), hpcs, 0);
[reg2, c2] = regress_eof(filtfilt(b1, a1, hc1), lpcs, 0);
[lat, lon, depth, lm] = getll('temp', lims);
default_global; FRAME = [0 360 -75 75];
clim = .25
figure(1); fo(1); clf;
sptalk(3,1,1);
  gcont(1e-8*reg1, 0.35);
  dc2(squeeze(reg1), 0.3, 1);
  color_shade(squeeze(c1.^2), clim, .7*[1 1 1]);
  title('< HPPC1, HC270 >')
sptalk(3,1,2);
  gcont(1e-8*reg2, 0.35);
  dc2(squeeze(reg2), 0.3, 1);
  color_shade(squeeze(c2.^2), clim, .7*[1 1 1]);
  title('< LPPC1, HC270 >')
  xlabel('Contour Interval:  35 \times 10^6 J m^-^2 std^-^1');
cd ~/Thesis/Talk
print -dps2 hc270_hp_lp_regression.ps