D0_massman function

Calculate D0