&star_job

      show_log_description_at_start = .false.
      !show_net_species_info = .true.

      load_saved_model = .true.
      saved_model_name = 'solar_calibration_input.mod'

      change_net = .true. ! switch nuclear reaction network
      new_net_name = 'pp_and_cno_extras.net'
      !new_net_name = 'basic.net'

      !save_fgong_when_terminate = .true.
      !save_fgong_filename = 'solar.fgong'
	   
      set_uniform_initial_composition = .true.
	   initial_zfracs = 3 ! GS98_zfracs = 3
      initial_h1 = 0.700262d0
      initial_h2 = 0
      initial_he3 = 0.279364d-4
      initial_he4 = 0.279364d0

	   set_rate_n14pg = 'Imbriani'
	   set_rate_c12ag = 'Kunz'

      chem_isotopes_filename = 'isotopes.data_real' ! this file is in chem_data in mesa_data_dir
      !chem_isotopes_filename = 'isotopes.data_approx' ! this file is in chem_data in mesa_data_dir

      kappa_file_prefix = 'OP'
      kappa_lowT_prefix = 'lowT_fa05_gs98' ! for lower temperatures.
      
      
      change_lnPgas_flag = .true.
      new_lnPgas_flag = .true.
      

/ ! end of star_job namelist


&controls

      D_norm_kappa = 1
      D_norm_alpha = 0 
      D_norm_atol = 1d-4
      D_norm_rtol = 1d-4
      
      use_lnE_for_eps_grav = .true.


      max_age = 4.57d9
      
      mixing_length_alpha = 2.02144520885243d0

      ! controls for output
      photostep = 50
      profile_interval = 50
      log_cnt = 1
      terminal_cnt = 2
      write_header_frequency = 10

      !calculate_Brunt_N2 = .true. 

      !report_hydro_solver_progress = .true. ! set true to see info about newton iterations
      report_ierr = .true. ! if true, produce terminal output when have some internal error

      !report_why_dt_limits = .true.
      !report_all_dt_limits = .true.
      
      report_hydro_dt_info = .false.
      
      !show_mesh_changes = .true.
      !mesh_dump_call_number = 79

      !trace_evolve = .true.
      
      ! hydro debugging
      !hydro_numerical_jacobian = .true.
      !hydro_check_everything = .true.
      !hydro_show_rcond = .true.
      !hydro_inspectB_flag = .true.
      !hydro_save_numjac_plot_data = .true.
      !hydro_dump_call_number = 1




      ! atmosphere
      which_atm_option = 'photosphere_tables'
      !which_atm_option = 'solar_Hopf_grey'   
      !which_atm_option = 'Eddington_grey'   
         
      ! element diffusion
      !show_diffusion_info = .true. ! terminal output for diffusion
      !diffusion_dump_call_number = 228

      do_element_diffusion = .true. ! determines whether or not we do diffusion
      diffusion_dt_limit = 7d11 ! no element diffusion if dt < this limit (in seconds)
      diffusion_atol = 1d-5
      diffusion_rtol = 1d-4
      diffusion_T_full_on = 1d3
      diffusion_T_full_off = 1d3
      
      !diffusion_min_dq_at_surface = 0.001
            ! treat at least this much at surface as a single cell for purposes of diffusion
      
      diffusion_calculates_ionization = .true.
      
!      diffusion_num_classes = 5 ! number of classes of species for diffusion calculations
!      diffusion_class_representative(1) = 'h1'
!      diffusion_class_representative(2) = 'he3'
!      diffusion_class_representative(3) = 'he4'
!      diffusion_class_representative(4) = 'o16'
!      diffusion_class_representative(5) = 'fe56'
   
      ! in ascending order.  species goes into 1st class with A_max >= species A
!      diffusion_class_A_max(1) = 2
!      diffusion_class_A_max(2) = 3
!      diffusion_class_A_max(3) = 4
!      diffusion_class_A_max(4) = 16
!      diffusion_class_A_max(5) = 10000
      
      diffusion_num_classes = 4 ! number of classes of species for diffusion calculations
      diffusion_class_representative(1) = 'h1'
      diffusion_class_representative(2) = 'he4'
      diffusion_class_representative(3) = 'o16'
      diffusion_class_representative(4) = 'fe56'
!   
!      ! in ascending order.  species goes into 1st class with A_max >= species A
      diffusion_class_A_max(1) = 2
      diffusion_class_A_max(2) = 4
      diffusion_class_A_max(3) = 16
      diffusion_class_A_max(4) = 10000
      
!      diffusion_num_classes = 3 ! number of classes of species for diffusion calculations
!      diffusion_class_representative(1) = 'h1'
!      diffusion_class_representative(2) = 'he4'
!      diffusion_class_representative(3) = 'fe56'
!   
!      ! in ascending order.  species goes into 1st class with A_max >= species A
!      diffusion_class_A_max(1) = 2
!      diffusion_class_A_max(2) = 4
!      diffusion_class_A_max(3) = 10000
 
 
      ! timesteps
      max_years_for_timestep = 1d7           
      varcontrol_target = 1d-4

      ! overshoot
      !overshoot_f_below_nonburn = 0.014


      ! timestep
      !dHe3_div_He3_limit_min_He3 = 1d-5
         ! dHe3_div_He3 limits only apply where xa(he3,k) >= this limit
      !dHe3_div_He3_limit = 1d-2
      !dHe3_decreases_only = .false.

      !dHe_div_He_limit_min_He = 1d-5
         ! dHe_div_He limits only apply where xa(he4,k) >= this limit
      !dHe_div_He_limit = 1d-3
      !dHe_decreases_only = .false.
      

      ! mesh adjustment
      mesh_delta_coeff = 0.8
   	
   	P_function_weight = 25
   	T_function1_weight = 75
      
      xtra_coef_czb_full_on = 1
      xtra_coef_czb_full_off = 1

      xtra_coef_a_l_nb_czb = 0.45 ! above lower nonburn convective boundary
      xtra_dist_a_l_nb_czb = 20 ! above lower nonburn convective boundary

      xtra_coef_b_l_nb_czb = 0.45 ! below lower nonburn convective boundary
      xtra_dist_b_l_nb_czb = 20 ! below lower nonburn convective boundary

         
      xa_function_species(1) = 'he4'  ! name of nuclide as defined in chem_def
		xa_function_weight(1) = 130
		xa_function_param(1) = 1d-2
         
      xa_function_species(2) = 'he3'  ! name of nuclide as defined in chem_def
		xa_function_weight(2) = 60
		xa_function_param(2) = 1d-5

		


      ! opacity
      cubic_interpolation_in_Z = .true.

      trace_solar_neutrinos = .false. ! if true, then output neutrino flux info

      

/ ! end of controls namelist

