function s=average_icq5(fl1,fl2,flout)
s1=read_icq4_nomex(fl1);
s2=read_icq4_nomex(fl2);

s=s2;
f=fields(s1);
n=length(f);
for k=1:n
    if ~ischar(f{k}) ,
    str=['s.',f{k},'=(s1.',f{k},'+s2.',f{k},')/2;'];
    disp(str)
    eval(str);
end
end

write_icq4(flout,s)
