[Mesa-users] Opacity discontinuity

[Zhang, Zhaoxi M. (Michael)]

Hello MESA users,


I'm using OP opacities to model a white dwarf, and getting discontinuities in the opacity (see plot).  The opacity suddenly drops from 20 to 2 within a single grid cell before climbing back up, all within a very narrow log(1-q) range.  These discontinuities are wreaking havoc on the timestep, as the structure around the discontinuities is forced to change on a very small timescale.


Does anyone know what might be causing these discontinuities around T=6960 and rho=2.3e-6?  I don't get any errors from OP.  I'm diving into the tables, but thought that I should ask here before diving too deep.



Cheers,


Michael
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20180309/758e4658/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: opacity_discontinuity.png
Type: image/png
Size: 31554 bytes
Desc: opacity_discontinuity.png
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20180309/758e4658/attachment.png>