From ee1d409bad4d5321c3853dd8ff2f1a2b3f8df438 Mon Sep 17 00:00:00 2001 From: Luv Sharma Date: Fri, 21 Nov 2014 09:03:54 +0000 Subject: [PATCH] critical shear can be set to zero for family not in use. --- code/damage_anisoDuctile.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/code/damage_anisoDuctile.f90 b/code/damage_anisoDuctile.f90 index 8e04e119b..5ba3f87a9 100644 --- a/code/damage_anisoDuctile.f90 +++ b/code/damage_anisoDuctile.f90 @@ -205,7 +205,7 @@ subroutine damage_anisoDuctile_init(fileUnit) damage_anisoDuctile_totalNslip(instance) = sum(damage_anisoDuctile_Nslip(:,instance)) if (damage_anisoDuctile_aTol_damage(instance) >= 1.0e-3_pReal) & damage_anisoDuctile_aTol_damage(instance) = 1.0e-3_pReal ! default absolute tolerance 1e-3 - if (any(damage_anisoDuctile_critAccShear(:,instance) <= 0.0_pReal)) & + if (any(damage_anisoDuctile_critAccShear(:,instance) < 0.0_pReal)) & call IO_error(211_pInt,el=instance,ext_msg='critical_accshear ('//LOCAL_DAMAGE_anisoDuctile_LABEL//')') if (damage_anisoDuctile_N(instance) <= 0.0_pReal) & call IO_error(211_pInt,el=instance,ext_msg='rate_sensitivity_damage ('//LOCAL_DAMAGE_anisoDuctile_LABEL//')')