[Mesa-users] Segfault while using other_kap subroutine

[Koushik Sen]

Hi Rob,

I have been able to get the code running.

You'll notice that if you print net_iso(ife56) with your existing setup
> (its really helpful to send the starting model and any networks you've
> altered so i don't need to tweak your code to get it to run) prints 0 which
> the comments for net_iso in the my_other_kap routine contain:
>
>  > value is 0 if the iso is not in the current net
>

This was not the case as the write (*,*) line to print net_iso(ife56) was
giving a value of 10, it was in the net. (as I had the required files given
in the website
https://github.com/orlox/mesa_input_data/tree/master/2016_binary_models,
but did not include in my mail, sorry!)

Also not sure why you my_other_kap is attempting to call my_other_kap
> inside the subroutine, thats going to be a recursive loop with no exit
> condition, maybe you meant to call kap_get to call mesa's implementation of
> the opacities?
>

This was the problem, along with some mismatch of arguments passed in the
kap_get() 10398 version and previous get_kap_type*() versions of MESA.
Thank you so much for your insights and help.

Yours faithfully,
Koushik
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20181113/d8089e71/attachment.html>