function mm = mescd(y,yt)
%
%  mm = mescd(y,yt)    computed mescd.
%
err = max(max( abs(y-yt)./(1+abs(yt)) ));
mm  = max(0,-log10(err));
return
