close all

figure(1)

fid = fopen('obs_float/obs_float_test.dat', 'rb');
fseek(fid, 0, 'eof');
fileSize = ftell(fid);
frewind(fid);
data = fread(fid, fileSize, 'uint8');
numLines = sum(data == 10) + 1;
fclose(fid);
m=numLines-1;
n=m/2;
load obs_float/obs_float_test.dat;

x=obs_float_test (1:2:4,4);
y=obs_float_test (2:2:4,4);
x=x+1;
y=y+1;
plot(x(1),y(1),'+')
hold on;
plot(x(2),y(2),'>')
hold on;
plot(x(1:2),y(1:2),'+');
hold on;
plot(x(1:2),y(1:2),'-r')
hold on;

% x3=obs_float_test (n+1:2:m,4);
% y3=obs_float_test (n+2:2:m,4);
% % x3=x3+1;
% % y3=y3+1;
% plot(x3(1),y3(1),'+')
% hold on;
% plot(x3(n/2),y3(n/2),'>')
% hold on;
% plot(x3(1:n/2),y3(1:n/2),'+');
% hold on;
% plot(x3(1:n/2),y3(1:n/2),'-m')
% hold on;


fid = fopen('obs_float/rfloat_total.dat', 'rb');
fseek(fid, 0, 'eof');
fileSize = ftell(fid);
frewind(fid);
data = fread(fid, fileSize, 'uint8');
numLines = sum(data == 10) + 1;
fclose(fid);
p=numLines-1;
u=p/2;
load obs_float/rfloat_total.dat
x2=rfloat_total(1:4,4);
y2=rfloat_total(1:4,3);
ligne=1;
for h=1:3:4
    x2(h)=obs_float_test(ligne,4);
    y2(h)=obs_float_test(ligne+1,4);
    ligne=ligne+2;
end
x2=x2+1;
y2=y2+1;
[k,h]=size(x2);
plot(x2(1),y2(1),'+')
hold on;
plot(x2(k),y2(k),'>')
hold on;
plot(x2(1:k),y2(1:k),'+')
hold on;
plot(x2(1:k),y2(1:k),'-g')
hold on;

% rfloat_total(:,:)
% x4=rfloat_total(u+1:p,4)
% y4=rfloat_total(u+1:p,3)
% ligne=n+1;
% for h=1:3:u
%     x4(h)=obs_float_test(ligne,4);
%     y4(h)=obs_float_test(ligne+1,4);
%     ligne=ligne+2;
% end
% % x4=x4+1;
% % y4=y4+1;
% [k,h]=size(x4);
% plot(x4(1),y4(1),'+')
% hold on;
% plot(x4(k),y4(k),'>')
% hold on;
% plot(x4(1:k),y4(1:k),'+')
% hold on;
% plot(x4(1:k),y4(1:k),'-g')
% hold on;


load obs_float/uv_total.dat;
whos;
X1=uv_total(1:5046,end-3);
Y1=uv_total(1:5046,end-2);
X1=X1+1;
Y1=Y1+1;
%change si pas de temps change
U1=zeros(5046,1);
V1=zeros(5046,1);
e=1;
for r=1:1:2
    
    Unew=zeros(1,5046);
    Vnew=zeros(1,5046);
    for count=1:5046
        Unew(count)=uv_total(e,end-1);
        Vnew(count)=uv_total(e,end);
        e=e+1;
    end
    U1=U1+Unew';
    V1=V1+Vnew';
   
    
end

U1=U1./2;
V1=V1./2;
quiver(X1,Y1,U1,V1,2,'r','Autoscale','off');

% load obs_float/unext3d.dat;
% 
% unext3d=roundn(unext3d,-6)
% whos;
% X3=unext3d(:,end-2);
% Y3=unext3d(:,end-1);
% U3=unext3d(:,end);
% load obs_float/vnext3d.dat
% vnext3d=roundn(vnext3d,-6)
% V3=vnext3d(:,end);
% quiver(X3,Y3,U3,V3,2,'b','Autoscale','off');

load obs_float/uzero.dat;


whos;
X2=uzero(:,end-2);
Y2=uzero(:,end-1);
U2=uzero(:,end);
load obs_float/vzero.dat

V2=vzero(:,end);
quiver(X2,Y2,U2,V2,2,'b','Autoscale','off');



