
&star_job

         show_log_description_at_start = .false. 

         create_pre_main_sequence_model = .true.

         set_rate_c12ag = 'Buchmann'

         num_special_rate_factors = 1
         reaction_for_special_factor(:) = 'r_c12_ag_o16'
         special_rate_factor(1) = 1.1
         
         kappa_file_prefix = 'gs98'

         !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 = 30
         max_number_retries = 60
         
         
      okay_to_reduce_gradT_excess = .true.
         
            ! these are for calculation of efficiency boosted gradT
            gradT_excess_f1 = 1d-4
            gradT_excess_f2 = 1d-2
            gradT_excess_age_fraction = 0.9d0
            
            ! turn on full all the time
            gradT_excess_lambda1 = 1.0
            gradT_excess_beta1 = 0.5
            gradT_excess_lambda2 = 1.0
            gradT_excess_beta2 = 0.5
            gradT_excess_dlambda = 1
            gradT_excess_dbeta = 1

      
         initial_mass = 40
         initial_z = 1d-2
         
         varcontrol_target = 1d-3
         mesh_delta_coeff = 1.4
         
         xa_central_upper_limit_species(1) = 'ne20'
         xa_central_upper_limit(1) = 0.08
         
         max_model_number = 10000

         rate_factor_for_fe56ec_to_cr56 = 1d-4

         RGB_wind_scheme = 'Dutch'
         AGB_wind_scheme = ''
         RGB_to_AGB_wind_switch = 1d-4
         Dutch_wind_eta = 0.4

         log_center_density_limit = 12 
         log_center_temp_limit = 10.3d0 ! max T about 2d10
         
         xa_scale = 1d-5
         newton_itermin = 2

         !mass_change_full_on_dt = 1d7 ! (seconds)
         !mass_change_full_off_dt = 1d6 ! (seconds)

         mixing_length_alpha = 2
         screening_mode = 'extended' ! empty string means no screening

         use_Ledoux_criterion = .true.
         alpha_semiconvection = 0.1

         overshoot_f_above_nonburn = 0.01
         overshoot_f_below_nonburn = 0.01
         
         min_timestep_limit = 1d-12 ! (seconds)

         use_artificial_viscosity = .false.
         
         tol_correction_high_T_limit = 1d9 ! if center T >= this limit, switch to following
         tol_correction_norm_high_T = 3d-4
         tol_max_correction_high_T = 3d-2
         
         delta_lgL_He_limit = -1

         dX_nuc_drop_max_A_limit = 52
         dX_nuc_drop_min_X_limit = 1d-4
         dX_nuc_drop_limit = 5d-3
         dX_nuc_drop_hard_limit = 1d99
         
         delta_lgTeff_limit = 0.005  
         delta_lgL_limit = 0.05
         delta_lgRho_cntr_limit = 0.02  

         mesh_dlog_pp_dlogP_extra = 0.15
         mesh_dlog_cno_dlogP_extra = 0.15
      
         mesh_dlog_3alf_dlogP_extra = 0.15
         mesh_dlog_burn_c_dlogP_extra = 0.15
         
         mesh_dlog_burn_n_dlogP_extra = 0.15
         mesh_dlog_burn_o_dlogP_extra = 0.15
      
         mesh_dlog_burn_ne_dlogP_extra = 0.15
         
         mesh_dlog_burn_na_dlogP_extra = 0.15
         mesh_dlog_burn_mg_dlogP_extra = 0.15
         
         mesh_dlog_cc_dlogP_extra = 0.15
         mesh_dlog_co_dlogP_extra = 0.15
         
         mesh_dlog_oo_dlogP_extra = 0.15
         
         mesh_dlog_burn_si_dlogP_extra = 0.15
         mesh_dlog_burn_s_dlogP_extra = 0.15
         mesh_dlog_burn_ar_dlogP_extra = 0.15
         mesh_dlog_burn_ca_dlogP_extra = 0.15
         mesh_dlog_burn_ti_dlogP_extra = 0.15
         mesh_dlog_burn_cr_dlogP_extra = 0.15
         mesh_dlog_burn_fe_dlogP_extra = 0.15
         
         mesh_dlog_pnhe4_dlogP_extra = 0.15
         mesh_dlog_other_dlogP_extra = 0.15
         mesh_dlog_photo_dlogP_extra = 1

         T_mix_limit = 1d4


         photostep = 100
         profile_interval = 100
         history_interval = 10
         terminal_interval = 10
         write_header_frequency = 10


      ! FOR DEBUGGING
      !max_model_number = 2818

      !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 = 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 = 3327

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