figure
subplot(2,2,1)
imagesc(squeeze(q(:,:,15))');colorbar;axis square
text(20,20,'q')
text(0,(0-Ny/60-15),descriptor)
text(1.75*Nx,(0-Ny/60-15),'timestep')
text(2.1*Nx,(0-Ny/60-15),it)
subplot(2,2,2)
imagesc(squeeze(qr(:,:,15))');colorbar;axis square
text(20,20,'qr')
subplot(2,2,3)
imagesc(squeeze(qf(:,:,15))');colorbar;axis square
text(20,20,'qf')
subplot(2,2,4)
imagesc(squeeze(qs(:,:,15))');colorbar;axis square
text(20,20,'qs')