From 0033175ab2111ad3bf2246330f203992e9e82a6f Mon Sep 17 00:00:00 2001
From: Dave Whipp <dwhipp@dal.ca>
Date: Tue, 8 Dec 2009 20:59:14 +0000
Subject: [PATCH] Modified e2d output to be sqrt of 2nd invariant
---
src/do_leaf_measurements.f90 | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/src/do_leaf_measurements.f90 b/src/do_leaf_measurements.f90
index acd2c16a..6c0fd08f 100644
--- a/src/do_leaf_measurements.f90
+++ b/src/do_leaf_measurements.f90
@@ -162,7 +162,9 @@ do i=1+iproc,nleaves,nproc
J2d=second_invariant(exx,eyy,ezz,exy,eyz,ezx)
J3d=third_invariant (exx,eyy,ezz,exy,eyz,ezx)
- e2dtemp(i)=J2d
+ !e2dtemp(i)=J2d
+ ! Set e2d as effective stress (sqrt of 2nd invariant) - dwhipp 12/09
+ e2dtemp(i)=sqrt(J2d)
if (J2d/=0.d0) lodetemp(i)=lode_angle (J2d,J3d)/pi
--
GitLab