[mesa-users] How does one set iekinr in ADIPLS within MESA?

[Warrick Ball]

Hi,

I'm doing some asteroseismology involving mode inertia (or energies) and 
I'd like them to be normalized by their total photospheric displacements. 
In ADIPLS, this is done setting iekinr=1, rather than the default value of 
0 (which corresponds to normalizing by the vertical displacement at the 
surface).

As a simple test, I tried setting iekinr to 1 in the adipls.c.in and 
adipls.c.pruned.in files included in the star/astero/work directory, and 
then trying a run.  In adipls.c.in, it's found on line 38:

37:  istdpr,nout,nprcen,irsord,iekinr
38: 9,10,,20,,,,,,,,,     @

so I replaced the second line with

38: 9,10,,20,1,,,,,,,,     @

But then, when I compile and run, the first block of output, which shows 
the ADIPLS parameters, has this line

   istdpr,nout,nprcen,irsord,iekinr?
            6          50           0           0           0

So I presume iekinr isn't being set.

I tried grepping for "iekinr" in star/astero/src/*.f and 
star/astero/work/src/*.f but got nothing.  I'll start trawling through the 
code to see if I can see where the input files are read or linked or 
something.  Even if I have to hardcode iekinr for the time being, that's 
okay.  But I currently can't find any control at all!

Cheers,
Warrick



------------
Warrick Ball
Postdoc, Institut für Astrophysik Göttingen
wball at astro.physik.uni-goettingen.de
+49 (0) 551 39 5069