From ad6467638fb5d0b58dacc99aa7129d0ba0af619e Mon Sep 17 00:00:00 2001
From: Pratheek Shanthraj
Date: Thu, 14 Aug 2014 09:59:16 +0000
Subject: [PATCH] minor change to strain energy calculation
---
code/damage_local.f90 | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/code/damage_local.f90 b/code/damage_local.f90
index 72f528298..33fb5f03b 100644
--- a/code/damage_local.f90
+++ b/code/damage_local.f90
@@ -294,7 +294,7 @@ subroutine damage_local_dotState(Tstar_v, Fe, Lp, ipc, ip, el)
if (pressure < 0.0_pReal) stress = stress - pressure*math_I3
trialDamage = min(1.0_pReal, &
(math_trace33(lattice_surfaceEnergy33(1:3,1:3,phase))/3.0_pReal)/ &
- (abs(sum(strain*stress)) + &
+ (sum(abs(strain*stress)) + &
damageState(phase)%state(1,constituent)))
damageState(phase)%dotState(1,constituent) = &