diff --git a/src/constitutive.f90 b/src/constitutive.f90
index 667a23127..36400fca7 100644
--- a/src/constitutive.f90
+++ b/src/constitutive.f90
@@ -202,11 +202,17 @@ module constitutive
       real(pReal) :: T
     end function constitutive_thermal_T
 
+
     module subroutine constitutive_mech_setF(F,co,ip,el)
       real(pReal), dimension(3,3), intent(in) :: F
       integer, intent(in) :: co, ip, el
     end subroutine constitutive_mech_setF
 
+    module subroutine constitutive_thermal_setT(T,co,ip,el)
+      real(pReal), intent(in) :: T
+      integer, intent(in) :: co, ip, el
+    end subroutine constitutive_thermal_setT
+
 ! == cleaned:end ===================================================================================
 
     module function crystallite_stress(dt,co,ip,el) result(converged_)
@@ -414,6 +420,7 @@ module constitutive
     constitutive_restartRead, &
     integrateSourceState, &
     constitutive_mech_setF, &
+    constitutive_thermal_setT, &
     constitutive_mech_getP, &
     constitutive_mech_getLp, &
     constitutive_mech_getF, &
diff --git a/src/constitutive_thermal.f90 b/src/constitutive_thermal.f90
index 1a05b983f..01d517124 100644
--- a/src/constitutive_thermal.f90
+++ b/src/constitutive_thermal.f90
@@ -139,4 +139,11 @@ module function constitutive_thermal_T(co,ip,el) result(T)
 end function constitutive_thermal_T
 
 
+! setter for homogenization
+module subroutine constitutive_thermal_setT(T,co,ip,el)
+  real(pReal), intent(in) :: T
+  integer, intent(in) :: co, ip, el
+end subroutine constitutive_thermal_setT
+
+
 end submodule constitutive_thermal
diff --git a/src/homogenization.f90 b/src/homogenization.f90
index 686bb9885..df7369096 100644
--- a/src/homogenization.f90
+++ b/src/homogenization.f90
@@ -27,6 +27,8 @@ module homogenization
 
 !--------------------------------------------------------------------------------------------------
 ! General variables for the homogenization at a  material point
+  real(pReal),   dimension(:),         allocatable, public :: &
+    homogenization_T
   real(pReal),   dimension(:,:,:),     allocatable, public :: &
     homogenization_F0, &                                                                            !< def grad of IP at start of FE increment
     homogenization_F                                                                                !< def grad of IP to be reached at end of FE increment
@@ -56,6 +58,9 @@ module homogenization
         num_homog                                                                                   !< pointer to mechanical homogenization numerics data
     end subroutine mech_init
 
+    module subroutine thermal_init
+    end subroutine thermal_init
+
     module subroutine mech_partition(subF,ip,el)
       real(pReal), intent(in), dimension(3,3) :: &
         subF
diff --git a/src/homogenization_thermal.f90 b/src/homogenization_thermal.f90
new file mode 100644
index 000000000..59e7357b6
--- /dev/null
+++ b/src/homogenization_thermal.f90
@@ -0,0 +1,37 @@
+!--------------------------------------------------------------------------------------------------
+!> @author Martin Diehl, KU Leuven
+!--------------------------------------------------------------------------------------------------
+submodule(homogenization) homogenization_thermal
+
+
+contains
+
+!--------------------------------------------------------------------------------------------------
+!> @brief Allocate variables and set parameters.
+!--------------------------------------------------------------------------------------------------
+module subroutine thermal_init()
+
+  print'(/,a)',   ' <<<+-  homogenization_thermal init  -+>>>'
+
+  allocate(homogenization_T(discretization_nIPs*discretization_Nelems), source=0.0_pReal)
+
+end subroutine thermal_init
+
+
+!--------------------------------------------------------------------------------------------------
+!> @brief Partition T onto the individual constituents.
+!--------------------------------------------------------------------------------------------------
+module subroutine thermal_partition(T,ip,el)
+
+  real(pReal), intent(in) :: T
+  integer,     intent(in) :: &
+    ip, &                                                                                           !< integration point
+    el                                                                                              !< element number
+
+
+  call constitutive_thermal_setT(T,1,ip,el)
+
+end subroutine thermal_partition
+
+
+end submodule homogenization_thermal