[mesa-users] questions about the "Eulerian" scheme for eps_grav

[Josiah Schwab]

Hi Masanobu,

> Please check the attached file “run_star_extras.f”. I succeeded in
> compiling, but just after running (./rn), I got the following message
> and the calculation terminated:
> ----------------------------------------------------------------------------------------------------------------------
> Program received signal SIGSEGV: Segmentation fault - invalid memory
> reference.
>
... snipped 10 lines ...
> ./rn: line 2: 21670 Segmentation fault: 11 ./star_vlm
> ----------------------------------------------------------------------------------------------------------------------
> I apologise if my modifications were not appropriate. I have not
> solved this problem yet.

I only had a minute to glance at your file, but you need to make sure to
get the star pointer (s) associated with a given id.

That is, before you can use "s" within a subroutine, you need to do 

    ierr = 0
    call star_ptr(id, s, ierr)
    if (ierr /= 0) return

If you take a look at the stock run_star_extras.f in r7184, you'll see
this idiom is now the first step of most of the extras_* subroutines.

Hope that helps,
Josiah