source: ETALON/CLIO/CLIO_scripts/makemean.m @ 751

function [Xpos,Ypos, Yerr]=makemean(X,Y);
% Xpos=sort(unique(X));
Xpos=3e+6:5e+4:3.4e+6;

TotM=mean(Y(~isnan(Y)));
for c=1:size(Xpos,2)
    diap=(X==Xpos(c))'&(abs(Y-TotM)<30);
    Ypos(c)=mean(Y(diap));
    Yerr(c)=rms(Y(diap)-Ypos(c));
 
end
%  Ypos(isnan( Ypos))=[];
% Yerr(isnan( Ypos))=[];