add Mb Me as Mfactor_min and Mfactor_max varying along axis, at z=0, Mb, z=zmax Me