% Last modified: 06/09/2012 by Pete

% plot rho
NodeRho = NodeRho ;
pcolor( XRho, YRho, rho_tmp(NodeRho(2:end-1,2:end-1)) )
colormap(pmkmp(256,'CubicL'))
shading interp, caxis([0.0,1.0]), freezeColors

% interpolate velocities to pressure grid coordinates and plot
% streamlines
u_intrp = (1/2) * ( u_tmp(Nodeu(2:Ny+1,1:Nx)) ...
    + u_tmp(Nodeu(2:Ny+1,2:Nx+1)) ) ;
v_intrp = (1/2) * ( v_tmp(Nodev(2:Ny+1,2:Nx+1)) ...
    + v_tmp(Nodev(1:Ny,2:Nx+1)) ) ;
h = streamslice( XP, YP, u_intrp, v_intrp, 'noarrows' ) ;
set(h,'color','k') ;

% set subplot properties
title( titlestr, 'FontWeight', 'bold' )
axis image, box on, xlim([0,a]), ylim([0,b])
set( gca, 'XTick', (0:3:a), 'YTick', (0:3:b) )

