we need porosity/water_content values on the RHS of the equation. As our mass matrix assembly also contains the porosity values.

This fix will make the curve right for KIN_REACT_GIA module. The fix for REACT_GIA module will come later.
following the revious commit, adding the porosity/water content values to the RHS of the equation in the REACT_GIA module.