
! 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

      read_extra_star_job_inlist1 = .true.
      extra_star_job_inlist1_name = '../inlist_test_suite'
      
      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

      ! check for retries and backups as part of test_suite
      ! you can/should delete this for use outside of test_suite
      max_number_backups = 5
      max_number_retries = 6




      prune_bad_cz_min_Hp_height = 0 ! lower limit on radial extent of cz
      remove_mixing_glitches = .true. ! if true, then okay to remove gaps and singletons
      
               ! the following controls are for different kinds of "glitches" that can be removed
            
               clip_D_limit = 0 ! zero mixing diffusion coeffs that are smaller than this         

               okay_to_remove_mixing_singleton = .true.
         
               min_convective_gap = -1 ! 0.1 
                  ! close gap between convective regions if smaller than this (< 0 means skip this)
                  ! gap measured radially in units of pressure scale height
         
               min_thermohaline_gap = 0
                  ! close gap between thermohaline mixing regions if smaller than this (< 0 means skip this)
                  ! gap measured radially in units of pressure scale height
               
               min_thermohaline_dropout = 0
               max_dropout_gradL_sub_grada = 1d-3
                  ! if find radiative region embedded in thermohaline,
                  ! and max(gradL - grada) in region is everywhere < max_dropout_gradL_sub_grada
                  ! and region height is < min_thermohaline_dropout
                  ! then convert the region to thermohaline
         
               min_semiconvection_gap = 0 
                  ! close gap between semiconvective mixing regions if smaller than this (< 0 means skip this)
                  ! gap measured radially in units of pressure scale height
         
               remove_embedded_semiconvection = .false.
                  ! if have a semiconvection region bounded on each side by convection,
                  ! convert it to be convective too.
         

      initial_mass =  0.003
      initial_z = 0.02d0
      
      max_model_number = 12
      max_years_for_timestep = 0.1
      
      mesh_delta_coeff = 0.7

      which_atm_option = 'tau_10_tables'

      photostep = 50
      profile_interval = 100
      history_interval = 10
      terminal_cnt = 10
      write_header_frequency = 10

/ ! end of controls namelist
