% Plots for Lock Exchange
clf
subplot(221)
surf(XP,YP,P(NodeP(2:app.Ny+1,2:app.Nx+1)),'LineStyle','none') 
xlim([min(xp) max(xp)]);ylim([min(yp) max(yp)]);
TT(4)=title(sprintf('Pressure time=%2.4g',T));%caxis([-1 1]*3);
view(2),shading interp
axis equal
axis off

subplot(222)
pcolor(XP,YP,rho(NodeP(2:app.Ny+1,2:app.Nx+1))) ; colorbar;
xlim([min(xp) max(xp)]);ylim([min(yp) max(yp)]);
TT(4)=title(sprintf('Rho (app.Nx=%g, app.Ny=%g)',app.Nx,app.Ny));%caxis([-1 1]*5);
view(2),shading interp
axis equal
axis off

subplot(223)
pcolor(xv,yv,v(Nodev(2:app.Ny,2:app.Nx+1))) ; colorbar
xlim([min(xv) max(xv)]);ylim([min(yv) max(yv)]);
TT(4)=title(sprintf('V velocity',app.Nx,app.Ny));%caxis([-1 1]*2);
view(2),shading interp
axis equal
axis off

subplot(224)
pcolor(xu,yu,u(Nodeu(2:app.Ny+1,2:app.Nx))) 
xlim([min(xu) max(xu)]);ylim([min(yu) max(yu)]);
TT(4)=title(sprintf('U velocity',app.Nx,app.Ny));%caxis([-1 1]*2);
view(2),shading interp
axis equal
axis off
colormap(mk_colmap)
GenericSaveScript;