

&star_job

      show_log_description_at_start = .true.
      list_net_reactions = .true.
      
      create_pre_main_sequence_model = .true.

      change_initial_net = .true.      
      new_net_name = 'postAGB.net' ! just the superAGB nuclear network + some relevant isotopes
      
      initial_zfracs = 6
      kappa_file_prefix = 'a09'
      
      pgstar_flag = .false.

/ !end of star_job namelist


&controls

! The initial mass of the model is somewhat arbitrary (see paper)
      
      initial_mass = 2.5
      initial_z = 0.02d0

! The following lines are to evolve the models in three steps:

! Step 1:
! Evolve the star until the desired helium core mass is reached
!      he_core_mass_limit = 0.55

! Step 2:
! Once the desired helium core mass is reached, change the stellar wind in run_star_extras to 1d-4 Msun/yr until an envelope mass of 0.02 Msun is left 
!      envelope_mass_limit = 0.02

! Step 3:
! We have now entered to post-AGB phase. We can continue the evolution while changing the parameters for accretion until we reach a temperature of around 50000 K
!      Teff_upper_limit = 5d4



      use_Ledoux_criterion = .true.
      alpha_semiconvection = 4d-2
      thermohaline_coeff = 2  
    


!     These lines increase the resolution in areas with strong nuclear burning

      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


! T_mix_limit is used to mix all layers that are convective down to 80000 K. See Sect. 6.1.3 in Oomen et al. (2019) for more information.
      T_mix_limit = 8d4

      
      use_other_adjust_mdot = .true.

      photo_interval = 1000
      profile_interval = 1000
      history_interval = 10
      terminal_interval = 10
      write_header_frequency = 100


! Reducing the gradT excess is only used for massive post-AGB stars of around 0.65 Msun to avoid convergence issues. See appendix C in Miller Bertolami (2016) for more information.

      okay_to_reduce_gradT_excess = .false.

            ! these are for calculation of efficiency boosted gradT
            gradT_excess_f1 = 1d-3
            gradT_excess_f2 = 1d-2
            gradT_excess_age_fraction = 0.9d0
            
            ! turn on full all the time
            gradT_excess_lambda1 = -1.0


/ ! end of controls namelist
