scan_2d_write_plot Subroutine

public subroutine scan_2d_write_plot(iscan, outvar, sweep_1_vals, sweep_2_vals)

Arguments

Type IntentOptional AttributesName
integer, intent(inout) :: iscan
real(kind=dp), intent(in), dimension(:,:):: outvar
real(kind=dp), intent(in), dimension(:):: sweep_1_vals
real(kind=dp), intent(in), dimension(:):: sweep_2_vals

Contents

Source Code


Source Code

  subroutine scan_2d_write_plot(iscan, outvar, sweep_1_vals, sweep_2_vals)
    use constants, only: nplot
    use global_variables, only: icase, xlabel, xlabel_2
    implicit none

    integer, intent(inout) :: iscan
    real(dp), dimension(:,:), intent(in) :: outvar
    real(dp), dimension(:), intent(in) :: sweep_1_vals, sweep_2_vals

    integer :: ivar
    character(len=48) :: tlabel
    character(len=25), dimension(noutvars), save :: plabel

    plabel( 1) = 'Ifail____________________'
    plabel( 2) = 'Sqsumsq__________________'
    plabel( 3) = 'Electric_cost_(mil/kwh)__'
    plabel( 4) = 'Capital_cost_(mil/kwh)___'
    plabel( 5) = 'Fuel_cost_(mil/kwh)______'
    plabel( 6) = 'Operations_cost_(mil/kwh)'
    plabel( 7) = 'Capital_cost_(millions)__'
    plabel( 8) = 'Core_costs_(millions)____'
    plabel( 9) = 'Direct_cost_(billions)___'
    plabel(10) = 'Major_Radius_(m)_________'
    plabel(11) = 'Aspect_Ratio_____________'
    plabel(12) = 'Plasma_Current_(MA)______'
    plabel(13) = 'B_Toroidal_Axis_(T)______'
    plabel(14) = 'B_total_on_axis_(T)______'
    plabel(15) = 'Safety_Factor____________'
    plabel(16) = 'qlim_(zero_if_ishape=0)__'
    plabel(17) = 'Beta_____________________'
    plabel(18) = 'Beta_Limit_______________'
    plabel(19) = 'Epsilon_Beta_Poloidal____'
    plabel(20) = 'Dens.weight_Te_(10keV)___'
    plabel(21) = 'Average_Dens_(10^20/m^3)_'
    plabel(22) = 'H-fact_Iter_Power________'
    plabel(23) = 'H-fact_Iter_Offset_______'
    plabel(24) = 'Fusion_Power_(MW)________'
    plabel(25) = 'nb_Fusion_Power_(MW)_____'
    plabel(26) = 'Wall_Load_(MW/m^2)_______'
    plabel(27) = 'Injection_Power_(MW)_____'
    plabel(28) = 'Inject_Pwr_Wall_Plug_(MW)'
    plabel(29) = 'Heating_Power_(MW)_______'
    plabel(30) = 'Current_Drive_(MW)_______'
    plabel(31) = 'Big_Q____________________'
    plabel(32) = 'Bootstrap_Fraction_______'
    plabel(33) = 'Neutral_Beam_Energy_(MeV)'
    plabel(34) = 'Divertor_Heat_(MW/m^2)___'
    plabel(35) = 'TF_coil_Power_(MW)_______'
    plabel(36) = 'TF_coil_weight_(kg)______'
    plabel(37) = 'vM_stress_in_TF_cond_(Pa)'
    plabel(38) = 'J_TF_inboard_leg_(MA/m^2)'
    plabel(39) = 'Centrepost_max_T_(TART)__'
    plabel(40) = 'Res_TF_inbrd_leg_Pwr_(MW)'
    plabel(41) = 'Coolant_Fraction_Ctr.____'
    plabel(42) = 'C/P_coolant_radius_(m)___'
    plabel(43) = 'C/P_coolant_velocity(m/s)'
    plabel(44) = 'C/P_pump_power_(MW)______'
    plabel(45) = 'PF_coil_Power_(MW)_______'
    plabel(46) = 'PF_coil_weight_(kg)______'
    plabel(47) = 'Gross_Elect_Pwr_(MW)_____'
    plabel(48) = 'Net_electric_Pwr_(MW)____'
    plabel(49) = 'Recirculating_Fraction___'
    plabel(50) = 'Psep/R___________________'
    plabel(51) = '' !OBSOLETE
    plabel(52) = 'Tot._radiation_power_(MW)'
    plabel(53) = 'First_wall_peak_temp_(K)_'
    plabel(54) = 'Cu_frac_TFC_conductor____'
    plabel(55) = 'Winding_pack_area_TFC(m2)'
    plabel(56) = 'Conductor_area_TFC_(m2)__'
    plabel(57) = 'Area_TF_inboard_leg_(m2)_'
    plabel(58) = 'Taup/taueff_lower_limit__'
    plabel(59) = 'Plasma_temp_at_sep_[keV]_'
    plabel(60) = 'SOL_density_at_OMP_______'
    plabel(61) = 'Power_through__separatrix'
    plabel(62) = 'neomp/nesep______________'
    plabel(63) = 'qtargettotal_____________'
    plabel(64) = 'Total_pressure_at_target_'
    plabel(65) = 'Temperature_at_target____'
    plabel(66) = 'Helium_fraction__________'
    plabel(67) = 'Momentum_loss_factor_____'
    plabel(68) = 'totalpowerlost_[W]_______'
    plabel(69) = 'H__concentration_________'
    plabel(70) = 'He_concentration_________'
    plabel(71) = 'Be_concentration_________'
    plabel(72) = 'C__concentration_________'
    plabel(73) = 'N__concentration_________'
    plabel(74) = 'O__concentration_________'
    plabel(75) = 'Ne_concentration_________'
    plabel(76) = 'Si_concentration_________'
    plabel(77) = 'Ar_concentration_________'
    plabel(78) = 'Fe_concentration_________'
    plabel(79) = 'Ni_concentration_________'
    plabel(80) = 'Kr_concentration_________'
    plabel(81) = 'Xe_concentration_________'
    plabel(82) = 'W__concentration_________'
    plabel(83) = 'teped____________________'
    plabel(84) = 'Max_field_on_TF_coil_____'

    tlabel = icase

  end subroutine scan_2d_write_plot