Confusing Xi_r results upon testing solar Model S
Posted: Wed Oct 09, 2013 1:11 am
I've been testing gyre on the standard solar Model S. My objective is to compute the derivatives of the mode displacement eigenfunctions. To give fair warning, this is my first foray into asteroseismology so it is quite possible that I have something seriously wrong in my understanding. With that caveat in mind, I have computed the l = 0-3 modes in the 1-10 mHz range using model S supplied in fgong format from Jorgen Christensen-Dalsgaard's website. When I look at, for example, the mode output for the first l=0 mode in this range, we see that at a fractional radius of about 1, the xi_r displacement rises to be on the order of 100. The documentation on the output variables page on the gyre wiki indicates that that this xi_r displacement perturbation is in units of the stellar radius. I can think of 3 possibilities for this problem: either I don't understand what Xi_r represents, which is entirely possible; the units in the documentation are incorrect; or there is something wrong with gyre, either on my computer or in general. I've attached the input file and an example mode output. If anybody has an idea as to what has occurred here, the information would be much appreciated.
Cheers,
Jack Muir
Cheers,
Jack Muir