Documentation of LP_HC_evolution_30s30n

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

Help text

  Look at movie:

Cross-Reference Information

Listing of script LP_HC_evolution_30s30n

clear
loadpcs;
lims = [100 300 -35 35];
tim = 101:1000;
lags = -8:2:6;
hc1 = getheat(lims, 1:3, tim);
hc2 = getheat(lims, 4:7, tim);
[b, a] = butter(9, 2/9);
hc1 = filtfilt(b, a, detrend(hc1));
hc2 = filtfilt(b, a, detrend(hc2));
[reg1, c1] = regress_eof(hc1, lpcs, lags);
[reg2, c2] = regress_eof(hc2, lpcs, 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:5; % length(lags);
  sptalk(6,2,2*i-1);
    gcont(1e-6*reg1(i,:,:), 10);
    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-0);
    gcont(1e-6*reg2(i,:,:), 20);
    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}< LPPC1, HC0-80 >');
sptalk(6,2,2);
  title('\fontsize{10}< LPPC1, HC80-270 >');
sptalk(6,2,11);
  xlabel('\fontsize{10}Contour Interval:  10 \times 10^6 J m^-^2 std^-^1');
sptalk(6,2,12);
  xlabel('\fontsize{10}Contour Interval:  20 \times 10^6 J m^-^2 std^-^1');
cd ~/Thesis/Talk
print -dps2 LP9_EQ_30s30n_hc_reg.ps

clear

loadpcs;

lims = [100 300 -65 65];
tim = 101:1000;
lags = -12:12;

hc1 = getheat(lims, 1:3, tim);
hc2 = getheat(lims, 4:7, tim);

[b, a] = butter(9, 2/9);
hc1 = filtfilt(b, a, detrend(hc1));
hc2 = filtfilt(b, a, detrend(hc2));

[reg1, c1] = regress_eof(hc1, lpcs, lags);
[reg2, c2] = regress_eof(hc2, lpcs, lags);

[lat, lon, depth, lm] = getll('temp', lims);

default_global; FRAME = [105 290 -60 60];

figure(2); fl(1); clf;

M = moviein(length(lags));

clev = [-100:10:100];

for i = 1:length(lags);
  subplot(1,1,1);
    mcolor(1e-6*reg1(i,:,:), clev, 'giso');
    drawmapm
    caxis([-80 80]);
    M(:,i) = getframe;
end

M2 = moviein(length(lags));

default_global; FRAME = [105 290 -60 60];

clev = [-220:20:220];
figure(2); fl(1); clf;
for i = 1:length(lags);
  sptalk(1,1,1);
    mcolor(1e-6*reg2(i,:,:), clev, 'giso');
    drawmapm
    caxis([-180 180]);
    M2(:,i) = getframe;
end

findex = [1 1 1 1:length(lags) length(lags)*[1 1 1]];
findex = [findex findex findex findex];

gotop = ['movie(M, findex, 4);'];
gobot = ['movie(M2, findex, 4);'];

eval(gotop);
eval(gobot);

cd ~/Thesis/Talk
%save movieheat.mat M M2 findex gotop gobot


%  Just play the movie:

cd ~/Thesis/Talk
load movieheat.mat

eval(gotop);
eval(gobot);