source: altifloat/matlab_toolbox/nemedofthreejtlag3.m @ 160

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

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?

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

x5=third_observation (1:2:m5,4);
y5=third_observation (2:2:m5,4);
x=x+1;
y=y+1;
plot(x5(1),y5(1),'+')
hold on;
plot(x5(n5),y5(n5),'>')
hold on;
plot(x5(1:n5),y5(1:n5),'+');
hold on;
plot(x5(1:n5),y5(1:n5),'-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:u,4);
y2=rfloat_total(1:u,3);
ligne=1;
for h=1:3:u
    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;

%%% To show the same path but forwarded with the bck that is beeing used 
fid = fopen('obs_float/generated_fwd.dat', 'rb');
fseek(fid, 0, 'eof');
fileSize = ftell(fid);
frewind(fid);
data = fread(fid, fileSize, 'uint8');
numLines = sum(data == 10) + 1;
fclose(fid);
p3=numLines-1;
u3=p3/3;

load obs_float/generated_fwd.dat
x11=generated_fwd (1:2:u3,4);
y11=generated_fwd (2:2:u3,4);
x=x+1;
y=y+1;
plot(x11(1),y11(1),'+')
hold on;
plot(x11(u3/2),y11(u3/2),'>')
hold on;
plot(x11(1:u3/2),y11(1:u3/2),'+');
hold on;
plot(x11(1:u3/2),y11(1:u3/2),'-b')
hold on;

x13=generated_fwd (u3+1:2:2*u3,4);
y13=generated_fwd (u3+2:2:2*u3,4);
x3=x3+1;
y3=y3+1;
plot(x13(1),y13(1),'+')
hold on;
plot(x13(u3/2),y13(u3/2),'>')
hold on;
plot(x13(1:u3/2),y13(1:u3/2),'+');
hold on;
plot(x13(1:u3/2),y13(1:u3/2),'-b')
hold on;

x15=generated_fwd (2*u3+1:2:p3,4);
y15=generated_fwd (2*u3+2:2:p3,4);
x3=x3+1;
y3=y3+1;
plot(x15(1),y15(1),'+')
hold on;
plot(x15(u3/2),y15(u3/2),'>')
hold on;
plot(x15(1:u3/2),y15(1:u3/2),'+');
hold on;
plot(x15(1:u3/2),y15(1:u3/2),'-b')
hold on;

%%%% To show the same path but forwarded with the mean that is beeing used 
fid = fopen('obs_float/forwardedbymean.dat', 'rb');
fseek(fid, 0, 'eof');
fileSize = ftell(fid);
frewind(fid);
data = fread(fid, fileSize, 'uint8');
numLines = sum(data == 10) + 1;
fclose(fid);
p4=numLines-1;
u4=p4/3;
% 
load obs_float/forwardedbymean.dat
x11=forwardedbymean (1:2:u4,4);
y11=forwardedbymean (2:2:u4,4);
% x=x+1;
% y=y+1;
plot(x11(1),y11(1),'+')
hold on;
plot(x11(u4/2),y11(u4/2),'>')
hold on;
plot(x11(1:u4/2),y11(1:u4/2),'+');
hold on;
plot(x11(1:u4/2),y11(1:u4/2),'-r')
hold on;

x13=forwardedbymean (u4+1:2:2*u4,4);
y13=forwardedbymean (u4+2:2:2*u4,4);
% x3=x3+1;
% y3=y3+1;
plot(x13(1),y13(1),'+')
hold on;
plot(x13(u4/2),y13(u4/2),'>')
hold on;
plot(x13(1:u4/2),y13(1:u4/2),'+');
hold on;
plot(x13(1:u4/2),y13(1:u4/2),'-r')
hold on;

x15=forwardedbymean (2*u4+1:2:p4,4);
y15=forwardedbymean (2*u4+2:2:p4,4);
% x3=x3+1;
% y3=y3+1;
plot(x15(1),y15(1),'+')
hold on;
plot(x15(u4/2),y15(u4/2),'>')
hold on;
plot(x15(1:u4/2),y15(1:u4/2),'+');
hold on;
plot(x15(1:u4/2),y15(1:u4/2),'-r')
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:117
    
    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./117;
V1=V1./117;
%%% Because if the jtlag is 3 we have a set of zeros each 3 time steps in
%%% the uv_total.
U1=U1;
V1=V1;
% U1=roundn(U1,-6);
% V1=roundn(V1,-6);

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,'g','Autoscale','off');

