plotting_jetdata_uncertainty.m

%Make video for uncertainty
clear;
clc;
basedir='/scratch/shannon/a/bhattac3/';

filedir='Planar_Uncertainty_work/Results/final_plots/03_21_2017/WS64/matfiles/';

load([basedir,filedir,'Jetdata.mat']);

v = VideoWriter('Jetdata64.avi');
v.FrameRate=10;
%%
% Create 13-level blue/red diverging color map:
level = 13; n = ceil(level/2);
cmap1 = [linspace(1, 1, n); linspace(0, 1, n); linspace(0, 1, n)]';
cmap2 = [linspace(1, 0, n); linspace(1, 0, n); linspace(1, 1, n)]';
cmap = [cmap1; cmap2(2:end, :)];
%%
open(v);
Nim=size(err_up,3);

for i=1:5:Nim
    
    figure(1);
    subplot(1,2,1);
    velmag=sqrt(Up(:,:,i).^2+Vp(:,:,i).^2);
    imagesc(Xp(:),Yp(:),velmag);caxis([0 3.5]);
    hold on;
    quiver(Xp,Yp,Up(:,:,i),Vp(:,:,i),1,'k');
    hold off;
    xlim([320 384]);
    axis off;
    
    if i<=Nim-20
        subplot(1,2,2);
        UMC=sqrt(UMCx(:,:,i).^2+UMCy(:,:,i).^2);
        imagesc(Xp(:),Yp(:),UMC);colormap('jet');caxis([0 0.12]);
        xlim([320 384]);
        axis off;
    else
        subplot(1,2,2);
        UMC=sqrt(mean(UMCx,3).^2+mean(UMCy,3).^2);
        imagesc(Xp(:),Yp(:),UMC);colormap('jet');caxis([0 0.12]);
        xlim([320 384]);
        axis off;
    end
    Frame1=getframe(gcf);
    pause(0.01);
    writeVideo(v,Frame1)
end
close(v);