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

[Warrick Ball]

Hi again,

Okay, so I've determined two things, and the problem is solved!  (Or, 
rather, it never actually existed.)

First, that block of output isn't helpful.  When I run ADIPLS on it's own, 
with iekinr definitely set to 1 (based on output), I still get the lead 
block that says iekinr=0.  This seems to be based on the initialization of 
parameters in adipls/adipack.c/adipls/adipls.c.d.f.  iekinr is set to zero 
around line 876, and printed to screen around line 1176, before there's a 
chance for anything to modify it.  So the displayed value is the raw 
initialization, before the input parameters are read.

Second, by setting iekinr=1 in adipls.c.pruned.in, I seem to be getting 
the correct results.  That is, the mode energies are being normalized at 
the photosphere.  I checked this by doing a comparison of MESA's output 
with a manual run of ADIPLS on the FGONG for the best-fit model along one 
track.

To summarize, I don't think there was a problem in the first place! 
iekinr *is* controlled by setting it in adipls.c.pruned.in.  ADIPLS just 
doesn't tell you what it's actual value is!

Sorry to have bothered the forum with this...

Cheers,
Warrick


On Thu, 15 May 2014, Warrick Ball wrote:

> 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


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