diff --git a/code/constitutive.f90 b/code/constitutive.f90
index e5608434e..41ab511b8 100644
--- a/code/constitutive.f90
+++ b/code/constitutive.f90
@@ -771,6 +771,7 @@ function constitutive_getConductionThermal(ipc, ip, el)
    pReal
  use material, only: &
    mappingHomogenization, &
+   material_phase, &
    fieldThermal, &
    field_thermal_type, &
    FIELD_DAMAGE_LOCAL_ID, &