
&star_job

      load_saved_model = .true.
      saved_model_name = '/home/wangzhht/project/He_1Z/He_2.5.mod'

      show_log_description_at_start = .false.

      show_timestep_limit_counts_when_terminate = .true.

      write_profile_when_terminate = .true.
      filename_for_profile_when_terminate = 'LOGS1/profile.data'


      !set initial age and model number to 0 again
      set_initial_age = .true.
      initial_age = 0.0
      set_initial_model_number = .true.
      initial_model_number = 0

      !OPACITY, NETWORK, RATE, etc.
      !network
      change_net = .true.
      new_net_name = 'approx21.net'
      show_net_species_info = .true.

      new_rotation_flag = .true.
      change_rotation_flag = .true.
      set_initial_surface_rotation_v = .true.
      new_surface_rotation_v = 0.0d0
      num_steps_to_relax_rotation = 50

      !use jina
      set_rates_preference = .true.
      new_rates_preference = 2

/ ! end of star_job namelist

&eos
      ! eos options
      ! see eos/defaults/eos.defaults

/ ! end of eos namelist

&kap
      ! kap options
      ! see kap/defaults/kap.defaults
      use_Type2_opacities = .true.
      Zbase = 1.42d-3

/ ! end of kap namelist


&controls
     x_ctrl(1) = 1.00d0 ! 1Eddinton limit

     use_other_wind=.false.
     !zams_filename = './zams_z0.0142m2_y0.2703.data' 
     initial_z = 1.42d-3
     opacity_max = 0.5
    
     warn_when_large_rel_run_E_err = 1d99

 
     use_gold_tolerances = .false.
     use_dedt_form_of_energy_eqn = .false.


      !OUTPUT OPTIONS
      photo_directory = 'photos1'
      log_directory = 'LOGS1'
      history_interval = 1
      terminal_interval = 10
      write_header_frequency = 100
      photo_digits = 5
      photo_interval = 0


      star_history_dbl_format = '(1pes32.16e3, 1x)'
      star_history_int_format = '(i32, 1x)'
      star_history_txt_format = '(a32, 1x)'

      star_history_name = 'history.data'
      write_profiles_flag = .false.

      !WHEN TO STOP
      gamma_center_limit = 1.0d1 !crystallization for WD

      !PHYSICS
      !which_atm_option = 'simple_photosphere'
      !which_atm_off_table_option = 'simple_photosphere'
      Pextra_factor = 2.0d0

      !CONVECTION + MIXING
      do_conv_premix = .true.
      conv_premix_avoid_increase = .true.

      mixing_length_alpha = 1.93d0
      mlt_option = 'Mihalas'
      mlt_use_rotation_correction = .true.
      okay_to_reduce_gradT_excess = .true.

      use_ledoux_criterion = .true.
      alpha_semiconvection = 1d0
      thermohaline_coeff = 1d0
      mixing_D_limit_for_log = 1d-50

      !H core overshoot MIST+Brott11: implemented in run_star_extra.f

      !ROTATION
      fitted_fp_ft_i_rot = .true.
      w_div_wcrit_max = 0.9d0


      am_nu_ST_factor = 1.0d0 !MANOS: 0:noTS, 1:TS
      D_visc_factor = 0.0d0
      D_ST_factor = 0.0d0
      D_SH_factor = 0.0d0 !this is different from MIST2
      D_GSF_factor = 1.0d0
      D_ES_factor = 1.0d0
      D_SSI_factor = 1.0d0
      D_DSI_factor = 1.0d0
      am_D_mix_factor = 0.0333333d0 ! f_c in Heger+2000
      am_gradmu_factor = 0.05d0 !f_mu = 0.05 in Heger+2000

      !MASS LOSS
      cool_wind_full_on_T = 1.0d0
      hot_wind_full_on_T = 2.0d0
      hot_wind_scheme = 'Dutch'
      cool_wind_RGB_scheme ='Dutch'
      cool_wind_AGB_scheme = 'Dutch'
      Dutch_scaling_factor = 0.667d0
!      max_wind = 1.0d-3


      !CORE MASS DEFINITION
      he_core_boundary_h1_fraction = 1d-2
      c_core_boundary_he4_fraction = 1d-1
      o_core_boundary_c12_fraction = 1d-1
      min_boundary_fraction = 1d-1

      !step scheme , alpha_ov = 0.1
      overshoot_scheme(1) = 'step'
      overshoot_zone_type(1) = 'any'
      overshoot_zone_loc(1) = 'any'
      overshoot_bdy_loc(1) = 'top'
      overshoot_f(1) = 0.11
      overshoot_f0(1) = 0.01

      !MESH AND TIMESTEP PARAMETERS
      mesh_delta_coeff = 1.0d0
      max_allowed_nz = 20000
      
      varcontrol_target = 1d-3
      min_allowed_varcontrol_target = 1d-4  !!terminate if varcontrol_target < min_allowed_varcontrol_target
      time_delta_coeff = 1d0  !!smaller forces smaller timesteps giving better time resolution. 
      !!the new control for temporal resolutions
      max_timestep_factor = 1.5d0
      

      solver_iters_timestep_limit = 20

      !stop when the center mass fraction of h1 drops below this limit
      !xa_central_lower_limit_species(1) = 'h1'
      !xa_central_lower_limit(1) = 0d0
      !xa_central_lower_limit_species(2) = 'he4'
      !xa_central_lower_limit(2) = 1d-3
      !xa_surface_lower_limit_species(2) = 'he4'
      !xa_surface_lower_limit(2) = 1d-3
/ ! end of controls namelist
