
&star_job
      
      show_log_description_at_start = .false. 

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

      eos_file_prefix = 'mesa'
      kappa_file_prefix = 'gs98'
      
      change_initial_net = .true.      
      new_net_name = 'co_burn_plus.net'
         
      change_v_flag = .true.
      new_v_flag = .true.

      change_rotation_flag = .true.
      new_rotation_flag = .true.
      new_omega_div_omega_crit = 0.1
      set_omega_div_omega_crit_step_limit = 10

      set_initial_age = .true.
      initial_age = 0 ! in years

      set_initial_model_number = .true.
      initial_model_number = 0

      set_initial_dt = .true.
      years_for_initial_dt = 10

      !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 = 100
         max_number_retries = 500
         max_model_number = 4000

         
         num_trace_history_values = 2
         trace_history_value_name(1) = 'surf_avg_v_rot'
         trace_history_value_name(2) = 'surf_avg_v_div_v_crit'
         


      initial_mass = 10
      initial_z = 0.02d0 
      
      use_Type2_opacities = .true.
      Zbase = 0.02d0

      varcontrol_target = 1d-3
      dX_nuc_drop_limit = 1d-2

      mesh_delta_coeff = 1.5
         
      okay_to_reduce_gradT_excess = .true.
      gradT_excess_lambda1 = -1 ! turn on full all the time

      log_L_lower_limit = 2

      mixing_length_alpha = 2
      use_Ledoux_criterion = .true.
      alpha_semiconvection = 0.01
      thermohaline_coeff = 2  

      am_nu_visc_factor = 0
      am_D_mix_factor = 0.0333333333333333d0
      D_DSI_factor = 0
      D_SH_factor = 1
      D_SSI_factor = 1
      D_ES_factor = 1
      D_GSF_factor = 1
      D_ST_factor = 1

      photostep = 50
      profile_interval = 50
      history_interval = 1
      terminal_cnt = 10
      write_header_frequency = 10
      
      max_age = 10d9
      max_years_for_timestep = 1d9
      
      RGB_wind_scheme = ''
      AGB_wind_scheme = 'Blocker'
      RGB_to_AGB_wind_switch = 1d-4
      Reimers_wind_eta = 0.1d0  
      Blocker_wind_eta = 10d0  

      which_atm_option = 'grey_and_kap'
      Pextra_factor = -1

      overshoot_f_below_burn_z_shell = 0.1    
      overshoot_f0_below_burn_z_shell = 0.01    

      ! 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
      !hydro_show_correction_info = .true.
      
      !max_years_for_timestep = 3.67628942044319d-05

      !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 = 1339 ! 1st model number to trace
      trace_newton_bcyclic_stephi = 1339 ! last model number to trace
      
      trace_newton_bcyclic_iterlo = 0 ! 1st newton iter to trace
      trace_newton_bcyclic_iterhi = 111 ! 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 = 0

/ ! end of controls namelist



&pgstar
         pause = .true.
         Profile_Panels4_win_flag = .true. ! Abundance, Power, Mixing_Ds, and Dynamo    
         History_Panels1_win_flag = .true. ! logT, logRho, log_L, log_Teff
         Grid8_win_flag = .true. ! Summary_Burn, Abundance, HR, TRho, TRho_Profile, Text_Summary1

/ ! end of pgstar namelist
