From e20271b89a2e4182f6c3c6a2f76051ef3ef107cc Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Thu, 8 Jan 2015 13:26:00 +0000 Subject: [PATCH] fix for temperature in spectral solver --- code/DAMASK_spectral_utilities.f90 | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/code/DAMASK_spectral_utilities.f90 b/code/DAMASK_spectral_utilities.f90 index febcfd843..f2347ffb8 100644 --- a/code/DAMASK_spectral_utilities.f90 +++ b/code/DAMASK_spectral_utilities.f90 @@ -804,6 +804,8 @@ subroutine utilities_constitutiveResponse(F_lastInc,F,temperature,timeinc,& materialpoint_F, & materialpoint_P, & materialpoint_dPdF + use thermal_isothermal, only: & + thermal_isothermal_temperature implicit none real(pReal), intent(in) :: temperature !< temperature (no field) @@ -837,7 +839,7 @@ subroutine utilities_constitutiveResponse(F_lastInc,F,temperature,timeinc,& call CPFEM_general(CPFEM_COLLECT,F_lastInc(1:3,1:3,1,1,1),F(1:3,1:3,1,1,1), & temperature,timeinc,1_pInt,1_pInt) - + thermal_isothermal_temperature(:) = temperature materialpoint_F = reshape(F,[3,3,1,product(grid)]) call debug_reset()