&star_job  ! HD49385   

      show_log_description_at_start = .false.

      create_pre_main_sequence_model = .true.    

      astero_just_call_my_extras_check_model = .true.

      pgstar_flag = .true.

      pre_ms_T_c = 2e5
     
!      history_columns_file = 'defaults/history_columns_ds_comp.list'
!      profile_columns_file = 'defaults/profile_columns_ds_comp.list'

/ ! end of star_job namelist


&controls  
      
   ! Gold tolerances

   use_gold_tolerances = .true.
   use_dedt_form_of_energy_eqn = .true.

     mesh_delta_coeff = 0.5 !0.75
     !  max_delta_x_for_merge = 0.05d0

      log_directory = 'LOGS/SMH100_approx_clean'
      star_history_name = 'history_SMH100_approx_clean.data'
      
      super_eddington_scaling_factor = 0
      super_eddington_wind_Ledd_factor = 0
  
      initial_mass =  5.d+0
      initial_z = 0.00d0
      star_mass_max_limit = 1d4
      
      mass_change =  1d-3 ! rate of accretion (Msun/year)
      !mass_change = 1d-1 ! rate of accretion (Msun/year)
      
      write_pulse_data_with_profile = .true.
      pulse_data_format = 'GYRE'
      !add_atmosphere_to_pulse_info = .false.
      !add_center_point_to_pulse_info = .false.
      
      profile_interval = 1
      history_interval = 1
      terminal_interval = 10
      write_header_frequency = 25

      max_num_profile_models = 99999
      
      ! move k_below_const_q inward from surface until q(k) <= max_q.
      ! then continue moving inward until reach logT(k) >= max_logT or q(k) <= min_q.
      max_logT_for_k_below_const_q = 6  !5
      max_q_for_k_below_const_q = 0.9 !0.96
      min_q_for_k_below_const_q = 0.9 !0.95

      use_other_energy_implicit = .true.
      okay_to_reduce_gradT_excess = .true.
      !eps_nuc_factor = 0.0 


     ! Testing different parameters to suppress superadiabacity even better (!!):

         gradT_excess_f1 = 1d-3
         gradT_excess_f2 = 1d-1
         gradT_excess_age_fraction = 0.9d0 ! fraction of old to mix with new to get next
      
         ! in some situations you might want to force alfa = 1.
         ! you can do that by setting gradT_excess_lambda1 < 0.
         ! following are for the normal calculation of gradT_excess_alfa

         gradT_excess_lambda1 = 1.0
         gradT_excess_beta1 = 0.35
        ! 
         gradT_excess_lambda2 = 0.5
         gradT_excess_beta2 = 0.25
         
         gradT_excess_dlambda = 0.3
         gradT_excess_dbeta = 0.3
         
              

      mixing_length_alpha = 1.9  ! 2
         ! the mixing length is this parameter times a local pressure scale height.
         ! to increase R vs. L, decrease alpha


      

      atm_option = 'table'
      atm_table = 'photosphere'

      ! create_atm_max_step_size =  0.1 ! 0.05  in units of log10_tau
      

      !when_to_stop_rtol = 1d-13		!TEST
      !when_to_stop_atol = 1d-13		!TEST

      min_timestep_limit = 1d-15 ! 1d-6 (seconds)
         ! stop if need timestep smaller than this limit

      max_number_retries = -1 ! ignore if < 0

      !min_timestep_factor = 0.80d0
      !max_timestep_factor = 1.20d0
      !max_timestep_factor = 0.80d0
      !max_timestep= -1
      !delta_lgRho_cntr_limit = 0.0004	 
      !delta_lgT_cntr_limit = 0.00008

      
/ ! end of controls namelist


&pgstar
         
      ! top level controls
         
         show_TRho_Profile_legend = .true.
         Profile_Panels2_xaxis_name = 'mass'
         Grid1_win_flag = .true.
         Grid1_win_width = 9
      
      !Kipp_win_flag = .true.
      !Kipp_show_burn = .false.
      !Kipp_file_flag = .false.
      !Kipp_show_luminosities = .true.

      rti_win_flag = .true.

      rti_file_flag = .true.


/ ! end of pgstar namelist
