From 5e39d38a86489becbb9c2fb08bffdff117c305a0 Mon Sep 17 00:00:00 2001
From: zhangc43 <zhangc43@egr.msu.edu>
Date: Mon, 25 Apr 2016 16:02:56 -0400
Subject: [PATCH] update formula for strain path modifier

---
 src/plastic_phenoplus.f90 | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/plastic_phenoplus.f90 b/src/plastic_phenoplus.f90
index 1646f5f41..047abe33b 100644
--- a/src/plastic_phenoplus.f90
+++ b/src/plastic_phenoplus.f90
@@ -980,8 +980,8 @@ subroutine plastic_phenoplus_microstructure(orientation,ipc,ip,el,F0,Fe,Fp,Tstar
           ENDIF
         ENDDO LOOPMYNEIGHBORS
 
-      plasticState(ph)%state(index_kappa+me_slip, of) = 1.0_pReal &
-                                                        + 0.15_pReal*(1.0_pReal + ERF(10.0_pReal*(maxval(ne_mprimes)-0.8_pReal)))
+      plasticState(ph)%state(index_kappa+me_slip, of) = 1.2_pReal &
+                                                        + 0.2_pReal*ERF(4.0_pReal * maxval(ne_mprimes) - 4.0_pReal)
       ENDDO LOOPMYSLIP
 
     ENDIF