Put this routine back in, since it is needed to compute element stiffness matrix.
Put this routine back in, since it is needed to compute element
stiffness matrix.


geodynamics/lithomop,geodynamics/lithomop,geodynamics/lithomop,geodynamics/lithomop