[mesa-users] How to use other_D_mix?

[Ehsan Moravveji]

Dear MESA users,
I am trying to plug in an ad-hoc linear diffusive mixing coefficient,
through using other_D_mix; however, I have little progress, since I don't
get what I expect!
So, typical, right? :-)

In brief, I use v.7385, and set use_other_D_mix to .true. in my inlist.
Also, I set
         s% other_D_mix => Change_D_mix
in extras_controls(). In the latter procedure, I have the following simple
do loop:

         do k = 1, nz
            ok_to_mix = .false.     ! reset to default at each mesh
            ok_to_mix = ((s% mixing_type(k) > softened_convective_mixing)
.or. &
                         (s% mixing_type(k) == 0))
            if (.not. ok_to_mix) cycle

            s% D_mix(k) = D_0 + v_diff * s% r(k)
            s% D_mix_non_rotation(k) = s% D_mix(k)
            print*, k, s% D_mix(k)
         enddo

With the print statement, I just make sure that this procedure works, and
some mesh points, k, that are not convective, overshooting and softening
are correctly selected to tweak s% D_mix(k). The printed value of D_mix
that I want to set is also OK. However, when I plot the profile of
log_D_mix and log_D_mix_non_rotation, I cannot see the envisaged increasing
linear trend in radiative parts of the model. A plot is attached.

Digging and "grep"ing in star/private/mix_info.f90 also did not take me any
further.
So, my question is: Which of the star_info variables are allowed to be
changed through other_D_mix procedure?
How can I see what I expect? ;-)

Kind regards
Ehsan
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150301/a923a07b/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: log_D_mix.png
Type: image/png
Size: 47336 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150301/a923a07b/attachment.png>