% 1 floater 2 pts

close all

figure(1)

ls;
load obs_float/obs_float_test.dat;
whos;
obs_float_test (1:4,:)
x=obs_float_test (1:2:350,4)
y=obs_float_test (2:2:350,4)
x=x+1
y=y+1
plot(x(1),y(1),'+')
hold on;
plot(x(175),y(175),'>')
hold on;
plot(x(1:175),y(1:175),'+');
hold on;
plot(x(1:175),y(1:175),'-r')
hold on;


x3=obs_float_test (351:2:700,4)
y3=obs_float_test (352:2:700,4)
x3=x3+1
y3=y3+1
plot(x3(1),y3(1),'+')
hold on;
plot(x3(175),y3(175),'>')
hold on;
plot(x3(1:175),y3(1:175),'+');
hold on;
plot(x3(1:175),y3(1:175),'-r')
hold on;

x5=obs_float_test (701:2:1050,4)
y5=obs_float_test (702:2:1050,4)
x5=x5+1
y5=y5+1
plot(x5(1),y5(1),'+')
hold on;
plot(x5(175),y5(175),'>')
hold on;
plot(x5(1:175),y5(1:175),'+');
hold on;
plot(x5(1:175),y5(1:175),'-r')
hold on;

x7=obs_float_test (1051:2:1400,4)
y7=obs_float_test (1052:2:1400,4)
x7=x7+1
y7=y7+1
plot(x7(1),y7(1),'+')
hold on;
plot(x7(175),y7(175),'>')
hold on;
plot(x7(1:175),y7(1:175),'+');
hold on;
plot(x7(1:175),y7(1:175),'-r')
hold on;

x9=obs_float_test (1401:2:1750,4)
y9=obs_float_test (1402:2:1750,4)
x9=x9+1
y9=y9+1
plot(x9(1),y9(1),'+')
hold on;
plot(x9(175),y9(175),'>')
hold on;
plot(x9(1:175),y9(1:175),'+');
hold on;
plot(x9(1:175),y9(1:175),'-r')
hold on;

load obs_float/rfloat_total.dat
rfloat_total(:,:)
x2=rfloat_total(1:145,4)
y2=rfloat_total(1:145,3)
ligne=1;
for h=1:36:144
    x2(h)=obs_float_test(ligne,4);
    y2(h)=obs_float_test(ligne+1,4);
    ligne=ligne+72;
end
x2=x2+1
y2=y2+1

[k,h]=size(x2);
plot(x2(1),y2(1),'+')
hold on;
plot(x2(k-1),y2(k-1),'>')
hold on;
plot(x2(1:k-1),y2(1:k-1),'+')
hold on;
plot(x2(1:k-1),y2(1:k-1),'-g')
hold on;

rfloat_total(:,:)
x4=rfloat_total(176:320,4)
y4=rfloat_total(176:320,3)
ligne=351;
for h=1:36:144
    x4(h)=obs_float_test(ligne,4);
    y4(h)=obs_float_test(ligne+1,4);
    ligne=ligne+72;
end
x4=x4+1
y4=y4+1
[k,h]=size(x4);
plot(x4(1),y4(1),'+')
hold on;
plot(x4(k-1),y4(k-1),'>')
hold on;
plot(x4(1:k-1),y4(1:k-1),'+')
hold on;
plot(x4(1:k-1),y4(1:k-1),'-g')
hold on;

rfloat_total(:,:)
x6=rfloat_total(351:495,4)
y6=rfloat_total(351:495,3)
ligne=701;
for h=1:36:144
    x6(h)=obs_float_test(ligne,4);
    y6(h)=obs_float_test(ligne+1,4);
    ligne=ligne+72;
end
x6=x6+1
y6=y6+1
[k,h]=size(x6);
plot(x6(1),y6(1),'+')
hold on;
plot(x6(k-1),y6(k-1),'>')
hold on;
plot(x6(1:k-1),y6(1:k-1),'+')
hold on
plot(x6(1:k-1),y6(1:k-1),'-g')
hold on;

