

&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>>
      
      write_profile_when_terminate = .true.
      !filename_for_profile_when_terminate = 'final_profile.data'
      
      change_initial_net = .true.      
      new_net_name = 'basic.net'

      set_initial_model_number = .true.
      initial_model_number = 0

      set_initial_age = .true.
      initial_age = <<initial_age>>  ! in years or 6e7

      set_initial_dt = .true.
      years_for_initial_dt =  <<initial_age>> ! in years

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

      set_tau_factor = .true. ! change tau_factor without reconverging.
      set_to_this_tau_factor = 100
      
      change_lnPgas_flag = .true.
      new_lnPgas_flag = .true.

      !pgstar_flag = .true.


/ ! end of star_job namelist



&controls

       
	! surface heating
	!column_depth_for_irradiation = <<irrad_col>> ! 3.d2 cm^2/g
        !irradiation_flux = <<flux_dayside>>   !555501654.562 ! 1.d9 erg/cm^2/s ! day side flux!!! Bill puts in 1/4 by hand
        min_timestep_limit = 1d5

 	mesh_delta_coeff = 0.5
	!varcontrol_target = 1d-3
	!dt_years_for_steps_before_max_age = 1d7
      
	use_Type2_opacities = .false.

	!initial_mass = 0.001
	!initial_z = 0.02d0

	mixing_length_alpha = 1.89d0 ! based on solar model
	MLT_option = 'Henyey'
      

	! when to stop
	max_age = <<maxage>> 
	!max_model_number = 3000 ! don't let it go on and on if run stuck
	Teff_lower_limit = 10
        logQ_limit=20000000000
	!max_years_for_timestep = 1d9
      


	T_mix_limit = 0

	! output to files and terminal
	!photostep = 100
	profile_interval = 450
	max_num_profile_models=2
	history_interval = 1
	!terminal_cnt = 10
	write_header_frequency = 10


   
        ! controls for atmospheric escape
        use_other_adjust_mdot= <<knob>>
        x_ctrl(1)= 0.3 ! mass loss efficiency, higher for low mass planets
        x_ctrl(2)= <<orbital_distance>> ! orbital seperation in AU

        ! controls for deep heating
        use_other_energy = .true.
        x_ctrl(3) = 1d10 ! core heat capacity  eps=erg/g/K
        x_ctrl(4) = 1.0 ! fraction of core mass in "chondrite-like" rocky material
	!x_ctrl(3) = <<Teq>> 
	!x_ctrl(4) =  2    !Gausss, dipole B field


      use_other_atm= .false.
      x_integer_ctrl(51)=2
      x_integer_ctrl(52)= 40
      x_ctrl(15)= <<kappa_v>> 
      x_ctrl(16)= 1     !bound fact  set at 1 for tau=2/3
      x_ctrl(20)= 1e-4
      x_ctrl(17)= 10
      x_ctrl(18)= 3000

      x_ctrl(19)=  <<Teq>> ! planetary equilibrium temperature

      
      which_atm_option = 'grey_irradiated'

      atm_grey_irradiated_atol = 1d-4
      atm_grey_irradiated_rtol = 1d-4
      atm_grey_irradiated_T_eq = <<Teq>> ! equilibrium temperature based on irradiation
         ! e.g., might use T_eq = Teff_star*(R_star/(2*Distance_to_star))^(1/2)
      atm_grey_irradiated_kap_v = 4d-3 ! opacity for irradiation
      atm_grey_irradiated_P_surf = 1d6 ! 1 bar in cgs units
      atm_grey_irradiated_max_tries = 50 ! limit on iterations
      trace_atm_grey_irradiated = .false.

      !for planetesimal absorbtion
      
      use_accreted_material_j = .true.
      use_other_adjust_mdot = .true.  !in run stars extras

      ! ratio of accreted am to keplerian
      x_ctrl(5) = 1.0
      ! log10 accretion rate in msun/secyer
      x_ctrl(6) = -3
      ! log10 max timestep in yr
      x_ctrl(7) = 7.0



      star_history_name = "<<historyName>>"

      profile_data_prefix = "<<historyName>>" 


      profile_data_suffix = '.data'


/ ! end of controls namelist


&pgstar
         



         
         


/ ! end of pgstar namelist
