! inlist_core

&star_job



      show_log_description_at_start = .false.

      load_saved_model = .true.
      saved_model_name = <<loadfile>>

      save_model_when_terminate = .true.
      save_model_filename = <<smwtfname>>
      
      !set_initial_age = .true.
      !initial_age = 0.0 ! in years
      !set_initial_dt = .true.
      !years_for_initial_dt = 1.d-5 ! in years

      ! controls for adding core of a certain density (alternatively, could fix radius)
      relax_core = .true.
      new_core_mass = <<new_core_mass>> ! 10 Mearth in Msun units
      core_avg_rho = <<core_avg_rho>> ! g/cm^3
      core_avg_eps = 0 ! ergs/g/sec
      dlg_core_mass_per_step = 0.01     ! decrease if convergence problems
      relax_core_years_for_dt = 0.2     ! increase(?) if convergence problems

      relax_R_center = .true.
      !relax_initial_R_center = .true.
      new_R_center = <<core_radius>>  ! in cm
      dlgR_per_step = 3d-3
      relax_R_center_dt = 3.1558149984d1 ! secyer*1d-6

      eos_file_prefix = 'mesa'
      kappa_file_prefix = 'gs98'
      kappa_lowT_prefix = 'lowT_Freedman11'

      change_lnPgas_flag = .true.
      new_lnPgas_flag = .true.
      
      !pgstar_flag = .true.

/ !end of star_job

&controls

        ! surface heating
        !column_depth_for_irradiation = 100.0
        !irradiation_flux = 1.d9     ! day side flux!!! Bill puts in 1/4 by hand

        varcontrol_target = 1d-4
 	mesh_delta_coeff = 0.5

	! when to stop
	max_age = 2.d3    
        max_model_number = 2000         ! don't let it go on and on if run stuck

	! knobs
	!use_lnE_for_eps_grav = .false.
	T_mix_limit = 0

	! output to files and terminal
	!photostep = 50
	profile_interval = 100
	max_num_profile_models=1000
	history_interval = 2
	!terminal_cnt = 10
	write_header_frequency = 10
/


&pgstar




/ ! end of pgstar namelist
