subroutine ComputeDivergenceContinuitySolid implemented in nonlinear.F for the spatial discretization of the solid continuity eqution