
&star_job

      show_log_description_at_start = .false.
      
      !show_net_species_info = .true.
      !show_net_reactions_info = .true.

      !create_pre_main_sequence_model = .true.
      
      !save_model_number = 90
      !save_model_filename = 'test.mod'

      load_saved_model = .true.
      saved_model_name = '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'
      
      change_lnPgas_flag = .true.
      new_lnPgas_flag = .true.

      change_initial_net = .true.      
      new_net_name = 'o18_and_ne22.net'

!      pre_ms_relax_num_steps = 50
      
!      new_surface_rotation_v = 2 ! solar (km sec^1)
!      set_near_zams_surface_rotation_v_steps = 10 ! to turn on rotation when near zams
         ! if rotation_flag is false and L_nuc_burn_total >= L_phot and this control > 0
         ! then turn on rotation and set set_surf_rotation_v_step_limit to
         ! the current model_number plus this control
   
      set_initial_age = .true.
      initial_age = 0 ! in years

      set_initial_model_number = .true.
      initial_model_number = 0

      !pgstar_flag = .true.


/ ! end of star_job namelist

&controls
         
      max_model_number = 1

      initial_mass = 1.0
      initial_z = 0.02d0
      
      varcontrol_target = 2d-4
      mesh_delta_coeff = 0.7

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

      log_L_lower_limit = -0.5

      smooth_convective_bdy = .true.                  
      convective_bdy_weight = 1
      
      RGB_wind_scheme = 'Reimers'
      AGB_wind_scheme = 'Blocker'
      RGB_to_AGB_wind_switch = 1d-4
      Reimers_wind_eta = 0.7d0  
      Blocker_wind_eta = 0.7d0  

!      set_min_D_mix = .true.
      min_D_mix = 1d2 ! D_mix will be at least this large
      min_center_Ye_for_min_D_mix = 0.47 ! min_D_mix is only used when center Ye >= this
         ! i.e., when center Ye drops below this, min_D_mix is 0
      
      
      num_cells_for_smooth_brunt_B = 0
      !use_brunt_gradmuX_form = .true.
      
      !show_diffusion_info = .true. ! terminal output for diffusion
         !show_diffusion_substep_info = .true.         
         !show_diffusion_timing = .true.

      do_element_diffusion = .true. ! determines whether or not we do diffusion
      
      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
         
      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'


      ! FOR DEBUGGING

      !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
      
      !max_years_for_timestep = 1d6

      !report_why_dt_limits = .true.
      !report_all_dt_limits = .true.
      
      !show_mesh_changes = .true.
      !mesh_dump_call_number = 5189
      !okay_to_remesh = .false.
      
      !trace_evolve = .true.
      

      ! hydro debugging
      !hydro_check_everything = .true.
      !hydro_inspectB_flag = .true.
      
      !hydro_numerical_jacobian = .true.
      !hydro_save_numjac_plot_data = .true.
      !small_mtx_decsol = 'lapack'
      !large_mtx_decsol = 'lapack'
      !hydro_dump_call_number = 195



      !trace_newton_bcyclic_solve_input = .true. ! input is "B" j k iter B(j,k)
      !trace_newton_bcyclic_solve_output = .true. ! output is "X" j k iter X(j,k)
      
      !trace_newton_bcyclic_steplo = 1 ! 1st model number to trace
      !trace_newton_bcyclic_stephi = 1 ! last model number to trace
      
      !trace_newton_bcyclic_iterlo = 2 ! 1st newton iter to trace
      !trace_newton_bcyclic_iterhi = 2 ! last newton iter to trace
      
      !trace_newton_bcyclic_nzlo = 1 ! 1st cell to trace
      !trace_newton_bcyclic_nzhi = 10000 ! last cell to trace; if < 0, then use nz as nzhi
      
      !trace_newton_bcyclic_jlo = 1 ! 1st var to trace
      !trace_newton_bcyclic_jhi = 100 ! last var to trace; if < 0, then use nvar as jhi
      
      !trace_k = 2


/ ! 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
