% Load in monthly coads climatology and display as a movie
% written Feb 96 for OC3030


reruns=1;                  % number of times movie is to play
fps=5;                     % frames per second

nframes = 12;              % number of frames in the movie
Frames = moviein(nframes); % initialize the matrix 'Frames'

title('SST Climatology')
colormap('jet')
map=jet;
% Load each month, contour it and save it in a frame:
   load sstjan.dat
   pcolor(sstjan); shading interp
   Frames(:,1)=getframe;

   load sstfeb.dat
  pcolor(sstfeb); shading interp
   Frames(:,2) = getframe;

   load sstmar.dat
   pcolor(sstmar); shading interp
   Frames(:,3)=getframe;

   load sstapr.dat
  pcolor (sstapr); shading interp
   Frames(:,4) = getframe;

   load sstmay.dat
   pcolor(sstmay); shading interp
   Frames(:,5)=getframe;

   load sstjun.dat
   pcolor(sstjun); shading interp
   Frames(:,6) = getframe;
 
   load sstjul.dat
   pcolor(sstjul); shading interp
   Frames(:,7)=getframe;

   load sstaug.dat
   pcolor(sstaug); shading interp
   Frames(:,8) = getframe;
 
   load sstsep.dat
   pcolor(sstsep); shading interp
   Frames(:,9)=getframe;

   load sstoct.dat
   pcolor(sstoct); shading interp
   Frames(:,10)= getframe;
 
   load sstnov.dat
   pcolor(sstnov); shading interp
   Frames(:,11)=getframe;

   load sstdec.dat
   pcolor(sstdec); shading interp
   Frames(:,12) = getframe;      


% Now play the movie: 
   movie(Frames,reruns,fps)

% Now save the movie as an mpeg file for use on the Web:
map=colormap         
%  mpgwrite(Frames,map,'sstmovie.mpg')