Global Index (short | long) | Local contents | Local Index (short | long)
Total movie
This script calls | |
---|---|
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);