From 3759024b42147524af94da94b101f41ed19c229f Mon Sep 17 00:00:00 2001 From: Christoph Kords Date: Sun, 19 May 2013 19:23:31 +0000 Subject: [PATCH] forgot to set neighboring_rho (also reverting changes made in rev 2384) --- code/constitutive_nonlocal.f90 | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/code/constitutive_nonlocal.f90 b/code/constitutive_nonlocal.f90 index 4f80c72bc..551ee1d51 100644 --- a/code/constitutive_nonlocal.f90 +++ b/code/constitutive_nonlocal.f90 @@ -2394,11 +2394,15 @@ if (.not. phase_localPlasticity(material_phase(g,ip,el))) then neighboring_v(1_pInt:ns,t) = state0(g,neighboring_ip,neighboring_el)%p((13_pInt+t)*ns+1_pInt:(14_pInt+t)*ns) neighboring_rhoSgl(1_pInt:ns,t) = max(state0(g,neighboring_ip,neighboring_el)%p((t-1_pInt)*ns+1_pInt:t*ns), 0.0_pReal) endforall + forall (t = 5_pInt:8_pInt) & + neighboring_rhoSgl(1_pInt:ns,t) = state0(g,neighboring_ip,neighboring_el)%p((t-1_pInt)*ns+1_pInt:t*ns) else forall (t = 1_pInt:4_pInt) neighboring_v(1_pInt:ns,t) = state(g,neighboring_ip,neighboring_el)%p((13_pInt+t)*ns+1_pInt:(14_pInt+t)*ns) neighboring_rhoSgl(1_pInt:ns,t) = max(state(g,neighboring_ip,neighboring_el)%p((t-1_pInt)*ns+1_pInt:t*ns), 0.0_pReal) endforall + forall (t = 5_pInt:8_pInt) & + neighboring_rhoSgl(1_pInt:ns,t) = state(g,neighboring_ip,neighboring_el)%p((t-1_pInt)*ns+1_pInt:t*ns) endif where (abs(neighboring_rhoSgl) * mesh_ipVolume(neighboring_ip,neighboring_el) ** 0.667_pReal & < constitutive_nonlocal_significantN(myInstance) &