Documentation of movie_hc_25yrlp

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

Cross-Reference Information

Listing of script movie_hc_25yrlp

clean
ct = getct;
ct = detrend(ct);
[b1, a1] = butter(9, 2/18);
[b2, a2] = butter(9, 2/40);
ct2 = filtfilt(b1, a1, ct)-filtfilt(b2, a2, ct);
lims = [90 305 -80 80];
tim = 101:1000;
varn = 'temp';
[hct, lat1, lon1] = getheat(lims, 1:3, 176:425, 'temp');
hcb = getheat(lims, 4:7, 176:425, 'temp');
hct = detrend(hct);
hcb = detrend(hcb);
hct = filtfilt(b1, a1, hct)-filtfilt(b2, a2, hct);
hcb = filtfilt(b1, a1, hcb)-filtfilt(b2, a2, hcb);
figure(2); fl(1); clf;
M = moviein(50);
clev = [-1:.1:1];
for i = 1:50;
  contourf(lon1, lat1, squeeze(1e-8*hct(2*i,:,:)), clev);
  shading flat
  caxis([-1 1]);
  M(:,i) = getframe;
end
M2 = moviein(50);
clev = [-2:.2:2];
for i = 1:50;
  contourf(lon1, lat1, squeeze(1e-8*hcb(2*i,:,:)), clev);
  shading flat
  caxis([-2 2]);
  M2(:,i) = getframe;
end