[Mesa-users] Segfault while using other_kap subroutine

[Rob Farmer]

Hi,
Without knowing what you did to make the basic2.net i would guess you are
missing fe56 in your nuclear network thus the line:

eff_metallicity = 0.017d0*xa(net_iso(ife56))/0.00124d0

in your run_star_extras.f90 will access an invalid array element. I would
try just using

change_net=.true.
new_net_name='approx21.net'

in your star_job inlist, to use a minimal nuclear network which contains
fe56.

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

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?

Rob




On Tue, 13 Nov 2018 at 14:34, Koushik Sen via Mesa-users <
mesa-users at lists.mesastar.org> wrote:

> Hi everyone,
>
> I am trying to use my own my_other_kap subroutine by editing the
> run_star_extras.f file as instructed in the tutorial in the MESA website
> but running into the following error.
>
> #####################################################
> ksen at aibn165:/vol/aibn165/data1/ksen/code/mesa/binary/1_project/5_single$
> ./rn
> DATE: 2018-11-13
> TIME: 14:16:07
>                                          version_number       10398
>  read inlist_project
>  read inlist_extra_controls
>  use_eosPTEH_for_low_density T
>  use_eosPTEH_for_high_Z T
> ./rn: line 9: 3198757 Segmentation fault      (core dumped) ./star
> DATE: 2018-11-13
> TIME: 14:16:08
> ######################################################
>
> Note: I am editing the files of 2016_bonary_models (
> https://github.com/orlox/mesa_input_data/tree/master/2016_binary_models)
> by Pablo Marchant. While compiling the code in the exact form as in the
> website, I had run into an error which Pablo told that the current version
> of MESA in incompatible with this other_kap subroutine. Subsequently, I
> removed the use_other_kap = .true. and the corresponding subroutine in the
> run_star_extras.f file and things were running. But now, I want to use my
> own subroutine my__other_kap instead of the opacity tables "BROTT_*.data"
> included in the website. I was able to compile the code using the
> my_other_kap subroutine as well by tweaking the pointers as required for
> the 10398 version.
>
> I also checked that this error is only happening when I am using the
> other_kap routine. I have included the inlist files and the
> run_star_extras.f subroutine with this mail.
>
> ksen at aibn165:/vol/aibn165/data1/ksen/code/mesa/binary/1_project/5_single$
> uname -a
> Linux aibn165 4.15.0-38-generic #41-Ubuntu SMP Wed Oct 10 10:59:38 UTC
> 2018 x86_64 x86_64 x86_64 GNU/Linux
> ksen at aibn165:/vol/aibn165/data1/ksen/code/mesa/binary/1_project/5_single$
> gfortran -v
> Using built-in specs.
> COLLECT_GCC=/vol/aibn165/data1/ksen/code/mesasdk/bin/gfortran.exec
>
> COLLECT_LTO_WRAPPER=/vol/aibn165/data1/ksen/code/mesasdk/bin/../libexec/gcc/x86_64-pc-linux-gnu/7.2.0/lto-wrapper
> Target: x86_64-pc-linux-gnu
> Configured with: /root/mesasdk-src/gcc/configure CC=gcc
> --build=x86_64-pc-linux-gnu --prefix=/root/mesasdk --with-gmp=/root/mesasdk
> --with-mpfr=/root/mesasdk --with-mpc=/root/mesasdk
> --enable-languages=c,c++,fortran --disable-multilib --disable-nls
> --disable-libsanitizer --enable-clocale=generic
> Thread model: posix
> gcc version 7.2.0 (GCC)
> ksen at aibn165:/vol/aibn165/data1/ksen/code/mesa/binary/1_project/5_single$
> echo $MESASDK_ROOT
> /vol/aibn165/data1/ksen/code/mesasdk
> ksen at aibn165:/vol/aibn165/data1/ksen/code/mesa/binary/1_project/5_single$
> echo $PATH
>
> /vol/aibn165/data1/ksen/code/mesasdk/bin:/usr/local/sbin:/usr/local/bin:/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/usr/local/games:/snap/bin
> ksen at aibn165:/vol/aibn165/data1/ksen/code/mesa/binary/1_project/5_single$
> echo $MESA_DIR
> /vol/aibn165/data1/ksen/code/mesa
>
> Awaiting your reply,
>
> Yours faithfully,
> Koushik Sen
> P.hD. Student
> Argelander Institute for Astronomy
>
>
>
>
>
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20181113/6f8157e0/attachment.html>