source: altifloat/matlab_toolbox/threefloaterslags.m @ 160

% 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:42,4)
y=obs_float_test (2:2:42,4)
x=x+1;
y=y+1;
plot(x(1),y(1),'+')
hold on;
plot(x(21),y(21),'>')
hold on;
plot(x(1:21),y(1:21),'+');
hold on;
plot(x(1:21),y(1:21),'-r')
hold on;


x3=obs_float_test (43:2:84,4)
y3=obs_float_test (44:2:84,4)
x3=x3+1;
y3=y3+1;
plot(x3(1),y3(1),'+')
hold on;
plot(x3(21),y3(21),'>')
hold on;
plot(x3(1:21),y3(1:21),'+');
hold on;
plot(x3(1:21),y3(1:21),'-r')
hold on;


x5=obs_float_test (85:2:126,4)
y5=obs_float_test (86:2:126,4)
x5=x5+1;
y5=y5+1;
plot(x5(1),y5(1),'+')
hold on;
plot(x5(21),y5(21),'>')
hold on;
plot(x5(1:21),y5(1:21),'+');
hold on;
plot(x5(1:21),y5(1:21),'-r')
hold on;


load obs_float/rfloat_total.dat
rfloat_total(:,:)
x2=rfloat_total(1:60,4)
y2=rfloat_total(1:60,3)
ligne=1;
for h=1:3:60
    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(63:122,4)
y4=rfloat_total(63:122,3)
ligne=43;
for h=1:3:60
    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;

rfloat_total(:,:)
x6=rfloat_total(125:184,4)
y6=rfloat_total(125:184,3)
ligne=85;
for h=1:3:60
    x6(h)=obs_float_test(ligne,4);
    y6(h)=obs_float_test(ligne+1,4);
    ligne=ligne+2;
end
x6=x6+1
y6=y6+1
[k,h]=size(x6);
plot(x6(1),y6(1),'+')
hold on;
plot(x6(k),y6(k),'>')
hold on;
plot(x6(1:k),y6(1:k),'+')
hold on
plot(x6(1:k),y6(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:55
    
    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/55;
V1=V1/55;
quiver(X1,Y1,U1,V1,2,'r','Autoscale','off');


load obs_float/uzero.dat;

uzero=roundn(uzero,-6)
whos;
X3=uzero(:,end-2);
Y3=uzero(:,end-1);
U3=uzero(:,end);
load obs_float/vzero.dat
vzero=roundn(vzero,-6)
V3=vzero(:,end);
quiver(X3,Y3,U3,V3,2,'c','Autoscale','off');