!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Inlist for MESA r24.08.1 to do rotating models.
! Please fill in the <> fields before running.
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

&star_job
      pgstar_flag = .true.

      show_log_description_at_start = .false.
      
      set_initial_model_number = .true.
      initial_model_number = 0

! PRE-MAIN SEQUENCE MODEL
     create_pre_main_sequence_model = .false.
     load_saved_model = .true.
     load_model_filename = 'model_ZAMS'
      
! OPACITY, NETWORK, RATE, etc.
      ! Nuclear network
      change_net = .true.
      new_net_name = 'pp_cno_extras_o18_ne22.net'

      ! Initial omega/omega_crit
      new_omega_div_omega_crit = 0.6d0
      set_near_zams_omega_div_omega_crit_steps = 20

/ ! end of star_job namelist

&eos

/ ! end of eos namelist

&kap

    ! opacities with GN93 abundances
    kap_file_prefix = 'gn93'
    kap_lowT_prefix = 'lowT_fa05_gn93'
    kap_CO_prefix = 'gn93_co'

    ! CO enhanced opacities
    use_Type2_opacities = .true.
    kap_Type2_full_off_X = 1d-3
    kap_Type2_full_on_X = 1d-6
    Zbase = 0.02d0


/ ! end of kap namelist

&controls

! WHEN TO STOP
      xa_central_lower_limit_species(1) = 'h1'
      xa_central_lower_limit(1) = 1d-3

! OUTPUT OPTIONS
      ! LOGS directory name
      log_directory = 'LOGS_0.6_DSI_ES_ST'
      
      history_interval       = 1
      terminal_interval      = 1
      write_header_frequency = 10
      photo_digits           = 5
      profile_interval       = 99999 ! Profiles are written out at certain Xc values.
      photo_interval         = 99999

      star_history_name = 'history.data'
      write_controls_info_with_profile = .false.
      write_profiles_flag = .true.
      max_num_profile_models = 250
      
! PHYSICS

       ! Viscosity
       set_uniform_am_nu_non_rot = .false. 
       uniform_am_nu_non_rot = 0d0
        
       ! Rotation
       am_D_mix_factor = 0d0               !0.033d0 commonly used, 0d0 default
       am_gradmu_factor = 0.05d0           !0.05d0 default
       am_nu_factor = 1                    ! 1 default
       am_nu_non_rotation_factor = 1       ! 1 default
       am_nu_omega_rot_factor = 1          ! 1 default
       am_nu_omega_non_rot_factor = 1      ! 1 default
       am_nu_j_rot_factor = 0              ! 0 default
       am_nu_j_non_rot_factor = 0          ! 0 default
       am_nu_visc_factor = 1               ! 1 default
       
       !! AM transport by rotationally-induced processes 
       D_DSI_factor = 1d0   ! dynamical shear instability
       D_SH_factor  = 0d0   ! Solberg-Hoiland
       D_SSI_factor = 0d0   ! secular shear instability
       D_ES_factor  = 1d0   ! Eddington-Sweet circulation
       D_GSF_factor = 0d0   ! Goldreich-Schubert-Fricke
       D_ST_factor  = 1d0   ! Spruit-Tayler dynamo


      ! Mixing
      set_min_D_mix = .true.
      min_D_mix = 1d3

      ! Atmosphere
      atm_option = 'T_tau'
      atm_T_tau_relation = 'Eddington'
      atm_T_tau_opacity = 'varying'

      ! Convection
      mixing_length_alpha = 1.713d0  ! Solar calibrated for Eddington_grey.
      mlt_option = 'Cox'
      
      ! Overshooting
      overshoot_scheme(1) = 'exponential'
      overshoot_zone_type(1) = 'burn_H'
      overshoot_zone_loc(1) = 'core'
      overshoot_bdy_loc(1) = 'top'
      overshoot_f(1) = 0.015d0
      overshoot_f0(1) = 0.005d0

      overshoot_scheme(2) = 'exponential'
      overshoot_zone_type(2) = 'any'
      overshoot_zone_loc(2) = 'shell'
      overshoot_bdy_loc(2) = 'any'
      overshoot_f(2) = 0.005d0
      overshoot_f0(2) = 0.005d0

! MESH AND TIMESTEP PARAMETERS
      mesh_delta_coeff  = 1.0
      varcontrol_target = 1d-3
      max_allowed_nz    = 2500
      time_delta_coeff = 0.5 ! Default (1.0) might be too high for the rotational evolution.

      !! This line was added to solve the rel_err energy warning in MESA r11532.
      energy_eqn_option = 'dedt'
      
      max_number_retries = 1000
      use_Ledoux_criterion = .false.
      
! OVERSHOOTING AND CONVECTION CONTROLS
    overshoot_D_min = 5d-2

    do_conv_premix = .true.



/ ! end of controls namelist

&pgstar


/ ! end of pgstar namelist
