adjustment.128x64x1/input/mkmovie.m

aea29c8517 Alis*0001 nx=128;ny=64;
                0002
                0003 x=[0.5:nx]/nx*360;
                0004 y=[0.5:ny]/ny*180-90;
                0005
                0006 frame=0;
                0007  for t=0:8:288;
                0008  frame=frame+1;
                0009
                0010  sf=sprintf('%10.10i',t);
                0011  P=rdmds(['Eta.' sf]);
                0012
                0013  Q=P;
                0014
                0015  plot(y,Q');
                0016
                0017 %pcol(x,y,Q');colorbar
                0018 %cax=caxis;
                0019 %caxis([-1 1]*max(abs(max(cax)),abs(min(cax))));colorbar
                0020
                0021  title(sprintf('Time = %4.1f hours',t*450/3600))
                0022  drawnow
                0023
                0024 %print('-djpeg60','-r60',sprintf('splash-%3.3i.jpg',frame));
                0025 end