% Plots for Rayleigh-Bernard

clf;
colormap(mk_colmap);

subplot(221)
pcolor(XP,YP,P(NodeP(2:app.Ny+1,2:app.Nx+1))) 
TT(4)=title(sprintf('Pressure (time=%0.3e)',T));
shading interp
colorbar;
axis equal
axis off
SetPlot;

subplot(222)
pcolor(XP,YP,rho(NodeP(2:app.Ny+1,2:app.Nx+1)))
TT(4)=title(sprintf('Rho (Nx=%g, Ny=%g)',app.Nx,app.Ny));
caxis([0 1]);
shading interp
colorbar;
axis equal
axis off
SetPlot;

subplot(223)
pcolor(XV,YV,v(Nodev(2:app.Ny,2:app.Nx+1)))
TT(4)=title(sprintf('V velocity'));
shading interp
colorbar;
axis equal
axis off
SetPlot;

subplot(224)
pcolor(XU,YU,u(Nodeu(2:app.Ny+1,2:app.Nx)))
TT(4)=title(sprintf('U velocity'));
shading interp
colorbar;
axis equal
axis off
SetPlot;

drawnow;
