% Plots for Lock Exchange

clf;

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

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

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

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

colormap(mk_colmap);

drawnow;