U=zeros(size(v_filtered{1,1}));
V=U;
L=numel(u_filtered);
F=U;
f=cell(L,1);
tic;
%%%%%%  Rifracting Indexes  %%%%%%
N_OIL = 1.403;
N_AIR = 1;
N_PERSPEX = 1.49; 

%%%%%%  Distances  %%%%%% [m]
H_TOT=1.71;
H_OIL = 3*10^-3; %??
H_PERSPEX = 4*10^-3; %??
H_AIR=25*10^-3;

alpha = 1 - N_AIR/N_OIL;
h_p = H_OIL + H_PERSPEX * (N_OIL/N_PERSPEX) + H_AIR * (N_PERSPEX/N_AIR);
h_star=((1/(alpha*h_p))-1/H_TOT)^-1;
X=x{1,1};
dx=X(1,2)-X(1,1); 
clear H_TOT H_OIL H_PERSPEX H_AIR alpha h_p X N_OIL N_AIR N_PERSPEX
for i=1:L;
    
    u=u_filtered{i,1}-mean(mean(u_filtered{i,1}));
    v=v_filtered{i,1}-mean(mean(v_filtered{i,1}));
    r=5; %how much to clear from each side of the original matrix
    u=u(r+1:(end-r),r+1:(end-r));
    v=v(r+1:(end-r),r+1:(end-r));
    u=-u/h_star;
    v=-v/h_star;
    
    
    F=intgrad2(u,v,dx,dx);
    F=F-(mean(mean(F)));
    f{i,1}=F;
    Axis=dx*(size(u_filtered{i,1},1)-1);
    [X,Y]=meshgrid([0:dx:Axis],[0:dx:Axis]);
    surf(X*1000,Y*1000,imgaussfilt(1000*F,0.6));
    
    set(gcf,'PaperPositionMode','auto')
    view(-30,80)
    xlim([0 Axis*1000]);ylim([0 Axis*1000]);zlim([-0.015 0.015]);%axis equal; 
    xlabel('mm','FontSize',16,'FontWeight','bold')
    ylabel('mm','FontSize',16,'FontWeight','bold')
    zlabel('mm','FontSize',16,'FontWeight','bold')
    
    colorbar
    caxis([-0.010 0.010])
    
    filename = [sprintf('%03d',i) '.jpg'];
    fullname = fullfile('directory to save at','video',filename);
    
    print(fullname,'-dtiffn');
end
toc
clear filename u v fullname