init_pfcoil_variables Subroutine

  subroutine init_pfcoil_variables
    !! Initialise module variables
    implicit none

    alfapf = 5.0D-10
    alstroh = 4.0D8
    i_cs_stress = 0
    a_cs_poloidal = 0.0D0
    a_cs_turn = 0.0D0
    awpoh = 0.0D0
    b_cs_peak_flat_top_end = 0.0D0
    b_cs_peak_pulse_start = 0.0D0
    b_pf_coil_peak = 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/)
    j_cs_pulse_start = 0.0D0
    j_cs_flat_top_end = 1.85D7
    c_pf_coil_turn = 0.0D0
    c_pf_coil_turn_peak_input = (/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/)
    c_pf_cs_coil_pulse_start_ma = 0.0D0
    c_pf_cs_coil_flat_top_ma = 0.0D0
    c_pf_cs_coil_pulse_end_ma = 0.0D0
    etapsu = 0.9D0
    f_j_cs_start_end_flat_top = 0.0D0
    f_j_cs_start_pulse_end_flat_top = 0.9D0
    fcuohsu = 0.7D0
    fcupfsu = 0.69D0
    fvs_cs_pf_total_ramp = 1.0
    i_pf_location = (/2,2,3,0,0,0,0,0,0,0/)
    i_pf_conductor = 0
    itr_sum = 0.0D0
    i_cs_superconductor = 1
    i_pf_superconductor = 1
    j_crit_str_cs = 0.0D0
    j_crit_str_pf = 0.0D0
    i_pf_current = 1
    i_sup_pf_shape = 0
    j_cs_conductor_critical_pulse_start = 0.0D0
    j_cs_conductor_critical_flat_top_end = 0.0D0
    jcableoh_bop = 0.0D0
    jcableoh_eof = 0.0D0
    n_pf_cs_plasma_circuits = 0
    n_pf_coils_in_group = (/1,1,2,0,0,0,0,0,0,0,0,0/)
    nfxfh = 7
    n_pf_coil_groups = 3
    n_cs_pf_coils = 0
    f_z_cs_tf_internal = 0.71D0
    f_a_cs_steel = 0.5D0
    pf_current_safety_factor = 1.0D0
    pfcaseth = 0.0D0
    rho_pf_coil = 2.5D-8
    rhopfbus = 3.93D-8
    m_pf_coil_max = 0.0D0
    r_pf_coil_outer_max = 0.0D0
    pfwpmw = 0.0D0
    p_cs_resistive_flat_top = 0.0D0
    p_pf_coil_resistive_total_flat_top = 0.0D0
    r_pf_coil_inner = 0.0D0
    r_pf_coil_outer = 0.0D0
    c_pf_cs_coils_peak_ma = 0.0D0
    j_pf_coil_wp_peak = (/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/)
    j_cs_critical_flat_top_end = 0.0D0
    j_cs_critical_pulse_start = 0.0D0
    j_pf_wp_critical = 0.0D0
    r_cs_middle = 0.0D0
    routr = 1.5D0
    r_pf_coil_middle = 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_shear_cs_peak = 0.0D0
    sigpfcalw = 500.0D0
    sigpfcf = 1.0D0
    ind_pf_cs_plasma_mutual = 0.0D0
    temp_cs_margin = 0.0D0
    n_pf_coil_turns = 0.0D0
    f_a_pf_coil_void = 0.3D0
    f_a_cs_void = 0.3D0
    vs_cs_pf_total_burn = 0.0D0
    vs_pf_coils_total_burn = 0.0D0
    vs_pf_coils_total_ramp = 0.0D0
    vs_pf_coils_total_pulse = 0.0D0
    vs_cs_total_pulse = 0.0D0
    vs_cs_burn = 0.0D0
    vs_cs_ramp = 0.0D0
    vs_cs_pf_total_ramp = 0.0D0
    vs_cs_pf_total_pulse = 0.0D0
    waves = 0.0D0
    m_pf_coil_conductor_total = 0.0D0
    m_pf_coil_structure_total = 0.0D0
    m_pf_coil_conductor = 0.0D0
    m_pf_coil_structure = 0.0D0
    z_pf_coil_upper = 0.0D0
    z_pf_coil_lower = 0.0D0
    z_pf_coil_middle = 0.0D0
    zref = (/3.6D0, 1.2D0, 2.5D0, &
      1.0D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0, 1.0D0/)
    b_cs_limit_max = 13.0
    fb_cs_limit_max = 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