init_pfcoil_variables Subroutine

  subroutine init_pfcoil_variables
    !! Initialise module variables
    implicit none

    alfapf = 5.0D-10
    alstroh = 4.0D8
    i_cs_stress = 0
    areaoh = 0.0D0
    a_oh_turn = 0.0D0
    awpoh = 0.0D0
    bmaxoh = 0.0D0
    bmaxoh0 = 0.0D0
    bpf = 0.0D0
    ccl0_ma = (/0.0D0, 0.0D0, 0.0D0, &
      0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0/)
    ccls_ma = (/0.0D0, 0.0D0, 0.0D0, &
      0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0, 0.0D0/)
    cohbop = 0.0D0
    coheof = 1.85D7
    cpt = 0.0D0
    cptdin = (/4.0D4, 4.0D4, 4.0D4, &
    4.0D4, 4.0D4, 4.0D4, 4.0D4, 4.0D4, 4.0D4, 4.0D4, &
    4.0D4, 4.0D4, 4.0D4, 4.0D4, 4.0D4, 4.0D4, 4.0D4, &
    4.0D4, 4.0D4, 4.0D4, 4.0D4, 4.0D4/)
    curpfb = 0.0D0
    curpff = 0.0D0
    curpfs = 0.0D0
    etapsu = 0.9D0
    fcohbof = 0.0D0
    fcohbop = 0.9D0
    fcuohsu = 0.7D0
    fcupfsu = 0.69D0
    fvssu = 1.0
    ipfloc = (/2,2,3,0,0,0,0,0,0,0/)
    ipfres = 0
    itr_sum = 0.0D0
    isumatoh = 1
    isumatpf = 1
    j_crit_str_cs = 0.0D0
    j_crit_str_pf = 0.0D0
    i_pf_current = 1
    i_sup_pf_shape = 0
    jscoh_bop = 0.0D0
    jscoh_eof = 0.0D0
    jcableoh_bop = 0.0D0
    jcableoh_eof = 0.0D0
    ncirt = 0
    ncls = (/1,1,2,0,0,0,0,0,0,0,0,0/)
    nfxfh = 7
    ngrp = 3
    nohc = 0
    ohhghf = 0.71D0
    oh_steel_frac = 0.5D0
    pf_current_safety_factor = 1.0D0
    pfcaseth = 0.0D0
    pfclres = 2.5D-8
    rhopfbus = 3.93D-8
    pfmmax = 0.0D0
    pfrmax = 0.0D0
    pfwpmw = 0.0D0
    powohres = 0.0D0
    powpfres = 0.0D0
    ra = 0.0D0
    rb = 0.0D0
    ric = 0.0D0
    rjconpf = (/3.0D7, 3.0D7, 3.0D7, &
      3.0D7, 3.0D7, 3.0D7, 3.0D7, 3.0D7, 3.0D7, 3.0D7, &
      3.0D7, 3.0D7, 3.0D7, 3.0D7, 3.0D7, 3.0D7, 3.0D7, &
      3.0D7, 3.0D7, 3.0D7, 3.0D7, 3.0D7/)
    rjohc = 0.0D0
    rjohc0 = 0.0D0
    rjpfalw = 0.0D0
    rohc = 0.0D0
    routr = 1.5D0
    rpf = 0.0D0
    rpf1 = 0.0D0
    rpf2 = -1.63D0
    rref = (/7.0D0, 7.0D0, 7.0D0, &
      7.0D0, 7.0D0, 7.0D0, 7.0D0, 7.0D0, 7.0D0, 7.0D0/)
    s_tresca_oh = 0.0D0
    sigpfcalw = 500.0D0
    sigpfcf = 1.0D0
    sxlg = 0.0D0
    tmargoh = 0.0D0
    turns = 0.0D0
    vf = 0.3D0
    vfohc = 0.3D0
    vsbn = 0.0D0
    vsefbn = 0.0D0
    vsefsu = 0.0D0
    vseft = 0.0D0
    vsoh = 0.0D0
    vsohbn = 0.0D0
    vsohsu = 0.0D0
    vssu = 0.0D0
    vstot = 0.0D0
    waves = 0.0D0
    whtpf = 0.0D0
    whtpfs = 0.0D0
    wtc = 0.0D0
    wts = 0.0D0
    zh = 0.0D0
    zl = 0.0D0
    zpf = 0.0D0
    zref = (/3.6D0, 1.2D0, 2.5D0, &
      1.0D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0/)
    bmaxcs_lim = 13.0
    fbmaxcs = 1.0
    ld_ratio_cst = 70/22
    l_cond_cst = 0.0D0
    d_cond_cst = 0.0D0
    r_in_cst = 0.0D0
    r_out_cst = 3.0D-3
  end subroutine init_pfcoil_variables