[mesa-users] unsmoothed element abundance for MESA version of 6208

[WU T.]

Hi everyone,

I used the version of 6208 to calculated the stellar evolution with 1Msun masses and 0.02 metallicity and encountered a problem that the abundance of hydrogen varying with masses shows an unsmoothed phenomenon. This unformal phenomenon is durning main sequence stars to red giant branch stars (the position of red bump). After the position of red bump, the unformal phenomenon disappear. I don know what does lead to this phenomenon and how to solve this problem? Who did encounter this problem and can give some advice for it?

I adopted the module of 1M_pre_ms_to_wd and the inlist file as follows:
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

&star_job

      show_log_description_at_start = .false.
     
      !show_net_species_info = .true.
      !show_net_reactions_info = .true.

      !create_pre_main_sequence_model = .true.

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

      kappa_file_prefix = 'gs98'
     
      change_lnPgas_flag = .true.
      new_lnPgas_flag = .true.

      change_initial_net = .true.     
      new_net_name = 'o18_and_ne22.net'

      pre_ms_relax_num_steps = 50
     
      new_surface_rotation_v = 2 ! solar (km sec^1)
      set_near_zams_surface_rotation_v_steps = 10 ! to turn on rotation when near zams

      !pgstar_flag = .true.


/ ! end of star_job namelist

&controls

      max_model_number = 20000

      initial_mass = 1.0
      initial_z = 0.02d0
     
      max_num_profile_models = 200
      photostep = 10
      profile_interval = 10
      history_interval = 1
      terminal_cnt = 10
      write_header_frequency = 10

      log_L_lower_limit = -1

      smooth_convective_bdy = .true.                 
      convective_bdy_weight = 1

      RGB_wind_scheme = 'Reimers'
      AGB_wind_scheme = 'Blocker'
      RGB_to_AGB_wind_switch = 1d-4
      Reimers_wind_eta = 0.0!0.7d0 
      Blocker_wind_eta = 0.0!0.7d0 
     
/ ! end of controls namelist
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

Best wiehes!

Tao Wu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20140626/1fbb2a6c/attachment.html>