rfloat_total(:,:)
x8=rfloat_total(526:670,4)
y8=rfloat_total(526:670,3)
ligne=1051;
for h=1:36:144
    x8(h)=obs_float_test(ligne,4);
    y8(h)=obs_float_test(ligne+1,4);
    ligne=ligne+72;
end
x8=x8+1
y8=y8+1
[k,h]=size(x8);
plot(x8(1),y8(1),'+')
hold on;
plot(x8(k-1),y8(k-1),'>')
hold on;
plot(x8(1:k-1),y8(1:k-1),'+')
hold on;
plot(x8(1:k-1),y8(1:k-1),'-g')
hold on;

rfloat_total(:,:)
x10=rfloat_total(701:845,4)
y10=rfloat_total(701:845,3)
ligne=1401;
for h=1:36:144
    x10(h)=obs_float_test(ligne,4);
    y10(h)=obs_float_test(ligne+1,4);
    ligne=ligne+72;
end
x10=x10+1
y10=y10+1
[k,h]=size(x10);
plot(x10(1),y10(1),'+')
hold on;
plot(x10(k-1),y10(k-1),'>')
hold on;
plot(x10(1:k-1),y10(1:k-1),'+')
hold on;
plot(x10(1:k-1),y10(1:k-1),'-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:140
%     
%     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./140;
% V1=V1./140;
% 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,'g','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:350,4)
% y=obs_float_test (2:2:350,4)
% x=x+1
% y=y+1
% plot(x(1),y(1),'+')
% hold on;
% plot(x(175),y(175),'>')
% hold on;
% plot(x(1:175),y(1:175),'+');
% hold on;
% plot(x(1:175),y(1:175),'-r')
% hold on;
% 
% 
% x3=obs_float_test (351:2:700,4)
% y3=obs_float_test (352:2:700,4)
% x3=x3+1
% y3=y3+1
% plot(x3(1),y3(1),'+')
% hold on;
% plot(x3(175),y3(175),'>')
% hold on;
% plot(x3(1:175),y3(1:175),'+');
% hold on;
% plot(x3(1:175),y3(1:175),'-r')
% hold on;
% 
% x5=obs_float_test (701:2:1050,4)
% y5=obs_float_test (702:2:1050,4)
% x5=x5+1
% y5=y5+1
% plot(x5(1),y5(1),'+')
% hold on;
% plot(x5(175),y5(175),'>')
% hold on;
% plot(x5(1:175),y5(1:175),'+');
% hold on;
% plot(x5(1:175),y5(1:175),'-r')
% hold on;
% 
% x7=obs_float_test (1051:2:1400,4)
% y7=obs_float_test (1052:2:1400,4)
% x7=x7+1
% y7=y7+1
% plot(x7(1),y7(1),'+')
% hold on;
% plot(x7(175),y7(175),'>')
% hold on;
% plot(x7(1:175),y7(1:175),'+');
% hold on;
% plot(x7(1:175),y7(1:175),'-r')
% hold on;
% 
% x9=obs_float_test (1401:2:1750,4)
% y9=obs_float_test (1402:2:1750,4)
% x9=x9+1
% y9=y9+1
% plot(x9(1),y9(1),'+')
% hold on;
% plot(x9(175),y9(175),'>')
% hold on;
% plot(x9(1:175),y9(1:175),'+');
% hold on;
% plot(x9(1:175),y9(1:175),'-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:350,4)
% y=obs_float_test (2:2:350,4)
% x=x+1
% y=y+1
% plot(x(1),y(1),'+')
% hold on;
% plot(x(175),y(175),'>')
% hold on;
% plot(x(1:175),y(1:175),'+');
% hold on;
% plot(x(1:175),y(1:175),'-r')
% hold on;
% 
% 
% x3=obs_float_test (351:2:700,4)
% y3=obs_float_test (352:2:700,4)
% x3=x3+1
% y3=y3+1
% plot(x3(1),y3(1),'+')
% hold on;
% plot(x3(175),y3(175),'>')
% hold on;
% plot(x3(1:175),y3(1:175),'+');
% hold on;
% plot(x3(1:175),y3(1:175),'-r')
% hold on;
% 
% x5=obs_float_test (701:2:1050,4)
% y5=obs_float_test (702:2:1050,4)
% x5=x5+1
% y5=y5+1
% plot(x5(1),y5(1),'+')
% hold on;
% plot(x5(175),y5(175),'>')
% hold on;
% plot(x5(1:175),y5(1:175),'+');
% hold on;
% plot(x5(1:175),y5(1:175),'-r')
% hold on;
% 
% x7=obs_float_test (1051:2:1400,4)
% y7=obs_float_test (1052:2:1400,4)
% x7=x7+1
% y7=y7+1
% plot(x7(1),y7(1),'+')
% hold on;
% plot(x7(175),y7(175),'>')
% hold on;
% plot(x7(1:175),y7(1:175),'+');
% hold on;
% plot(x7(1:175),y7(1:175),'-r')
% hold on;
% 
% x9=obs_float_test (1401:2:1750,4)
% y9=obs_float_test (1402:2:1750,4)
% x9=x9+1
% y9=y9+1
% plot(x9(1),y9(1),'+')
% hold on;
% plot(x9(175),y9(175),'>')
% hold on;
% plot(x9(1:175),y9(1:175),'+');
% hold on;
% plot(x9(1:175),y9(1:175),'-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:350,4)
% y=obs_float_test (2:2:350,4)
% x=x+1
% y=y+1
% plot(x(1),y(1),'+')
% hold on;
% plot(x(175),y(175),'>')
% hold on;
% plot(x(1:175),y(1:175),'+');
% hold on;
% plot(x(1:175),y(1:175),'-r')
% hold on;
% 
% 
% x3=obs_float_test (351:2:700,4)
% y3=obs_float_test (352:2:700,4)
% x3=x3+1
% y3=y3+1
% plot(x3(1),y3(1),'+')
% hold on;
% plot(x3(175),y3(175),'>')
% hold on;
% plot(x3(1:175),y3(1:175),'+');
% hold on;
% plot(x3(1:175),y3(1:175),'-r')
% hold on;
% 
% x5=obs_float_test (701:2:1050,4)
% y5=obs_float_test (702:2:1050,4)
% x5=x5+1
% y5=y5+1
% plot(x5(1),y5(1),'+')
% hold on;
% plot(x5(175),y5(175),'>')
% hold on;
% plot(x5(1:175),y5(1:175),'+');
% hold on;
% plot(x5(1:175),y5(1:175),'-r')
% hold on;
% 
% x7=obs_float_test (1051:2:1400,4)
% y7=obs_float_test (1052:2:1400,4)
% x7=x7+1
% y7=y7+1
% plot(x7(1),y7(1),'+')
% hold on;
% plot(x7(175),y7(175),'>')
% hold on;
% plot(x7(1:175),y7(1:175),'+');
% hold on;
% plot(x7(1:175),y7(1:175),'-r')
% hold on;
% 
% x9=obs_float_test (1401:2:1750,4)
% y9=obs_float_test (1402:2:1750,4)
% x9=x9+1
% y9=y9+1
% plot(x9(1),y9(1),'+')
% hold on;
% plot(x9(175),y9(175),'>')
% hold on;
% plot(x9(1:175),y9(1:175),'+');
% hold on;
% plot(x9(1:175),y9(1:175),'-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:350,4)
% y=obs_float_test (2:2:350,4)
% x=x+1
% y=y+1
% plot(x(1),y(1),'+')
% hold on;
% plot(x(175),y(175),'>')
% hold on;
% plot(x(1:175),y(1:175),'+');
% hold on;
% plot(x(1:175),y(1:175),'-r')
% hold on;
% 
% 
% x3=obs_float_test (351:2:700,4)
% y3=obs_float_test (352:2:700,4)
% x3=x3+1
% y3=y3+1
% plot(x3(1),y3(1),'+')
% hold on;
% plot(x3(175),y3(175),'>')
% hold on;
% plot(x3(1:175),y3(1:175),'+');
% hold on;
% plot(x3(1:175),y3(1:175),'-r')
% hold on;
% 
% x5=obs_float_test (701:2:1050,4)
% y5=obs_float_test (702:2:1050,4)
% x5=x5+1
% y5=y5+1
% plot(x5(1),y5(1),'+')
% hold on;
% plot(x5(175),y5(175),'>')
% hold on;
% plot(x5(1:175),y5(1:175),'+');
% hold on;
% plot(x5(1:175),y5(1:175),'-r')
% hold on;
% 
% x7=obs_float_test (1051:2:1400,4)
% y7=obs_float_test (1052:2:1400,4)
% x7=x7+1
% y7=y7+1
% plot(x7(1),y7(1),'+')
% hold on;
% plot(x7(175),y7(175),'>')
% hold on;
% plot(x7(1:175),y7(1:175),'+');
% hold on;
% plot(x7(1:175),y7(1:175),'-r')
% hold on;
% 
% x9=obs_float_test (1401:2:1750,4)
% y9=obs_float_test (1402:2:1750,4)
% x9=x9+1
% y9=y9+1
% plot(x9(1),y9(1),'+')
% hold on;
% plot(x9(175),y9(175),'>')
% hold on;
% plot(x9(1:175),y9(1:175),'+');
% hold on;
% plot(x9(1:175),y9(1:175),'-r')
% hold on;
% figure (6)
% imagesc(F1)
% 
% axis([71 77 33 37])
% hold on 
% obs_float_test (1:4,:)
% obs_float_test (1:4,:)
% x=obs_float_test (1:2:350,4)
% y=obs_float_test (2:2:350,4)
% x=x+1
% y=y+1
% plot(x(1),y(1),'+')
% hold on;
% plot(x(175),y(175),'>')
% hold on;
% plot(x(1:175),y(1:175),'+');
% hold on;
% plot(x(1:175),y(1:175),'-r')
% hold on;
% 
% 
% x3=obs_float_test (351:2:700,4)
% y3=obs_float_test (352:2:700,4)
% x3=x3+1
% y3=y3+1
% plot(x3(1),y3(1),'+')
% hold on;
% plot(x3(175),y3(175),'>')
% hold on;
% plot(x3(1:175),y3(1:175),'+');
% hold on;
% plot(x3(1:175),y3(1:175),'-r')
% hold on;
% 
% x5=obs_float_test (701:2:1050,4)
% y5=obs_float_test (702:2:1050,4)
% x5=x5+1
% y5=y5+1
% plot(x5(1),y5(1),'+')
% hold on;
% plot(x5(175),y5(175),'>')
% hold on;
% plot(x5(1:175),y5(1:175),'+');
% hold on;
% plot(x5(1:175),y5(1:175),'-r')
% hold on;
% 
% x7=obs_float_test (1051:2:1400,4)
% y7=obs_float_test (1052:2:1400,4)
% x7=x7+1
% y7=y7+1
% plot(x7(1),y7(1),'+')
% hold on;
% plot(x7(175),y7(175),'>')
% hold on;
% plot(x7(1:175),y7(1:175),'+');
% hold on;
% plot(x7(1:175),y7(1:175),'-r')
% hold on;
% 
% x9=obs_float_test (1401:2:1750,4)
% y9=obs_float_test (1402:2:1750,4)
% x9=x9+1
% y9=y9+1
% plot(x9(1),y9(1),'+')
% hold on;
% plot(x9(175),y9(175),'>')
% hold on;
% plot(x9(1:175),y9(1:175),'+');
% hold on;
% plot(x9(1:175),y9(1:175),'-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