[mesa-users] kap subroutine returns NaNs

[Jared Brooks]

Problem solved!

If calling subroutines from the kap module (without running mesa/star), you
need to load the opacity tables by calling the kap_init routine first. Here
is what I used, with default values:

call kap_init('gn93','gn93_co','lowT_fa05_gs98',4.1d0,4.0d0,3.
8d0,.true.,'',info)

Here is a link to this explanation:
http://mesastar.org/teaching-materials/using-kap-module/

-Jared

On Wed, Oct 1, 2014 at 7:00 PM, Jared Brooks <jcbrooks at physics.ucsb.edu>
wrote:

> Hey mesa-users,
>
> I'm trying to use the subroutine kap_get_elect_cond_opacity available in
> mesa/kap/public/kap_lib.f, and the three output values, kap, dlnkap_dlnRho,
> dlnkap_dlnT, consistently return NaNs, while ierr=0. I'm having trouble
> tracking this down because there is no error message, and I quickly got
> lost following all the calculations this subroutine uses from kap/private.
>
> I'm using version 5596 here.
>
> -Jared
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20141002/bda91446/attachment.html>