From f845b9dc87d42ac8c97e96622a9132e7fc75d453 Mon Sep 17 00:00:00 2001 From: Christoph Kords Date: Fri, 6 Mar 2009 12:29:47 +0000 Subject: [PATCH] corrected indices for calculation of dLp_dTstar --- trunk/constitutive_j2.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/trunk/constitutive_j2.f90 b/trunk/constitutive_j2.f90 index f7ed4ce83..0963c9883 100644 --- a/trunk/constitutive_j2.f90 +++ b/trunk/constitutive_j2.f90 @@ -303,7 +303,7 @@ subroutine constitutive_j2_LpAndItsTangent(Lp,dLp_dTstar,Tstar_v,Temperature,sta forall (k=1:3,l=1:3,m=1:3,n=1:3) & dLp_dTstar3333(k,l,m,n) = Tstar33(k,l)*Tstar33(m,n) * (constitutive_j2_n(matID)-1.0_pReal)/squarenorm_Tstar forall (k=1:3,l=1:3) & - dLp_dTstar3333(k,k,l,l) = dLp_dTstar3333(k,k,l,l) + 1.0_pReal + dLp_dTstar3333(k,l,k,l) = dLp_dTstar3333(k,l,k,l) + 1.0_pReal dLp_dTstar = math_Plain3333to99(factor * dLp_dTstar3333) return