
! inlist_sdB


&star_job

      show_log_description_at_start = .false. 

      !save_model_when_terminate = .true.
      !save_model_filename = 'test_1m2_envelope.mod'

      load_saved_model = .true.
      saved_model_name = 'test_1m2_envelope.mod'
      !saved_model_name = 'sdb_from_1M.mod'
      !saved_model_name = 'he_ignition_peak.mod'
      !saved_model_name = '1M_he_ignition.mod'

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

      eos_file_prefix = 'mesa'
      kappa_file_prefix = 'OP_gs98'
   
      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 = 1d4 ! in years

      !set_tau_factor = .true.  
      !set_to_this_tau_factor = 150

      change_rotation_flag = .true.
      new_rotation_flag = .false.

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

      !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_retries = 40
      !max_number_backups = 20
      
      !max_model_number = 1111

      initial_mass = 1
      initial_z = 0.02d0
      
      mixing_length_alpha = 2
      use_Henyey_MLT = .true.

      varcontrol_target = 5d-4 ! 1d-5
      mesh_delta_coeff = 0.75 ! 1.5
      
      !min_dq = 1d-10
      !max_surface_cell_dq = 1d-14
      
      photostep = 10
      profile_interval = 10
      history_interval = 1
      terminal_cnt = 1
      write_header_frequency = 10      
      
      !dH_div_H_limit_min_H = 1d99 ! disable this
      !dHe_div_He_limit_min_He = 1d99 ! disable this

      !prune_bad_cz_min_Hp_height = 0.1 ! lower limit on radial extent of cz
      
      !overshoot_f_above_nonburn = 0.01
      !overshoot_f_below_nonburn = 0.01
      
      !max_years_for_timestep = 1d4
      
      ! element diffusion
         
         show_diffusion_info = .true.
         !show_diffusion_substep_info = .true.         
         !show_diffusion_timing = .true.
         
         
!      surface_mix_log_D0 = 4
!      surface_mix_logRho_lo_full_on = 0.1d0
!      surface_mix_logRho_hi_full_off = 0.5d0
      
      do_element_diffusion = .true. ! determines whether or not we do diffusion
      
      diffusion_calculates_ionization = .true.


      diffusion_num_classes = 4 ! number of classes of species for diffusion calculations
      diffusion_class_representative(1) = 'h1'
      diffusion_class_representative(2) = 'he4'
      diffusion_class_representative(3) = 'o16'
      diffusion_class_representative(4) = 'fe56'

      ! in ascending order.  species goes into 1st class with A_max >= species A
      diffusion_class_A_max(1) = 2
      diffusion_class_A_max(2) = 4
      diffusion_class_A_max(3) = 16
      diffusion_class_A_max(4) = 10000

      diffusion_nsmooth_v_advection = 0
      diffusion_nsmooth_x_in_fixup = 2
      diffusion_max_logT_for_smooth_x = 99
      
      !diffusion_use_isolve = .true.
      diffusion_rtol_for_isolve = 1d-4
      diffusion_atol_for_isolve = 1d-5
      diffusion_maxsteps_for_isolve = 1000
      diffusion_isolve_solver = 'ros2_solver'
         
!      envelope_mass_limit = 1d-2 ! 4d-4
      when_to_stop_rtol = 1d-3
      when_to_stop_atol = 1d-5
         
!      RGB_wind_scheme = 'Reimers'
!      AGB_wind_scheme = 'Blocker'
      RGB_to_AGB_wind_switch = 0.75d0
      Reimers_wind_eta = 1d3 ! 0.1d0  
      Blocker_wind_eta = 0.2d0  

      num_cells_for_smooth_brunt_B = 10
         
      use_Type2_opacities = .true.
   
      ! switch to Type1 if X too large
      kap_Type2_full_off_X = 0.71d0 ! Type2 full off for X >= this (< 0 means use kap default)
      kap_Type2_full_on_X = 0.69d0 ! Type2 full on for X <= this (< 0 means use kap default)


      ! 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
      
      !max_years_for_timestep = 1d6

      !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




/ ! end of controls namelist



&pgstar
         
         

/ ! end of pgstar namelist
