The GYRE Forums

Hi Rich,

It is as you describe — the damping rates are very ragged (quite large variations of either sign). I can live with the second-order results for now. Thanks for your help!

Hi Rich,

Here's what I get with the patch applied for COLLOC_GL2 (and also Magnus). These are the absolute values of the damping rates (there's a sign change at the large discontinuity). Still seeing issues with the work integrals (both mine and GYRE's) for the 4th-order difference schemes.

Am I correct in understanding that the ZERO_R boundary condition is not intended to be used with evolutionary models that already specify a central point? I get errors with GYRE 5.2 in either case (just different errors). With the central point included, I obtain an assertion error: ASSERT 'cx%pt_i%...

Hi Rich, Sorry for the delay. The model in question can be found here: http://hyad.es/sendfile/payloads/Patched.GYRE Namelist follows: &model model_type = 'EVOL' file = '../Patched.GYRE' file_format = 'MESA' ! file_format = 'FGONG' ! data_format = '(1P5E16.9,x)' / &constants G_GRAVITY = 6.67408d-8 /...

Hi Rich, Thanks for the response. 2. and 3. I've checked my integrals against the values of W returned from the nonadiabatic calculation and they appear to be correct-ish — but the values of W also appear to not be consistent with the imaginary eigenvalue components: http://hyad.es/work_plot_2.png T...

Hi all, Has anyone tried determining damping rates from the eigenfunction work integrals? In particular, following Unno+ 1987 and/or Belkacem+ 2012 A&A 540, L7, we can cast the damping rate in terms of a work integral over the nonadiabatic eigenfunctions: http://hyad.es/work.png However, when I actu...

Certainly, please see attached.

When computing frequencies for some red giant models, I receive a SIGFPE signal (Floating-point exception - erroneous arithmetic operation). Any suggestions on how I should go around diagnosing why this is happening? For reference, my input namelists are as follows: &model model_type = 'EVOL' file =...

For what it's worth, I'm also running into memory issues from remeshing the grid even for computing radial p-modes (so no g-mode mixing expected). I can confirm that increasing the minimum frequency of the search interval relieves the memory pressure (I use Δν/10 as the lower bound of the interval, ...

Hi Rich, I'm running GYRE 5.2 (which came with MESA). So far I've observed this behaviour mostly on NFS mounts on my local HPC cluster, and also occasionally on department machines, but not on my laptop (although that might be because I don't run GYRE on my laptop very often). All of these run some ...