From d40e4a3800a14c80c110523cd6934cac81243a8f Mon Sep 17 00:00:00 2001 From: Vitesh Shah Date: Tue, 30 Jul 2019 10:53:19 +0200 Subject: [PATCH] SFE at 0K can be negative as well --- src/plastic_dislotwin.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/plastic_dislotwin.f90 b/src/plastic_dislotwin.f90 index 3553630dc..57be24480 100644 --- a/src/plastic_dislotwin.f90 +++ b/src/plastic_dislotwin.f90 @@ -807,7 +807,7 @@ subroutine plastic_dislotwin_dotState(Mp,T,instance,of) else sigma_cl = norm2(matmul(Mp,prm%n0_sl(1:3,i))) ! ToDo: MD: correct? - if (prm%SFE_0K > 0.0_pReal) then ! ToDo: MD: I'm not really sure if this is correct. Maybe Gamma(0K) = 0 + if (prm%SFE_0K == 0.0_pReal) then ! ToDo: MD: I'm not really sure if this is correct. Maybe Gamma(0K) = 0 b_d = 24.0_pReal*PI*(1.0_pReal - prm%nu)/(2.0_pReal + prm%nu)* Gamma/(prm%mu*prm%b_sl(i)) else b_d = 1.0_pReal