diff --git a/source/fortran/initial.f90 b/source/fortran/initial.f90 index 2adc158e93..ad86dc1e20 100755 --- a/source/fortran/initial.f90 +++ b/source/fortran/initial.f90 @@ -347,6 +347,12 @@ 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) should be removed from the input file, or set to the default value 0.1D0.' + stop 1 + end if + ! Impurity fractions do imp = 1,nimp impurity_arr_frac(imp) = fimp(imp)