

&star_job

      show_log_description_at_start = .false. 

      load_saved_model = .true.
      saved_model_name = 'wd_diff_test.mod'

      save_model_when_terminate = .true.
      save_model_filename = 'final.mod'
      
      write_profile_when_terminate = .true.
      filename_for_profile_when_terminate = 'final_profile.data'

      kappa_file_prefix = 'gs98' 
      eos_file_prefix = 'mesa'  

      set_tau_factor = .true.  
      set_to_this_tau_factor = 15

      eos_file_prefix = 'mesa'  

      !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 = 100
         max_model_number = 1000

      min_dq = 1d-20
      
      mesh_delta_coeff = 0.75
      varcontrol_target = 2d-4
      
      initial_mass =  0.20
      initial_z = 0.01d0
      
      use_Type2_opacities = .true.
      Zbase = 0.01d0

      which_atm_option = 'grey_and_kap' ! use this for getting pulsation data including atmosphere
      !which_atm_option = 'Paczynski_grey' ! use this for getting pulsation data including atmosphere
      !Paczynski_atm_R_surf_errtol = 3d-4
      !create_atm_max_step_size = 0.05 ! in units of log10_tau
      
      
      !max_age = 1d8
      !max_model_number = 111

      !when_to_stop_rtol = 1d-3
      !when_to_stop_atol = 1d-3



      eps_nuc_factor = 1 
      
      Teff_lower_limit = 5d3 ! 1d4 !12000

      photostep = 10
      profile_interval = 10
      history_interval = 1
      terminal_cnt = 1
      write_header_frequency = 10

      
      use_Ledoux_criterion = .true.
      thermohaline_coeff = 1000  

      MLT_option = 'ML2'
      mixing_length_alpha = 1.0

      ! element diffusion
                     
      show_diffusion_info = .true. ! terminal output for diffusion
      !show_diffusion_substep_info = .true. ! terminal output for diffusion
      !show_diffusion_timing = .true.

      do_element_diffusion = .true.

      diffusion_num_classes = 8 ! 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) = 'c12'
      diffusion_class_representative(5) = 'n14'
      diffusion_class_representative(6) = 'o16'
      diffusion_class_representative(7) = 'ne20'
      diffusion_class_representative(8) = '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) = 12
      diffusion_class_A_max(5) = 14
      diffusion_class_A_max(6) = 16
      diffusion_class_A_max(7) = 20
      diffusion_class_A_max(8) = 10000
         
      diffusion_use_isolve = .true.
      diffusion_rtol_for_isolve = 1d-4
      diffusion_atol_for_isolve = 1d-5
      diffusion_maxsteps_for_isolve = 1000
      diffusion_isolve_solver = 'ros2_solver'
         
      
      
      tol_residual_norm1 = 1d-5
      tol_max_residual1 = 1d-2
      tol_residual_norm2 = 1d99
      tol_max_residual2 = 1d99
      iter_for_resid_tol2 = 3
      iter_for_resid_tol3 = 12


/ ! end of controls namelist


&pgstar
         


/ ! end of pgstar namelist
