From ea5ac59b696899e3dec152732048d64a2e2c14a2 Mon Sep 17 00:00:00 2001 From: ym1906 Date: Wed, 14 Feb 2024 09:56:15 +0000 Subject: [PATCH 1/2] Added check to see if fimp(2) deviates from default value, if it does write a warning message to use ralpne instead. --- source/fortran/initial.f90 | 5 +++++ 1 file changed, 5 insertions(+) diff --git a/source/fortran/initial.f90 b/source/fortran/initial.f90 index 2adc158e93..f6ded44fca 100755 --- a/source/fortran/initial.f90 +++ b/source/fortran/initial.f90 @@ -347,6 +347,11 @@ subroutine check trithtmw = 0.0D0 end if + if (fimp(2) .ne. 0.1D0) then + write(*,*)'The thermal alpha/electron density ratio should be controlled using ralpne (itv 109) and not fimp(2).' + write(*,*)'fimp(2) will be ignored.' + end if + ! Impurity fractions do imp = 1,nimp impurity_arr_frac(imp) = fimp(imp) From 9882f8bb94043a205bbefe13ee91d63c2133aa15 Mon Sep 17 00:00:00 2001 From: ym1906 Date: Wed, 14 Feb 2024 15:42:59 +0000 Subject: [PATCH 2/2] Added stop to enforce the message, mofidified message --- source/fortran/initial.f90 | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/source/fortran/initial.f90 b/source/fortran/initial.f90 index f6ded44fca..ad86dc1e20 100755 --- a/source/fortran/initial.f90 +++ b/source/fortran/initial.f90 @@ -349,7 +349,8 @@ subroutine check if (fimp(2) .ne. 0.1D0) then write(*,*)'The thermal alpha/electron density ratio should be controlled using ralpne (itv 109) and not fimp(2).' - write(*,*)'fimp(2) will be ignored.' + write(*,*)'fimp(2) should be removed from the input file, or set to the default value 0.1D0.' + stop 1 end if ! Impurity fractions