From 175bab4073d9a223e8b0f995d87256637dadc9c6 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Thu, 24 May 2018 00:25:29 +0200 Subject: [PATCH] 1e-6 seems to be a more reasoable absolute tolerance for accumulated shear --- 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 cb14b4cc2..23906bee3 100644 --- a/src/plastic_dislotwin.f90 +++ b/src/plastic_dislotwin.f90 @@ -1138,7 +1138,7 @@ subroutine plastic_dislotwin_init(fileUnit) endIndex=endIndex+ns state(instance)%accshear_slip=>plasticState(phase)%state(startIndex:endIndex,:) dotState(instance)%accshear_slip=>plasticState(phase)%dotState(startIndex:endIndex,:) - plasticState(phase)%aTolState(startIndex:endIndex) = 1e6_pReal + plasticState(phase)%aTolState(startIndex:endIndex) = 1e-6_pReal startIndex=endIndex+1 endIndex=endIndex+nt @@ -1150,7 +1150,7 @@ subroutine plastic_dislotwin_init(fileUnit) endIndex=endIndex+nt state(instance)%accshear_twin=>plasticState(phase)%state(startIndex:endIndex,:) dotState(instance)%accshear_twin=>plasticState(phase)%dotState(startIndex:endIndex,:) - plasticState(phase)%aTolState(startIndex:endIndex) = 1e6_pReal + plasticState(phase)%aTolState(startIndex:endIndex) = 1e-6_pReal startIndex=endIndex+1 endIndex=endIndex+nr