constraint_eqn_080 Subroutine

public subroutine constraint_eqn_080(tmp_cc, tmp_con, tmp_err, tmp_symbol, tmp_units)

Equation for pdivt lower limit author: J Morris, Culham Science Centre args : output structure : residual error; constraint value; residual error in physical units; output string; units string Lower limit pdivt #=# physics #=#=# fpdivlim, pdivt Logic change during pre-factoring: err, symbol, units will be assigned only if present. fpdivlim : input : F-value for lower limit on pdivt (cons. 80, itvar 153) pdivtlim : input : Minimum power crossing separatrix pdivt [MW] pdivt : input : Power crossing separatrix [MW]

Arguments

Type IntentOptional AttributesName
real(kind=dp), intent(out) :: tmp_cc
real(kind=dp), intent(out) :: tmp_con
real(kind=dp), intent(out) :: tmp_err
character(len=1), intent(out) :: tmp_symbol
character(len=10), intent(out) :: tmp_units

Contents

Source Code


Source Code

   subroutine constraint_eqn_080(tmp_cc, tmp_con, tmp_err, tmp_symbol, tmp_units)
      !! Equation for pdivt lower limit
      !! author: J Morris, Culham Science Centre
      !! args : output structure : residual error; constraint value; residual error in physical units;
      !! output string; units string
      !! Lower limit pdivt
      !! #=# physics
      !! #=#=# fpdivlim, pdivt
      !! Logic change during pre-factoring: err, symbol, units will be assigned only if present.
      !! fpdivlim : input : F-value for lower limit on pdivt (cons. 80, itvar 153)
      !! pdivtlim : input : Minimum power crossing separatrix pdivt [MW]
      !! pdivt : input : Power crossing separatrix [MW]
      use physics_variables, only: fpdivlim, pdivt
      use constraint_variables, only : pdivtlim
      implicit none

            real(dp), intent(out) :: tmp_cc
      real(dp), intent(out) :: tmp_con
      real(dp), intent(out) :: tmp_err
      character(len=1), intent(out) :: tmp_symbol
      character(len=10), intent(out) :: tmp_units
      tmp_cc     = 1.0D0 - fpdivlim * pdivt / pdivtlim
      tmp_con    = pdivtlim
      tmp_err    = pdivt * tmp_cc
      tmp_symbol = '>'
      tmp_units  = 'MW'

   end subroutine constraint_eqn_080