init_physics_variables Subroutine

public subroutine init_physics_variables()

Initialise module variables

Arguments

None

Contents


Source Code

  subroutine init_physics_variables
    !! Initialise module variables
    implicit none

    abeam = 0.0D0
    afuel = 0.0D0
    aion = 0.0D0
    alphaj = 1.0D0
    alphan = 0.25D0
    alphap = 0.0D0
    alpharate = 0.0D0
    alphat = 0.5D0
    aspect = 2.907D0
    beamfus0 = 1.0D0
    beta = 0.042D0
    betaft = 0.0D0
    betalim = 0.0D0
    betalim_lower = 0.0D0
    betanb = 0.0D0
    betap = 0.0D0
    normalised_total_beta = 0.0D0
    betbm0 = 1.5D0
    bp = 0.0D0
    bt = 5.68D0
    btot = 0.0D0
    burnup = 0.0D0
    burnup_in = 0.0D0
    bvert = 0.0D0
    c_beta = 0.5D0
    csawth = 1.0D0
    cvol = 1.0D0
    cwrmax = 1.35D0
    dene = 9.8D19
    deni = 0.0D0
    dlamee = 0.0D0
    dlamie = 0.0D0
    dlimit = 0.0D0
    dnalp = 0.0D0
    dnbeam = 0.0D0
    dnbeam2 = 0.0D0
    dnbeta = 3.5D0
    dnelimt = 0.0D0
    dnitot = 0.0D0
    dnla = 0.0D0
    dnprot = 0.0D0
    dntau = 0.0D0
    dnz = 0.0D0
    epbetmax = 1.38D0
    eps = 0.34399724802D0
    aux_current_fraction = 0.0D0
    inductive_current_fraction = 0.0D0
    falpe = 0.0D0
    falpha = 0.95D0
    falpi = 0.0D0
    fdeut = 0.5D0
    ftar = 1.0D0
    ffwal = 0.92D0
    fgwped = 0.85D0
    fgwsep = 0.50D0
    fhe3 = 0.0D0
    figmer = 0.0D0
    fkzohm = 1.0D0
    fplhsep = 1.0D0
    fpdivlim = 1.0D0
    fne0 = 1.0D0
    ftrit = 0.5D0
    fusionrate = 0.0D0
    fvsbrnni = 1.0D0
    gamma = 0.4D0
    gammaft = 0.0D0
    hfac = 0.0D0
    hfact = 1.0D0
    taumax = 10.0D0
    i_bootstrap_current = 3
    iculbl = 0
    i_plasma_current = 4
    i_diamagnetic_current = 0
    idensl = 7
    idivrt = 2
    ifalphap = 1
    ignite = 0
    iinvqd = 1
    ipedestal = 1
    i_pfirsch_schluter_current = 0
    neped = 4.0D19
    nesep = 3.0D19
    alpha_crit = 0.0D0
    nesep_crit = 0.0D0
    plasma_res_factor = 1.0D0
    rhopedn = 1.0D0
    rhopedt = 1.0D0
    tbeta = 2.0D0
    teped = 1.0D0
    tesep = 0.1D0
    iprofile = 1
    iradloss = 1
    isc = 34
    iscrp = 1
    ishape = 0
    itart = 0
    itartpf = 0
    iwalld = 1
    kappa = 1.792D0
    kappa95 = 1.6D0
    kappaa = 0.0D0
    kappaa_IPB = 0.d0
    ne0 = 0.0D0
    ni0 = 0.0D0
    m_s_limit = 0.3D0
    p0 = 0.0D0
    palppv = 0.0D0
    palpepv = 0.0D0
    palpfwmw = 0.0D0
    palpipv = 0.0D0
    palpmw = 0.0D0
    palpnb = 0.0D0
    pchargemw = 0.0D0
    pchargepv = 0.0D0
    pcoef = 0.0D0
    pinnerzoneradmw = 0.0D0
    pcoreradpv = 0.0D0
    pdd = 0.0D0
    pdhe3 = 0.0D0
    pdivt = 0.0D0
    pdivl = 0.0D0
    pdivu = 0.0D0
    pdivmax = 0.0D0
    pdt = 0.0D0
    pouterzoneradmw = 0.0D0
    pedgeradpv = 0.0D0
    pfuscmw = 0.0D0
    phiint = 0.0D0
    photon_wall = 0.0D0
    piepv = 0.0D0
    plasma_current = 0.0D0
    pneutmw = 0.0D0
    pneutpv = 0.0D0
    pohmmw = 0.0D0
    pohmpv = 0.0D0
    powerht = 0.0D0
    powfmw = 0.0D0
    pperim = 0.0D0
    pradmw = 0.0D0
    pradpv = 0.0D0
    pradsolmw = 0.0D0
    protonrate = 0.0D0
    psolradmw = 0.0D0
    psyncpv = 0.0D0
    ilhthresh = 19
    plhthresh = 0.0D0
    pthrmw = 0.0D0
    ptremw = 0.0D0
    ptrepv = 0.0D0
    ptrimw = 0.0D0
    pscalingmw = 0.0D0
    ptripv = 0.0D0
    q = 3.0D0
    q0 = 1.0D0
    q95 = 0.0D0
    qfuel = 0.0D0
    tauratio = 1.0D0
    qlim = 0.0D0
    qstar = 0.0D0
    rad_fraction_sol = 0.8D0
    rad_fraction_total = 0.0D0
    ralpne = 0.10D0
    protium = 0.0D0
    rli = 0.9D0
    rlp = 0.0D0
    rmajor = 8.14D0
    rminor = 0.0D0
    rnbeam = 0.005D0
    rncne = 0.0D0
    rndfuel = 0.0D0
    rnfene = 0.0D0
    rnone = 0.0D0
    rpfac = 0.0D0
    rplas = 0.0D0
    res_time = 0.0D0
    sarea = 0.0D0
    sareao = 0.0D0
    sf = 0.0D0
    i_single_null = 1
    ssync = 0.6D0
    tauee = 0.0D0
    tauee_in = 0.0D0
    taueff = 0.0D0
    tauei = 0.0D0
    taup = 0.0D0
    te = 12.9D0
    te0 = 0.0D0
    ten = 0.0D0
    ti = 12.9D0
    ti0 = 0.0D0
    tin = 0.0D0
    tratio = 1.0D0
    triang = 0.36D0
    triang95 = 0.24D0
    vol = 0.0D0
    vsbrn = 0.0D0
    vshift = 0.0D0
    vsind = 0.0D0
    vsres = 0.0D0
    vsstt = 0.0D0
    wallmw = 0.0D0
    wtgpd = 0.0D0
    xarea = 0.0D0
    zeff = 0.0D0
    zeffai = 0.0D0
  end subroutine init_physics_variables