[Mesa-users] Non-adiabatic GYRE and input format?

Mon Nov 18 20:00:42 UTC 2024

Dear all,

I have been trying to compute oscillations with GYRE for a MESA model 
that I saved in the OSC format, to no avail. The short question is: can 
GYRE do non-adiabatic calculations from an OSC input file?

More details for the real aficionados: while the adiabatic calculation 
seems to go smoothly, the non-adiabatic calculation fails right as soon 
as it actually tries to compute complex eigenvalues. From the error 
message, it appears some variables are missing from the OSC files. I 
have attached my GYRE inlist, the output I get (both stdout and error), 
and my 5MB OSC model file can be found here 
https://drop.infini.fr/r/j7Ey8qmRM7#uX/Rewom6a73q4D4EHmVbq063lJ53waucJ6gasUs4aI=. 

I am using GYRE 7.2.1, the model was generated with MESA 23.05.1 on a 
machine running Ubuntu 22.04.1 LTS.

I also noted that all the non-adiabatic test cases use the MESA input 
(the one that's called GYRE in MESA), is that the only one that works 
for the nonad case? I quickly checked the documentation of both formats, 
they seem to contain almost the exact same info except that OSC has 
(eps_nuc + eps_grav) while GYRE/MESA has the two quantities as two 
separate variables, am I missing some subtlety here?

Thank you in advance for your help,
Best.
-- 
Giovanni.
-------------- next part --------------
gyre [7.2.1]
------------

OpenMP Threads   : 20
Input filename   : gyre_ad.in

Model Init
----------

Reading from OSC file
   File name /home/gio/Work/patrick/test_track/profile180.data.OSC
   File version 2000
   Read 4936 points
   No need to add central point

Mode Search
-----------

Mode parameters
   l : 0
   m : 0

Building frequency grid (REAL axis)
   added scan interval :  0.3000E+01 ->  0.1000E+02 (50 points, LINEAR)

Building spatial grid
   Scaffold grid from model
   Refined 0 subinterval(s) in iteration 1
   Final grid has 1 segment(s) and 4936 point(s):
      Segment 1 : x range 0.0000 -> 1.0000 (1 -> 4936)

Starting search (adiabatic)

Evaluating discriminant
  Time elapsed :     0.741 s

Root Solving
   l    m    n_pg    n_p    n_g       Re(omega)       Im(omega)        chi n_iter
   0    0       1      1      0  0.36981957E+01  0.00000000E+00 0.7788E-14      7
   0    0       2      2      0  0.46454449E+01  0.00000000E+00 0.2989E-14      7
   0    0       3      3      0  0.57064555E+01  0.00000000E+00 0.6508E-15      6
   0    0       4      4      0  0.68023094E+01  0.00000000E+00 0.8751E-14      7
   0    0       5      5      0  0.78820586E+01  0.00000000E+00 0.8386E-14      8
   0    0       6      6      0  0.89145219E+01  0.00000000E+00 0.4528E-13      6
   0    0       7      7      0  0.99621378E+01  0.00000000E+00 0.6488E-14      6
  Time elapsed :      0.043 s

Model lacks the following structure coefficients:
Model lacks the following structure coefficients:
Model lacks the following structure coefficients:

           1           2           3           4           5           8           7           6          10          11          13          14          15          16          19          20
 ABORT at line 148 of /home/gio/Software/gyre-7.2.1/src/model/model_util_m.fypp
 model check failed
ERROR STOP 