&star_job

      load_saved_model = .true.
      saved_model_name = 'shock_part4.mod'

      save_model_when_terminate = .true.
      save_model_filename = 'shock_part5.mod'

      limit_initial_dt = .true.
      seconds_for_initial_dt = 1d-1

      save_stella_data_when_terminate = .true.
      save_stella_data_filename = 'mesa'
      stella_show_headers = .false.
      stella_min_surf_logRho = -9d0
      stella_min_velocity = 5d7
      stella_skip_inner_dm = 0.01d0 !default is 0.01, in Msun
      stella_num_points = 600 ! number of zones in STELLA (set this to 700 when including CSM)
      stella_nz_extra = 0 !set to 0 for no csm, 30 is default (set this to 100 when including CSM)
      stella_mdot_years_for_wind = 0.01d0 !wind mass-loss duration (vary this when making a CSM grid)
      stella_mdot_for_wind = 0.30d0 !wind mass-loss rate (vary this when making a CSM grid)
      stella_v_wind = 10d0 !wind velocity

      !pgstar_flag = .false.

/ ! end of star_job namelist

&controls

      !max_number_backups = 50
      !max_number_retries = 250
      !max_model_number = 20000
         
      x_integer_ctrl(1) = 5 ! part number
      
      ! x_ctrl(16) = 0.05 ! stop when shock reaches this mass below surface
      
      ! use these if want to go to breakout
      !log_L_upper_limit = 8.8 ! 10
      x_ctrl(16) = -1 ! if < 0, run until hit some limit
      max_model_number = 1

      x_integer_ctrl(3) = 3 ! number of boxcar passes
      x_ctrl(47) = 0.8  ! width of boxcar (Msun)
      
      ! extent of mixing region is set to entire model
      x_ctrl(45) = 0.0
      x_ctrl(46) = 90.0

      ! set final Ni mass
      x_ctrl(12) = 0.019d0
      x_ctrl(35) = -1d99 ! this means rescale Ni in place
      x_ctrl(36) = -1d99
      x_ctrl(37) = 5d7 ! will discard center cells with u < this (cm/s), so don't count for Ni56.

      x_ctrl(55) = 5d7 ! for velocity cut from Jared's

! mesh controls
      split_merge_amr_okay_to_split_nz = .false.
      split_merge_amr_log_zoning = .false.
      split_merge_amr_nz_baseline = 1600

! RTI

! equation controls
      use_compression_outer_BC = .true.

! timestep controls
      delta_lgL_limit = 0.01d0

! solver controls

   use_gold_tolerances = .true.

! output controls
      num_trace_history_values = 9
      trace_history_value_name(1) = 'rel_E_err'
      trace_history_value_name(2) = 'log_rel_run_E_err'
      trace_history_value_name(3) = 'total_mass ni56'
      trace_history_value_name(4) = 'total_ni_co_56'
      trace_history_value_name(5) = 'm_center'
      trace_history_value_name(6) = 'shock_mass'
      trace_history_value_name(7) = 'total_energy'
      trace_history_value_name(8) = 'shock_radius'
      trace_history_value_name(9) = 'photosphere_r'

      !photo_interval = 1
      !profile_interval = 1
      !history_interval = 1
      !terminal_interval = 1


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

      !fill_arrays_with_NaNs = .true.
      !stop_for_NaNs = .true.

      !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_sizequ_flag = .true.
      
      !hydro_get_a_numerical_partial = 1d-3
      !hydro_test_partials_k = 798
      
      !hydro_numerical_jacobian = .true.
      !hydro_save_numjac_plot_data = .true.
      !hydro_dump_call_number = 617
      !hydro_dump_iter_number = 0
      !hydro_epsder_struct = 1d-8
      !hydro_epsder_chem = 1d-6

      !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_matrix_input = .true.
      !trace_newton_bcyclic_matrix_output = .true.
      
      !trace_newton_bcyclic_steplo = 1 ! 1st model number to trace
      !trace_newton_bcyclic_stephi = 1 ! last model number to trace
      
      !trace_newton_bcyclic_iterlo = 2 ! 1st newton iter to trace
      !trace_newton_bcyclic_iterhi = 2 ! 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 = 0


/ ! end of controls namelist

&pgstar
         
!pause = .true.
         
Abundance_xmin = -101
Mixing_xmin = -101
Profile_Panels1_xmin = -101d0
Profile_Panels2_xmin = -101d0

Profile_Panels1_yaxis_name(2) = 'v_div_csound'
Profile_Panels1_ymax(2) = -101
Profile_Panels1_ymin(2) = 0 ! -101
Profile_Panels1_other_yaxis_name(2) = 'csound'
Profile_Panels1_other_ymax(2) = -101
Profile_Panels1_other_ymin(2) = -101

Profile_Panels2_yaxis_name(2) = 'luminosity'
Profile_Panels2_ymin(2) = -101
Profile_Panels2_other_yaxis_name(2) = 'log_dr_div_cs'
Profile_Panels2_other_ymin(2) = -101

         Text_Summary1_name(5,3) = 'v_div_csound_surf'


      Grid2_plot_name(4) = 'History_Track1'        
         History_Track1_title = 'Breakout Light Curve'
         History_Track1_xname = 'star_age_hr'
         History_Track1_yname = 'log_L'
         History_Track1_xaxis_label = 'hours'
         History_Track1_yaxis_label = 'log L/L\d\(2281)'
         History_Track1_reverse_xaxis = .false.
         History_Track1_reverse_yaxis = .false.
         History_Track1_log_xaxis = .false. ! show log10 of abs value
         History_Track1_log_yaxis = .false. ! show log10 of abs value

/ ! end of pgstar namelist
