From 2505c02c75c1f2af3112b18410945200ee223ab8 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Mon, 11 Mar 2019 23:55:04 +0100 Subject: [PATCH] bugfix --- src/plastic_nonlocal.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/plastic_nonlocal.f90 b/src/plastic_nonlocal.f90 index 3dbb5adc2..f35950620 100644 --- a/src/plastic_nonlocal.f90 +++ b/src/plastic_nonlocal.f90 @@ -991,8 +991,6 @@ forall (s = 1_pInt:ns) & !*** coefficients are corrected for the line tension effect !*** (see Kubin,Devincre,Hoc; 2008; Modeling dislocation storage rates and mean free paths in face-centered cubic crystals) -myInteractionMatrix = 0.0_pReal -myInteractionMatrix(1:ns,1:ns) = prm%interactionSlipSlip(1:ns,1:ns) if (lattice_structure(ph) == LATTICE_bcc_ID .or. lattice_structure(ph) == LATTICE_fcc_ID) then ! only fcc and bcc do s = 1_pInt,ns myRhoForest = max(rhoForest(s),prm%significantRho) @@ -1000,8 +998,10 @@ if (lattice_structure(ph) == LATTICE_bcc_ID .or. lattice_structure(ph) == LATTI + prm%linetensionEffect & * log(0.35_pReal * prm%burgers(s) * sqrt(myRhoForest)) & / log(0.35_pReal * prm%burgers(s) * 1e6_pReal)) ** 2.0_pReal - myInteractionMatrix(s,1:ns) = correction * myInteractionMatrix(1:ns,s) + myInteractionMatrix(1:ns,s) = correction * prm%interactionSlipSlip(1:ns,s) enddo +else + myInteractionMatrix = prm%interactionSlipSlip endif forall (s = 1_pInt:ns) & dst%tau_threshold(s,of) = prm%mu * prm%burgers(s) &