Possibly found source of problems

The density, rho, becomes negative in the first step. The radiative
pressure must therefore become greater than the total pressure for this
to happen. Will look for the error later.