
&one_zone

      !show_net_reactions_info = .true.
      
      mesa_data_dir = '../../data'   	
      !net_name = 'pp_and_cno_extras.net'
      !net_name = 'rp_mg.net'
      !net_name = 'rp_si.net'
      !net_name = 'rp_s.net'
      !net_name = 'rp_ar.net'
      !net_name = 'rp_ca.net'
      !net_name = 'rp_ca38_x.net'
      !net_name = 'rp_ti.net'
      !net_name = 'rp_cr.net'
      !net_name = 'rp_fe.net'
      !net_name = 'rp_ni.net'
      !net_name = 'rp_zn.net'
      !net_name = 'rp_ge.net'
      !net_name = 'rp.net'
      !net_name = 'ns_h.net'
      !net_name = 'ns_h_smaller.net'
      net_name = 'arnett.net'
      
      !net_name = 'explosive_burning_big.net'

   	!burn_rho = 100
   	!burn_temp = 19e6
   	!data_heading_line = 'Solar burn at T = 19e6, Rho = 100'
   	!burn_tend = 3.16d17 ! 10 Gyr

   	!burn_rho = 100
   	!burn_temp = 4e8
   	!data_heading_line = 'hot cno solar burn at T = 4e8, Rho = 100'
   	!burn_tend = 5e4 ! seconds

   	!burn_rho = 4.2d5
   	!burn_temp = 9.98d8

   	burn_rho = 1e7
   	burn_temp = 5e9

      !rattab_logT_lower_bound = 8
      !rattab_logT_upper_bound = 10

   	data_heading_line = 'T9 = 5, rho = 1e7'
   	burn_tend = 10 ! seconds 
      
      min_for_show_peak_abundances = 1d-14 ! output peak value if it is > this

      data_filename = 'one_zone_burn.data'
      !data_output_min_t = 1d-6
      
      !save_final_abundances = .true.
      !final_abundances_filename = 'x_at_1_second.data'
      
      !read_initial_abundances = .true.
      !initial_abundances_filename = 'x_at_1_second.data'
      
      !uniform_Xinit = .true. ! for debugging

      num_isos_for_Xinit = 2
      names_of_isos_for_Xinit(1) = 'c12'
      values_for_Xinit(1) = 0.5d0
      names_of_isos_for_Xinit(2) = 'o16'
      values_for_Xinit(2) = 0.5d0


      num_names_of_isos_to_show = 16
      names_of_isos_to_show(1) = 'he4'
      names_of_isos_to_show(2) = 'c12'
      names_of_isos_to_show(3) = 'o16'
      names_of_isos_to_show(4) = 'ne20'
      names_of_isos_to_show(5) = 'mg24'
      names_of_isos_to_show(6) = 'si28'
      names_of_isos_to_show(7) = 's32'
      names_of_isos_to_show(8) = 'ar36'
      names_of_isos_to_show(9) = 'ca40'
      names_of_isos_to_show(10) = 'ti44'
      names_of_isos_to_show(11) = 'cr48'
      names_of_isos_to_show(12) = 'fe52'
      names_of_isos_to_show(13) = 'ni56'
      names_of_isos_to_show(14) = 'zn60'
      names_of_isos_to_show(15) = 'ge64'
      names_of_isos_to_show(16) = 'se68'      
      
   	
      burn_rtol = 1d-8
      burn_atol = 1d-9
      which_solver = 'seulex_solver'
      !which_solver = 'sodex_solver'

   	
      !burn_rtol = 1d-3
      !burn_atol = 1d-4
      !which_solver = 'ros2_solver'
      
      
      
      decsol_switch = 20
         ! if current number of species <= switch,
            ! then use small_mtx_decsol,
            ! else use large_mtx_decsol.
      small_mtx_decsol = 'lapack'
      large_mtx_decsol = 'klu'
      
      cache_suffix = '' ! use this only for reactions from rate_tables
      
      screening_mode = 2 ! extended_screening
      theta_e_for_graboske_et_al = 1
      which_rates_choice = 2
         ! rates_NACRE_if_available = 1
         ! rates_JR05_if_available = 2 << for jina reaclib


/ 
