Changed initial conditions to included totalmass0. Output is density0.