From fede1dcd09e80e0a61baeb1b82f832b753668174 Mon Sep 17 00:00:00 2001
From: Martin Diehl <mail@martin-diehl.net>
Date: Sun, 24 Jan 2021 11:16:17 +0100
Subject: [PATCH] dot_T is needed

---
 src/constitutive_thermal.f90 | 7 +++----
 1 file changed, 3 insertions(+), 4 deletions(-)

diff --git a/src/constitutive_thermal.f90 b/src/constitutive_thermal.f90
index 2746ff889..610b0b4a6 100644
--- a/src/constitutive_thermal.f90
+++ b/src/constitutive_thermal.f90
@@ -10,9 +10,9 @@ submodule(constitutive) constitutive_thermal
   end enum
 
   type :: tDataContainer
-    real(pReal), dimension(:), allocatable :: T
+    real(pReal), dimension(:), allocatable :: T, dot_T
   end type tDataContainer
-  integer(kind(THERMAL_UNDEFINED_ID)),     dimension(:,:), allocatable :: &
+  integer(kind(THERMAL_UNDEFINED_ID)),  dimension(:,:), allocatable :: &
     thermal_source
 
   type(tDataContainer), dimension(:), allocatable :: current
@@ -94,6 +94,7 @@ module subroutine thermal_init(phases)
     Nconstituents = count(material_phaseAt == ph) * discretization_nIPs
 
     allocate(current(ph)%T(Nconstituents),source=300.0_pReal)
+    allocate(current(ph)%dot_T(Nconstituents),source=0.0_pReal)
     phase => phases%get(ph)
     if(phase%contains('thermal')) then
       thermal => phase%get('thermal')
@@ -147,13 +148,11 @@ module subroutine constitutive_thermal_getRate(TDot, ip, el)
   integer :: &
     ph, &
     homog, &
-    instance, &
     me, &
     so, &
     co
 
    homog  = material_homogenizationAt(el)
-   instance = thermal_typeInstance(homog)
 
   TDot = 0.0_pReal
   do co = 1, homogenization_Nconstituents(homog)