
&star_job

      show_log_description_at_start = .false.

      load_saved_model = .true.
      saved_model_name = 'solar_calibration_input.mod'
      
      save_model_when_terminate = .true.
      save_model_filename = 'solar.mod'

      change_net = .true. ! switch nuclear reaction network
      new_net_name = 'pp_and_cno_extras.net'
	   
      set_uniform_initial_composition = .true.
	   initial_zfracs = 3 ! GS98_zfracs = 3
      initial_h1 = 0.69936748978639d0
      initial_h2 = 0
      initial_he3 = 0.00002803008078d0
      initial_he4 = 0.28027277776185d0
   ! initial_z = 0.02033399508500
	   
	   set_rate_n14pg = 'Imbriani'
	   set_rate_c12ag = 'Kunz'

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


/ ! end of star_job namelist


&controls

      ! check for retries and backups as part of test_suite
      ! you can/should delete this for use outside of test_suite
         max_number_backups = 10
         max_number_retries = 80

         max_model_number = 2000
         
         
      max_age = 4.57d9
      
      mixing_length_alpha = 1.91054420292493d0
      
      overshoot_f_below_nonburn = 0.01242157889075d0

      ! controls for output
      photostep = 50
      profile_interval = 50
      history_interval = 10
      terminal_cnt = 10
      write_header_frequency = 10


      ! atmosphere
      which_atm_option = 'photosphere_tables'
      
      ! atomic diffusion
      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_T_full_on = 1d3
      diffusion_T_full_off = 1d3
      
      diffusion_calculates_ionization = .true.

      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
 
      ! timesteps
      max_years_for_timestep = 1d7           
      varcontrol_target = 1d-4

      ! 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
      
      use_Type2_opacities = .false.
      cubic_interpolation_in_Z = .true.


/ ! end of controls namelist



&pgstar
         
      ! top level controls

         !pause = .true. 
            ! if true, the code waits for user to enter a RETURN on the command line
      

/ ! end of pgstar namelist
