Added positivity limiting/mass filling, put initial conditions into its own subroutine.