% Last modified: 06/09/2012 by Pete

% interpolate velocities to pressure grid coordinates, plot velocity
% magnitudes, 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)) ) ;
pcolor( XP, YP, ( u_intrp.^2 + v_intrp.^2 ).^(1/2) )
shading interp, colormap(hot), caxis([0,2])
h = streamslice( XP, YP, u_intrp, v_intrp, 'noarrows' ) ;
set( h, 'color', 'k' ) ;

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

