diff --git a/documentation/proc-pages/physics-models/plasma_beta/plasma_beta.md b/documentation/proc-pages/physics-models/plasma_beta/plasma_beta.md index c8e8541b..3e262b1f 100644 --- a/documentation/proc-pages/physics-models/plasma_beta/plasma_beta.md +++ b/documentation/proc-pages/physics-models/plasma_beta/plasma_beta.md @@ -191,6 +191,8 @@ The value of `beta_poloidal_max` can be set to the desired maximum poloidal beta This constraint can be activated by stating `icc = 84` in the input file. +The value of `beta_max` can be set to the desired minimum total beta. The scaling value `fbeta_min` can be varied also. + [^1]: N.A. Uckan and ITER Physics Group, 'ITER Physics Design Guidelines: 1989', [^2]: D.J. Ward, 'PROCESS Fast Alpha Pressure', Work File Note F/PL/PJK/PROCESS/CODE/050 diff --git a/source/fortran/constraint_equations.f90 b/source/fortran/constraint_equations.f90 index b35ea75f..c7fe3c76 100755 --- a/source/fortran/constraint_equations.f90 +++ b/source/fortran/constraint_equations.f90 @@ -3137,9 +3137,8 @@ subroutine constraint_eqn_084(tmp_cc, tmp_con, tmp_err, tmp_symbol, tmp_units) !! fbeta_min : input real : f-value for constraint beta-beta_fast_alpha > beta_min !! beta_min : input real : Lower limit for beta !! beta : input real : plasma beta - !! beta_fast_alpha : input real : Alpha particle beta - use physics_variables, only: beta_min, beta, beta_fast_alpha + use physics_variables, only: beta_min, beta use constraint_variables, only: fbeta_min implicit none real(dp), intent(out) :: tmp_cc @@ -3149,9 +3148,9 @@ subroutine constraint_eqn_084(tmp_cc, tmp_con, tmp_err, tmp_symbol, tmp_units) character(len=10), intent(out) :: tmp_units - tmp_cc = 1.0D0 - fbeta_min * (beta-beta_fast_alpha)/beta_min + tmp_cc = 1.0D0 - fbeta_min * (beta)/beta_min tmp_con = beta_min * (1.0D0 - tmp_cc) - tmp_err = (beta-beta_fast_alpha) * tmp_cc + tmp_err = (beta) * tmp_cc tmp_symbol = '>' tmp_units = ''