Back to home page

MITgcm

 
 

    


Warning, /verification/adjustment.cs-32x32x1/input/mkmovie.m is written in an unsupported language. File is not indexed.

view on githubraw file Latest commit aea29c85 on 2001-05-29 14:02:00 UTC
aea29c8517 Alis*0001 load GRID
                0002 
                0003 nx=32;olx=5;
                0004 j=0;
                0005 
                0006 cax=[1 1];
                0007 frame=0;
                0008 %for t=48:48:192;
                0009  for t=0:8:288;
                0010 %for t=0:8:64;
                0011  frame=frame+1;
                0012 
                0013  sf=sprintf('%10.10i',t);
                0014  P=read_tiles(['Eta.' sf],nx*[1 1]);
                0015 %T=read_tiles(['T.' sf],nx*[1 1]);
                0016 %S=read_tiles(['S.' sf],nx*[1 1]);
                0017 %W=read_tiles(['W.' sf],nx*[1 1]);
                0018 %GT=read_tiles(['GT.' sf],nx*[1 1]);
                0019 
                0020  Q=P;
                0021 
                0022 %displaytiles(Q);
                0023 
                0024 %plotcube(XG,YG,ZG,Q)
                0025 %axis square
                0026 %axis off
                0027 %view(130+0*t/256*360,20+0*sin(t/256*pi));
                0028 %axis([-1 1 -1 1 -1 1]*sqrt(2/3))
                0029 %cax=caxis;
                0030 %caxis([-1 1]*max(abs(max(cax)),abs(min(cax))))
                0031 
                0032  merccube(lonG,latG,Q);colorbar
                0033  cax=caxis;
                0034  caxis([-1 1]*max(abs(max(cax)),abs(min(cax))));colorbar
                0035 
                0036  title(sprintf('Time = %4.1f hours',t*450/3600))
                0037  drawnow
                0038 
                0039 %print('-djpeg60','-r60',sprintf('splash-%3.3i.jpg',frame));
                0040 end