reinke_tsep Function

  function reinke_tsep(bt, flh, qstar, rmajor, eps, fgw, kappa, lhat)

    !! Function for calculating upstream temperature(keV) in Reinke model
    !! author: H Lux, CCFE/UKAEA
    !! bt      : input real : toroidal field on axis (T)
    !! flh     : input real : fraction of Psep/P_LH
    !! qstar   : input real : safety factor similar to q95 (see #707)
    !! rmajor  : input real : major radius (m)
    !! eps     : input real : inverse aspect ratio
    !! fgw     : input real : ratio of volume averaged density to n_GW
    !! kappa   : input real : elongation
    !! lhat    : input real : connection length factor
    !! This function calculates the upstream temperature in the
    !! divertor/SoL model used for the Reinke citerion.
    !! Issue #707
    !! M.L. Reinke 2017 Nucl. Fusion 57 034004

    implicit none
    real(dp) :: reinke_tsep
    real(dp) :: bt, flh, qstar, rmajor, eps, fgw, kappa, lhat
    real(dp), parameter :: kappa_0 = 2D3 !Stangeby W/m/eV^(7/2)

    reinke_tsep = bt**0.72 * flh**0.29 * fgw**0.21 * qstar**0.08 * rmajor**0.33
    !reinke_tsep = bt**0.72 * flh**0.2857 * fgw**0.2057 * qstar**0.08 * rmajor**0.3314

    reinke_tsep = reinke_tsep * eps**0.15 * (1.D0 + kappa**2.)**0.34
    !reinke_tsep = reinke_tsep * eps**0.1486 * (1.D0 + kappa**2.)**0.34

    reinke_tsep = reinke_tsep * lhat**0.29 * kappa_0 **(-0.29) * 0.285D0
    !reinke_tsep = reinke_tsep * lhat**0.2857 * kappa_0 **(-0.2857) * 0.285D0

  end function reinke_tsep