Sanity checks for SFE

This commit is contained in:
Su Leen Wong 2015-07-08 15:37:16 +00:00
parent df5cdb6b01
commit 7a4415b28a
1 changed files with 4 additions and 0 deletions

View File

@ -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) &