Plotting model output
Download this model output to your home directory:
   model.mat

Information on model variables:
theq(time,y,x): equivalent potential temperature at the surface (K)
clw_column(time,y,x): column-integrated cloud water (mm)
x, y: model dimensions (m)

clear

load model.mat

figure
for ti = 1:12,		% loop through time
  clf			% clear the current figure
  pcolor(x/1e3, y/1e3, squeeze(theq(ti,:,:)));	
			% plot Theta_e
  shading interp;	% interpolate the color map
  colorbar 
  hold on		% this holds the current plot, so next line with overlay instead of overwrite
  contour(x/1e3, y/1e3, squeeze(clw_column(ti,:,:)), [1 1], 'w-');
			% draw the 1 mm contour of clw_column 
  pause			% wait for the user to continue
end

caxis([335 355])	% change the color axis
colorbar

axis([200 700 50 250]);	% change the domain axis