From 7a4415b28a21dd32aa740ebaa7f0a1d3c3864024 Mon Sep 17 00:00:00 2001 From: Su Leen Wong Date: Wed, 8 Jul 2015 15:37:16 +0000 Subject: [PATCH] Sanity checks for SFE --- code/plastic_dislotwin.f90 | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/code/plastic_dislotwin.f90 b/code/plastic_dislotwin.f90 index 844a884a2..ae99f791f 100644 --- a/code/plastic_dislotwin.f90 +++ b/code/plastic_dislotwin.f90 @@ -693,6 +693,10 @@ subroutine plastic_dislotwin_init(fileUnit) call IO_error(211_pInt,el=instance,ext_msg='aTolTwinFrac ('//PLASTICITY_DISLOTWIN_label//')') endif if (sum(plastic_dislotwin_Ntrans(:,instance)) > 0_pInt) then + if (abs(plastic_dislotwin_SFE_0K(instance)) <= tiny(0.0_pReal) .and. & + abs(plastic_dislotwin_dSFE_dT(instance)) <= tiny(0.0_pReal) .and. & + lattice_structure(phase) == LATTICE_fcc_ID) & + call IO_error(211_pInt,el=instance,ext_msg='SFE0K ('//PLASTICITY_DISLOTWIN_label//')') if (plastic_dislotwin_aTolTransFrac(instance) <= 0.0_pReal) & call IO_error(211_pInt,el=instance,ext_msg='aTolTransFrac ('//PLASTICITY_DISLOTWIN_label//')') if (plastic_dislotwin_Cdwp(instance) < 0.0_pReal) &