figure
subplot(2,2,1)
imagesc(meanq(:,:,2));colorbar;axis square
text(0,(0-Ny/60-15),descriptor)
text(1.75*Nx,(0-Ny/60-15),'timestep')
text(2.1*Nx,(0-Ny/60-15),int2str(eval(itstart)))
text(2.4*Nx,(0-Ny/60-15),'to')
text(2.5*Nx,(0-Ny/60-15),int2str(eval(itend)))
%title('mean q')
subplot(2,2,2)
imagesc(meanqr(:,:,2));colorbar;axis square
%title('mean qr')
subplot(2,2,3)
imagesc(meanqf(:,:,2));colorbar;axis square
%title('mean qf')
subplot(2,2,4)
imagesc(meanqs(:,:,2));colorbar;axis square
%title('mean qs')