
&star_job

      ! set up initial profile
      show_log_description_at_start = .false.

      ! set up output file
      save_model_when_terminate = .true.
      save_model_filename = 'ready_to_bounce.mod'
 
      ! set up end program record     
      write_profile_when_terminate = .true.
      filename_for_profile_when_terminate = 'final_profile.data'

      ! initial profile
      !create_pre_main_sequence_model = .true.
      load_saved_model = .true.
      saved_model_name = 'M140_Z0004_Omega0.mod'

      ! no rotation needed
      !new_rotation_flag = .true.
      !change_rotation_flag = .false.
 
      ! use simpler network
      change_net = .true.      
      new_net_name = 'approx21_plus_co56.net'

      ! with dynamical velocity
      change_v_flag = .true.

      ! relax to the required mass 
      relax_initial_mass = .true.
      new_mass = 60
      !lg_max_abs_mdot = -2

      ! relax to the required metallicity
      !relax_initial_Z = .true. 
      !new_Z = 0.01D0

      ! relax to a given composition
      file_for_uniform_xa = 'xini2.dat'
      set_uniform_initial_xa_from_file = .true.

      !pgstar_flag = .true.

/ ! end of star_job namelist

&controls

      ! for testing 
      !MLT_option = 'none'	! default mute this command
      !mixing_length_alpha = 2e-5  ! default 1.5

      ! 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 = 500
      max_model_number = 10000

      ! default output
      photo_directory = 'photos'
      log_directory = 'LOGS'

      log_center_temp_limit = 9.4
      !min_timestep_limit = 1d-0
 
      ! some default muted statements
      !no_wind_if_no_rotation = .true.
   
      ! initial mass setting
      initial_mass = 140
      initial_z = 4d-3
      Zbase = 4d-3

      ! default
      use_Type2_opacities = .true.
      
      ! muted to slightly delay the termination time
      !non_fe_core_infall_limit = 6d7 ! cm/sec
      
      ! for producing a consistent model with the H-star
      hot_wind_scheme = 'Dutch'
      Dutch_scaling_factor = 0.8
      hot_wind_full_on_T = 1d0
      hot_wind_full_off_T = 0d0

      ! more on atmosphere
      !which_atm_option = 'grey_and_kap'
      !Pextra_factor = -1

      ! output parameters
      photostep = 50
      profile_interval = 50   ! 20
      history_interval = 10 ! 10
      terminal_cnt = 10
      write_header_frequency = 10

      ! timestep controls      
      dX_nuc_drop_limit = 2d-3

      ! mesh controls     
      okay_to_remesh = .true.     
      mesh_delta_coeff = 1.0 ! 0.6 for production
      mesh_delta_coeff_for_highT  = 1.0 ! 0.6 for production
      !min_dq = 1d-4
      max_center_cell_dq = 1d-4
      max_surface_cell_dq = 1d-4

      ! original but muted inputs
      !mach1_minus_dr_factor = 1d-2                                  
      !mach1_plus_dr_factor = 0d0
      !tol_correction_norm = 1d-4
      !tol_max_correction = 1d-3

      ! stablize
      corr_coeff_limit = 1d0 ! else dies from failed attempts to use tiny corrections

      ! equation controls (muted because it is quasi hydrostatic)
      !use_ODE_var_eqn_pairing = .true.
      !use_dvdt_form_of_momentum_eqn = .true.
      !use_dPrad_dm_form_of_T_gradient_eqn = .true.
      !use_dedt_form_of_energy_eqn = .true.
      !use_momentum_outer_BC = .true.      
      !use_ODE_form_of_density_eqn = .true.     ! no diffusion of rho unless this is true

      ! 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

      !report_why_dt_limits = .true.
      !report_all_dt_limits = .true.
      !report_hydro_dt_info = .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 = 195


/ ! end of controls namelist

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

/ ! end of pgstar namelist
