
&star_job
      show_log_description_at_start = .false.

      load_saved_model = .true.
      saved_model_name = 'ready_to_remove_center.mod'

      save_model_when_terminate = .true.
      save_model_filename = 'ready_to_shock.mod'

      remove_initial_center_by_mass_msun = 2.403

      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 = 1d-13

      !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 = 10
         max_number_retries = 20
      
      x_ctrl(1) = -1

      initial_mass = 7
      initial_z = 0.001d0
      Zbase = 0.02d0

      max_model_number = 20

      velocity_logT_lower_bound = 8

      fe_core_infall_limit = 1d9 ! cm/sec

      R_function_weight = 150
      P_function_weight = 40
      T_function1_weight = 110

		xa_function_species(1) = 'h1'
		xa_function_weight(1) = 0 ! 10
		xa_function_param(1) = 1d-3

		xa_function_species(2) = 'he4'
		xa_function_weight(2) = 0 ! 50
		xa_function_param(2) = 1d-3

! No Wind
      RGB_wind_scheme = 'Dutch'
      AGB_wind_scheme = 'Dutch'
      RGB_to_AGB_wind_switch = 1d-4
      Dutch_wind_eta = 0.0d0 ! 0.8

! mixing controls
      mixing_length_alpha = 1.89
      MLT_option = 'Henyey'
      T_mix_limit = 0
      min_T_for_acceleration_limited_conv_velocity = 0
      okay_to_reduce_gradT_excess = .false.

! mesh
      okay_to_remesh = .false.

      photostep = 10000
      profile_interval = 10000
      profile_dbl_format = '(1pes27.16e3, 1x)'  ! '(1pes40.16e3, 1x)'
      profile_int_format = '(i27, 1x)'          ! '(i40, 1x)'
      profile_txt_format = '(a27, 1x)'          ! '(a40, 1x)'

! output
      history_interval = 1
      terminal_interval = 10
      write_header_frequency = 10

!      num_trace_history_values = 0
!      trace_history_value_name(1) = 'split_mixing_choice'
!      trace_history_value_name(2) = 'max_T_lgT'
!      trace_history_value_name(3) = 'max_T_mass'
!      trace_history_value_name(4) = 'burn_max_num_iters'


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


      !dX_nuc_drop_limit = 1d-3


/ ! end of controls namelist

&pgstar


   !pause = .true.

         !Profile_Panels2_win_flag = .true.

      Abundance_xaxis_name = 'radius'
      Abundance_xmin = -101d0 ! only used if /= -101d0
Abundance_xmax = 0.0023 ! only used if /= -101d0

      Profile_Panels1_xaxis_name = 'radius'
      Profile_Panels1_xmin = -101d0 ! only used if /= -101d0
Profile_Panels1_xmax = 0.0023 ! -101d0 ! only used if /= -101d0

         Profile_Panels1_yaxis_name(3) = 'v_div_csound'
         !Profile_Panels1_yaxis_log(3) = .true. ! show log10 of abs value
         Profile_Panels1_ymin(3) = -101d0 ! only used if /= -101d0
         Profile_Panels1_ymax(3) = -101d0 ! only used if /= -101d0
         Profile_Panels1_ymargin(3) = 0.1

         !History_Panels1_win_flag = .true.
         History_Panels1_win_width = 9
         History_Panels1_yaxis_log(3) = .true. ! show log10 of abs value
         History_Panels1_other_yaxis_name(3) = 'star_age'
         History_Panels1_other_yaxis_log(3) = .true. ! show log10 of abs value

!Profile_Panels1_xmax = 2.2 ! -101d0 ! only used if /= -101d0
!Abundance_xmax = 2.2 ! -101 ! only used if /= -101d0


      Grid2_win_flag = .true.

! file output
!Grid2_file_flag = .true.
   Grid2_file_dir = 'png'
   Grid2_file_prefix = 'grid_'
   Grid2_file_cnt = 1 ! 1 ! output when mod(model_number,Grid2_file_cnt)==0
   Grid2_file_width = -1 ! negative means use same value as for window
   Grid2_file_aspect_ratio = -1 ! negative means use same value as for window



/ ! end of pgstar namelist


         Profile_Panels1_num_panels = 3
         Profile_Panels1_yaxis_name(1) = 'velocity'
         Profile_Panels1_ymin(1) = -101d0 ! only used if /= -101d0
         Profile_Panels1_ymax(1) = -101d0 ! only used if /= -101d0
         Profile_Panels1_ymargin(1) = 0

         Profile_Panels1_other_yaxis_name(1) = 'dlnd_dt'
         Profile_Panels1_other_ymin(1) = -1d-3 ! only used if /= -101d0
         Profile_Panels1_other_ymax(1) = 1d-3 ! only used if /= -101d0
         Profile_Panels1_other_ymargin(1) = 0

         Profile_Panels1_yaxis_name(2) = 'div_v'
         Profile_Panels1_ymin(2) = -5d-3  ! only used if /= -101d0
         Profile_Panels1_ymax(2) = 5d-3   ! only used if /= -101d0
         Profile_Panels1_ymargin(2) = 0

         Profile_Panels1_other_yaxis_name(2) = 'd_v_div_r_dr'
         Profile_Panels1_other_ymin(2) = -1d-11 ! only used if /= -101d0
         Profile_Panels1_other_ymax(2) = 1d-11 ! only used if /= -101d0
         Profile_Panels1_other_ymargin(2) = 0

         Profile_Panels1_yaxis_name(3) = 'eps_visc'
         Profile_Panels1_ymin(3) = 0 ! only used if /= -101d0
         Profile_Panels1_ymax(3) = 1d9 ! only used if /= -101d0
         Profile_Panels1_ymargin(3) = 0

         Profile_Panels1_other_yaxis_name(3) = 'dvdt_visc'
         Profile_Panels1_other_ymin(3) = -0.1 ! only used if /= -101d0
         Profile_Panels1_other_ymax(3) = 0.1 ! only used if /= -101d0
         Profile_Panels1_other_ymargin(3) = 0
