From d2ae859b1c982bdeca73d3fdf75fa9c9e5e1c379 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Sun, 17 Jan 2021 15:19:36 +0100 Subject: [PATCH] use new structure --- src/DAMASK_marc.f90 | 2 +- src/grid/grid_thermal_spectral.f90 | 3 +++ 2 files changed, 4 insertions(+), 1 deletion(-) diff --git a/src/DAMASK_marc.f90 b/src/DAMASK_marc.f90 index 0f9d37ddb..9464595b1 100644 --- a/src/DAMASK_marc.f90 +++ b/src/DAMASK_marc.f90 @@ -365,7 +365,7 @@ subroutine flux(f,ts,n,time) f f(2) = 0.0_pReal - call thermal_conduction_getSource(f(1), ts(3), n(3),mesh_FEM2DAMASK_elem(n(1))) + call thermal_conduction_getSource(f(1), n(3),mesh_FEM2DAMASK_elem(n(1))) end subroutine flux diff --git a/src/grid/grid_thermal_spectral.f90 b/src/grid/grid_thermal_spectral.f90 index 9ed0b24c1..d164f11e7 100644 --- a/src/grid/grid_thermal_spectral.f90 +++ b/src/grid/grid_thermal_spectral.f90 @@ -16,6 +16,7 @@ module grid_thermal_spectral use spectral_utilities use discretization_grid use thermal_conduction + use homogenization use YAML_types use config use material @@ -190,6 +191,7 @@ function grid_thermal_spectral_solution(timeinc) result(solution) call thermal_conduction_putTemperatureAndItsRate(T_current(i,j,k), & (T_current(i,j,k)-T_lastInc(i,j,k))/params%timeinc, & 1,ce) + homogenization_T(ce) = T_current(i,j,k) enddo; enddo; enddo call VecMin(solution_vec,devNull,T_min,ierr); CHKERRQ(ierr) @@ -231,6 +233,7 @@ subroutine grid_thermal_spectral_forward(cutBack) (T_current(i,j,k) - & T_lastInc(i,j,k))/params%timeinc, & 1,ce) + homogenization_T(ce) = T_current(i,j,k) enddo; enddo; enddo else T_lastInc = T_current