diff --git a/src/constitutive.f90 b/src/constitutive.f90
index 2ec947f2b..6a73cc071 100644
--- a/src/constitutive.f90
+++ b/src/constitutive.f90
@@ -274,12 +274,10 @@ module constitutive
         dPhiDot_dPhi
     end subroutine constitutive_damage_getRateAndItsTangents
 
-    module subroutine constitutive_thermal_getRate(TDot, T,ip,el)
+    module subroutine constitutive_thermal_getRate(TDot, ip,el)
       integer, intent(in) :: &
         ip, &                                                                                       !< integration point number
         el                                                                                          !< element number
-      real(pReal), intent(in) :: &
-        T
       real(pReal), intent(out) :: &
         TDot
     end subroutine constitutive_thermal_getRate
diff --git a/src/constitutive_thermal.f90 b/src/constitutive_thermal.f90
index cabd03c1e..8a18f030e 100644
--- a/src/constitutive_thermal.f90
+++ b/src/constitutive_thermal.f90
@@ -127,13 +127,11 @@ end subroutine thermal_init
 !----------------------------------------------------------------------------------------------
 !< @brief calculates thermal dissipation rate
 !----------------------------------------------------------------------------------------------
-module subroutine constitutive_thermal_getRate(TDot, T, ip, el)
+module subroutine constitutive_thermal_getRate(TDot, ip, el)
 
   integer, intent(in) :: &
     ip, &                                                                                           !< integration point number
     el                                                                                              !< element number
-  real(pReal), intent(in) :: &
-    T                                                                                             !< plastic velocity gradient
   real(pReal), intent(out) :: &
     TDot
 
diff --git a/src/grid/grid_thermal_spectral.f90 b/src/grid/grid_thermal_spectral.f90
index c0a1b0f40..9ed0b24c1 100644
--- a/src/grid/grid_thermal_spectral.f90
+++ b/src/grid/grid_thermal_spectral.f90
@@ -278,7 +278,7 @@ subroutine formResidual(in,x_scal,f_scal,dummy,ierr)
   ce = 0
   do k = 1, grid3;  do j = 1, grid(2);  do i = 1,grid(1)
     ce = ce + 1
-    call thermal_conduction_getSource(Tdot, T_current(i,j,k), 1, ce)
+    call thermal_conduction_getSource(Tdot, 1,ce)
     scalarField_real(i,j,k) = params%timeinc*(scalarField_real(i,j,k) + Tdot) &
                             + thermal_conduction_getMassDensity (1,ce)* &
                               thermal_conduction_getSpecificHeat(1,ce)*(T_lastInc(i,j,k)  - &
diff --git a/src/thermal_conduction.f90 b/src/thermal_conduction.f90
index 79fe0d6cd..1711c5512 100644
--- a/src/thermal_conduction.f90
+++ b/src/thermal_conduction.f90
@@ -91,13 +91,11 @@ end subroutine thermal_conduction_init
 !--------------------------------------------------------------------------------------------------
 !> @brief return heat generation rate
 !--------------------------------------------------------------------------------------------------
-subroutine thermal_conduction_getSource(Tdot, T,ip,el)
+subroutine thermal_conduction_getSource(Tdot, ip,el)
 
   integer, intent(in) :: &
     ip, &                                                                                           !< integration point number
     el                                                                                              !< element number
-  real(pReal), intent(in) :: &
-    T
   real(pReal), intent(out) :: &
     Tdot
 
@@ -105,7 +103,7 @@ subroutine thermal_conduction_getSource(Tdot, T,ip,el)
     homog
 
   homog = material_homogenizationAt(el)
-  call constitutive_thermal_getRate(TDot, T,ip,el)
+  call constitutive_thermal_getRate(TDot, ip,el)
 
   Tdot = Tdot/real(homogenization_Nconstituents(homog),pReal)