function []=averagesFields(doInitFiles,doWriteMean,rep_out,suff_out,varargin);
global averagesFields_count;

%algorithm:
%	first call (doInitFiles=1 doWriteMean=0) -> init files, arrays and set counter to 0	
%	loop calls (doInitFiles=0 doWriteMean=0) -> increment counter and arrays
%	dump calls (doInitFiles=0 doWriteMean=1) -> write one record, re-initialize arrays and counter
%	final call (doInitFiles=-1 doWriteMean=0) -> close files

%input: varargin is a series of structures

if doInitFiles==1&doWriteMean==0;
averagesFields_count=1;
for icur=1:length(varargin);
eval(['global struct' num2str(icur) '; struct' num2str(icur) '=varargin{icur}; ']);
vars_struct=fieldnames(varargin{icur});
eval(['global fid_struct' num2str(icur) '; fid_struct' num2str(icur) '=zeros(1,length(vars_struct));']); 
for vcur=1:length(vars_struct);
tmp1=cell2mat(vars_struct(vcur)); tmp2=[rep_out tmp1 suff_out];
eval(['fid_struct' num2str(icur) '(vcur)=fopen(''' tmp2 ''',''w'',''b'');']);
end
end

elseif doInitFiles==0&doWriteMean==0;
averagesFields_count=averagesFields_count+1;
for icur=1:length(varargin);
eval(['global struct' num2str(icur) ';']);
vars_struct=fieldnames(varargin{icur});
for vcur=1:length(vars_struct);
tmp1=cell2mat(vars_struct(vcur)); 
eval(['struct' num2str(icur) '.' tmp1 '=struct' num2str(icur) '.' tmp1 '+varargin{icur}.' tmp1 ';']);
end
end

elseif doInitFiles==0&doWriteMean==1;
for icur=1:length(varargin);
eval(['global fid_struct' num2str(icur) ' struct' num2str(icur) ';']);
vars_struct=fieldnames(varargin{icur});
for vcur=1:length(vars_struct);
tmp1=cell2mat(vars_struct(vcur));
eval(['fwrite(fid_struct' num2str(icur) '(vcur),struct' num2str(icur) '.' tmp1 '/averagesFields_count,''float32'');']);
eval(['struct' num2str(icur) '.' tmp1 '(:)=0;']);
end
end
averagesFields_count=0;

elseif doInitFiles==-1&doWriteMean==0;
for icur=1:length(varargin);
eval(['global fid_struct' num2str(icur) ' struct' num2str(icur) ';']);
vars_struct=fieldnames(varargin{icur});
for vcur=1:length(vars_struct); eval(['fclose(fid_struct' num2str(icur) '(vcur));']); end;
eval(['clear struct' num2str(icur) ';']);
end
averagesFields_count=[];

else
fprintf('error in averagesFields \n');
averagesFields_count=[];
end