[Mesa-users] Lithium diffusion problem

[WAN AISHAH WAN HARUN STUDENT]

Hi everyone,
I have a problem in calculating my lithium abundance using this diffusion
control.
I used pp extras in my net reaction,so I include this in my &control list

       ! atomic diffusion
      do_element_diffusion = .true. ! determines whether or not we do
diffusion
      diffusion_dt_limit = 7d11 ! no element diffusion if dt < this limit
(in seconds)
      diffusion_T_full_on = 1d3
      diffusion_T_full_off = 1d3

      diffusion_calculates_ionization = .true.

       diffusion_num_classes = 11 ! number of classes of species for
diffusion calculations
      diffusion_class_representative(1) = 'h1'
      diffusion_class_representative(2) = 'he3'
      diffusion_class_representative(3) = 'he4'
      diffusion_class_representative(4) = 'li7'
       diffusion_class_representative(5) = 'be7'
      diffusion_class_representative(6) = 'b8'
      diffusion_class_representative(7) = 'c12'
      diffusion_class_representative(8) = 'n14'
      diffusion_class_representative(9) = 'o16'
      diffusion_class_representative(10) = 'ne20'
      diffusion_class_representative(11) = 'mg24'

      ! in ascending order.  species goes into 1st class with A_max >=
species A
      diffusion_class_A_max(1) = 2
      diffusion_class_A_max(2) = 3
      diffusion_class_A_max(3) = 4
      diffusion_class_A_max(4) = 7
      diffusion_class_A_max(5) = 7
      diffusion_class_A_max(6) = 8
      diffusion_class_A_max(7) = 12
      diffusion_class_A_max(8) = 14
      diffusion_class_A_max(9) = 16
      diffusion_class_A_max(10) = 20
      diffusion_class_A_max(11) = 24

      diffusion_use_isolve = .true.
      diffusion_rtol_for_isolve = 1d-4
      diffusion_atol_for_isolve = 1d-5
      diffusion_maxsteps_for_isolve = 1000
      diffusion_isolve_solver = 'ros2_solver'

I read the documentation,but I am still quite lost here on how to use the
diffusion control.
Results produced for He4 looks quite alright to me.However the lithium
abundance at surface and convection zone calculated in diffusion is higher
than the non-diffusion case(which is weird because li7 should be diffuse
towards the center too ).Please let me know where did I do wrong here.
Attached herewith are plotting results for he4 and li7  abundances vs
age(from PMS)


Thanks
Aisha
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190310/2ecd83dd/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: he4_diff.png
Type: image/png
Size: 56500 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190310/2ecd83dd/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: li7_surface.png
Type: image/png
Size: 31846 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190310/2ecd83dd/attachment-0001.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: li7_cz.png
Type: image/png
Size: 34935 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20190310/2ecd83dd/attachment-0002.png>