
&star_job

   save_model_when_terminate = .true.
   save_model_filename = 'final.mod'

   ! no H, so do blend higher up
   kappa_blend_logT_upper_bdy = 4.1d0
   kappa_blend_logT_lower_bdy = 4.0d0   

  write_profile_when_terminate = .true.
  filename_for_profile_when_terminate = 'final_prof.prof'

   pgstar_flag = .true.
   save_pgstar_files_when_terminate = .true.
   create_pre_main_sequence_model = .true.

   timescale_for_relax_composition = 1d5

   !change_initial_net = .true.
   new_net_name = "approx21.net"

   ! I re-implement this to turn on u_flag
   ! check run_star_extras
   gamma1_integral_for_v_flag = -1d50

   ! this are options for the relaxation of the model in-between pulses
   timescale_for_relax_entropy = 1d-15
   max_dt_for_relax_entropy = 5d-17
   num_timescales_for_relax_entropy = 30
   max_steps_to_relax_entropy = 10000

   change_initial_conv_vel_flag = .true.
   new_conv_vel_flag = .true.

   set_initial_cumulative_energy_error = .true.
   new_cumulative_energy_error = 0d0

   change_RTI_flag = .true.
   new_RTI_flag = .false.

   num_special_rate_factors=0
   reaction_for_special_factor(1)='r_c12_ag_o16'
   special_rate_factor(1)=1.0d0 ! Used for linear scaling of rates

   ! rate_tables_dir='../rate_tables' ! Used when using STARLIB or other custom rate library

   ! pre_ms_T_c =3.1e5 ! Needed for low Z models to get through create_pre_main_sequence

/ !end of star_job namelist


