
! inlist_cool_brown_dwarf

&star_job
      show_log_description_at_start = .false.
      
      create_pre_main_sequence_model = .true.
      pre_ms_T_c = 3d5

      save_model_when_terminate = .true.
      save_model_filename = 'final.mod'
      
      write_profile_when_terminate = .true.
      filename_for_profile_when_terminate = 'final_profile.data'

      change_initial_net = .true.      
      new_net_name = 'pp_extras.net'

      set_uniform_initial_composition = .true.
      initial_h1 = 0.70597d0 ! for Z = 0.019
      initial_h2 = 3d-5
      initial_he3 = 2d-5
      initial_he4 = 0.27498d0
      initial_zfracs = 3 ! GS98_zfracs = 3

      eos_file_prefix = 'mesa'
      kappa_file_prefix = 'gs98'

      show_log_description_at_start = .false. 

      set_initial_dt = .true.
      years_for_initial_dt = 1d-7
      
      change_lnPgas_flag = .true.
      new_lnPgas_flag = .true.
      
      !pgstar_flag = .true.


/ ! end of star_job namelist


&controls

      ! check for retries and backups as part of test_suite
      ! you can/should delete this for use outside of test_suite
         max_number_backups = 100
         max_number_retries = 200
            
         max_model_number = 4000
         
         
            
      initial_mass = 0.05
      initial_z = 0.019d0

      varcontrol_target = 3d-5
      !max_years_for_timestep = 1d7
      
      mixing_length_alpha = 1.89d0 ! based on solar model
      MLT_option = 'Henyey'
      
      mesh_delta_coeff = 0.5

      which_atm_option = 'tau_100_tables'

      photostep = 50
      profile_interval = 50
      history_interval = 1
      terminal_cnt = 1
      write_header_frequency = 10
      
      max_age = 10.1d9


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

      ! hydro debugging
      !hydro_check_everything = .true.
      !hydro_inspectB_flag = .true.
      
      !hydro_numerical_jacobian = .true.
      !hydro_save_numjac_plot_data = .true.
      !hydro_dump_call_number = 195
      !hydro_dump_iter_number = 5

      !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_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
         
      ! top level controls
         
         pgstar_sleep = 0

         !pause = .true. 
            ! if true, the code waits for user to enter a RETURN on the command line
      


/ ! end of pgstar namelist
