
! note: this inlist is done after 1st doing inlist_massive_defaults
! so can override defaults as desired.


! relax composition to desired uniform values

! for masses >= 8 Msun

! set "initial_mass" to the desired mass
! set "accretion_species" to the desired composition


&star_job

      relax_initial_to_xaccrete = .true.
      num_steps_to_relax_composition = 50

      save_model_when_terminate = .true.
      save_model_filename = 'final.mod'
      
      write_profile_when_terminate = .true.
      filename_for_profile_when_terminate = 'final_profile.data'

      change_initial_net = .true.      
      new_net_name = 'approx21_plus_co56.net'

      !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 = 150
         max_number_retries = 600
         
         max_model_number = 9000
         
         initial_mass = 12
         initial_z = 0.02d0
         Zbase = 0.02d0
         
         
         

         varcontrol_target = 1d-3
         dX_nuc_drop_limit = 2d-2 ! 5d-3

         delta_lgT_cntr_limit = 0.01  
            ! limit for magnitude of change in log10 temperature at center   
         delta_lgT_cntr_hard_limit = 0.03

         delta_lgRho_cntr_limit = 0.05  
            ! limit for magnitude of change in log10 density at center   
         delta_lgRho_cntr_hard_limit = 0.15




         RGB_wind_scheme = ''
         AGB_wind_scheme = ''

         accrete_same_as_surface = .false. 
         accrete_given_mass_fractions = .true. 

         ! PURE HE
         ! set target uniform composition for relax_initial_to_xaccrete
!         num_accretion_species = 1
!         accretion_species_id(1) = 'he4'
!         accretion_species_xa(1) = 1d0


! 12C = −0.035 × ln(MC/O) + 0.253  (KEPLER)
         
         ! C/O MIX
         ! set target uniform composition for relax_initial_to_xaccrete
         num_accretion_species = 2
         accretion_species_id(1) = 'c12'
         accretion_species_xa(1) = 0.1582d0
         accretion_species_id(2) = 'o16'
         accretion_species_xa(2) = 0.8418d0


         photostep = 10
         profile_interval = 10
         history_interval = 1
         terminal_cnt = 2
         write_header_frequency = 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
      
      !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 = 2005
      !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 = 5612

/ ! end of controls namelist


&pgstar
   !clear_history = .true.
         
!pause = .true.
         
Profile_Panels1_xmax = 2.1 ! -101d0 ! only used if /= -101d0
Abundance_xmax = 2.1 ! -101 ! only used if /= -101d0

/ ! end of pgstar namelist
