From 62f254c9de63e8b0015e3bf34ed4555e1e368354 Mon Sep 17 00:00:00 2001 From: Haiming Zhang Date: Mon, 8 Dec 2014 09:09:54 +0000 Subject: [PATCH] age results and write restart data for vacancy concentration --- code/CPFEM.f90 | 2 ++ 1 file changed, 2 insertions(+) diff --git a/code/CPFEM.f90 b/code/CPFEM.f90 index 820a29c66..4a33314f5 100644 --- a/code/CPFEM.f90 +++ b/code/CPFEM.f90 @@ -304,6 +304,7 @@ subroutine CPFEM_general(mode, ffn, ffn1, temperature, dt, elFE, ip) damageState, & homogState, & thermalState, & + vacancyState,& mappingConstitutive, & material_phase, & phase_plasticity, & @@ -405,6 +406,7 @@ subroutine CPFEM_general(mode, ffn, ffn1, temperature, dt, elFE, ip) forall ( i = 1:size(plasticState)) plasticState(i)%state0= plasticState(i)%state ! copy state in this lenghty way because: A component cannot be an array if the encompassing structure is an array forall ( i = 1:size(damageState)) damageState(i)%state0 = damageState(i)%state ! copy state in this lenghty way because: A component cannot be an array if the encompassing structure is an array forall ( i = 1:size(thermalState)) thermalState(i)%state0= thermalState(i)%state ! copy state in this lenghty way because: A component cannot be an array if the encompassing structure is an array + forall ( i = 1:size(vacancyState)) vacancyState(i)%state0= vacancyState(i)%state ! copy state in this lenghty way because: A component cannot be an array if the encompassing structure is an array if (iand(debug_level(debug_CPFEM), debug_levelBasic) /= 0_pInt) then write(6,'(a)') '<< CPFEM >> aging states' if (debug_e <= mesh_NcpElems .and. debug_i <= mesh_maxNips) then