Documentation of HP_HC_evolution_30s30n

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

Cross-Reference Information

Listing of script HP_HC_evolution_30s30n

clear
loadpcs;
lims = [100 300 -35 35];
tim = 101:1000;
lags = -2:2;
hc1 = getheat(lims, 1:3, tim);
hc2 = getheat(lims, 4:7, tim);
[b, a] = butter(9, 2/10);
hc1 = hc1 - filtfilt(b, a, detrend(hc1));
hc2 = hc2 - filtfilt(b, a, detrend(hc2));
[reg1, c1] = regress_eof(hc1, hpcs, lags);
[reg2, c2] = regress_eof(hc2, hpcs, lags);
[lat, lon, depth, lm] = getll('temp', lims);
default_global; FRAME = [105 290 -30 30];
clim = 0.25;
figure(1); fo; clf;
for i = 1:length(lags);
  sptalk(6,2,2*i-1);
    gcont(1e-6*reg1(i,:,:), 15);
    dc2(lm, 0.3, 1);
    color_shade(squeeze(c1(i,:,:).^2), clim, .7*[1 1 1]);
    ylabel(['Lag = ' num2str(lags(i))]);
    set(gca, 'fontsize', 9);
  sptalk(6,2,2*i);
    gcont(1e-6*reg2(i,:,:), 30);
    dc2(lm, 0.3, 1);
    color_shade(squeeze(c2(i,:,:).^2), clim, .7*[1 1 1]);
    set(gca, 'fontsize', 9);
end
sptalk(6,2,1); 
  title('\fontsize{10}< HPPC1, HC0-80 >');
sptalk(6,2,2);
  title('\fontsize{10}< HPPC1, HC80-270 >');
sptalk(6,2,9);
  xlabel('\fontsize{10}Contour Interval:  15 \times 10^6 J m^-^2 std^-^1');
sptalk(6,2,10);
  xlabel('\fontsize{10}Contour Interval:  30 \times 10^6 J m^-^2 std^-^1');
cd ~/Thesis/Talk
print -dps2 HP10_EQ_30s30n_hc_reg.ps