ls;
load obs_float_test.dat;
whos;
obs_float_test (:,:)
x=obs_float_test (1:2:24,4)
y=obs_float_test (2:2:24,4)
x=x+1
y=y+1

   plot(x(1),y(1),'+')
   hold on;
   plot(x(12),y(12),'>')
   hold on;
   plot(x(1:12),y(1:12),'+')
   hold on;
   plot(x(1:12),y(1:12),'-b')
   hold on;1


load r_float_optim25float-filtre15-inc_16.dat;
whos;
r_float_optim25float_filtre15_inc_16(:,:)
x1=r_float_optim25float_filtre15_inc_16(1:8:24,4)
y1=r_float_optim25float_filtre15_inc_16(5:8:24,4)
x1=x1+1
y1=y1+1
plot(x1(1),y1(1),'+')
hold on;
plot(x1(3),y1(3),'>')
hold on;
plot(x1(1:3),y1(1:3),'+')
plot(x1(1:3),y1(1:3),'g')
hold on;
ls;
r_float_optim25float_filtre15_inc_16(:,:)
x2=r_float_optim25float_filtre15_inc_16(2:8:24,4)
y2=r_float_optim25float_filtre15_inc_16(6:8:24,4)
x2=x2+1
y2=y2+1
plot(x2(1),y2(1),'+')
hold on;
plot(x2(3),y2(3),'>')
hold on;
plot(x2(1:3),y2(1:3),'+')
plot(x2(1:3),y2(1:3),'g')
hold on;
ls;
r_float_optim25float_filtre15_inc_16(:,:)
x3=r_float_optim25float_filtre15_inc_16(3:8:24,4)
y3=r_float_optim25float_filtre15_inc_16(7:8:24,4)
x3=x3+1
y3=y3+1
plot(x3(1),y3(1),'+')
hold on;
plot(x3(3),y3(3),'>')
hold on;
plot(x3(1:3),y3(1:3),'+')
plot(x3(1:3),y3(1:3),'g')
hold on;
ls;
r_float_optim25float_filtre15_inc_16(:,:)
x4=r_float_optim25float_filtre15_inc_16(4:8:24,4)
y4=r_float_optim25float_filtre15_inc_16(8:8:24,4)
x4=x4+1
y4=y4+1
plot(x4(1),y4(1),'+')
hold on;
plot(x4(3),y4(3),'>')
hold on;
plot(x4(1:3),y4(1:3),'+')
plot(x4(1:3),y4(1:3),'g')
hold on;
ls;


load uv_total.dat
    
X=uv_total(1:5046,end-3);
Y=uv_total(1:5046,end-2);
U1=uv_total(1:5046,end-1);
V1=uv_total(1:5046,end);
U2=uv_total(5047:10092,end-1);
V2=uv_total(5047:10092,end);
U3=uv_total(10093:15138,end-1);
V3=uv_total(10093:15138,end);
    
U=(U1+U2+U3)/3;
V=(V1+V2+V3)/3;

quiver(X,Y,U,V)




% for k=1:16
% load uzero.dat;
% whos;
% X=uzero(:,end-2);
% Y=uzero(:,end-1);
% U=uzero(:,end);
% ls;
% load vzero.dat;
% whos;
% V=vzero(:,end);
% figure(k);
% quiver(X,Y,U,V);
% end