Bug/problem reports for any of the GYRE executables (gyre_ad, gyre_nad, etc)
-
warrick
- Posts: 84
- Joined: Wed Aug 28, 2013 2:47 am
Post
by warrick » Mon Mar 03, 2014 5:26 am
Hi,
I'm trying to do a calculation that involves the "energy integral" defined by equations (3.139) in Aerts, JCD & Kurtz, and (4.47) in JCD's lecture notes:
Is this number available as default output in GYRE? I tried using 'E' in the summary output, but its value always seems to be one, at least in adiabatic calculations. I checked the definition in the source (since it's so easy when it's written so well!) and it
looks like the same thing.
W
-
rhtownsend
- Site Admin
- Posts: 400
- Joined: Sun Mar 31, 2013 4:22 pm
Post
by rhtownsend » Mon Mar 03, 2014 8:21 pm
Hi Warrick --
In fact, all eigenfunctions are normalized to have E = 1 --- hence your results!
What you might want to consider is the normalized inertia E_norm -- see eqn. 3.140 of Aerts+2010.
cheers,
RIch
-
warrick
- Posts: 84
- Joined: Wed Aug 28, 2013 2:47 am
Post
by warrick » Tue Mar 04, 2014 1:22 am
Rich,
That's what I figured. It seems a bit strange for this formula to be based on the norm of the eigenvector, given that it isn't uniquely normalized... I'll dig around in the notes to see what's going on. Is the value also always 1 in the non-adiabatic calculation?
W
-
rhtownsend
- Site Admin
- Posts: 400
- Joined: Sun Mar 31, 2013 4:22 pm
Post
by rhtownsend » Tue Mar 04, 2014 9:24 am
Yes, it's always 1 -- see the mode_t_ routine in src/common/gyre_mode.fpp
cheers,
Rich
