function data=setnans(data,fax)
% data=SETNANS(data,fax)
%
% INPUT:
%
% data     The data matrix
% fax      faxth fraction of max(abs(data)) below which data are replaced
%          by NaN [default: 100]
% 
% OUTPUT:
%
% data    The new data matrix
%
% Last modified by fjsimons-at-alum.mit.edu, 03/15/2011

defval('fax',100)
% Check for -Inf or +Inf and get rid of those 
isi=isinf(data);
if sum(isi(:))>0
  fprintf(1,'%i infinite values replaced by 0 and thus by NaN\n',...
          sum(isi(:)))
  data(isinf(data))=0;
end
data(abs(data)<max(abs(data(:)))/fax)=NaN;