troubles of oscillation calculations for fast rotators

[wutao]

Dear Rich and guys,
I'm Tao, concentrating on the research of SPB asteroseismology and trying to calculate the theoretical modes with GYRE code for a fast rotating SPB star. For the modes of (l,m) = (1,-1) and (1,0), they can be regularly calculated, while, for the modes of (l,m) = (1,+1), the code does not work. Would you give some bits of advice for this? Thanks in advance.

The MESA file and the input file of GYRE ('gyre.in') are attached. The screen output is as follows.

Thanks again!
cheers.
Tao

screen output:


gyre [master]
-------------

OpenMP Threads : 30
Input filename : ../../gyre.in_master-c0.log

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

Reading from MESA file
File name profile636.data.GYRE
File version 1.01
Read 2181 points
No need to add central point

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

Mode parameters
l : 1
m : 0

Building frequency grid (REAL axis)
added scan interval : 0.1233E+00 -> 0.2589E+01 (848 points, INVERSE)

Building spatial grid
Refined 362 subinterval(s) in iteration 1
Refined 210 subinterval(s) in iteration 2
Refined 87 subinterval(s) in iteration 3
Refined 16 subinterval(s) in iteration 4
Refined 0 subinterval(s) in iteration 5
Final grid has 1 segment(s) and 2856 point(s):
Segment 1 : x range 0.0000 -> 1.0051 (1 -> 2856)

Starting search (adiabatic)

Evaluating discriminant
Time elapsed : 3.784 s

