module stellarator_variables !! author: S. Muldrew (UKAEA), F. Warmer, J. Lion (IPP Greifswald) !! !! Module containing global variables relating to the stellarator model !! !!### References !! !! - A general stellarator version of the Systems code PROCESS, J. Lion et al, Nucl. Fus. (2021), https://doi.org/10.1088/1741-4326/ac2dbf !! - Stellarator Divertor Model for the Systems Code PROCESS, F. Warmer, 21/06/2013 #ifndef dp use, intrinsic :: iso_fortran_env, only: dp=>real64 #endif implicit none public integer :: istell !! Switch for stellarator option (set via `device.dat`): !! !! - =0 use tokamak model !! - =1 use stellarator model: Helias5 !! - =2 use stellarator model: Helias4 !! - =3 use stellarator model: Helias3 !! - =4 use stellarator model: Wendelstein 7-X with 50 Coils !! - =5 use stellarator model: Wendelstein 7-X with 30 Coils !! - =6 use stellarator model: Use stella_conf.json file (any modulear stellarator, see documentation) real(dp) :: bmn !! relative radial field perturbation real(dp) :: f_asym !! divertor heat load peaking factor real(dp) :: f_rad !! radiated power fraction in SOL real(dp) :: f_w !! island size fraction factor real(dp) :: fdivwet !! wetted fraction of the divertor area real(dp) :: flpitch !! field line pitch (rad) real(dp) :: hportamax !! maximum available area for horizontal ports (m2) real(dp) :: hportpmax !! maximum available poloidal extent for horizontal ports (m) real(dp) :: hporttmax !! maximum available toroidal extent for horizontal ports (m) real(dp) :: iotabar !! rotational transform (reciprocal of tokamak q) for stellarator confinement time scaling laws integer :: isthtr !! Switch for stellarator auxiliary heating method: !! !! - = 1electron cyclotron resonance heating !! - = 2lower hybrid heating !! - = 3neutral beam injection integer :: m_res !! poloidal resonance number (1) real(dp) :: max_gyrotron_frequency !! Maximal available gyrotron frequency (input parameter) (Hz) integer :: n_res !! toroidal resonance number (1) real(dp) :: shear !! magnetic shear, derivative of iotabar (1) real(dp) :: te0_ecrh_achievable !! maximal central electron temperature as achievable by the ECRH, input. (keV) real(dp) :: vportamax !! maximum available area for vertical ports (m2) real(dp) :: vportpmax !! maximum available poloidal extent for vertical ports (m) real(dp) :: vporttmax !! maximum available toroidal extent for vertical ports (m) real(dp) :: powerht_constraint real(dp) :: powerscaling_constraint contains subroutine init_stellarator_variables !! Initialise module variables implicit none istell = 0 bmn = 1.0D-3 f_asym = 1.0D0 f_rad = 0.85D0 f_w = 0.5D0 fdivwet = 0.333333333333333D0 flpitch = 1.0D-3 hportamax = 0.0D0 hportpmax = 0.0D0 hporttmax = 0.0D0 iotabar = 1.0D0 isthtr = 3 m_res = 5 n_res = 5 shear = 0.5D0 vportamax = 0.0D0 vportpmax = 0.0D0 vporttmax = 0.0D0 end subroutine init_stellarator_variables end module stellarator_variables