% To run this script, type 
% 	PlotEquiv
% in the Matlab command window

clear

load model.mat

figure
colormap(hsv)   %default is 'jet'
for ti = 1:12,
  clf
  pcolor(x/1e3, y/1e3, squeeze(theq(ti,:,:)));
  shading interp;
  caxis([335 355])
  colorbar 
  hold on
  contour(x/1e3, y/1e3, squeeze(clw_column(ti,:,:)), [1 1], 'w-');
  axis([200 700 50 250]);
  pause
end


