! inlist_massive_defaults


&star_job
      show_log_description_at_start = .false.

      change_net = .true.      
      new_net_name = 'approx26.net'
      
!      reaction_for_special_factor(1) = 'r1212'
!      special_rate_factor(1) = 8
!      num_special_rate_factors = 1
      
      ! rate limits for reactions involving a neutron (and optionally protons)
      T_lo_neut = 2d9 ! below this, limit doesn't apply
      T_hi_neut= 5d9 ! above this, limit in full force
      min_factor_neut = 1d-9 ! use this factor for T >= T_hi_neut
         ! factor gradually goes to 1.0 as T drops from T_hi_neut to T_lo_neut

      ! rate limits for reactions involving a proton (but no neutrons)
      T_lo_prot = 2d9 ! below this, limit doesn't apply
      T_hi_prot = 5d9 ! above this, limit in full force
      min_factor_prot = 1d-12 ! use this factor for T >= T_hi_prot
         ! factor gradually goes to 1.0 as T drops from T_hi_prot to T_lo_prot

      ! rate limit factors for high temperatures
      turn_on_T_limits = .true.
      turn_off_T_limits = .false.

      
      kappa_file_prefix = 'gs98'

      change_v_flag = .true.
      new_v_flag = .true.
      
      !new_rotation_flag = .true.
      !change_rotation_flag = .false.  ! rotation off until near zams
      !new_omega_div_omega_crit = 0.5  
      !near_zams_relax_omega_div_omega_crit = .true.
      !num_steps_to_relax_rotation = 50 ! use this many steps to change value

      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
         
         !min_T_for_time_averaged_conv_velocity = 3e9
            
         ! high center T limit to avoid negative mass fractions
         sig_min_factor_for_high_Tcenter = 0.01
          ! inactive when >= 1d0
            ! if Tcenter >= Tcenter_min_for_sig_min_factor_full_on,
            ! then okay to reduce sig by as much as this factor
            ! as needed to prevent causing negative abundances
         Tcenter_min_for_sig_min_factor_full_on = 2d9
            ! if Tcenter >= this, factor = sig_min_factor_for_neg_abundances
            ! this should be > Tcenter_max_for_sig_min_factor_full_off.
         Tcenter_max_for_sig_min_factor_full_off = 1d9
            ! if Tcenter <= this, factor = 1, so has no effect
            ! this should be < Tcenter_min_for_sig_min_factor_full_on.
         ! for T > full_off and < full_on, factor changes linearly with Tcenter

         velocity_logT_lower_bound = 8
         max_dt_yrs_for_velocity_logT_lower_bound = 0.1

         mesh_delta_coeff = 1.5
         mesh_delta_coeff_for_highT = 2.5
         logT_max_for_standard_mesh_delta_coeff = 9.0 
         logT_min_for_highT_mesh_delta_coeff = 9.5 

         varcontrol_target = 7d-4         
         dX_nuc_drop_limit = 7d-3
         delta_Ye_highT_limit = 1d-4

         okay_to_reduce_gradT_excess = .true. 
         
            gradT_excess_f1 = 1d-4
            gradT_excess_f2 = 1d-2
            gradT_excess_age_fraction = 0.9d0

            gradT_excess_lambda1 = 1.0
            gradT_excess_beta1 = 0.35
            
            gradT_excess_lambda2 = 0.5
            gradT_excess_beta2 = 0.25
            
            gradT_excess_dlambda = 0.1
            gradT_excess_dbeta = 0.1


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

         include_composition_in_eps_grav = .true.

         use_Type2_opacities = .true.
         Zbase = 0.02

         xa_scale = 1d-5
         newton_itermin = 2

         mixing_length_alpha = 1.5
         use_Henyey_MLT = .true.

         use_Ledoux_criterion = .true.
         alpha_semiconvection = 0.01
         thermo_haline_coeff = 2
         thermohaline_option = 'Kippenhahn'

         overshoot_f_above_nonburn = 0.001
         overshoot_f_above_burn_h = 0.001
         overshoot_f_above_burn_he = 0.001
         overshoot_f_above_burn_z = 0.001
         
         overshoot_f_below_nonburn = 0.001
         overshoot_f_below_burn_h = 0.001
         overshoot_f_below_burn_he = 0.001
         overshoot_f_below_burn_z = 0.001

         set_min_D_mix = .true.
         min_D_mix = 1d3 ! D_mix will be at least this large
         min_center_Ye_for_min_D_mix = 0

         screening_mode = 'extended'
         
         max_iter_for_resid_tol1 = 3
         tol_residual_norm1 = 1d-5
         tol_max_residual1 = 1d-2
         
         max_iter_for_resid_tol2 = 12
         tol_residual_norm2 = 1d99
         tol_max_residual2 = 1d99
         
         min_timestep_limit = 1d-12 ! (seconds)

         delta_lgL_He_limit = 0.1 ! limit for magnitude of change in lgL_He

         dX_nuc_drop_max_A_limit = 52
         dX_nuc_drop_min_X_limit = 1d-4
         dX_nuc_drop_hard_limit = 1d99
         
         delta_lgTeff_limit = 0.5  
         delta_lgL_limit = 0.5
         delta_lgRho_cntr_limit = 0.02  

!         mesh_dlog_pp_dlogP_extra = 0.4
!         mesh_dlog_cno_dlogP_extra = 0.4
!      
!         mesh_dlog_burn_n_dlogP_extra = 0.4
!         mesh_dlog_3alf_dlogP_extra = 0.4
!         mesh_dlog_burn_c_dlogP_extra = 0.4

			xa_function_species(1) = 'h1'
			xa_function_weight(1) = 70
			xa_function_param(1) = 1d-3

			xa_function_species(2) = 'he4'
			xa_function_weight(2) = 70
			xa_function_param(2) = 1d-3

         T_mix_limit = 0


/ ! end of controls namelist


&pgstar


/ ! end of pgstar namelist
