From 5ee73dee728a9f156d150b27576446b5f05ceb11 Mon Sep 17 00:00:00 2001 From: Christoph Kords Date: Mon, 21 Mar 2011 13:06:11 +0000 Subject: [PATCH] age dotState before resetting it to zero, otherwise the statedamping does not work. --- code/crystallite.f90 | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/code/crystallite.f90 b/code/crystallite.f90 index ec157c4f3..86846db1d 100644 --- a/code/crystallite.f90 +++ b/code/crystallite.f90 @@ -2183,10 +2183,10 @@ endif call constitutive_microstructure(crystallite_Temperature(g,i,e), crystallite_Tstar_v(1:6,g,i,e), & crystallite_Fe, crystallite_Fp, g, i, e) ! update dependent state variables to be consistent with basic states endif + constitutive_previousDotState2(g,i,e)%p = constitutive_previousDotState(g,i,e)%p ! age previous dotState + constitutive_previousDotState(g,i,e)%p = constitutive_dotState(g,i,e)%p ! age previous dotState constitutive_dotState(g,i,e)%p = 0.0_pReal ! reset dotState to zero - constitutive_previousDotState(g,i,e)%p = constitutive_dotState(g,i,e)%p - constitutive_previousDotState2(g,i,e)%p = constitutive_previousDotState(g,i,e)%p - enddo; enddo; enddo + enddo; enddo; enddo !$OMP ENDDO @@ -2295,9 +2295,9 @@ do while (any(crystallite_todo) .and. NiterationState < nState ) call constitutive_microstructure(crystallite_Temperature(g,i,e), crystallite_Tstar_v(1:6,g,i,e), & crystallite_Fe, crystallite_Fp, g, i, e) ! update dependent state variables to be consistent with basic states endif + constitutive_previousDotState2(g,i,e)%p = constitutive_previousDotState(g,i,e)%p ! age previous dotState + constitutive_previousDotState(g,i,e)%p = constitutive_dotState(g,i,e)%p ! age previous dotState constitutive_dotState(g,i,e)%p = 0.0_pReal ! reset dotState to zero - constitutive_previousDotState(g,i,e)%p = constitutive_dotState(g,i,e)%p - constitutive_previousDotState2(g,i,e)%p = constitutive_previousDotState(g,i,e)%p enddo; enddo; enddo !$OMP ENDDO