source: altifloat/matlab_toolbox/twinexp2floaters.m @ 199

% 1 floater 2 pts
ls
load obs_float_test.dat
x1=obs_float_test (1:2:40,4)
y1=obs_float_test (2:2:40,4)
x1=x1+1;
y1=y1+1;
plot(x1(1),y1(1),'+')
hold on
plot(x1(20),y1(20),'>')
hold on
plot(x1(1:20),y1(1:20),'-r')
hold on
x2=obs_float_test (41:2:80,4)
y2=obs_float_test (42:2:80,4)
x2=x2+1;
y2=y2+1;
plot(x2(1),y2(1),'+')
hold on
plot(x2(20),y2(20),'>')
hold on
plot(x2(1:20),y2(1:20),'-r')
hold on
load rfloat_total.dat;
whos;
rfloat_total(:,:)
x3=rfloat_total(1:20,4)
y3=rfloat_total(1:20,3)
x3=x3+1;
y3=y3+1;
plot(x3(2),y3(2),'+')
hold on
plot(x3(19),y3(19),'>')
hold on
plot(x3(2:19),y3(2:19),'-g')
hold on
x4=rfloat_total(21:40,4)
y4=rfloat_total(21:40,3)
x4=x4+1;
y4=y4+1;
plot(x4(2),y4(2),'+')
hold on
plot(x4(19),y4(19),'>')
hold on
plot(x4(2:19),y4(2:19),'-g')
hold on
ls;


load uv_total.dat
X=uv_total(:,end-3);
Y=uv_total(:,end-2);
X=X+1;
Y=Y+1;
U=uv_total(:,end-1);
V=uv_total(:,end);
quiver(X,Y,U,V,'r');

load vel_bck_all.dat;
whos;
X=vel_bck_all(:,end-3);
Y=vel_bck_all(:,end-2);
U=vel_bck_all(:,end-1);
V=vel_bck_all(:,end);
quiver(X,Y,U,V,'b');

load vel_next_all.dat;
whos;
X=vel_next_all(:,end-3);
Y=vel_next_all(:,end-2);
U=vel_next_all(:,end-1);
V=vel_next_all(:,end);
quiver(X,Y,U,V,'g');