Modified the calculation of PIVIT_MAT for porous media. Now it is only calculated once or if the mesh changes. Also in ASSEMB_FORCE_CTY added an if to avoid doing some costly loops that were not adding anything for porous media. Test cases pass.