load obs_float/uzero.dat;

uzero=roundn(uzero,-6)
whos;
X2=uzero(:,end-2);
Y2=uzero(:,end-1);
U2=uzero(:,end);
load obs_float/vzero.dat
vzero=roundn(vzero,-6)
V2=vzero(:,end);
quiver(X2,Y2,U2,V2,2,'b','Autoscale','off');



% 
% 
%  figure(2)
% load obs_float/vel_bck_all.dat
% vel_bck_all=roundn(vel_bck_all,-6);
% 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/uv_total.dat
% 
% % steps=50
% % counterr=1;
% % 
% % Cu=0;
% % Cv=0;
% % for ll=1:steps
% %     Iu=0;
% %     Iv=0;
% %     it=counterr:(counterr+5045);
% % U1=roundn(U1,-6);
% % V1=roundn(V1,-6);
% 
%     Cu=reshape(U1,58,87);
%     Cv=reshape(V1,58,87);
% %     counterr=counterr+5046;
% %     Cu=Cu+Iu;
% %     Cv=Cv+Iv;
% % end
%     
% % Cu=Cu/50;
% % Cv=Cv/50;
% F1=(Tu-Cu).^2+(Tv-Cv).^2;
% imagesc(F1)
% axis xy
% colorbar
% 
% % 
% %  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
% % % caxis([0 15*10^-15])
% % hold on 
% % x=obs_float_test (1:2:n,4);
% % y=obs_float_test (2:2:n,4);
% % x=x+1;
% % y=y+1;
% % plot(x(1),y(1),'+')
% % hold on;
% % plot(x(n/2),y(n/2),'>')
% % hold on;
% % plot(x(1:n/2),y(1:n/2),'+');
% % hold on;
% % plot(x(1:n/2),y(1:n/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;
% % 
% % figure(3)
% % 
% % load obs_float/uv_total.dat
% % 
% % steps=100
% % 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/100;
% % Cv=Cv/100;
% % 
% % F1=(Tu-Cu).^2+(Tv-Cv).^2;
% % imagesc(F1)
% % axis xy
% % colorbar
% % hold on 
% % x=obs_float_test (1:2:n,4);
% % y=obs_float_test (2:2:n,4);
% % x=x+1;
% % y=y+1;
% % plot(x(1),y(1),'+')
% % hold on;
% % plot(x(n/2),y(n/2),'>')
% % hold on;
% % plot(x(1:n/2),y(1:n/2),'+');
% % hold on;
% % plot(x(1:n/2),y(1:n/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;
% % 
% % % F2=(Bu-Cu).^2+(Bv-Cv).^2;
% % % figure(4)
% % % imagesc(F2)
% % % axis xy 
% % % colorbar
% % % hold on
% % x=obs_float_test (1:2:n,4);
% % y=obs_float_test (2:2:n,4);
% % x=x+1;
% % y=y+1;
% % plot(x(1),y(1),'+')
% % hold on;
% % plot(x(n/2),y(n/2),'>')
% % hold on;
% % plot(x(1:n/2),y(1:n/2),'+');
% % hold on;
% % plot(x(1:n/2),y(1:n/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;
% % 
% % % figure (5)
% % % imagesc(F)
% % % axis([71 77 33 37])
% % % hold on
% % x=obs_float_test (1:2:n,4);
% % y=obs_float_test (2:2:n,4);
% % x=x+1;
% % y=y+1;
% % plot(x(1),y(1),'+')
% % hold on;
% % plot(x(n/2),y(n/2),'>')
% % hold on;
% % plot(x(1:n/2),y(1:n/2),'+');
% % hold on;
% % plot(x(1:n/2),y(1:n/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;
% % 
% % % 
% % % figure (6)
% % % imagesc(F1)
% % % 
% % % axis([71 77 33 37])
% % % hold on 
% % x=obs_float_test (1:2:n,4);
% % y=obs_float_test (2:2:n,4);
% % x=x+1;
% % y=y+1;
% % plot(x(1),y(1),'+')
% % hold on;
% % plot(x(n/2),y(n/2),'>')
% % hold on;
% % plot(x(1:n/2),y(1:n/2),'+');
% % hold on;
% % plot(x(1:n/2),y(1:n/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;
% % 
% % % 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