% 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
   contourf(sstjan)
   Frames(:,1)=getframe;

   load sstfeb.dat
   contourf(sstfeb)
   Frames(:,2) = getframe;

   load sstmar.dat
   contourf(sstmar)
   Frames(:,3)=getframe;

   load sstapr.dat
   contourf(sstapr)
   Frames(:,4) = getframe;

   load sstmay.dat
   contourf(sstmay)
   Frames(:,5)=getframe;

   load sstjun.dat
   contourf(sstjun)
   Frames(:,6) = getframe;
 
   load sstjul.dat
   contourf(sstjul)
   Frames(:,7)=getframe;

   load sstaug.dat
   contourf(sstaug)
   Frames(:,8) = getframe;
 
   load sstsep.dat
   contourf(sstsep)
   Frames(:,9)=getframe;

   load sstoct.dat
   contourf(sstoct)
   Frames(:,10)= getframe;
 
   load sstnov.dat
   contourf(sstnov)
   Frames(:,11)=getframe;

   load sstdec.dat
   contourf(sstdec)
   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')