The GYRE Forums

Hi Rich, thank you so much. Before I try them out, I just want to see if I understand it correctly. So, the write_pulse_data_with_profile = .true. gives output of the .data.GYRE file keep_surface_point_for_pulse_data = .true. means it keeps the point where tau=1d-3 (as default for the integration of...

Hi Rich,

I have attached the main inlist which I am using for cooling the model. For the graph attached before the x-axis is the log of pressure from the center (on the left) to the surface (on the right).

Thank you for all the help

Hi Rich, I am using a .GYRE model from MESA. I have calculated a running work integral and plotted it (picture attached). I have the region marked where I want more points so that it is well resolved. I don't know how to approach that. Is there a way to do it? (I have tried playing around with wosc,...

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 va...

Hi Rich,

thank you so much. I will look into them.

Anwesha

Hi, Where are the code files located in gyre? I want to see the formula used for the work integrals (non-adiabatic, quasi-adiabatic). I am just learning how to do analytical calculations using raw data and would like to match the calculations with those outputted as data (Like using raw values to ca...

Hi Rich,

thank you so much for the clarification. So, if I want the eigen frequencies in 'rad/s' to compare with lamb and Brunt Vaisala (both of which I have found in rad/s) that would simply mean I will be using sigma (without having any additional factor of 2pi)?

Hi, I am a bit confused with the units of the non-adiabatic exp(-i*sigma*t) term mode period = 2*pi/sigma (from documentation) --> so here sigma is the angular frequency having unit (1/rad) now the output omega by gyre is dimensionless. sigma = omega*sqrt(GM/R^3) --> sqrt(GM/R^3) has unit of (1/s) -...