Root Solving
l m n_pg n_p n_g Re(omega) Im(omega) chi n_iter
1 0 -60 0 60 0.22813292E+00 0.00000000E+00 0.6637E-13 6
1 0 -59 0 59 0.23022875E+00 0.00000000E+00 0.3411E-13 6
1 0 -58 0 58 0.23225546E+00 0.00000000E+00 0.6026E-13 6
1 0 -57 0 57 0.23451345E+00 0.00000000E+00 0.3830E-13 6
1 0 -56 0 56 0.23663492E+00 0.00000000E+00 0.9930E-13 6
1 0 -55 0 55 0.23901986E+00 0.00000000E+00 0.2157E-13 6
1 0 -54 0 54 0.24131713E+00 0.00000000E+00 0.4790E-13 6
1 0 -53 0 53 0.24380399E+00 0.00000000E+00 0.1530E-12 6
1 0 -52 0 52 0.24632654E+00 0.00000000E+00 0.1059E-12 6
1 0 -51 0 51 0.24888319E+00 0.00000000E+00 0.7883E-13 6
1 0 -50 0 50 0.25167016E+00 0.00000000E+00 0.1924E-12 6
1 0 -49 0 49 0.25430286E+00 0.00000000E+00 0.4592E-13 6
1 0 -48 0 48 0.25736765E+00 0.00000000E+00 0.5687E-13 7
1 0 -47 0 47 0.26010408E+00 0.00000000E+00 0.8283E-13 6
1 0 -46 0 46 0.26339658E+00 0.00000000E+00 0.9785E-13 6
1 0 -45 0 45 0.26637712E+00 0.00000000E+00 0.5207E-13 6
1 0 -44 0 44 0.26986943E+00 0.00000000E+00 0.3092E-13 6
1 0 -43 0 43 0.27320734E+00 0.00000000E+00 0.7150E-13 5
1 0 -42 0 42 0.27681405E+00 0.00000000E+00 0.9061E-14 6
1 0 -41 0 41 0.28062849E+00 0.00000000E+00 0.1504E-13 6
1 0 -40 0 40 0.28438593E+00 0.00000000E+00 0.8598E-14 6
1 0 -39 0 39 0.28874905E+00 0.00000000E+00 0.1335E-12 6
1 0 -38 0 38 0.29262992E+00 0.00000000E+00 0.3755E-13 6
1 0 -37 0 37 0.29753958E+00 0.00000000E+00 0.2203E-13 6
1 0 -36 0 36 0.30178291E+00 0.00000000E+00 0.8516E-13 5
1 0 -35 0 35 0.30728836E+00 0.00000000E+00 0.6452E-13 7
1 0 -34 0 34 0.31204065E+00 0.00000000E+00 0.6112E-13 6
1 0 -33 0 33 0.31795522E+00 0.00000000E+00 0.1215E-12 6
1 0 -32 0 32 0.32352116E+00 0.00000000E+00 0.1322E-12 5
1 0 -31 0 31 0.32991976E+00 0.00000000E+00 0.1523E-12 5
1 0 -30 0 30 0.33644218E+00 0.00000000E+00 0.1186E-12 7
1 0 -29 0 29 0.34311672E+00 0.00000000E+00 0.7121E-13 6
1 0 -28 0 28 0.35089886E+00 0.00000000E+00 0.3945E-13 6
1 0 -27 0 27 0.35824325E+00 0.00000000E+00 0.5873E-13 7
1 0 -26 0 26 0.36755942E+00 0.00000000E+00 0.1787E-13 7
1 0 -25 0 25 0.37549876E+00 0.00000000E+00 0.1102E-12 7
1 0 -24 0 24 0.38656923E+00 0.00000000E+00 0.2075E-13 6
1 0 -23 0 23 0.39597300E+00 0.00000000E+00 0.1139E-12 7
1 0 -22 0 22 0.40914175E+00 0.00000000E+00 0.2009E-13 7
1 0 -21 0 21 0.42035285E+00 0.00000000E+00 0.9758E-13 7
1 0 -20 0 20 0.43672386E+00 0.00000000E+00 0.1849E-12 6
1 0 -19 0 19 0.45122932E+00 0.00000000E+00 0.4059E-13 6
1 0 -18 0 18 0.47039216E+00 0.00000000E+00 0.2238E-13 7
1 0 -17 0 17 0.48822096E+00 0.00000000E+00 0.5377E-14 6
1 0 -16 0 16 0.51213095E+00 0.00000000E+00 0.1819E-12 8
1 0 -15 0 15 0.53741609E+00 0.00000000E+00 0.5022E-14 6
1 0 -14 0 14 0.56743062E+00 0.00000000E+00 0.7775E-13 7
1 0 -13 0 13 0.60394679E+00 0.00000000E+00 0.7461E-13 6
1 0 -12 0 12 0.64090467E+00 0.00000000E+00 0.1750E-12 8
1 0 -11 0 11 0.69913932E+00 0.00000000E+00 0.1407E-12 6
1 0 -10 0 10 0.74631554E+00 0.00000000E+00 0.3328E-13 8
1 0 -9 0 9 0.83662623E+00 0.00000000E+00 0.2387E-13 8
1 0 -8 0 8 0.90131944E+00 0.00000000E+00 0.1331E-12 7
1 0 -7 0 7 0.10457359E+01 0.00000000E+00 0.4086E-13 7
1 0 -6 0 6 0.11587942E+01 0.00000000E+00 0.8458E-14 7
1 0 -5 0 5 0.13675556E+01 0.00000000E+00 0.3228E-12 11
1 0 -5 0 5 0.14313213E+01 0.00000000E+00 0.1621E-12 11
Time elapsed : 6.249 s

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

Mode parameters
l : 1
m : -1

Building frequency grid (REAL axis)
added scan interval : 0.8220E-01 -> 0.1726E+01 (848 points, INVERSE)

Building spatial grid
Refined 42 subinterval(s) in iteration 1
Refined 10 subinterval(s) in iteration 2
Refined 0 subinterval(s) in iteration 3
Final grid has 1 segment(s) and 2233 point(s):
Segment 1 : x range 0.0000 -> 1.0051 (1 -> 2233)

Starting search (adiabatic)

Evaluating discriminant
Time elapsed : 3.288 s

