

&star_job

      show_log_description_at_start = .false.

      load_saved_model = .true.
      saved_model_name = "planet_evolve_5_ME_0.15.mod"

      save_model_when_terminate = .true.
      save_model_filename = "planet_evolve2_5_ME_0.15.mod"
      
      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 = 100000000.0 ! in years   

      set_initial_dt = .true.    
      years_for_initial_dt = 1d3 ! in years make changes depending on mass and flux recieved

      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 = 49.1126209662 ! 3.d2 cm^2/g
        !irradiation_flux = 131931642.958   !555501654.562 ! 1.d9 erg/cm^2/s ! day side flux!!! Bill puts in 1/4 by hand


 	mesh_delta_coeff = 0.5
	!varcontrol_target = 1d-2
	!dt_years_for_steps_before_max_age = 1d6
      
	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 = 10000000000.0 
	!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
        !min_timestep_limit = 1e4
      


	T_mix_limit = 0

	! output to files and terminal
	photo_interval  = 20000
	profile_interval = 2000
	max_num_profile_models=2
	history_interval = 1
	!terminal_cnt = 10
	write_header_frequency = 10


   
        ! controls for atmospheric escape
        use_other_adjust_mdot= .true.
        x_ctrl(1)= 0.1 ! mass loss efficiency
        x_ctrl(2)= 0.1 ! 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) = 862.198735264 
	!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)= 0.0203613649674 
      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)=  862.198735264 ! 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 = 862.198735264 ! 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.

      star_history_name = "0.1_5_0.15_6.0"

      profile_data_prefix = "0.1_5_0.15_6.0" 


      profile_data_suffix = '.data'

      !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) = 8.0
      !knob for planetesimal injection
      x_ctrl(8) = 1.0





/ ! end of controls namelist


&pgstar
         


         
         


/ ! end of pgstar namelist
