! inlist_agb


&star_job
      show_log_description_at_start = .false.

      mesa_data_dir = '../../../data' ! for eos and opacity tables
      
      !create_pre_main_sequence_model = .true.

      !save_model_number = 20
      save_model_filename = '2m_z1m2_agb.mod'

      load_saved_model = .true.
      saved_model_name = '2m_z1m2_agb.mod'

      change_net = .true.      
      new_net_name = 'basic.net'

      !set_initial_dt = .true.
      years_for_initial_dt = 1d-8 ! in years


      log_columns_file = 'log_columns.list' 


      !change_lnPgas_flag = .true.
      !new_lnPgas_flag = .false.
      !new_lnPgas_flag = .true.


/ ! end of star_job namelist

&controls

      D_norm_kappa = 1
      D_norm_alpha = 2
      D_norm_atol = 1d-3
      D_norm_rtol = 1d-3
   
      use_lnE_for_eps_grav = .true.


      initial_mass = 2.0
      initial_z = 0.01d0

      ! check for retries and backups as part of test_suite
      max_number_backups = 10
      max_number_retries = 80
      
      xa_surface_upper_limit_species(1) = 'c12'
      xa_surface_upper_limit(1) = 0.0013

      max_model_number = 8000

      photostep = 100
      profile_interval = 100
      log_cnt = 100
      terminal_cnt = 10
      write_header_frequency = 10
         
      !min_xa_hard_limit = -1d-10 

      !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

      report_why_dt_limits = .false.
      report_all_dt_limits = .false.
      
      !show_mesh_changes = .true.
      !mesh_dump_call_number = 8875
      !okay_to_remesh = .false.
      

      trace_evolve = .false.
      
      ! 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
      
 
      mesh_dlog_pp_dlogP_extra = 0.25
      mesh_dlog_cno_dlogP_extra = 0.25
      
      mesh_dlog_3alf_dlogP_extra = 0.225
      mesh_dlog_burn_c_dlogP_extra = 0.225
      mesh_dlog_burn_n_dlogP_extra = 0.225
      mesh_dlog_burn_o_dlogP_extra = 0.225
      
      
      mesh_logX_species(1) = 'h1'
      mesh_logX_min_for_extra(1) = -6
      mesh_dlogX_dlogP_extra(1) = 0.25
      
      mesh_logX_species(1) = 'he4'
      mesh_logX_min_for_extra(1) = -6
      mesh_dlogX_dlogP_extra(1) = 0.25

         
         
      ! multiply mesh_delta_coeff in overshooting regions by this factor
      ! value < 1 gives increased resolution
      xtra_coef_os_above_nonburn = 0.5
      xtra_coef_os_below_nonburn = 0.25
   
      xtra_coef_os_above_burn_h = 0.25
      xtra_coef_os_below_burn_h = 0.15
   
      xtra_coef_os_above_burn_he = 0.25
      xtra_coef_os_below_burn_he = 0.25
   
      xtra_coef_os_above_burn_z = 0.25
      xtra_coef_os_below_burn_z = 0.25

      ! continue to increase resolution for this distance beyond
      ! the edge of the overshooting region, with distance measured in units
      ! of the pressure scale height at the edge of the overshooting region.
      ! this applies to both edges of the overshooting region.
      xtra_dist_os_above_nonburn = 0.1
      xtra_dist_os_below_nonburn = 0.2
   
      xtra_dist_os_above_burn_h = 0.2
      xtra_dist_os_below_burn_h = 0.2
   
      xtra_dist_os_above_burn_he = 0.2
      xtra_dist_os_below_burn_he = 0.2
   
      xtra_dist_os_above_burn_z = 0.2
      xtra_dist_os_below_burn_z = 0.2


      overshoot_f_above_nonburn = 0.014
      overshoot_f_below_nonburn = 0.014
      overshoot_f_above_burn_h = 0.014
      overshoot_f_below_burn_h = 0.014
      overshoot_f_above_burn_he = 0.014
      overshoot_f_below_burn_he = 0.014    
      overshoot_f_above_burn_z = 0.014
      overshoot_f_below_burn_z = 0.014    
      
      overshoot_below_noburn_factor = 10
              
      
      RGB_wind_scheme = ''
      AGB_wind_scheme = ''
      RGB_to_AGB_wind_switch = 1d-4
      Reimers_wind_eta = 0.5d0  
      Blocker_wind_eta = 0.1d0  
      
      use_CO_enhanced_opacities = .true.
      base_Z = 0.01
      base_fC = 0.17123643187d0 ! base value of XC/Z
      base_fO = 0.4812175d0 ! base value of XO/Z


/ ! end of controls namelist