% 
% 
%  figure(2)
% load obs_float/vel_bck_all.dat
% it=find(vel_bck_all(:,1)==1);
% Tu=reshape(vel_bck_all(it,4),58,87);
% Tv=reshape(vel_bck_all(it,5),58,87);
% 
% load obs_float/vel_next_all.dat
% it=find(vel_next_all(:,1)==1);
% Bu=reshape(vel_next_all(it,4),58,87);
% Bv=reshape(vel_next_all(it,5),58,87);
% 
% F=(Bu-Tu).^2+(Bv-Tv).^2
% imagesc(F)
% axis xy
% colorbar
% hold on 
% obs_float_test (1:4,:)
% x=obs_float_test (1:2:214,4)
% y=obs_float_test (2:2:214,4)
% x=x+1
% y=y+1
% plot(x(1),y(1),'+')
% hold on;
% plot(x(40),y(40),'>')
% hold on;
% plot(x(1:40),y(1:40),'+');
% hold on;
% plot(x(1:40),y(1:40),'-r')
% hold on;
% 
% 
% x3=obs_float_test (289:2:502,4)
% y3=obs_float_test (290:2:502,4)
% x3=x3+1
% y3=y3+1
% plot(x3(1),y3(1),'+')
% hold on;
% plot(x3(40),y3(40),'>')
% hold on;
% plot(x3(1:40),y3(1:40),'+');
% hold on;
% plot(x3(1:40),y3(1:40),'-r')
% hold on;
% 
% x5=obs_float_test (577:2:790,4)
% y5=obs_float_test (578:2:790,4)
% x5=x5+1
% y5=y5+1
% plot(x5(1),y5(1),'+')
% hold on;
% plot(x5(34),y5(34),'>')
% hold on;
% plot(x5(1:34),y5(1:34),'+');
% hold on;
% plot(x5(1:34),y5(1:34),'-r')
% hold on;
% 
% x7=obs_float_test (865:2:348,4)
% y7=obs_float_test (866:2:348,4)
% x7=x7+1
% y7=y7+1
% plot(x7(1),y7(1),'+')
% hold on;
% plot(x7(34),y7(34),'>')
% hold on;
% plot(x7(1:34),y7(1:34),'+');
% hold on;
% plot(x7(1:34),y7(1:34),'-r')
% hold on;
% 
% figure(3)
% 
% load obs_float/uv_total.dat
% 
% steps=35
% counterr=1;
% 
% Cu=0;
% Cv=0;
% for ll=1:steps
%     Iu=0;
%     Iv=0;
%     it=counterr:(counterr+5045);
%     Iu=reshape(uv_total(it,4),58,87);
%     Iv=reshape(uv_total(it,5),58,87);
%     counterr=counterr+5046;
%     Cu=Cu+Iu;
%     Cv=Cv+Iv;
% end
%     
% Cu=Cu/35;
% Cv=Cv/35;
% 
% F1=(Tu-Cu).^2+(Tv-Cv).^2;
% imagesc(F1)
% axis xy
% colorbar
% hold on 
% obs_float_test (1:4,:)
% x=obs_float_test (1:2:214,4)
% y=obs_float_test (2:2:214,4)
% x=x+1
% y=y+1
% plot(x(1),y(1),'+')
% hold on;
% plot(x(34),y(34),'>')
% hold on;
% plot(x(1:34),y(1:34),'+');
% hold on;
% plot(x(1:34),y(1:34),'-r')
% hold on;
% 
% 
% x3=obs_float_test (289:2:502,4)
% y3=obs_float_test (290:2:502,4)
% x3=x3+1
% y3=y3+1
% plot(x3(1),y3(1),'+')
% hold on;
% plot(x3(34),y3(34),'>')
% hold on;
% plot(x3(1:34),y3(1:34),'+');
% hold on;
% plot(x3(1:34),y3(1:34),'-r')
% hold on;
% 
% x5=obs_float_test (577:2:790,4)
% y5=obs_float_test (578:2:790,4)
% x5=x5+1
% y5=y5+1
% plot(x5(1),y5(1),'+')
% hold on;
% plot(x5(34),y5(34),'>')
% hold on;
% plot(x5(1:34),y5(1:34),'+');
% hold on;
% plot(x5(1:34),y5(1:34),'-r')
% hold on;
% 
% x7=obs_float_test (865:2:348,4)
% y7=obs_float_test (866:2:348,4)
% x7=x7+1
% y7=y7+1
% plot(x7(1),y7(1),'+')
% hold on;
% plot(x7(34),y7(34),'>')
% hold on;
% plot(x7(1:34),y7(1:34),'+');
% hold on;
% plot(x7(1:34),y7(1:34),'-r')
% hold on;
% F2=(Bu-Cu).^2+(Bv-Cv).^2;
% figure(4)
% imagesc(F2)
% axis xy 
% colorbar
% hold on
% obs_float_test (1:4,:)
% x=obs_float_test (1:2:214,4)
% y=obs_float_test (2:2:214,4)
% x=x+1
% y=y+1
% plot(x(1),y(1),'+')
% hold on;
% plot(x(34),y(34),'>')
% hold on;
% plot(x(1:34),y(1:34),'+');
% hold on;
% plot(x(1:34),y(1:34),'-r')
% hold on;
% 
% 
% x3=obs_float_test (289:2:502,4)
% y3=obs_float_test (290:2:502,4)
% x3=x3+1
% y3=y3+1
% plot(x3(1),y3(1),'+')
% hold on;
% plot(x3(34),y3(34),'>')
% hold on;
% plot(x3(1:34),y3(1:34),'+');
% hold on;
% plot(x3(1:34),y3(1:34),'-r')
% hold on;
% 
% x5=obs_float_test (577:2:790,4)
% y5=obs_float_test (578:2:790,4)
% x5=x5+1
% y5=y5+1
% plot(x5(1),y5(1),'+')
% hold on;
% plot(x5(34),y5(34),'>')
% hold on;
% plot(x5(1:34),y5(1:34),'+');
% hold on;
% plot(x5(1:34),y5(1:34),'-r')
% hold on;
% 
% x7=obs_float_test (865:2:348,4)
% y7=obs_float_test (866:2:348,4)
% x7=x7+1
% y7=y7+1
% plot(x7(1),y7(1),'+')
% hold on;
% plot(x7(34),y7(34),'>')
% hold on;
% plot(x7(1:34),y7(1:34),'+');
% hold on;
% plot(x7(1:34),y7(1:34),'-r')
% hold on;
% 
% figure (5)
% imagesc(F)
% axis([71 77 33 37])
% hold on
% obs_float_test (1:4,:)
% x=obs_float_test (1:2:214,4)
% y=obs_float_test (2:2:214,4)
% x=x+1
% y=y+1
% plot(x(1),y(1),'+')
% hold on;
% plot(x(34),y(34),'>')
% hold on;
% plot(x(1:34),y(1:34),'+');
% hold on;
% plot(x(1:34),y(1:34),'-r')
% hold on;
% 
% 
% x3=obs_float_test (289:2:502,4)
% y3=obs_float_test (290:2:502,4)
% x3=x3+1
% y3=y3+1
% plot(x3(1),y3(1),'+')
% hold on;
% plot(x3(34),y3(34),'>')
% hold on;
% plot(x3(1:34),y3(1:34),'+');
% hold on;
% plot(x3(1:34),y3(1:34),'-r')
% hold on;
% 
% x5=obs_float_test (577:2:790,4)
% y5=obs_float_test (578:2:790,4)
% x5=x5+1
% y5=y5+1
% plot(x5(1),y5(1),'+')
% hold on;
% plot(x5(34),y5(34),'>')
% hold on;
% plot(x5(1:34),y5(1:34),'+');
% hold on;
% plot(x5(1:34),y5(1:34),'-r')
% hold on;
% 
% x7=obs_float_test (865:2:348,4)
% y7=obs_float_test (866:2:348,4)
% x7=x7+1
% y7=y7+1
% plot(x7(1),y7(1),'+')
% hold on;
% plot(x7(34),y7(34),'>')
% hold on;
% plot(x7(1:34),y7(1:34),'+');
% hold on;
% plot(x7(1:34),y7(1:34),'-r')
% hold on;
% 
% figure (6)
% imagesc(F1)
% 
% axis([71 77 33 37])
% hold on 
% obs_float_test (1:4,:)
% x=obs_float_test (1:2:214,4)
% y=obs_float_test (2:2:214,4)
% x=x+1
% y=y+1
% plot(x(1),y(1),'+')
% hold on;
% plot(x(34),y(34),'>')
% hold on;
% plot(x(1:34),y(1:34),'+');
% hold on;
% plot(x(1:34),y(1:34),'-r')
% hold on;
% 
% 
% x3=obs_float_test (289:2:502,4)
% y3=obs_float_test (290:2:502,4)
% x3=x3+1
% y3=y3+1
% plot(x3(1),y3(1),'+')
% hold on;
% plot(x3(34),y3(34),'>')
% hold on;
% plot(x3(1:34),y3(1:34),'+');
% hold on;
% plot(x3(1:34),y3(1:34),'-r')
% hold on;
% 
% x5=obs_float_test (577:2:790,4)
% y5=obs_float_test (578:2:790,4)
% x5=x5+1
% y5=y5+1
% plot(x5(1),y5(1),'+')
% hold on;
% plot(x5(34),y5(34),'>')
% hold on;
% plot(x5(1:34),y5(1:34),'+');
% hold on;
% plot(x5(1:34),y5(1:34),'-r')
% hold on;
% 
% x7=obs_float_test (865:2:348,4)
% y7=obs_float_test (866:2:348,4)
% x7=x7+1
% y7=y7+1
% plot(x7(1),y7(1),'+')
% hold on;
% plot(x7(34),y7(34),'>')
% hold on;
% plot(x7(1:34),y7(1:34),'+');
% hold on;
% plot(x7(1:34),y7(1:34),'-r')
% hold on;
% sum=0;
% count=0;
% for r=71:76
%     for t=33:37
%         sum=sum+F(t,r);
%         count=count+1;
%     end
% end
% sum=sum/count
% 
% sum3=0;
% count3=0;
% for r=71:76
%     for t=33:37
%         sum3=sum3+F1(t,r);
%         count3=count3+1;
%     end
% end
% sum3=sum3/count3
% sum=sqrt(sum/count);
% sum3=sqrt(sum3/count3);
% errorbefore=sum
% errorafter=sum3