Root Solving
l m n_pg n_p n_g Re(omega) Im(omega) chi n_iter
1 -1 -60 0 60 0.92214651E-01 0.00000000E+00 0.3015E-12 6
1 -1 -59 0 59 0.95031456E-01 0.00000000E+00 0.3821E-14 7
1 -1 -58 0 58 0.97725078E-01 0.00000000E+00 0.9343E-12 6
1 -1 -57 0 57 0.10071088E+00 0.00000000E+00 0.2378E-12 7
1 -1 -56 0 56 0.10350327E+00 0.00000000E+00 0.1963E-12 9
1 -1 -55 0 55 0.10659929E+00 0.00000000E+00 0.4848E-12 7
1 -1 -54 0 54 0.10958853E+00 0.00000000E+00 0.2125E-12 5
1 -1 -53 0 53 0.11275562E+00 0.00000000E+00 0.2034E-12 7
1 -1 -52 0 52 0.11598688E+00 0.00000000E+00 0.1203E-12 5
1 -1 -51 0 51 0.11918280E+00 0.00000000E+00 0.2521E-13 5
1 -1 -50 0 50 0.12268083E+00 0.00000000E+00 0.7769E-13 7
1 -1 -49 0 49 0.12591477E+00 0.00000000E+00 0.6287E-13 6
1 -1 -48 0 48 0.12966951E+00 0.00000000E+00 0.2393E-12 8
1 -1 -47 0 47 0.13298144E+00 0.00000000E+00 0.3091E-13 6
1 -1 -46 0 46 0.13690964E+00 0.00000000E+00 0.2274E-12 6
1 -1 -45 0 45 0.14045928E+00 0.00000000E+00 0.1012E-13 6
1 -1 -44 0 44 0.14451453E+00 0.00000000E+00 0.3656E-12 6
1 -1 -43 0 43 0.14840992E+00 0.00000000E+00 0.1717E-12 5
1 -1 -42 0 42 0.15249021E+00 0.00000000E+00 0.6040E-13 7
1 -1 -41 0 41 0.15683316E+00 0.00000000E+00 0.6995E-13 6
1 -1 -40 0 40 0.16098359E+00 0.00000000E+00 0.7904E-13 6
1 -1 -39 0 39 0.16581487E+00 0.00000000E+00 0.2505E-13 6
1 -1 -38 0 38 0.17001387E+00 0.00000000E+00 0.6294E-13 6
1 -1 -37 0 37 0.17530064E+00 0.00000000E+00 0.1416E-12 6
1 -1 -36 0 36 0.17980988E+00 0.00000000E+00 0.8570E-14 7
1 -1 -35 0 35 0.18557451E+00 0.00000000E+00 0.8712E-13 6
1 -1 -34 0 34 0.19052970E+00 0.00000000E+00 0.1603E-12 7
1 -1 -33 0 33 0.19657210E+00 0.00000000E+00 0.8814E-13 6
1 -1 -32 0 32 0.20226706E+00 0.00000000E+00 0.2892E-13 5
1 -1 -31 0 31 0.20866155E+00 0.00000000E+00 0.1306E-12 7
1 -1 -30 0 30 0.21520749E+00 0.00000000E+00 0.1218E-14 7
1 -1 -29 0 29 0.22177260E+00 0.00000000E+00 0.1346E-12 6
1 -1 -28 0 28 0.22946517E+00 0.00000000E+00 0.1261E-13 7
1 -1 -27 0 27 0.23661094E+00 0.00000000E+00 0.3487E-13 7
1 -1 -26 0 26 0.24568646E+00 0.00000000E+00 0.1402E-12 7
1 -1 -25 0 25 0.25336890E+00 0.00000000E+00 0.8218E-13 6
1 -1 -24 0 24 0.26410317E+00 0.00000000E+00 0.3378E-13 6
1 -1 -23 0 23 0.27318320E+00 0.00000000E+00 0.5011E-13 6
1 -1 -22 0 22 0.28590360E+00 0.00000000E+00 0.1539E-12 6
1 -1 -21 0 21 0.29673085E+00 0.00000000E+00 0.1631E-13 6
1 -1 -20 0 20 0.31246583E+00 0.00000000E+00 0.1610E-13 6
1 -1 -19 0 19 0.32652316E+00 0.00000000E+00 0.2873E-14 6
1 -1 -18 0 18 0.34487795E+00 0.00000000E+00 0.7591E-13 7
1 -1 -17 0 17 0.36227412E+00 0.00000000E+00 0.1245E-12 8
1 -1 -16 0 16 0.38518773E+00 0.00000000E+00 0.6000E-13 6
1 -1 -15 0 15 0.41004188E+00 0.00000000E+00 0.1261E-13 8
1 -1 -14 0 14 0.43887495E+00 0.00000000E+00 0.1583E-14 7
1 -1 -13 0 13 0.47496557E+00 0.00000000E+00 0.1888E-13 8
1 -1 -12 0 12 0.51060500E+00 0.00000000E+00 0.1684E-13 6
1 -1 -11 0 11 0.56823079E+00 0.00000000E+00 0.1298E-12 8
1 -1 -10 0 10 0.61418365E+00 0.00000000E+00 0.7195E-13 6
1 -1 -9 0 9 0.70356133E+00 0.00000000E+00 0.8211E-13 7
1 -1 -8 0 8 0.76730487E+00 0.00000000E+00 0.1266E-12 7
1 -1 -7 0 7 0.91024555E+00 0.00000000E+00 0.4517E-13 6
1 -1 -6 0 6 0.10231281E+01 0.00000000E+00 0.1350E-12 7
1 -1 -5 0 5 0.12633996E+01 0.00000000E+00 0.2716E-13 8
Time elapsed : 1.719 s

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

