Documentation of movies


Global Index (all files) (short | long) | Local Index (files in subdir) (short | long)


Help text

  Total movie

Cross-Reference Information

This script calls

Listing of script movies


clear
filin = 'nocor2.cdf'
filin = 'fplane3.cdf';
filin = 'mlbeta.cdf';
filin = 'eqbeta.cdf';
filin = 'sal.cdf';
cd /home/disk/hayes2/dvimont/ocean/mid_lat/dam
nc = netcdf(filin, 'nowrite');
  h = nc{'h'}(:,:,:,:);
nc = close(nc);
whos
figure(2); fl;
n = 50;
M = moviein(n);
lims = 1:72;
for i = 1:n
  mesh(squeeze(h(4*i,1,lims,lims)-15000));
  axis([1 72 1 72 -1 1]);
  view([-37.5, 70])
  caxis([-1 1])
  M(i) = getframe;
end
cd /home/disk/hayes2/dvimont/ocean/mid_lat/dam
nc = netcdf(filin, 'nowrite');
  u = nc{'u'}(:,:,:,:);
  v = nc{'v'}(:,:,:,:);
nc = close(nc);
figure(2); fl(1);
n = 50;
M = moviein(n);
for i = 1:n
  quiver(1:3:72, 1:3:72, squeeze(u(i,1,1:3:72,1:3:72)), ...
         squeeze(v(i,1,1:3:72,1:3:72)), -0.5, 'r');
     axis([1 72 1 72]); 
  M(i) = getframe(gca);
end
cd /home/disk/hayes2/dvimont/ocean/mid_lat/dam
nc = netcdf(filin, 'nowrite');
  q = nc{'q'}(:,:,:,:);
nc = close(nc);
qm = mean(mean(mean(mean(q))))
figure(1); fl;
n = 20;
M = moviein(n);
for i = 1:n
  mesh(squeeze(1e12*(q(i,1,:,:)-qm)));
  axis([1 72 1 72 -1 1]);
  caxis([-0.4 0.4])
  M(i) = getframe;
end

clear
filin = 'nocor2.cdf'
filin = 'fplane3.cdf';
filin = 'mlbeta.cdf';
filin = 'eqbeta.cdf';
filin = 'sal.cdf';


cd /home/disk/hayes2/dvimont/ocean/mid_lat/dam
nc = netcdf(filin, 'nowrite');
  h = nc{'h'}(:,:,:,:);
  u = nc{'u'}(:,:,:,:);
  v = nc{'v'}(:,:,:,:);
nc = close(nc);

figure(1); fo(1);
subplot(2,1,1); set(gca, 'units', 'pixel')
     pos1 = get(gca, 'Position');
subplot(2,1,2); set(gca, 'units', 'pixel')
     pos2 = get(gca, 'Position');
     rect = [pos1(1) pos2(2) pos1(2)+pos1(4)-pos2(2) pos1(4)];


n = 50;

ind = 1;
etot = 0.5*sum(sum(sum(sum((h(ind,:,:,:)-15000).^2)))) + ...
       0.5*sum(sum(sum(sum(u(ind,:,:,:).^2+v(ind,:,:,:).^2))));

lims = 1:72;
M = moviein(n);
clear pe ke
cd ~/matlab/Models/figs2

n = 50;
for i = 1:n
  subplot(2,1,1);
%  subplot('position', [0.1300 0.63 0.7750 0.33])
    mesh(squeeze(h(i,1,lims,lims)-15000));
    axis([1 72 1 72 -1 1]);
    view([-37.5, 70])
    caxis([-1 1])
  subplot(2,1,2); %set(gca, 'Position', [0.13 0.11 0.5270 0.3439])
%  subplot('position', [0.1300 0.28 0.7750 0.33])
    quiver(1:3:72, 1:3:72, squeeze(u(i,1,1:3:72,1:3:72)), ...
           squeeze(v(i,1,1:3:72,1:3:72)), 1.4, 'r');
    hold on;
    contour(squeeze(h(i,1,lims,lims)-15000), [-1:.2:-.2 .2:.2:1]);
    hold off;
    caxis([-1 1])
    axis([1 72 1 72]); %axis square
    eval(['print -dpsc fig' num2str(i) '.ps;']);
%    pe(i) = 9.8*0.5*sum(sum(sum(sum((h(i,:,:,:)-15000).^2))));
%    ke(i) = 15000*0.5*sum(sum(sum(sum(u(i,:,:,:).^2+v(i,:,:,:).^2))));
%  subplot('position', [0.1300 .04 .7750 .20])
%    plot(1:i, pe, '-b', 1:i, ke, '--b');
%    axis([1 n 0 25000]);
%  M(i) = getframe(gcf);
end

movie(gcf, M, 4, 8)


cd ~/matlab/Models
%save movies.mat fplane_movie no_cor_movie 

load movies.mat
movie(gcf, fplane_movie, 4, 8);
movie(gcf, no_cor_movie, 4, 8);