[mesa-users] Bug in MLT?

[Jim Fuller]

Hi all,

I may have found a bug in MLT, or maybe I just don't understand how the
code works. I am using mesa8845. In star/private/mlt_info.f90, there is a
statement:

         if (s% v_flag) then
            if (s% v(k) > s% csound_face(k)) cv_old = 0d0
         end if

However, I think this should say something like

         if (s% v_flag) then
            if (s% v(k)/s% csound_face(k) > s% max_v_div_cs_for_convection)
cv_old = 0d0
         end if

Am I wrong? I find artificial jumps in convection velocities in regions
with large velocities if use the upper statement, but the bottom statement
appears to fix it. Any feedback would be helpful. Best,

-Jim
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20160913/25720c19/attachment.html>