%Plot script for Sudden Expansion Test Case
clf
colormap(mk_colmap)

%Quiver Plot
tmp1 = 0.5*(u(Nodeu(2:app.Ny+1,1:app.Nx))+u(Nodeu(2:app.Ny+1,2:app.Nx+1)));
tmp2 = 0.5*(v(Nodev(2:app.Ny+1,2:app.Nx+1))+v(Nodev(1:app.Ny,2:app.Nx+1)));
pcolor(XP, YP, sqrt( tmp1.^2 + tmp2.^2 ))
shading interp; 
hold on
streamline( XP, YP, tmp1, tmp2, [a/app.Nx a/app.Nx a/app.Nx a/4 a/4],[b/2 5*b/8 3*b/8 b/4 3*b/4])
% surf(XP,YP,sqrt(tmp1.^2+tmp2.^2),'FaceAlpha',0.2,'EdgeAlpha',0.0), colorbar

hold off
xlim([0 a]), ylim([0 b])
TT(1)=xlabel('x(m)');
TT(2)=ylabel('y(m)');
TT(3)=zlabel('T(^oC)');
TT(4)=title(sprintf('\nVelocity, time =%03.03g',T));
formatplot([],gca,TT);view(2)

drawnow;
GenericSaveScript;
