set(0,'DefaultLineLineWidth',1.5)
set(0,'DefaultLineMarkerSize',15)
set(0,'DefaultAxesFontSize',20)
set(0,'DefaultFigureColor',[1,1,1])
set(0,'DefaultTextFontSize',18)
set(0,'DefaultTextFontName','Times-Roman')
set(0,'DefaultAxesFontName','Times-Roman')
load(['final_data_run1.mat'])
[Xq, Yq] = meshgrid(Xi,Yi);
fig = figure('Renderer', 'painters', 'Position', [1300 10 947 900])
quiver(Yq,-Xq,Vq,-Uq, 1);
yline(0);xline(0); hold on
l = title({'test_line'}, 'Interpreter', 'none');
title(['Averaged Velocities and Streamlines (',runname,')'])
rectangle('Position',[xtip,ytip,width,height],'FaceColor','k')
rectangle('Position',[xtip,ytip-gap,width,height],'FaceColor','k')
startx = (min(min(Yq))+.25)*ones(size(starty));
hgridinterp = streamline(Yq,-Xq,Vq,-Uq,startx,starty); hold on;
set(hgridinterp,'color','red','linewidth',3);