diff --git a/src/mesh/FEM_utilities.f90 b/src/mesh/FEM_utilities.f90
index c0e37001e..c152a2533 100644
--- a/src/mesh/FEM_utilities.f90
+++ b/src/mesh/FEM_utilities.f90
@@ -192,7 +192,6 @@ subroutine utilities_constitutiveResponse(timeinc,P_av,forwardData)
    debug_reset, &
    debug_info
  use math, only: &
-   math_transpose33, &
    math_rotate_forward33, &
    math_det33
  use FEsolving, only: &