The GYRE Forums

Hi Ehsan, I just wrote a long reply but I got timed out and lost my text. I'm not gonna input it again, so I try to keep it short. I computed frequencies for the raw MESA model (*.data file, not *.data.GYRE file) and a version that I interpolated myself (going from 1600->5000 points). This interpola...

Hi Ehsan, There is not much you can do about this, my guess is that this is all due to numerics. Recall that a MESA model is only in hydrostatic equilibrium according to the very strict definition you use. If you change the definition, e.g. from Langrangian to Eulerian, you're not in HSE. If you cha...

Hi Rich, ah sorry, I was looking at the OSC format for my column description... Only columns 37 to 40 are unused, but 37 and 38 are unofficially reserved for nabla and nablarad, that leaves only 39 and 40 for the other variables... that's a no go.. I will just stick to the GSM format then. I don't k...

Hi Rich, I'm trying to run gyre_nad on an FGONG file, but I get Root Solving l n_pg n_p n_g Re(omega) Im(omega) chi n_iter n ABORT at line 794 <gyre_coeffs_evol:tau_thm_1>: tau_thm is undefined STOP Program aborted after the adiabatic calculations went fine. I'm not familiar with the source code at ...

Hi Rich, I get a similar error (see below) when trying to compute the frequencies in the attached model / input file "gyre.in.tar.gz". Honestly I would not exactly mind having the error, but it also freezes up my terminal, so I have to manually interrupt the process before proceeding. I'm trying to ...

Hi Rich, I think there is a strong case for a uniform scaling of the adiabatic vs nonadiabatic modes. I don't have a real opinion on what the scaling needs to be exactly, as long as it is fairly easy to compare them. Having them scaled so that the total kinetic energy is unity could be a good idea, ...

Hi Rich, is it possible that the scaling of the adiabatic and nonadiabatic eigenfunctions is different? If I run the attached gyre_nad.in and model with "gyre_nad" and "gyre_ad" I get a factor of ~6.32 x 10^7 difference. I tried looking into the code where the difference could be, but I couldn't fin...

Thanks! And thanks for the GSM format, to which I'm switching now. Concering this GSM: I see in the documentation on https://bitbucket.org/rhdtownsend/gyre/src/tip/doc/gsm-format.pdf 1. that "nabla" is defined as dlnp/dlnT. Shouldn't this read dlnT/dlnp? 2. that in gyre_gsm_file.f90, "c_p" is no lon...

Hi Rich, when I try: $:> gyre_ad < gyre_nad.in (yes, I use gyre_ad but call my input file gyre_nad, sorry about that) I get the following output: gyre_ad [hg] OpenMP Threads : 4 MPI Processors : 1 Initialization ============== At line 164 of file gyre_b3_file.f90 Fortran runtime error: Array bound m...