source: altifloat/matlab_toolbox/takecorrectionasbck.m~ @ 199

function [ output_args ] = create_bck( totaltime )

totaltime=totaltime-1;
outfile='obs_float/ucorr.dat';
load obs_float/uv_total.dat;
load obs_float/uzero.dat

umean(:,1)=uzero(:,1);
umean(:,2)=uzero(:,2);
umean(:,3)=uzero(:,3);
U1=zeros(5046,1);
V1=zeros(5046,1);
e=1;
for r=1:1:totaltime
    
    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./totaltime;
V1=V1./totaltime;
%%% Because if the jtlag is 3 we have a set of zeros each 3 time steps in
%%% the uv_total.
U1=U1;
V1=V1;
umean(:,4)=U1;


fid=fopen(outfile,'w');
fprintf(fid,'%d %d %d %f\n',umean');
fclose(fid);

outfile='obs_float/vcorr.dat';
load obs_float/uv_total.dat;
load obs_float/uzero.dat

vmean(:,1)=uzero(:,1);
vmean(:,2)=uzero(:,2);
vmean(:,3)=uzero(:,3);
U1=zeros(5046,1);
V1=zeros(5046,1);
e=1;
for r=1:1:totaltime
    
    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./totaltime;
V1=V1./totaltime;
%%% Because if the jtlag is 3 we have a set of zeros each 3 time steps in
%%% the uv_total.
U1=U1;
V1=V1;
vmean(:,4)=V1;


fid=fopen(outfile,'w');
fprintf(fid,'%d %d %d %f\n',vmean');
fclose(fid);
end