
&star_job

      show_log_description_at_start = .false.

      create_RSP_model = .true.

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

      kappa_file_prefix = 'a09'
      kappa_lowT_prefix = 'lowT_fa05_a09p'
      kappa_CO_prefix = 'a09_co'
      initial_zfracs = 6

      color_num_files=3

      color_file_names(1) = 'lcb98cor.dat'
      color_num_colors(1) = 11

      color_file_names(2)='blackbody_johnson.dat'
      color_num_colors(2)=5

      color_file_names(3)='custom_color_file.dat' !Kepler, TESS, Gaia passbands
      color_num_colors(3)=13

      set_initial_age = .true.
      initial_age = 0

      set_initial_model_number = .true.
      initial_model_number = 0

      profile_starting_model = .true.
      
      set_initial_cumulative_energy_error = .true.
      new_cumulative_energy_error = 0d0

/ ! end of star_job namelist


&controls  

      !### radial stellar pulsations (RSP)
         ! inspired by Radec Smolec's program
         ! if you use RSP, please cite the following in addition to mesa:
         ! R.Smolec and P.Moskalik, ACTA ASTRONOMICA, 58(2008), pp. 193-232.
         
         ! must set mass, Teff, L, X, and Z.
         RSP_mass = -1 ! (Msun)
         RSP_Teff = -1 ! (K)
         RSP_L = -1 ! (Lsun)
         RSP_X = -1 ! hydrogen mass fraction
         RSP_Z = -1 ! metals mass fraction

         ! Parameters of the convection model.
         ! Note that RSP_alfap, RSP_alfas, RSP_alfac, RSP_alfad and RSP_gammar
         ! are expressed in the units of standard values.
         ! Standard values are the ones for which static version of the Kuhfuss
         ! model reduces to standard MLT.
         ! See Table 1 in Smolec & Moskalik 2008
         ! http://adsabs.harvard.edu/abs/2008AcA....58..193S
         ! for standard values and the description of the parameters.
         
         RSP_alfa =      1.5d0         ! mixing length; alfa = 0 gives a purely radiative model.
         RSP_alfac =     1.0d0         ! convective flux; Lc ~ RSP_alfac
         RSP_alfas =     1.0d0         ! turbulent source; Lc ~ 1/ALFAS; PII ~ RSP_alfas
         RSP_alfad =     1.0d0         ! turbulent dissipation; damp ~ RSP_alfad
         RSP_alfap =     1.0d0         ! turbulent pressure; Pt ~ alfap     
         RSP_alfat =     0.01d0         ! turbulent flux; Lt ~ RSP_alfat; overshooting.
         RSP_alfam =     0.40d0        ! eddy viscosity; Chi & Eq ~ RSP_alfam
         RSP_gammar =    0.0d0         ! radiative losses; dampR ~ RSP_gammar          
                  
         ! time weighting for end-of-step vs start-of-step values in equations.
         ! 1 corresponds to fully implicit scheme - stable, but can have large numerical damping.
         ! 0.5 corresponds to trapezoidal rule integration - gives least numerical damping.
         ! do not use values less than 0.5.  strongly recommend 0.5 for theta and thetat.
         ! don't mess with any of these unless you know what you are doing or like to watch the code crash.
         
         RSP_theta =      0.5d0        ! Pgas and Prad
         RSP_thetat =     0.5d0        ! Pturb
         RSP_thetae =     0.5d0        ! erad in terms using f_Edd
         RSP_thetaq =     1.0d0        ! avQ
         RSP_thetau =     1.0d0        ! Eq and Uq
         RSP_wtr =        0.6667d0     ! Lr
         RSP_wtc =        0.6667d0     ! Lc
         RSP_wtt =        0.6667d0     ! Lt
         RSP_gam =        1.0d0        ! Et src_snk
         
         ! controls for building the initial model
         RSP_nz = 200 ! total number of zones in initial model
         RSP_nz_outer = 60 ! number of zones in outer region of initial model
         RSP_T_anchor = 11d3 ! approx temperature at base of outer region
         RSP_T_inner = 2d6 ! T at inner boundary of initial model
         
         RSP_max_outer_dm_tries = 100 ! give up if fail to find outer dm in this many attempts
         RSP_max_inner_scale_tries = 100 ! give up if fail to find inner dm scale factor in this many attempts
         RSP_T_anchor_tolerance = 1d-8 
            ! allowed relative difference between T at base of outer region and T_anchor
            ! if fail trying to create initial model, try increasing this to 1d-6 or more
         RSP_T_inner_tolerance = 1d-8 
            ! allowed relative difference between T at inner boundary and T_inner
            ! if fail trying to create initial model, try increasing this to 1d-6 or more
         
         RSP_relax_initial_model = .true.
         RSP_relax_alfap_before_alfat = .true. ! else reverse the order         
         RSP_relax_max_tries = 1000   
         RSP_relax_dm_tolerance = 1d-6
         
         ! Initial kick makes use of the scaled linear velocity eigenvector
         ! of a given mode or of the linear combination of the eigenvectors
         ! for the fundamental mode and first two radial overtones. 
         ! The surface velocity is set to RSP_kick_vsurf_km_per_sec and
         ! the mode content is set by RSP_fraction_1st_overtone and RSP_fraction_2nd_overtone
         RSP_kick_vsurf_km_per_sec = 0.1d0 ! can be negative
         RSP_fraction_1st_overtone = 0d0
         RSP_fraction_2nd_overtone = 0d0
         ! fraction from fundamental = 1d0 - (1st + 2nd)
         ! Note: This is important for models in which two or more modes are linearly unstable.
         ! Appropriate setting may help to arrive at the desired mode, since the final pulsation
         ! state may depend on initial conditions set by the three parameters above.
         ! Integration of the same model with different initial kicks is a way to study
         ! the nonlinear mode selection -
         ! for an example see Fig. 1 in Smolec & Moskalik 2010
         ! http://adsabs.harvard.edu/abs/2010A%26A...524A..40S .
         
         ! random initial velocity profile.  added to any kick from eigenvector.
         RSP_Avel = 0d0 ! kms. linear in mesh points from 0 at inner boundary to this at surface
         RSP_Arnd = 0d0 ! kms. random fluctuation at each mesh point.
         
         ! period controls
         RSP_target_steps_per_cycle = 600
         RSP_save_profile_at_end_of_each_period = .true.
         RSP_min_max_R_for_periods = -1 ! skip if < 0; Rsun units.
         RSP_min_PERIOD_div_PERIODLIN = 0.5d0
         RSP_mode_for_setting_PERIODLIN = 1 ! 0, 1, or 2 for fundamental, 1st, and 2nd overtones.
         
         ! when to stop
         RSP_max_num_periods = 4000 ! ignore if < 0
	 !max_model_number = 1
         RSP_GREKM_avg_abs_frac_new = 0.1d0 ! fraction of new for updating avg at each cycle.
         RSP_GREKM_avg_abs_limit = -1 ! 0.01d0 ! ignore if < 0; else stop when GREKM_avg_abs < this limit.
         
         ! timestep limiting
         RSP_initial_dt_factor = 1d-2 ! set initial timestep to this times linear period/target_steps_per_cycle
            ! start with smaller timestep to give time for initial model to adjust
         RSP_Uq_threshold_for_dt_limit = 100d0
         RSP_v_div_cs_threshold_for_dt_limit = 0.8d0
         RSP_max_dt_times_min_dr_div_cs = 5d0 ! limit dt by this
            ! i.e., make dt <= this times min sound crossing time dr/cs
            !      for cells with abs(Uq) > threshold
         RSP_max_dt = -1 ! seconds
         RSP_report_limit_dt = .false.

         ! artificial viscosity controls
         ! for the equations see: Appendix C in Stellingwerf 1975
         ! http://adsabs.harvard.edu/abs/1975ApJ...195..441S.
         ! In principle, for not too-non-adiabatic convective models artificial viscosity is not
         ! needed or should be very small. Hence a large cut-off parameter below (in purely
         ! radiative models the default value for cut-off was 0.01)
         RSP_cq = 4.0d0 ! viscosity parameter (viscosity pressure proportional to cq)
         RSP_zsh = 0.1d0 ! "turn-on" compression in units of sound speed.
            ! zsh > 0 delays onset of artificial viscosity
            ! can eliminate most/all interior dissipation while still providing for extreme cases.
            ! using this parameter the dependence of limiting amplitude on cq is very weak.
         RSP_Qvisc_quadratic = 0d0 ! 1d-3 ! for Tscharnuter & Winkler form of artificial viscosity
            ! as described in section 4.2 of mesa3, 2015.
            ! RSP_Qvisc_quadratic is analogous to shock_spread_quadratic
            ! if switch to this form, set RSP_cq = 0 to shut off the Neumann & Richtmyer form.
            ! note that this form also uses RSP_zsh to delay onset of artficial viscosity
         
         ! surface pressure.  provides outer boundary condition for momentum equation.
         RSP_use_Prad_for_Psurf = .false.
         RSP_use_atm_grey_with_kap_for_Psurf = .false.
         RSP_tau_surf_for_atm_grey_with_kap = 3d-3 ! for atm_grey_with_kap
         RSP_fixed_Psurf = .true.
         RSP_Psurf = 0d0 ! ignore if < 0.  else use as surface pressure.
         set_RSP_Psurf_to_multiple_of_initial_P1 = -1 ! ignore if < 0. 
            ! set RSP_Psurf to this times initial surface cell pressure
         
         ! solver controls
         RSP_tol_max_corr = 1d-8
         RSP_tol_max_resid = 1d-6
         RSP_max_iters_per_step = 100
         RSP_max_retries_per_step = 8
         max_number_retries = -1 ! default for rsp is to ignore total num retries
         RSP_report_undercorrections = .false.
         RSP_nz_div_IBOTOM = 30d0 ! set IBOTOM = RSP_nz/RSP_nz_div_IBOTOM
         RSP_erad_eqn_dt_term_limit = 1d7 ! to avoid numeric issues when heat exchange time scale << dt
         RSP_Fr_eqn_always_use_d_erad_dm = .false. ! even for cases where f_Edd /= 1/3
         RSP_min_tau_for_turbulent_flux = 2d2 ! optionally turn off turbulent flux near surface 
         
         ! output data for work integrals during a particular period
         RSP_work_period = -1
         RSP_work_filename = 'work.data'
         
         ! output data for 3d map.  format same as for gnuplot pm3d 
         RSP_write_map = .false.
         RSP_map_columns_filename = 'map_columns.list' ! subset of profile columns list
            ! items listed in your map columns must also appear in your profile columns
         RSP_map_filename = 'map.data'
         RSP_map_first_period = -1
         RSP_map_last_period = -1
         RSP_map_zone_interval = 2
         RSP_map_history_filename = 'map_history.data'
         
         ! rsp hooks
         use_other_RSP_starting_velocities = .false. ! e.g., GYRE
         use_other_RSP_eddington_factor = .false.
         use_other_RSP_build_model = .false.

         ! for special tests
         ! can set ALFA = 0 for pure radiative with no turbulence or convection.
         ! can set gamma_law_hydro > 0 and it will be used by rsp for eos.
         ! can set zero_gravity = .true.
         ! can set opacity to be constant times density.
         RSP_kap_density_factor = -1 ! ignore if < 0
            ! else set opacity to this times density
         
         ! rsp misc
         RSP_efl0 =      1.0d2        
         RSP_nmodes = 3 ! number of modes LINA will calculate for initial model
         RSP_trace_RSP_build_model = .false.
         RSP_use_nlin_constants = .false.

         ! -------- following are not ready for general use --------
         
         ! rad hydro - not ready
         RSP_use_diffusion_limit = .true.
         RSP_use_T_form_of_erad_eqn = .true.
         RSP_use_T_form_of_Fr_eqn = .true.
         RSP_accel_eqn_use_Prad_instead_of_Fr_term = .true.
         
         ! f_Edd - not ready
         RSP_nz_div_number_of_core_rays_for_f_Edd = 5d0 ! adjust number of core rays to nz
         RSP_max_tau_for_variable_f_Edd = 1d3 ! for tau > this, f_Edd = 1/3
         
         ! parameters for radiative energy equation - not ready
         RSP_kapE_factor = 1d0 ! absorption mean opacity is approximated by this factor times Rosseland
         RSP_kapP_factor = 1d0 ! Planck mean opacity is approximated by this factor times Rosseland
                  
         ! mesh adjustment - not ready
         RSP_remesh_interval = 1 ! remesh when mod(model_number,this) == 0
         ! note: the following default weights are not appropriate for all situations.
         ! check the test cases for examples.
         RSP_mesh_fcn_weight_lnd = 0d0
         RSP_mesh_fcn_weight_lnT = 0d0
         RSP_mesh_fcn_weight_lnP = 0d0
         RSP_mesh_fcn_weight_lnE = 1d0
         RSP_mesh_fcn_weight_lnm = 4d0
         RSP_mesh_fcn_weight_lnkap = 1d0
         RSP_mesh_fcn_weight_avg_charge_H = 0d0
         RSP_mesh_fcn_weight_avg_charge_He = 0d0
         
         RSP_max_dq = 1d-1
         RSP_min_dq = 1d-8
         RSP_max_dq_jump = 1.5d0 ! max ratio of adjacent dq's
         RSP_max_periods_for_remeshing = 0 ! 0 means no remeshing; ignore this param if < 0
         ! note that rsp remeshing is off if okay_to_remesh is false
         
         RSP_remesh_test_partials = .false.
         trace_RSP_remesh = .false.

