From 3526c0a757ac71b1486d90aa058018f0bbaa4c28 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Sat, 15 Sep 2018 07:49:14 +0200 Subject: [PATCH] index missing, test failed --- src/plastic_dislotwin.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/plastic_dislotwin.f90 b/src/plastic_dislotwin.f90 index 6b9e5fa79..e79ebd644 100644 --- a/src/plastic_dislotwin.f90 +++ b/src/plastic_dislotwin.f90 @@ -1451,7 +1451,7 @@ subroutine kinetics_slip(prm,stt,mse,of,S,temperature,gdot_slip,dgdot_dtau_slip) integer(pInt) :: i do i = 1_pInt, prm%totalNslip - tau = math_mul33xx33(S,prm%Schmid_slip(1:3,1:3,i)) + tau(i) = math_mul33xx33(S,prm%Schmid_slip(1:3,1:3,i)) enddo significantStress: where((abs(tau)-mse%threshold_stress_slip(:,of)) > tol_math_check) @@ -1512,7 +1512,7 @@ subroutine kinetics_twin(prm,stt,mse,of,S,temperature,gdot_slip,gdot_twin,dgdot_ integer(pInt) :: i,s1,s2 do i = 1_pInt, prm%totalNtwin - tau = math_mul33xx33(S,prm%Schmid_twin(1:3,1:3,i)) + tau(i) = math_mul33xx33(S,prm%Schmid_twin(1:3,1:3,i)) isFCC: if (prm%isFCC) then s1=prm%fcc_twinNucleationSlipPair(1,i) s2=prm%fcc_twinNucleationSlipPair(2,i)