! inlist_axion_cooling


&star_job

         show_log_description_at_start = .false. 

         load_saved_model = .true.
         saved_model_name = 'axion_test.mod'

         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.
         
         kappa_file_prefix = 'gs98'


/ ! 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 = 50
         max_number_retries = 100
         
         max_model_number = 50000
      
         smooth_convective_bdy = .true.                  
         convective_bdy_weight = 1

         use_other_neu = .true.

         initial_mass = 10.5
         initial_z = 2d-2


         photostep = 50
         profile_interval = 50
         history_interval = 1
         terminal_interval = 10
         write_header_frequency = 10

         
         num_trace_history_values = 0
         trace_history_value_name(1) = 'power_neutrinos'

      prune_bad_cz_min_Hp_height = 0 ! lower limit on radial extent of cz
      remove_mixing_glitches = .true. ! if true, then okay to remove gaps and singletons
      
               ! the following controls are for different kinds of "glitches" that can be removed
            
               clip_D_limit = 0 ! zero mixing diffusion coeffs that are smaller than this         

               okay_to_remove_mixing_singleton = .true.
         
               min_convective_gap = -1 ! 0.1 
                  ! close gap between convective regions if smaller than this (< 0 means skip this)
                  ! gap measured radially in units of pressure scale height
         
               min_thermohaline_gap = 0
                  ! close gap between thermohaline mixing regions if smaller than this (< 0 means skip this)
                  ! gap measured radially in units of pressure scale height
               
               min_thermohaline_dropout = 0
               max_dropout_gradL_sub_grada = 1d-3
                  ! if find radiative region embedded in thermohaline,
                  ! and max(gradL - grada) in region is everywhere < max_dropout_gradL_sub_grada
                  ! and region height is < min_thermohaline_dropout
                  ! then convert the region to thermohaline
         
               min_semiconvection_gap = 0 
                  ! close gap between semiconvective mixing regions if smaller than this (< 0 means skip this)
                  ! gap measured radially in units of pressure scale height
         
               remove_embedded_semiconvection = .false.
                  ! if have a semiconvection region bounded on each side by convection,
                  ! convert it to be convective too.
         
         
         ! DISABLE THESE
         !sum_xa_tolerance = 1d99
         !min_xa_hard_limit = -1d99 
         
         mesh_delta_coeff = 2         

         RGB_wind_scheme = 'Dutch'
         AGB_wind_scheme = 'Dutch'
         RGB_to_AGB_wind_switch = 1d-4
         Dutch_wind_eta = 0.8

         log_center_density_limit = 12 
         log_center_temp_limit = 10.3d0 ! max T about 2d10
         
         xa_scale = 1d-5
         newton_itermin = 2

         !mass_change_full_on_dt = 1d7 ! (seconds)
         !mass_change_full_off_dt = 1d6 ! (seconds)

         mixing_length_alpha = 1.6  ! 1.0 for M > 30

         xa_central_upper_limit_species(1) = 'c12'
         xa_central_upper_limit(1) = 0.2

         max_num_profile_models = 10000
         
         !xa_central_upper_limit_species(1) = 'cr56'
         !xa_central_upper_limit(1) = 0.80
         
         tol_residual_norm1 = 1d-5
         tol_max_residual1 = 1d-2
         tol_residual_norm2 = 1d99
         tol_max_residual2 = 1d99
         iter_for_resid_tol2 = 3
         iter_for_resid_tol3 = 12

         !correction_xa_limit = 1d-4
         
         min_timestep_limit = 1d-12 ! (seconds)

         !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 = .true.
         !report_all_dt_limits = .true.
         
         !show_mesh_changes = .true.
         !mesh_dump_call_number = 1598
         
         !okay_to_remesh = .false.
         
         !trace_evolve = .true.
         
         ! hydro debugging
         !hydro_numerical_jacobian = .true.
         !hydro_check_everything = .true.
         !hydro_show_rcond = .true.
         !hydro_inspectB_flag = .true.
         !hydro_save_numjac_plot_data = .true.
         !hydro_dump_call_number = 1

         use_artificial_viscosity = .false.
         
         delta_lgL_He_limit = -1



         !dX_nuc_max_A_limit = 52
         !dX_limit_min_X = 1d-3
         !dX_limit = 2d-3         
         !dX_nuc_limit = 2d-3
         !dX_nuc_factor = 5d2
         !dX_nuc_min_yrs_for_dt = 5d-10


         dX_nuc_drop_max_A_limit = 52
         dX_nuc_drop_min_X_limit = 1d-4
         dX_nuc_drop_limit = 5d-3
         dX_nuc_drop_hard_limit = 1d99
         
         
         
         
         
         delta_lgTeff_limit = 0.005  
         delta_lgL_limit = 0.05
         delta_lgRho_cntr_limit = 0.02  


         mesh_dlog_pp_dlogP_extra = 0.15
         mesh_dlog_cno_dlogP_extra = 0.15
      
         mesh_dlog_3alf_dlogP_extra = 0.15
         mesh_dlog_burn_c_dlogP_extra = 0.15
         
         mesh_dlog_burn_n_dlogP_extra = 0.15
         mesh_dlog_burn_o_dlogP_extra = 0.15
      
         mesh_dlog_burn_ne_dlogP_extra = 0.15
         
         mesh_dlog_burn_na_dlogP_extra = 0.15
         mesh_dlog_burn_mg_dlogP_extra = 0.15
         
         mesh_dlog_cc_dlogP_extra = 0.15
         mesh_dlog_co_dlogP_extra = 0.15
         
         mesh_dlog_oo_dlogP_extra = 0.15
         
         mesh_dlog_burn_si_dlogP_extra = 0.15
         mesh_dlog_burn_s_dlogP_extra = 0.15
         mesh_dlog_burn_ar_dlogP_extra = 0.15
         mesh_dlog_burn_ca_dlogP_extra = 0.15
         mesh_dlog_burn_ti_dlogP_extra = 0.15
         mesh_dlog_burn_cr_dlogP_extra = 0.15
         mesh_dlog_burn_fe_dlogP_extra = 0.15
         
         mesh_dlog_pnhe4_dlogP_extra = 0.15
         mesh_dlog_other_dlogP_extra = 0.15
         mesh_dlog_photo_dlogP_extra = 1
         
         
         xa_function_species(1) = 'ne20'
   		xa_function_weight(1) = 20
   		xa_function_param(1) = 1d-2
         
         xa_function_species(1) = 'si28'
   		xa_function_weight(1) = 20
   		xa_function_param(1) = 1d-2
         
         xa_function_species(1) = 's32'
   		xa_function_weight(1) = 20
   		xa_function_param(1) = 1d-2
         
         xa_function_species(1) = 'fe54'
   		xa_function_weight(1) = 20
   		xa_function_param(1) = 1d-2
         
         xa_function_species(1) = 'fe56'
   		xa_function_weight(1) = 20
   		xa_function_param(1) = 1d-2

 
         overshoot_f_above_nonburn_core = 0.01
         overshoot_f_above_nonburn_shell = 0.01
         overshoot_f_below_nonburn_shell = 0.01
         overshoot_f_above_burn_h_core = 0.01
         overshoot_f_above_burn_h_shell = 0.01
         overshoot_f_below_burn_h_shell = 0.01
 
         overshoot_f0_above_nonburn_core = 0.005
         overshoot_f0_above_nonburn_shell = 0.005
         overshoot_f0_below_nonburn_shell = 0.005
         overshoot_f0_above_burn_h_core = 0.005
         overshoot_f0_above_burn_h_shell = 0.005
         overshoot_f0_below_burn_h_shell = 0.005


         min_overshoot_q = 1d-3 

         T_mix_limit = 1d4


/ ! end of controls namelist




&pgstar
         
      ! top level controls

         !pause = .true. 
            ! if true, the code waits for user to enter a RETURN on the command line
      


/ ! end of pgstar namelist