! general 
      set_rho_to_dm_div_dV = .false.

      use_eosDT2 = .true.
      use_eosELM = .true.

      ! check for retries and backups as part of test_suite
         !max_number_backups = 0
         !max_number_retries = 0
         !max_model_number = 4328
      
      mix_factor = 0d0
      dxdt_nuc_factor = 0d0      
      non_nuc_neu_factor = 0d0
      eps_nuc_factor = 0d0
      
! eos controls
      use_eosPTEH_for_low_density = .true.
      use_eosPTEH_for_high_Z = .true.
      use_eosELM = .true.
      use_eosDT2 = .true.
      
! opacity controls
      cubic_interpolation_in_X = .false.
      cubic_interpolation_in_Z = .false.
      include_electron_conduction = .true.
      use_Zbase_for_Type1_blend = .true.
      use_Type2_opacities = .false.
      kap_Type2_full_off_X  =  0.71d0
      kap_Type2_full_on_X  =  0.70d0
      kap_Type2_full_off_dZ  =  1d-3 
      kap_Type2_full_on_dZ  =  1d-2

! output controls

   ! for cases in which you have a run that has reached steady pulses
   ! and you want to look at plots of just a few periods to see the details,
   ! i suggest the following method. interrupt the run soon after
   ! it makes a photo.  remove or delete LOGS/history.data to force
   ! creation of a new one when restart.  edit the inlist to set
   ! history_interval to 1 to get maximum time resolution.  
   ! restart the run and let it go for as many periods as desired.   

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

      terminal_show_age_in_years = .false.
      terminal_show_age_in_days = .true.

      photo_interval = 1000     
      profile_interval = 1000          
      history_interval = 1            
      terminal_interval = 4000     
      max_num_profile_models = 100     
      log_directory='LOGS'
      photo_directory='photos'

      !photo_interval = 10
      !profile_interval = 10
      !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
      !stop_for_bad_nums = .true.
      !trace_evolve = .true.

      ! hydro debugging
      !hydro_check_everything = .true.
      !hydro_inspectB_flag = .true.
      !hydro_sizequ_flag = .true.
      
      ! for making logs, uncomment hydro_dump_call_number plus the following
      ! to make residual logs, uncomment hydro_sizequ_flag
      ! to make correction logs, uncomment hydro_inspectB_flag
      ! to make jacobian logs, uncomment hydro_numerical_jacobian, hydro_save_numjac_plot_data
      ! to do dfridr test, uncomment hydro_get_a_numerical_partial, hydro_test_partials_k,
      !     hydro_numerical_jacobian, hydro_save_numjac_plot_data, hydro_dump_iter_number
         
      !hydro_get_a_numerical_partial = 1d-7
      !hydro_test_partials_k = 26
      !hydro_numerical_jacobian = .true.
      !hydro_save_numjac_plot_data = .true.
      !hydro_dump_call_number = 449
      !hydro_dump_iter_number = 6
      
      !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.

      !report_why_dt_limits = .true.
      !report_all_dt_limits = .true.
      !report_hydro_dt_info = .true.
      !report_dX_nuc_drop_dt_limits = .true.
      !report_bad_negative_xa = .true.
      
      !show_mesh_changes = .true.
      !mesh_dump_call_number = 95
      

      !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

         pgstar_xaxis_label_scale = 1.0
         pgstar_left_yaxis_label_scale = 1.0
         pgstar_right_yaxis_label_scale = 1.0

         pgstar_xaxis_label_disp = 2.5
         pgstar_left_yaxis_label_disp = 4.2 ! 3.4
         pgstar_right_yaxis_label_disp = 4.4


         Grid2_win_width = 10
         Grid2_win_aspect_ratio = 0.75 ! aspect_ratio = height/width
         Grid2_num_cols = 3 ! 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) = 'Profile_Panels1'
         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.02 ! fraction of full window width for padding on left
         Grid2_plot_pad_right(1) = 0.09 ! fraction of full window width for padding on right
         Grid2_plot_pad_top(1) = -0.02 ! fraction of full window height for padding at top
         Grid2_plot_pad_bot(1) = -0.10 ! fraction of full window height for padding at bottom
         Grid2_txt_scale_factor(1) = 0.7 ! multiply txt_scale for subplot by this
         
         Grid2_plot_name(2) = 'History_Panels1'
         Grid2_plot_row(2) = 4 ! number from 1 at top
         Grid2_plot_rowspan(2) = 3 ! plot spans this number of rows
         Grid2_plot_col(2) =  1 ! number from 1 at left
         Grid2_plot_colspan(2) = 2 ! plot spans this number of columns 
         Grid2_plot_pad_left(2) = -0.02 ! fraction of full window width for padding on left
         Grid2_plot_pad_right(2) = 0.09 ! fraction of full window width for padding on right
         Grid2_plot_pad_top(2) = 0.18 ! fraction of full window height for padding at top
         Grid2_plot_pad_bot(2) = -0.02 ! fraction of full window height for padding at bottom
         Grid2_txt_scale_factor(2) = 0.7 ! multiply txt_scale for subplot by this
         
         Grid2_plot_name(3) = 'Text_Summary1'
         Grid2_plot_row(3) = 7 ! number from 1 at top
         Grid2_plot_rowspan(3) = 2 ! plot spans this number of rows
         Grid2_plot_col(3) =  1 ! number from 1 at left
         Grid2_plot_colspan(3) = 3 ! plot spans this number of columns 
         Grid2_plot_pad_left(3) = -0.08 ! fraction of full window width for padding on left
         Grid2_plot_pad_right(3) = 0.0 ! fraction of full window width for padding on right
         Grid2_plot_pad_top(3) = 0.12 ! fraction of full window height for padding at top
         Grid2_plot_pad_bot(3) = -0.04 ! fraction of full window height for padding at bottom
         Grid2_txt_scale_factor(3) = 0.2 ! multiply txt_scale for subplot by this
         
         Grid2_plot_name(4) = 'logL_R'
         Grid2_plot_row(4) = 5 ! number from 1 at top
         Grid2_plot_rowspan(4) = 2 ! plot spans this number of rows
         Grid2_plot_col(4) =  3 ! number from 1 at left
         Grid2_plot_colspan(4) = 1 ! plot spans this number of columns 
         Grid2_plot_pad_left(4) = 0.07 ! fraction of full window width for padding on left
         Grid2_plot_pad_right(4) = 0.0 ! fraction of full window width for padding on right
         Grid2_plot_pad_top(4) = 0.06 ! fraction of full window height for padding at top
         Grid2_plot_pad_bot(4) = -0.02 ! fraction of full window height for padding at bottom
         Grid2_txt_scale_factor(4) = 0.7 ! multiply txt_scale for subplot by this
         logL_R_title = ''      
         
         Grid2_plot_name(5) = 'logL_v'
         Grid2_plot_row(5) = 3 ! number from 1 at top
         Grid2_plot_rowspan(5) = 2 ! plot spans this number of rows
         Grid2_plot_col(5) = 3 ! number from 1 at left
         Grid2_plot_colspan(5) = 1 ! plot spans this number of columns 
         Grid2_plot_pad_left(5) = 0.07 ! fraction of full window width for padding on left
         Grid2_plot_pad_right(5) = 0.0 ! fraction of full window width for padding on right
         Grid2_plot_pad_top(5) = 0.02 ! fraction of full window height for padding at top
         Grid2_plot_pad_bot(5) = 0.02 ! fraction of full window height for padding at bottom
         Grid2_txt_scale_factor(5) = 0.7 ! multiply txt_scale for subplot by this
         logL_v_title = ''      
         
         Grid2_plot_name(6) = 'logL_Teff'
         Grid2_plot_row(6) = 1 ! number from 1 at top
         Grid2_plot_rowspan(6) = 2 ! plot spans this number of rows
         Grid2_plot_col(6) = 3 ! number from 1 at left
         Grid2_plot_colspan(6) = 1 ! plot spans this number of columns 
         Grid2_plot_pad_left(6) = 0.07 ! fraction of full window width for padding on left
         Grid2_plot_pad_right(6) = 0.0 ! fraction of full window width for padding on right
         Grid2_plot_pad_top(6) = -0.02 ! fraction of full window height for padding at top
         Grid2_plot_pad_bot(6) = 0.06 ! fraction of full window height for padding at bottom
         Grid2_txt_scale_factor(6) = 0.7 ! multiply txt_scale for subplot by this
         logL_Teff_title = ''      

   !   Grid2_file_flag = .true.
         file_digits = 7
         Grid2_file_dir = 'png'
         Grid2_file_prefix = 'grid1'
         Grid2_file_interval = 10 ! output when mod(model_number,Grid2_file_interval)==0


         History_Panels1_title = ''      
         History_Panels1_num_panels = 2

         !History_Panels1_xaxis_name='model_number'
         History_Panels1_xaxis_name='star_age_day'
         History_Panels1_max_width = 50 ! only used if > 0.  causes xmin to move with xmax.
         
         History_Panels1_yaxis_name(1) = 'effective_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) = 'log_L' 
         History_Panels1_other_yaxis_reversed(1) = .false.
         History_Panels1_other_ymin(1) = -101d0 ! only used if /= -101d0
         History_Panels1_other_ymax(1) = -101d0 ! only used if /= -101d0        
         History_Panels1_other_dymin(1) = -1 
         
         History_Panels1_yaxis_name(2) = 'radius'  !  'star_mass' 
         History_Panels1_yaxis_reversed(2) = .false.
         History_Panels1_ymin(2) = -101d0 ! only used if /= -101d0
         History_Panels1_ymax(2) = -101d0 ! only used if /= -101d0        
         History_Panels1_dymin(2) = -1 
         
         History_Panels1_other_yaxis_name(2) = 'v_surf'
         History_Panels1_other_yaxis_reversed(2) = .false.
         History_Panels1_other_ymin(2) = -101d0 ! only used if /= -101d0
         History_Panels1_other_ymax(2) = -101d0 ! only used if /= -101d0        
         History_Panels1_other_dymin(2) = -1 
         
         History_Panels1_yaxis_name(3) = 'r_center'  !  'star_mass' 
         History_Panels1_yaxis_reversed(3) = .false.
         History_Panels1_ymin(3) = -101d0 ! only used if /= -101d0
         History_Panels1_ymax(3) = -101d0 ! only used if /= -101d0        
         History_Panels1_dymin(3) = -1 
         
         History_Panels1_other_yaxis_name(3) = 'v_center'
         History_Panels1_other_yaxis_reversed(3) = .false.
         History_Panels1_other_ymin(3) = -101d0 ! only used if /= -101d0
         History_Panels1_other_ymax(3) = -101d0 ! only used if /= -101d0        
         History_Panels1_other_dymin(3) = -1 


      Profile_Panels1_win_flag = .false.
         Profile_Panels1_title = ''      
         
         Profile_Panels1_xaxis_name = 'logtau'
         Profile_Panels1_xaxis_reversed = .true.
         Profile_Panels1_xmin = -101D0
         Profile_Panels1_xmax = 6.0 ! -101D0
         Profile_Panels1_show_mix_regions_on_xaxis = .true.
         
         Profile_Panels1_num_panels = 4 
                 
         Profile_Panels1_yaxis_name(1) = 'v_div_cs'         
         Profile_Panels1_ymax(1) = -101D0 ! 1
         Profile_Panels1_ymin(1) = -101D0 ! -1
         Profile_Panels1_other_yaxis_name(1) = 'vel_km_per_s'
         Profile_Panels1_other_ymax(1) = -101D0 !  2e6 
         Profile_Panels1_other_ymin(1) = -101D0 ! -2e6
                  
         Profile_Panels1_yaxis_name(2) = 'entropy'         
         Profile_Panels1_other_yaxis_name(2) = 'opacity' 
         
         Profile_Panels1_yaxis_name(3) = 'logRho'         
         Profile_Panels1_ymax(3) = -101D0 ! 5.7    
         Profile_Panels1_ymin(3) = -101D0 ! 4.2     
         Profile_Panels1_other_yaxis_name(3) = 'logT' 
         Profile_Panels1_other_ymax(3) = -101D0 ! 6.45 
         Profile_Panels1_other_ymin(3) = -101D0 ! 4
                  
         Profile_Panels1_yaxis_name(4) = 'logL'  ! 'zone' !        
         Profile_Panels1_ymax(4) = -101D0 ! 3.7  
         Profile_Panels1_ymin(4) = -101D0 ! 3.1 
         Profile_Panels1_other_yaxis_name(4) = 'lum_div_Ledd' 
         Profile_Panels1_other_ymax(4) = -101D0 ! 
         Profile_Panels1_other_ymin(4) = -101D0 ! 
         

         Profile_Panels1_file_flag = .false.
         Profile_Panels1_file_dir = 'png'
         Profile_Panels1_file_prefix = 'profile_panels1_'
         Profile_Panels1_file_interval = 1 ! output when mod(model_number,Profile_Panels1_file_interval)==0

         Text_Summary1_num_rows = 4 ! <= 20
         Text_Summary1_num_cols = 4 ! <= 20
         Text_Summary1_name(:,:) = ''
         
         Text_Summary1_name(1,1) = 'model_number'
         Text_Summary1_name(2,1) = 'star_age_day'
         Text_Summary1_name(3,1) = 'time_step_sec'
         Text_Summary1_name(4,1) = 'rsp_num_periods'
         Text_Summary1_name(5,1) = ''
         
         Text_Summary1_name(1,2) = 'v_surf_km_s'
         Text_Summary1_name(2,2) = 'radius'
         Text_Summary1_name(3,2) = 'log_R'
         Text_Summary1_name(4,2) = 'rsp_period_in_days'
         Text_Summary1_name(5,2) = ''
         
         Text_Summary1_name(1,3) = 'effective_T'
         Text_Summary1_name(2,3) = 'log_Teff'
         Text_Summary1_name(3,3) = 'luminosity'
         Text_Summary1_name(4,3) = 'log_L'
         Text_Summary1_name(5,3) = ''
         
         Text_Summary1_name(1,4) = 'star_mass'
         Text_Summary1_name(2,4) = 'num_zones'
         Text_Summary1_name(3,4) = 'num_retries'
         Text_Summary1_name(4,4) = 'num_backups'
         Text_Summary1_name(5,4) = ''

/ ! end of pgstar namelist
