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