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