Documentation of movie_ceofs

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

Help text

  Try a movie of this

Cross-Reference Information

Listing of script movie_ceofs

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Data
load ceof_L1-10.mat
ctlim = [120 285 -65 65];
get_global
default_global;  FRAME = ctlim;
nkp = length(kp);
num = 1;
clims = [-1:.025:1];
figure(1);
for i = 1:4;
  wgt = conj(exp(j * (i-1) * pi/4) * ones(nkp, 1));
  pat = real(wgt .* lds(:,num));
  subplot(2,2,i);
     tem = NaN * ones(nlat*nlon, 1);
     tem(kp, 1) = pat;
     tem = reshape(tem, nlat, nlon);
     if ismap(gca); clma; end;
     mcont(tem, clims, 'giso');
     title(['CEOF1 Phase = ' num2str((i-1)*(180/4))]);
end

figure(1);  clf;  orient landscape;
num = 2;
clims = [-.26:.0125:.26];
nfrm = 16;

clear M
M = moviein(nfrm);

for i = 1:nfrm;

  wgt = conj(exp(j * (i-1) * 2 * pi/nfrm) * ones(nkp, 1));
  pat = real(wgt .* lds(:,num));

     tem = NaN * ones(nlat*nlon, 1);
     tem(kp, 1) = pat;
     tem = reshape(tem, nlat, nlon);
     tem(40,1) = 0.26; tem(40,2) = -0.26;

%     if ismap(gca); clma; end;
     gshade(tem, clims); dc;
     title(['CEOF1 Phase = ' num2str((i-1)*(360/nfrm))]);

  M(:,i) = getframe;

end

movie(M, 8, 8);

cd /home/disk/tao/dvimont/matlab/CSIRO/Data
save movie_ceof.mat M XAX YAX FRAME

figure(1)
set(gcf,'Position',[2 250 560 722])
load /home/disk/tao/dvimont/matlab/CSIRO/Data/movie_ceof.mat
n_repeat = 4;
frames_per_sec = 8;
movie(M, n_repeat, frames_per_sec)

%  Plot the temporal phase of the time series

cd /home/disk/tao/dvimont/matlab/CSIRO/Data
%load ceof_L1-10.mat
load ceof_L1-10_try3.mat

num = 1;
tphase = atan2(imag(pcs(:,num)),real(pcs(:,num))) * 180/pi;

for m = 1:4
     subplot(4,1,m)
     tind = (ntim/4)*(m-1)+[1:(ntim/4)];
     plot(tind, tphase(tind), '-k')
     axis([min(tind-1) max(tind+1) -185 185]);
     set(gca, 'YTick', [-180:45:180], 'XTick', [0:25:ntim]);
     grid
end

%  Movie patterns through actual phases

tphase = atan2(imag(pcs(:,num)),real(pcs(:,num)));
tamp = conj(pcs(:,num)).*pcs(:,num);

inds = 125:225;
figure(1);  clf;  orient landscape;
num = 1;
clims = [-.5:.025:.5];
nfrm = length(inds);

clear M
M = moviein(length(inds));

for i = inds;

  wgt = conj(exp(j * tphase(i)) * ones(nkp, 1));
  pat = real(tamp(i) .* wgt .* lds(:,num));

     tem = NaN * ones(nlat*nlon, 1);
     tem(kp, 1) = pat;
     tem = reshape(tem, nlat, nlon);
     tem(40,1) = 0.5; tem(40,2) = -0.5;

%     if ismap(gca); clma; end;
     gshade(tem, clims); dc;
     title(['CEOF1 Phase = ' num2str((i-1)*(360/nfrm))]);

  M(:,(i-124)) = getframe;

end

movie(M, 3, 8)