Routine to initialise the variables relevant to stellarators author: P J Knight, CCFE, Culham Science Centre author: F Warmer, IPP Greifswald None This routine initialises the variables relevant to stellarators. Many of these may override the values set in routine initial. !
subroutine stinit
!! Routine to initialise the variables relevant to stellarators
!! author: P J Knight, CCFE, Culham Science Centre
!! author: F Warmer, IPP Greifswald
!! None
!! This routine initialises the variables relevant to stellarators.
!! Many of these may override the values set in routine
!! <A HREF="initial.html">initial</A>.
!! !
! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
use build_variables, only: gapoh, iohcl, ohcth, tfootfi
use current_drive_variables, only: irfcd
use pfcoil_variables, only: ohhghf
use physics_variables, only: aspect, dnbeta, kappa, kappa95, q, rmajor, &
triang, hfac, tauscl
use numerics, only: boundl, boundu
use stellarator_variables, only: istell
use tfcoil_variables, only: n_tf
use times_variables, only: tburn, tcycle, tdown, tdwell, t_fusion_ramp, tohs, &
tpulse, tqnch, tramp
use global_variables, only: icase
use constants, only: pi, rmu0, nout
implicit none
! Arguments
! Local variables
!real(dp) :: fsum
! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! This routine is called before (!!!) the input file. put everything that depends on the input file in stcaller
if (istell == 0) return
boundu(1) = 40.0D0 ! allow higher aspect ratio
! Numerics quantities
!boundl(1) = 5.0D0
!boundu(3) = 30.0D0
!boundu(29) = 20.0D0
! These lines switch off tokamak specifics (solenoid, pf coils, pulses etc.).
! Are they still up to date? (26/07/22 JL)
! Build quantities
ohcth = 0.0D0
iohcl = 0
ohhghf = 0.0D0
gapoh = 0.0D0
tfootfi = 1.0D0
! Physics quantities
dnbeta = 0.0D0
kappa95 = 1.0D0
triang = 0.0D0
q = 1.03D0
! Turn off current drive
irfcd = 0
! Times for different phases
tramp = 0.0D0
tohs = 0.0D0
tburn = 3.15576D7 ! one year
tqnch = 0.0D0
tpulse = tohs + t_fusion_ramp + tburn + tqnch
tdown = tramp + tohs + tqnch + tdwell
tcycle = tramp + tohs + t_fusion_ramp + tburn + tqnch + tdwell
end subroutine stinit