&controls

   relax_max_number_retries = -1
   relax_max_number_backups = -1
   photo_interval=100000
   okay_to_remove_mixing_singleton = .false.
   smooth_convective_bdy = .false.

   num_trace_history_values = 2
   trace_history_value_name(1) = 'log_rel_run_E_err'
   trace_history_value_name(2) = 'rel_E_err'

   warn_when_large_rel_run_E_err = 1d-2

   use_dedt_form_of_energy_eqn = .true.

   ! adjustements to the newton solver
   scale_max_correction = 1d0
   conv_vel_ignore_coeff = .true.
   conv_vel_D = 0d0
   relax_use_gold_tolerances = .false.
   use_gold_tolerances = .false.
   maxT_for_gold_tolerances = 5d9 ! this is to ease core-collapse (~logT=9.7)
                                       ! deactivated for test_suite as it only does one pulse
                                       ! RECOMMENDED 5d9
   gold_iter_for_resid_tol2 = 10
   gold_newton_iterations_limit = 15
   newton_iterations_limit = 50
   gold_iter_for_resid_tol3 = 10
   gold_tol_residual_norm3 = 1d-6
   gold_tol_max_residual3 = 1d-3

   use_eosDT2 = .true.
   use_eosELM = .true.

   max_model_number = 100000 ! RECOMMENDED -1
   max_number_retries = -1 ! RECOMMENDED 4000
   max_number_backups = -1 ! RECOMMENDED 2000
   relax_max_number_retries = -1
   relax_max_number_backups = -1
   read_extra_controls_inlist1 = .true.
   extra_controls_inlist1_name = 'inlist_extra'

    initial_z=1d-3
    initial_y=0.999d0
    initial_he3 = 0d0

   ! in principle this is the only thing that needs changing
   ! it is set in inlist_extra
   !initial_mass = 55
 
   ! We use 10% solar, with the definition of Asplund+ 2009
   use_Type2_opacities = .true.
   zbase=1d-3 ! Change with initial_z

   ! our wind implementation follows Brott+ 2011
   use_other_wind = .true.

   ! convection controls
   use_ledoux_criterion = .true.
   mixing_length_alpha = 2d0
   alpha_semiconvection = 1d0
   thermohaline_coeff = 0d0
   num_cells_for_smooth_gradL_composition_term = 0
   !for some reason this crashes the initial model creation
   !its turned on in run_star_extras instead
   !max_conv_vel_div_csound = 1d0
   max_abs_du_div_cs_for_convection = 0.03d0
   max_v_div_cs_for_convection = 1d2
   max_v_for_convection = 1d4

   ! eos controls
   logQ_min_limit = -1d99
   logQ_limit = 5d0

   ! overshoot controls
   ! we only include a bit of exponential overshooting to smooth things out
   overshoot_f0_above_nonburn_core = 0.005         
   overshoot_f0_above_burn_h_core = 0.005
   overshoot_f0_above_burn_he_core = 0.005
   overshoot_f0_above_burn_z_core = 0.005

   overshoot_f_above_nonburn_core = 0.01
   overshoot_f_above_burn_h_core = 0.01
   overshoot_f_above_burn_he_core = 0.01
   overshoot_f_above_burn_z_core = 0.01

   overshoot_f0_above_nonburn_shell = 0.005         
   overshoot_f0_above_burn_h_shell = 0.005
   overshoot_f0_above_burn_he_shell = 0.005
   overshoot_f0_above_burn_z_shell = 0.005

   overshoot_f_above_nonburn_shell = 0.01
   overshoot_f_above_burn_h_shell = 0.01
   overshoot_f_above_burn_he_shell = 0.01
   overshoot_f_above_burn_z_shell = 0.01

   overshoot_f0_below_nonburn_shell = 0.005         
   overshoot_f0_below_burn_h_shell = 0.005
   overshoot_f0_below_burn_he_shell = 0.005
   overshoot_f0_below_burn_z_shell = 0.005

   overshoot_f_below_nonburn_shell = 0.01
   overshoot_f_below_burn_h_shell = 0.01
   overshoot_f_below_burn_he_shell = 0.01
   overshoot_f_below_burn_z_shell = 0.01

   ! termination conditions
   ! make this a bit higher to avoid interrupting pair-instability collapse
   fe_core_infall_limit = 8d8

   ! timestep options
   varcontrol_target = 5d-4
   max_timestep_factor = 1.05d0
   min_timestep_factor = 0.9d0
   dX_nuc_drop_limit = 5d-2
   dX_nuc_drop_limit_at_high_T = 1d-2 ! for center logT > 9.45
   delta_Ye_highT_limit = 1d-3
   dX_nuc_drop_max_A_limit = 52
   dX_nuc_drop_min_X_limit = 1d-4
   dX_nuc_drop_hard_limit = 1d99
   delta_lgTeff_limit = 1d0
   delta_lgL_phot_limit = -1d0
   delta_lgL_limit = -1d0
   delta_lgL_He_limit = -1d0
   lgL_nuc_burn_min = 4d0

   ! extra controls for timestep
   delta_lg_star_mass_limit = 2d-3 ! RECOMMENDED 2d-3
   delta_lgRho_cntr_limit = 0.0025d0 ! RECOMMENDED 0.0025d0
   delta_lgRho_cntr_hard_limit = 0.005d0 ! RECOMMENDED 0.005d0
   dt_div_min_dr_div_cs_limit = 1d99
   min_timestep_limit = 1d-20 ! (seconds)
   relax_hard_limits_after_backup = .false.
   relax_hard_limits_after_retry = .false.

   ! mesh controls
   okay_to_remesh = .true.
   max_dq = 1d-3 ! RECOMMENDED 1d-3
   mesh_delta_coeff = 0.8d0 ! RECOMMENDED 0.8d0
   !! this one is turned on in run_star_extras
   !use_split_merge_amr = .true.
   split_merge_amr_log_zoning = .true.
   split_merge_amr_nz_baseline = 6000 ! RECOMMENDED 6000
   split_merge_amr_MaxLong = 1.25d0
   split_merge_amr_MaxShort = 2.5d0
   split_merge_amr_max_iters = 1000
   split_merge_amr_okay_to_split_nz = .false.

    ! Fixing the position of the Lagrangian region of the mesh helps
    ! convergence near the Eddington limit
    max_logT_for_k_below_const_q = 100
    max_q_for_k_below_const_q = 0.99
    min_q_for_k_below_const_q = 0.99
    max_logT_for_k_const_mass = 100
    max_q_for_k_const_mass = 0.98
    min_q_for_k_const_mass = 0.98
    fix_eps_grav_transition_to_grid = .true.

   !! equation controls
   !! these are turned on in run_star_extras with Riemann hydro
   !! this is done by reading inlist_hydro_on
   !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.

   history_interval = 1
   max_num_profile_models = 10000 ! RECOMMENDED 10000
   profile_interval = 100000
   warn_rates_for_high_temp = .false.

   ! Here begins the inferno of additional paramaters to control turning hydro on/off and
   ! relaxing the star to a lower mass
   ! 'Abandon hope all ye who enter here'

   ! in order for Riemann hydro to be turned on, the following two conditions need to be met
   x_ctrl(2) = 9d0 ! logT must be larger than this
   x_ctrl(3) = 0.01d0 ! the integrated gamma1 must be smaller than this
   
   ! or just turn it on if neutrino luminosity is above this limit
   ! independent of the other two
   x_ctrl(11) = 10d0 ! RECOMMENDED 10d0

   ! parameters controlling swtich between hydrostatic and Riemann hydro
   ! after the pulse starts, we estimate the dynamical timescale of the star
   ! (excluding some of the outer layers), and model the dynamical evolution
   ! for at least a fixed number of these timescales before turning off hydro
   x_ctrl(4) = 50d0 ! velocity in km/s at any zone for which we consider the system to be undergoing a pulse
                    ! time will start counting from this point before switching off hydro
                    ! NOTE: we only consider the velocities for q < x_ctrl(7)
   x_ctrl(5) = 0.9d0 ! dynamical timescale will be estimated for the inside of the
                     ! star, up to a mass fraction equal to this value of the CO core mass.
                     ! If there is no CO core yet (unlikely but just to be sure),
                     ! the He core mass is used. The timescale is calculated as 1/sqrt(G*<rho>),
                     ! where <rho>=M/(4 pi R^3/3) up to the specified q
   x_ctrl(6) = 50d0 ! number of dynamical timescales to evolve before turning off hydro
   x_ctrl(1) = 1d0 ! When its time to remove the outer layers, all mass above any
                     ! layer that reaches the escape velocity times this factor
                     ! will be removed
   x_ctrl(7) = 0.99d0 ! after the set number of timescales have passed, before turning off hydro
                    ! two more conditions must be met. This is to allow more time for additional
                    ! shocks to develop on the outermost layers, which can add additional heat there.
                    ! These conditions must be met at all points inside a given mass fraction of
                    ! what would remain of the star, and that fraction is given by this option.
                    ! For example, for x_ctrl(7) = 0.95, if a star with 55 Msun would eject 5 Msun
                    ! leaving a 50 Msun star, then the additional conditions need to be satisfied
                    ! in the inner 47.5 Msun of the star.
                    ! The additional conditions are as follows
   x_ctrl(8) = 20d0 ! absolute velocities must be below this in km/s
   x_ctrl(9) = 0.5d0 ! the ratio of the (absolute) velocity to the sound speed must be below this
   x_ctrl(16) = 11d0 ! do not relax if log10 Lneu/Lsun is bigger than this
   x_ctrl(20) = 10d0 ! do not relax if log10 Lnuc/Lsun is bigger than this

   x_integer_ctrl(1) = 100 ! conditions for relax must be satisfied for this number of consecutive steps
                           ! before the outer layers are removed


   x_ctrl(10) = 9.6d0 ! if temperature is above this, just turn on Riemann hydro and keep it on
   x_ctrl(17) = 8.5d0   ! After first pulse, if logL_neu is above this, turn on Riemann hydro. Helps when the outer
                      ! layers are very bloated to prevent them going crazy when the core evolves too fast.
                      ! RECOMMENDED 8.5d0

   x_ctrl(18) = 0.001d0 ! central mass fraction changes of He, C and O are limited to this (by adjusting timestep).
                       ! RECOMMENDED 0.001d0

   x_ctrl(12) = 2d4 ! if surface velocity exceeds this, switch to fixed_vsurf BC with the current vsurf

   x_logical_ctrl(1) = .true. ! Remove outer layers that exceed 1d4 Rsun, used to simulate post pulse evolution
                              ! for longer without worrying about the outer layers exceeding the range of the EOS

   x_ctrl(23) = 1d3 ! After a pulse begins, limit the timestep to this (in seconds). Ignored after outer layers expand
                    ! beyond 1d4 Rsun when x_logical_ctrl(1) is true
                    ! RECOMMENDED 1d3

   x_ctrl(19) = 8d0 ! turn off wind mass loss  if log10(Lneu/Lsun) exceeds this

   x_ctrl(21) = 0.025d0 ! log(L_nuc) timestep limit is set to this, hard limit to double this value.
                        ! This is turned off during relax and near core collapse
                        ! RECOMMENDED 0.025d0

   x_ctrl(22) = 15d0 ! Ignore the L_nuc limit of x_ctrl(21) if L_nuc exceeds this value, otherwise step to step
                     ! variations can lead to tiny timesteps

   x_logical_ctrl(2) = .false. ! If true, terminate at the onset of the second pulse
                              ! meant for use in the test_suite
                              ! RECOMMENDED .false.

   x_ctrl(24)=1.d0 ! Wind eta
   x_integer_ctrl(2)=1 ! wind choice 1 brott 2 NL
   x_integer_ctrl(3)=1 ! 1 NL 2 tramper
   x_integer_ctrl(4)=0 ! 1 used 1d-3 for Z in winds *NOT USED*

   use_other_neu=.true.
   x_ctrl(26)=3 ! Number of neutrino species
   x_ctrl(25)=0.d0 ! Sigmas from median neutrino rate
   x_ctrl(27)=0.23867d0 ! Weinberg angle

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

      ! hydro debugging
      !hydro_check_everything = .true.
      !hydro_inspectB_flag = .true.
      !hydro_sizequ_flag = .true.
      
      !hydro_get_a_numerical_partial = 1d-4
      !hydro_test_partials_k = 1
      !hydro_numerical_jacobian = .true.
      !hydro_save_numjac_plot_data = .true.
      !hydro_dump_call_number = 195
      !hydro_dump_iter_number = 5
      !hydro_dump_stage_number = -1
      !hydro_epsder_struct = 1d-6
      !hydro_epsder_chem = 1d-6
      !hydro_save_photo = .true. ! Saves a photo when hydro_call_number = hydro_dump_call_number -1

      !fill_arrays_with_NaNs = .true.
      !stop_for_NaNs = .true.
      
      !max_years_for_timestep = 3.67628942044319d-05

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

      !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
   pgstar_interval = 10
   !pause = .true.
   
   pgstar_age_disp = 2.5
   pgstar_model_disp = 2.5
   
   !### scale for axis labels
   pgstar_xaxis_label_scale = 1.3
   pgstar_left_yaxis_label_scale = 1.3
   pgstar_right_yaxis_label_scale = 1.3
   
   Grid2_win_flag = .true.
   !Grid2_win_flag = .false.
   
   Grid2_win_width = 15
   Grid2_win_aspect_ratio = 0.65 ! aspect_ratio = height/width
   
   Grid2_plot_name(4) = 'Mixing'
   
   Grid2_num_cols = 7 ! divide plotting region into this many equal width cols
   Grid2_num_rows = 8 ! divide plotting region into this many equal height rows
   
   Grid2_num_plots = 6 ! <= 10
   
   Grid2_plot_name(1) = 'TRho_Profile'
   Grid2_plot_row(1) = 1 ! number from 1 at top
   Grid2_plot_rowspan(1) = 3 ! plot spans this number of rows
   Grid2_plot_col(1) =  1 ! number from 1 at left
   Grid2_plot_colspan(1) = 2 ! plot spans this number of columns 
   Grid2_plot_pad_left(1) = -0.05 ! fraction of full window width for padding on left
   Grid2_plot_pad_right(1) = 0.01 ! fraction of full window width for padding on right
   Grid2_plot_pad_top(1) = 0.00 ! fraction of full window height for padding at top
   Grid2_plot_pad_bot(1) = 0.05 ! fraction of full window height for padding at bottom
   Grid2_txt_scale_factor(1) = 0.65 ! multiply txt_scale for subplot by this


   Grid2_plot_name(5) = 'Kipp'
   Grid2_plot_row(5) = 4 ! number from 1 at top
   Grid2_plot_rowspan(5) = 3 ! plot spans this number of rows
   Grid2_plot_col(5) =  1 ! number from 1 at left
   Grid2_plot_colspan(5) = 2 ! plot spans this number of columns 
   Grid2_plot_pad_left(5) = -0.05 ! fraction of full window width for padding on left
   Grid2_plot_pad_right(5) = 0.01 ! fraction of full window width for padding on right
   Grid2_plot_pad_top(5) = 0.03 ! fraction of full window height for padding at top
   Grid2_plot_pad_bot(5) = 0.0 ! fraction of full window height for padding at bottom
   Grid2_txt_scale_factor(5) = 0.65 ! multiply txt_scale for subplot by this
   Kipp_title = ''
   Kipp_show_mass_boundaries = .false.

   Grid2_plot_name(6) = 'History_Panels1'
   Grid2_plot_row(6) = 6 ! number from 1 at top
   Grid2_plot_rowspan(6) = 3 ! plot spans this number of rows
   Grid2_plot_col(6) =  6 ! number from 1 at left
   Grid2_plot_colspan(6) = 2 ! plot spans this number of columns 
   !Grid2_plot_pad_left(6) = 0.00 ! fraction of full window width for padding on left
   !Grid2_plot_pad_right(6) = 0.05 ! fraction of full window width for padding on right
   !Grid2_plot_pad_top(6) = 0.03 ! fraction of full window height for padding at top
   !Grid2_plot_pad_bot(6) = 0.0 ! fraction of full window height for padding at bottom
   !Grid2_txt_scale_factor(6) = 0.65 ! multiply txt_scale for subplot by this

   Grid2_plot_pad_left(6) = 0.05 ! fraction of full window width for padding on left
   Grid2_plot_pad_right(6) = 0.03 ! fraction of full window width for padding on right
   Grid2_plot_pad_top(6) = 0.0 ! fraction of full window height for padding at top
   Grid2_plot_pad_bot(6) = 0.0 ! fraction of full window height for padding at bottom
   Grid2_txt_scale_factor(6) = 0.65 ! multiply txt_scale for subplot by this

   History_Panels1_title = ''      
   History_Panels1_num_panels = 2

   History_Panels1_xaxis_name='model_number'
   History_Panels1_max_width = -1 ! only used if > 0.  causes xmin to move with xmax.
   
   History_Panels1_yaxis_name(1) = 'log_center_T' 
   History_Panels1_yaxis_reversed(1) = .false.
   History_Panels1_ymin(1) = -101d0 ! only used if /= -101d0
   History_Panels1_ymax(1) = -101d0 ! only used if /= -101d0        
   History_Panels1_dymin(1) = 1 
   
   History_Panels1_other_yaxis_name(1) = 'gamma_integral' 
   History_Panels1_other_yaxis_reversed(1) = .false.
   History_Panels1_other_ymin(1) = -0.06d0 ! only used if /= -101d0
   History_Panels1_other_ymax(1) = 0.1d0 ! only used if /= -101d0        
   History_Panels1_other_dymin(1) = 0.05d0
   
   History_Panels1_yaxis_name(2) = 'log_Lnuc' 
   History_Panels1_yaxis_reversed(2) = .false.
   History_Panels1_ymin(2) = 0d0 ! only used if /= -101d0
   History_Panels1_ymax(2) = 20d0 ! only used if /= -101d0        
   History_Panels1_dymin(2) = 1 
   
   History_Panels1_other_yaxis_name(2) = 'log_Lneu' 
   History_Panels1_other_yaxis_reversed(2) = .false.
   History_Panels1_other_ymin(2) = 0d0 ! only used if /= -101d0
   History_Panels1_other_ymax(2) = 20d0 ! only used if /= -101d0        
   History_Panels1_other_dymin(2) = 1 
   
   Grid2_plot_name(2) = 'Text_Summary1'
   Grid2_plot_row(2) = 7 ! number from 1 at top
   Grid2_plot_rowspan(2) = 2 ! plot spans this number of rows
   Grid2_plot_col(2) = 1 ! number from 1 at left
   Grid2_plot_colspan(2) = 4 ! plot spans this number of columns 
   Grid2_plot_pad_left(2) = -0.08 ! fraction of full window width for padding on left
   Grid2_plot_pad_right(2) = -0.10 ! fraction of full window width for padding on right
   Grid2_plot_pad_top(2) = 0.08 ! fraction of full window height for padding at top
   Grid2_plot_pad_bot(2) = -0.04 ! fraction of full window height for padding at bottom
   Grid2_txt_scale_factor(2) = 0.19 ! multiply txt_scale for subplot by this
         
   Text_Summary1_name(2,1) = 'star_age'
   Text_Summary1_name(3,1) = 'time_step_sec'
   Text_Summary1_name(8,1) = 'yr_since_coll'
   Text_Summary1_name(3,2) = 'co_core_mass'
   Text_Summary1_name(4,2) = 'M_below_vesc'
   Text_Summary1_name(5,2) = 'gamma_integral'
   Text_Summary1_name(6,2) = 'log_R_vesc_098'
   Text_Summary1_name(7,2) = 'log_R_vesc_095'
   Text_Summary1_name(8,2) = 'log_R_vesc_090'
   
   Grid2_plot_name(3) = 'Profile_Panels3'
   Profile_Panels3_title = 'Abundance-Power-Mixing'
   Profile_Panels3_num_panels = 3
   Profile_Panels3_yaxis_name(1) = 'Abundance'
   Profile_Panels3_yaxis_name(2) = 'Power'
   Profile_Panels3_yaxis_name(3) = 'Mixing'
   
   Profile_Panels3_xaxis_name = 'mass'
   Profile_Panels3_xaxis_reversed = .false.
   Profile_Panels3_xmin = -101d0 ! only used if /= -101d0
   Profile_Panels3_xmax = -101d0 ! 10 ! only used if /= -101d0
   
   Grid2_plot_row(3) = 1 ! number from 1 at top
   Grid2_plot_rowspan(3) = 6 ! plot spans this number of rows
   Grid2_plot_col(3) = 3 ! plot spans this number of columns 
   Grid2_plot_colspan(3) = 3 ! plot spans this number of columns 
   
   Grid2_plot_pad_left(3) = 0.09 ! fraction of full window width for padding on left
   Grid2_plot_pad_right(3) = 0.07 ! fraction of full window width for padding on right
   Grid2_plot_pad_top(3) = 0.0 ! fraction of full window height for padding at top
   Grid2_plot_pad_bot(3) = 0.0 ! fraction of full window height for padding at bottom
   Grid2_txt_scale_factor(3) = 0.65 ! multiply txt_scale for subplot by this

   Grid2_plot_name(4) = 'Profile_Panels1'
   Grid2_plot_row(4) = 1 ! number from 1 at top
   Grid2_plot_rowspan(4) = 5 ! plot spans this number of rows
   Grid2_plot_col(4) =  6 ! number from 1 at left
   Grid2_plot_colspan(4) = 2 ! plot spans this number of columns 
   Grid2_plot_pad_left(4) = 0.05 ! fraction of full window width for padding on left
   Grid2_plot_pad_right(4) = 0.03 ! fraction of full window width for padding on right
   Grid2_plot_pad_top(4) = 0.0 ! fraction of full window height for padding at top
   Grid2_plot_pad_bot(4) = 0.07 ! fraction of full window height for padding at bottom
   Grid2_txt_scale_factor(4) = 0.65 ! multiply txt_scale for subplot by this

      
   Abundance_line_txt_scale_factor = 1.1 ! relative to other text
   Abundance_legend_txt_scale_factor = 1.1 ! relative to other text
   Abundance_legend_max_cnt = 0
   Abundance_log_mass_frac_min = -3.5 ! only used if < 0
   
   Profile_Panels1_title = ''
   
   Profile_Panels1_xaxis_name = 'mass'
   
   Profile_Panels1_xmin = -101d0 ! only used if /= -101d0
   Profile_Panels1_xmax = -101d0 ! 100.1 ! -101d0 ! only used if /= -101d0
   Profile_Panels1_num_panels = 3
   
   Profile_Panels1_yaxis_name(1) = 'logRho'
   Profile_Panels1_other_yaxis_name(1) = 'gamma1'
   Profile_Panels1_other_ymax(1) = 1.41
   Profile_Panels1_other_ymin(1) = 1.28
   
   Profile_Panels1_yaxis_name(2) = 'vel_km_per_s'
   Profile_Panels1_dymin(2) = 50
   
   Profile_Panels1_other_yaxis_name(2) = 'v_div_csound' ! 'vel_km_per_s'
   Profile_Panels1_other_ymax(2) = 4.5 ! -101d0
   Profile_Panels1_other_ymin(2) = -101d0
   Profile_Panels1_other_ymargin(2) = 0.1 ! 10
   Profile_Panels1_other_dymin(2) = 0.5
   
   
   Profile_Panels1_yaxis_name(3) = 'logR_cm'
   Profile_Panels1_ymax(3) = -101d0
   Profile_Panels1_ymin(3) = 7
   Profile_Panels1_dymin(3) = 0.1
   
   Profile_Panels1_other_yaxis_name(3) = 'v_div_v_escape' ! 'mass'
   Profile_Panels1_other_ymax(3) = 5d0
   Profile_Panels1_other_ymin(3) = -2d0
   Profile_Panels1_other_dymin(3) = 0.15
   
   show_TRho_Profile_kap_regions = .false.
   show_TRho_Profile_gamma1_4_3rd = .true.
            TRho_Profile_xmin = -30
            !TRho_Profile_xmax = 10
            TRho_Profile_ymin = 1
            !TRho_Profile_ymax = 10
   
   Grid2_file_flag = .true.
   Grid2_file_dir="Grid1"
   Grid2_file_prefix = 'grid_'
   Grid2_file_interval=1000
   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

   !show_TRho_Profile_eos_regions = .true.

/ ! end of pgstar namelist
