
! inlist_rotation


&star_job
      show_log_description_at_start = .false.

      load_saved_model = .true.
      saved_model_name = 'he_burn.mod'

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

      eos_file_prefix = 'mesa'
      kappa_file_prefix = 'gs98'
               
!      change_v_flag = .true.
!      new_v_flag = .true.
!               
!      change_rotation_flag = .true.
!      new_rotation_flag = .true.
!      
!      set_surface_rotation_v = .true.
!      new_surface_rotation_v = 200 ! km/sec   << saved model is rotating 

      !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 = 0
            max_number_retries = 0


         initial_mass = 15
         initial_z = 2d-2

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

         am_nu_visc_factor = 0
         am_D_mix_factor = 0.0333333333333333d00
         D_DSI_factor = 0
         D_SH_factor = 1
         D_SSI_factor = 1
         D_ES_factor = 1
         D_GSF_factor = 1
         D_ST_factor = 1

         log_center_density_limit = 12 
         log_center_temp_limit = 10.3d0 ! max T about 2d10

         photostep = 50
         profile_interval = 50
         history_interval = 1
         terminal_cnt = 10
         write_header_frequency = 10
         
         max_model_number = 2001

         xa_central_lower_limit_species(1) = 'he4'
         xa_central_lower_limit(1) = 0.6                       
      
         use_Ledoux_criterion = .true.
   
         alpha_semiconvection = 4d-2
         thermohaline_coeff = 2  

         dX_mix_dist_limit = 1d-4 
         
         omega_function_weight = 20
         
         iter_for_resid_tol2 = 4
         iter_for_resid_tol3 = 7 


         mdot_omega_power = 0.43

         clip_D_limit = 10
         
         xa_scale = 1d-5
         newton_itermin = 2

         tol_correction_norm_alt = 1d-3
         tol_max_correction_alt = 1d0

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

         mixing_length_alpha = 1.6
         
         delta_lgL_He_limit = -1

         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

         min_overshoot_q = 1d-3 

         T_mix_limit = 1d4

      ! 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
      !hydro_show_correction_info = .true.
      
      !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 = 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 = 1

      !trace_newton_bcyclic_solve_input = .true. ! input is "B" j k iter B(j,k)
      !trace_newton_bcyclic_solve_output = .true. ! output is "X" j k iter X(j,k)
      
      trace_newton_bcyclic_steplo = 1 ! 1st model number to trace
      trace_newton_bcyclic_stephi = 1 ! last model number to trace
      
      trace_newton_bcyclic_iterlo = 0 ! 1st newton iter to trace
      trace_newton_bcyclic_iterhi = 1 ! last newton iter to trace
      
      trace_newton_bcyclic_nzlo = 1 ! 1st cell to trace
      trace_newton_bcyclic_nzhi = 10000 ! last cell to trace; if < 0, then use nz as nzhi
      
      trace_newton_bcyclic_jlo = 1 ! 1st var to trace
      trace_newton_bcyclic_jhi = 100 ! last var to trace; if < 0, then use nvar as jhi
      
      !trace_k = 1

/ ! end of controls namelist



&pgstar
         

/ ! end of pgstar namelist
