wrapper for homogenization of temperature rate

src/homogenization.f90

@@ -75,6 +75,10 @@ module homogenization
       integer,     intent(in) :: ce
     end subroutine thermal_partition
 
+    module subroutine thermal_homogenize(ip,el)
+      integer, intent(in) :: ip,el
+    end subroutine thermal_homogenize
+
     module subroutine mech_homogenize(dt,ip,el)
      real(pReal), intent(in) :: dt
      integer, intent(in) :: &
@@ -272,6 +276,7 @@ subroutine materialpoint_stressAndItsTangent(dt,FEsolving_execIP,FEsolving_execE
               print*, ' Integration point ', ip,' at element ', el, ' terminally ill'
             terminallyIll = .true.                                                                  ! ...and kills all others
          endif
+         call thermal_homogenize(ip,el)
         enddo
       enddo
     enddo

src/homogenization_thermal.f90

@@ -21,7 +21,7 @@ end subroutine thermal_init
 
 
 !--------------------------------------------------------------------------------------------------
-!> @brief Partition T onto the individual constituents.
+!> @brief Partition temperature onto the individual constituents.
 !--------------------------------------------------------------------------------------------------
 module subroutine thermal_partition(T,ce)
 
@@ -37,4 +37,16 @@ module subroutine thermal_partition(T,ce)
 end subroutine thermal_partition
 
 
+!--------------------------------------------------------------------------------------------------
+!> @brief Homogenize temperature rates
+!--------------------------------------------------------------------------------------------------
+module subroutine thermal_homogenize(ip,el)
+
+  integer, intent(in) :: ip,el
+
+  call constitutive_thermal_getRate(homogenization_dot_T((el-1)*discretization_nIPs+ip), ip,el)
+
+end subroutine thermal_homogenize
+
+
 end submodule homogenization_thermal