The GYRE Forums

Hi JJ --

I'll take a look at the crashing case in the coming days, to see if I can pinpoint why things are going awry.

cheers,

Rich

Hi JJ -- Further to Warrick's comments, it looks like both lag_T_eff and xi_r are zero at the stellar surface, causing the error. This shouldn't usually happen, but may be a consequence of using the MAGNUS_GL4 diff scheme. Could you try using MAGNUS_GL2 instead and see if that works better? (The hig...

Hi Gautam -- The ability to calculate the first-order effects of the Coriolis force has been in GYRE since at least v6.0. The domega_rot output item gives the (additive) correction that the user can apply (via post-processing) to the calculated dimensionless frequencies: omega -> omega + domega_rot....

I recommend using the 'keep_surface_point_for_pulse_data = .true.' option -- this will make sure that the GYRE model includes the outermost point of the MESA model, which in this case corresponds to tau=2/3.

Hi Anwesha --

What model atmosphere is being used for the MESA calculations -- i.e., what is the setting for atm_option? It would help if you could post your MESA inlist file.

Also, regarding your plot, what is the x-axis?

cheers,

Rich

Hi, I would like to calculate the optical depth of the surface point of the gyre model I am using T^4 = (Teff^4)*(2/3 + tau)*(3/4) for T I am using (addata['T'])[-1] for the surface point, Teff I am using the Teff of the input model. Now that leads me to a tau of ~ 32, whereas the input model has a...

Hi Rich, just another quick question, Do the underscore terms for 'c_thk' and 'c_thn' in the structure equation symbol list correspond to thick and thin optical depth respectively? If not, then what do they mean? what should be a tau value for it to be considered as thick and what should be a tau v...

Hi Anwesha --

Code files in GYRE reside in the $GYRE_DIR/src subdirectory. For the definition of the differential work (and work integrals), take a look at $GYRE_DIR/src/mode/gyre_wave.fpp

cheers,

Rich

Hi folks -- I'm pleased to announce the release of GYRE 7.1, which offers a treatment of turbulent convective damping and an important bug fix to the 'KAWALER' zeta_scheme. Source code: https://github.com/rhdtownsend/gyre/releases/download/v7.1/gyre-7.1.tar.gz Release notes: https://github.com/rhdto...

Hi Cayden -- Good question -- this is little buried in the docs, but the way to do this is using the alpha_frq parameter of the &tide namelist group. The inertial-frame forcing frequency used in gyre_tides is sigma = alpha_frq * k * Omega_orb ...where Omega_orb is the orbital frequency and k the har...