Mode parameters
l : 1
m : 1

Building frequency grid (REAL axis)
added scan interval : 0.2466E+00 -> 0.5178E+01 (848 points, INVERSE)

Building spatial grid
Refined 2180 subinterval(s) in iteration 1
Refined 4360 subinterval(s) in iteration 2
Refined 8720 subinterval(s) in iteration 3
Refined 17440 subinterval(s) in iteration 4
Refined 34880 subinterval(s) in iteration 5
Refined 69758 subinterval(s) in iteration 6
Refined 139475 subinterval(s) in iteration 7
Refined 278571 subinterval(s) in iteration 8
Refined 553824 subinterval(s) in iteration 9
Refined 1080563 subinterval(s) in iteration 10
Refined 2025533 subinterval(s) in iteration 11
Refined 3559384 subinterval(s) in iteration 12
Refined 6420509 subinterval(s) in iteration 13
Refined 9554998 subinterval(s) in iteration 14
Refined 15900037 subinterval(s) in iteration 15
Refined 25859650 subinterval(s) in iteration 16
Refined 25886720 subinterval(s) in iteration 17
Refined 0 subinterval(s) in iteration 18
Final grid has 1 segment(s) and 91398783 point(s):
Segment 1 : x range 0.0000 -> 1.0051 (1 -> 91398783)

[rhtownsend]

Hi Tao --

It looks as though the re-gridding algorithm is running into trouble for the l=1, m=1 modes -- GYRE is trying to create a grid with 91,398,783 points, which likely triggers the code to crash from lack of memory.

My guess is that the problem arises because of your chosen frequency range:

Code: Select all

&scan

 grid_type = 'INVERSE'
 freq_min =    4.00945351658092     
 freq_max =    84.1985238481993     
 freq_units = 'UHZ'
 n_freq =          848
 tag_list = 'l12'
 
/

These parameters specify the range of frequencies to scan in the inertial reference frame. When GYRE transforms to the co-rotating reference frame, I expect that the lower frequency bound (freq_min) becomes a very small number --- corresponding to very high-order g-modes. GYRE attempts to resolve the spatial oscillation of these modes by creating a extremely fine grid.

The solution is to specify the frequency range in the co-rotating frame of reference. To do this, add the following parameters to the &scan namelist group:

Code: Select all

grid_frame = 'COROT_I'
freq_frame = 'COROT_I'

Best wishes,

Rich