diff --git a/code/DAMASK_abaqus_exp.f b/code/DAMASK_abaqus_exp.f
index 37ee9db76..575a3308e 100644
--- a/code/DAMASK_abaqus_exp.f
+++ b/code/DAMASK_abaqus_exp.f
@@ -95,6 +95,13 @@ end module DAMASK_interface
 #include "homogenization_RGC.f90"
 #include "homogenization.f90"
 #include "CPFEM.f90"
+#include "constitutive_damage.f90"
+#include "damage_none.f90"
+#include "damage_gradient.f90"
+#include "constitutive_thermal.f90"
+#include "thermal_none.f90"
+#include "thermal_conduction.f90"
+#include "thermal_adiabatic.f90"
 
 subroutine vumat(nBlock, nDir, nshr, nStateV, nFieldV, nProps, lAnneal, &
                  stepTime, totalTime, dt, cmName, coordMp, charLength, &
diff --git a/code/DAMASK_abaqus_std.f b/code/DAMASK_abaqus_std.f
index 0c6435f30..08dda9088 100644
--- a/code/DAMASK_abaqus_std.f
+++ b/code/DAMASK_abaqus_std.f
@@ -95,6 +95,14 @@ end module DAMASK_interface
 #include "homogenization_RGC.f90"
 #include "homogenization.f90"
 #include "CPFEM.f90"
+#include "constitutive_damage.f90"
+#include "damage_none.f90"
+#include "damage_gradient.f90"
+#include "constitutive_thermal.f90"
+#include "thermal_none.f90"
+#include "thermal_conduction.f90"
+#include "thermal_adiabatic.f90"
+
  
 subroutine UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,&
                 RPL,DDSDDT,DRPLDE,DRPLDT,STRAN,DSTRAN,&
diff --git a/code/DAMASK_marc.f90 b/code/DAMASK_marc.f90
index e4b1c57e6..0b1f61e92 100644
--- a/code/DAMASK_marc.f90
+++ b/code/DAMASK_marc.f90
@@ -124,6 +124,13 @@ end module DAMASK_interface
 #include "homogenization_RGC.f90"
 #include "homogenization.f90"
 #include "CPFEM.f90"
+#include "constitutive_damage.f90"
+#include "damage_none.f90"
+#include "damage_gradient.f90"
+#include "constitutive_thermal.f90"
+#include "thermal_none.f90"
+#include "thermal_conduction.f90"
+#include "thermal_adiabatic.f90"
 
 
 !--------------------------------------------------------------------------------------------------