
! inlist_rc1

! load 0.003, then remove mass = desired rocky core mass.
!   e.g., if want rocky core = 10 M_earth = 3e-5 Msun
!   the change to mass = 3e-2 - 3e-5 = 0.00297

&star_job

      mesa_data_dir = '../../../data'
      
      show_log_description_at_start = .false. 
      
      load_saved_model = .true.
      saved_model_name = '../../../data/star_data/very_low_mass_COND_tau100_models/0.003Msun.mod'

      save_model_number = 11
      save_model_filename = 'rc1.mod'

      kappa_file_prefix = 'gs98'
      kappa_lowT_prefix = 'lowT_Freedman11'

      relax_mass_scale = .true.
      new_mass = 0.00297


/ ! end of star_job namelist

&controls
      
      max_model_number = 12
      max_years_for_timestep = 0.1

      initial_mass =  0.003
      initial_z = 0.02d0

      !newly_nonconvective_weight = 0
      !smooth_newly_nonconvective = .false.

      ! check for retries and backups as part of test_suite
      max_number_backups = 5
      max_number_retries = 6
      
      mesh_delta_coeff = 0.7

      which_atm_option = 'tau_10_tables'

      photostep = 100
      profile_interval = 100
      log_cnt = 10
      terminal_cnt = 10
      write_header_frequency = 10
      
      !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

      report_why_dt_limits = .false.
      report_all_dt_limits = .false.
      
      !show_mesh_changes = .true.
      !mesh_dump_call_number = 131
      
      !max_model_number = 21100
      !max_model_number = 11000
      !max_model_number = 18000 ! << increase HR xmax

      !trace_evolve = .true.
      
      ! hydro debugging
      !hydro_check_everything = .true.
      !hydro_inspectB_flag = .true.
      
      !hydro_numerical_jacobian = .true.
      !hydro_save_numjac_plot_data = .true.
      !small_mtx_decsol = 'lapack'
      !large_mtx_decsol = 'lapack'
      !hydro_dump_call_number = 195

/ ! end of controls namelist
