source: altifloat/matlab_toolbox/make_mean.m @ 160

function [ output_args ] = make_mean(  )
%UNTITLED Summary of this function goes here
%   Detailed explanation goes here

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:54
    
    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./54;
V1=V1./54;
%%% Because if the jtlag is 3 we have a set of zeros each 3 time steps in
%%% the uv_total.
U1=U1*1.5;
V1=V1*1.5;
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:81
    
    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./81;
V1=V1./81;
%%% Because if the jtlag is 3 we have a set of zeros each 3 time steps in
%%% the uv_total.
U1=U1*2.25;
V1=V1*2.25;
vmean(:,4)=V1;


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