Lower asymptotic period spacing than expected

General discussion of all things GYRE-related (e.g., results, talks, ideas, tips)
francisco
Posts: 8
Joined: Thu Mar 11, 2021 2:35 pm

Lower asymptotic period spacing than expected

Post by francisco » Tue May 04, 2021 3:49 pm

Hi all,
im new here, so hello !
I have a question that arises in my first steps working with gyre.
Im trying to find the asimptotic period spacing of red clump stars. As a first try, I
developed a model with mesa, with certain parameters and then run Gyre over an specific model.
When I look at the delta_g value computed by GYre its lower than what I would expect in principle, when looking that the period spacing.
Here I show the period spacing obtained from my computations and the what i want to obtain.
I know that im doing something wrong, because it also appears when i compute other models, but dont know what it is.
For your knowledge, I added the following lines in the MESA inlist in order to get a proper resolution of the model.
mesh_delta_coeff = 0.5
add_double_points_to_pulse_data = .true.
Thanks in advance


Francisco
Attachments
deltag.png
deltag.png (40.87 KiB) Viewed 2621 times

User avatar
rhtownsend
Site Admin
Posts: 397
Joined: Sun Mar 31, 2013 4:22 pm

Re: Lower asymptotic period spacing than expected

Post by rhtownsend » Tue May 04, 2021 7:57 pm

Hi Francisco --

Welcome to the forum! As a first suggestion, can you post the namelist input file you use for GYRE?

cheers,

Rich

francisco
Posts: 8
Joined: Thu Mar 11, 2021 2:35 pm

Re: Lower asymptotic period spacing than expected

Post by francisco » Wed May 05, 2021 8:06 am

Hi Rich,
thanks for your reply. Yesterday I posted the input list, dont know why is not here.
Here it is.

Thanks!
Attachments
gyre-template.in
(1.49 KiB) Downloaded 103 times

User avatar
rhtownsend
Site Admin
Posts: 397
Joined: Sun Mar 31, 2013 4:22 pm

Re: Lower asymptotic period spacing than expected

Post by rhtownsend » Wed May 05, 2021 8:26 am

Hi Francisco --

At first, I thought there may be an issue with your input file -- but that looks fine.

On further consideration, I think the result you are finding for the period spaceings looks correct. Your model star shows a period spacing with regular dips -- this is the symptom of mixed modes, where there are two wave-trapping cavities coupled together. Away from the dips, the period spacing tends to be larger than the mean period spacing given by the asymptotic formula. I think this occurs to compensate for the dips; I see it in my own models.

It might be worth calculating the average period spacing for all modes you plot, and see how that compares to the green line you've plotted.

cheers,

Rich

francisco
Posts: 8
Joined: Thu Mar 11, 2021 2:35 pm

Re: Lower asymptotic period spacing than expected

Post by francisco » Wed May 05, 2021 10:33 am

Hi again,
thanks for your feedback Rich.
Well, I computed the mean period spacing :

awk '$1==1 && NR>7 {print 1/(p/1000000)-1/($3/1000000)} {p=$3}' summary.txt | paste -sd+ | bc --> 83003.809

then, i divided this value by

awk '$1==1 && NR>6 {print 1/(p/1000000)-1/($3/1000000)} {p=$3}' summary.txt | wc -l --> 353

which are the number of l=1 modes in the summary.txt file

83003.80/353= 235
that seems more reasonable.

I was trying to reproduce the plot that Ive attached, and there the mean period spacing looks fine, with a mean period spacing of 246 s.
I know that the trapping properties, and the depth of the period spacing could impact the
mean period spacing, but i think there is other issue behind the problem presented.

Thanks again!
Attachments
Captura de pantalla de 2021-05-04 16-48-24.png
Captura de pantalla de 2021-05-04 16-48-24.png (26.88 KiB) Viewed 2604 times

User avatar
rhtownsend
Site Admin
Posts: 397
Joined: Sun Mar 31, 2013 4:22 pm

Re: Lower asymptotic period spacing than expected

Post by rhtownsend » Wed May 05, 2021 11:04 am

OK, thanks for checking that out!

But I still think this behavior is physical -- see, e.g., Miglio et al. (2008, MNRAS, 386, 1487), especially their Fig. 4.

cheers,

Rich

francisco
Posts: 8
Joined: Thu Mar 11, 2021 2:35 pm

Re: Lower asymptotic period spacing than expected

Post by francisco » Thu May 06, 2021 9:39 am

Hi Rich,
thanks for your reply. I will check that paper.

Best

Francisco

francisco
Posts: 8
Joined: Thu Mar 11, 2021 2:35 pm

Re: Lower asymptotic period spacing than expected

Post by francisco » Mon May 31, 2021 6:47 pm

Hi Rich,
I returned to this problem because I find this behaviour in other calculations.
Again I found Deltag values quite lower than expected. So I write the N on the MESA
file and proceed to compute "by hand" the intergral to obtain the value of the asimptotic period spacing with the
Tassoul formula.
Doing this, I found a more reasonable value for the deltag. In the figure attached I show the forward period spacing, the asymptotic "by hand" (blue) and the one from the Gyre (green). So now I suspect that im not reading well the MESA output with Gyre.
I dont think so but, could be a units problem? some imcompatibility with the MESA version (r12778)?
Thanks!
Attachments
profile.data.GYRE
(1.25 MiB) Downloaded 162 times

francisco
Posts: 8
Joined: Thu Mar 11, 2021 2:35 pm

Re: Lower asymptotic period spacing than expected

Post by francisco » Mon May 31, 2021 6:49 pm

Summary: (removed last ten modes)
1
freq_units
UHZ
1 2 3 4 5 6 7
l n_pg Re(freq) Im(freq) E_norm Delta_g Delta_p
1 -381 0.1000408359671362E+002 0.0000000000000000E+000 0.2348883678527399E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -380 0.1003092519544762E+002 0.0000000000000000E+000 0.2298268352902712E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -379 0.1005740053558879E+002 0.0000000000000000E+000 0.2559786460746876E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -378 0.1008489045650239E+002 0.0000000000000000E+000 0.2640566791186644E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -377 0.1011150169169086E+002 0.0000000000000000E+000 0.2698587675447434E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -376 0.1013894258828632E+002 0.0000000000000000E+000 0.2981772769813862E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -375 0.1016646457920399E+002 0.0000000000000000E+000 0.3017325896157959E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -374 0.1019402756933756E+002 0.0000000000000000E+000 0.3158626984144028E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -373 0.1022169052291665E+002 0.0000000000000000E+000 0.3237158952686296E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -372 0.1024943187552410E+002 0.0000000000000000E+000 0.3409841725277896E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -371 0.1027757399372484E+002 0.0000000000000000E+000 0.3553535460897765E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -370 0.1030595776101364E+002 0.0000000000000000E+000 0.3634591927936856E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -369 0.1033415788456931E+002 0.0000000000000000E+000 0.3617971173436948E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -368 0.1036237134586404E+002 0.0000000000000000E+000 0.3937921433340106E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -367 0.1039148349695676E+002 0.0000000000000000E+000 0.3961080598808728E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -366 0.1041991991919860E+002 0.0000000000000000E+000 0.3918543201952280E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -365 0.1044916575799062E+002 0.0000000000000000E+000 0.4213388367162311E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -364 0.1047822649945652E+002 0.0000000000000000E+000 0.4093530639051636E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -363 0.1050732117424367E+002 0.0000000000000000E+000 0.4406765214984947E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -362 0.1053719300732450E+002 0.0000000000000000E+000 0.4257427341281573E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -361 0.1056621803492536E+002 0.0000000000000000E+000 0.4383813777119263E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -360 0.1059662644494302E+002 0.0000000000000000E+000 0.4634849418956510E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -359 0.1062625912918370E+002 0.0000000000000000E+000 0.4288639469206295E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -358 0.1065608113887184E+002 0.0000000000000000E+000 0.4761933723176304E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -357 0.1068700558296665E+002 0.0000000000000000E+000 0.4660955863242302E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -356 0.1071690713736387E+002 0.0000000000000000E+000 0.4502373674676867E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -355 0.1074755687510903E+002 0.0000000000000000E+000 0.4729559981000134E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -354 0.1077841190703607E+002 0.0000000000000000E+000 0.4701123664870696E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -353 0.1080943736728626E+002 0.0000000000000000E+000 0.4805040376932042E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -352 0.1084064870042312E+002 0.0000000000000000E+000 0.4674428783138377E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -351 0.1087149568851835E+002 0.0000000000000000E+000 0.4608374870856435E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -350 0.1090292694537286E+002 0.0000000000000000E+000 0.4884034650989485E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -349 0.1093520532350784E+002 0.0000000000000000E+000 0.4817063892423477E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -348 0.1096668326696963E+002 0.0000000000000000E+000 0.4480490680177212E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -347 0.1099811308281723E+002 0.0000000000000000E+000 0.4796578874318862E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -346 0.1103070023904874E+002 0.0000000000000000E+000 0.4770411097390327E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -345 0.1106281714181796E+002 0.0000000000000000E+000 0.4566845171034120E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -344 0.1109554547035444E+002 0.0000000000000000E+000 0.4647376316701258E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -343 0.1112794032624595E+002 0.0000000000000000E+000 0.4413137763845201E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -342 0.1116020959510912E+002 0.0000000000000000E+000 0.4666505058503896E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -341 0.1119381180288202E+002 0.0000000000000000E+000 0.4576030721474191E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -340 0.1122662318413558E+002 0.0000000000000000E+000 0.4229627356412715E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -339 0.1125964549251721E+002 0.0000000000000000E+000 0.4421445171181710E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -338 0.1129315077450335E+002 0.0000000000000000E+000 0.4339798819612113E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -337 0.1132645670442988E+002 0.0000000000000000E+000 0.4234382987956833E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -336 0.1136026896760155E+002 0.0000000000000000E+000 0.4171221525090208E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -335 0.1139395352823714E+002 0.0000000000000000E+000 0.4025496771806195E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -334 0.1142770777517988E+002 0.0000000000000000E+000 0.4045986512452701E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -333 0.1146201699494306E+002 0.0000000000000000E+000 0.3973779913466477E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -332 0.1149615735710540E+002 0.0000000000000000E+000 0.3776486144562763E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -331 0.1153020547946156E+002 0.0000000000000000E+000 0.3784257696713230E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -330 0.1156497358402700E+002 0.0000000000000000E+000 0.3712976352242135E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -329 0.1159970566865161E+002 0.0000000000000000E+000 0.3462508623679910E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -328 0.1163441601650710E+002 0.0000000000000000E+000 0.3481231746095759E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -327 0.1166973633264109E+002 0.0000000000000000E+000 0.3441501431833636E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -326 0.1170503510742503E+002 0.0000000000000000E+000 0.3193237853671321E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -325 0.1174070434425088E+002 0.0000000000000000E+000 0.3075499160008905E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -324 0.1177680418744174E+002 0.0000000000000000E+000 0.3018600975963829E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -323 0.1181296073811503E+002 0.0000000000000000E+000 0.2979510781721443E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -322 0.1184953995308955E+002 0.0000000000000000E+000 0.2791010607854792E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -321 0.1188620190949230E+002 0.0000000000000000E+000 0.2572354586330788E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -320 0.1192340184138540E+002 0.0000000000000000E+000 0.2586739889824280E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -319 0.1196124418585827E+002 0.0000000000000000E+000 0.2503268372299452E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -318 0.1199859026557684E+002 0.0000000000000000E+000 0.2294916826364999E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -317 0.1203635522084154E+002 0.0000000000000000E+000 0.2204878888768735E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -316 0.1207493657558975E+002 0.0000000000000000E+000 0.2103617604921924E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -315 0.1211361184092655E+002 0.0000000000000000E+000 0.2019436584597769E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -314 0.1215243109941435E+002 0.0000000000000000E+000 0.1894555892680501E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -313 0.1219118852758867E+002 0.0000000000000000E+000 0.1765684920357227E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -312 0.1223068558733532E+002 0.0000000000000000E+000 0.1691413275602704E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -311 0.1227056420557933E+002 0.0000000000000000E+000 0.1576790981003573E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -310 0.1231040003506442E+002 0.0000000000000000E+000 0.1491894795178707E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -309 0.1235068022481015E+002 0.0000000000000000E+000 0.1390960136262996E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -308 0.1239109034735895E+002 0.0000000000000000E+000 0.1273852962860627E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -307 0.1243162343646188E+002 0.0000000000000000E+000 0.1199742936194069E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -306 0.1247246373125414E+002 0.0000000000000000E+000 0.1155494956988225E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -305 0.1251327884304651E+002 0.0000000000000000E+000 0.1070188139363765E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -304 0.1254690886389745E+002 0.0000000000000000E+000 0.1806688346294050E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -303 0.1256336158032102E+002 0.0000000000000000E+000 0.1397722843551913E+000 0.2138764005532469E+003 0.2409948417554577E+000
1 -302 0.1260022014487496E+002 0.0000000000000000E+000 0.9148752766291274E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -301 0.1264194288385454E+002 0.0000000000000000E+000 0.8012278676766445E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -300 0.1268384826239044E+002 0.0000000000000000E+000 0.7161395703923347E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -299 0.1272664640377496E+002 0.0000000000000000E+000 0.6344600828130101E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -298 0.1276949396767963E+002 0.0000000000000000E+000 0.5889909512643771E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -297 0.1281355113917977E+002 0.0000000000000000E+000 0.5407309342309641E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -296 0.1285670072667544E+002 0.0000000000000000E+000 0.4463995055122651E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -295 0.1290000766332233E+002 0.0000000000000000E+000 0.4187194061763468E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -294 0.1294518784120930E+002 0.0000000000000000E+000 0.3645139078526165E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -293 0.1298961299982005E+002 0.0000000000000000E+000 0.3087522509285260E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -292 0.1303441312697371E+002 0.0000000000000000E+000 0.2733120470574320E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -291 0.1307915267516003E+002 0.0000000000000000E+000 0.2343376399820351E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -290 0.1312481881758401E+002 0.0000000000000000E+000 0.1992247799269172E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -289 0.1317065628972489E+002 0.0000000000000000E+000 0.1596442933365440E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -288 0.1321648874591244E+002 0.0000000000000000E+000 0.1380324851489787E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -287 0.1326305006266300E+002 0.0000000000000000E+000 0.1122436372049076E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -286 0.1330906005225124E+002 0.0000000000000000E+000 0.8731297950338765E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -285 0.1335558258619916E+002 0.0000000000000000E+000 0.6893446205770012E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -284 0.1340253143052628E+002 0.0000000000000000E+000 0.5419336205987168E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -283 0.1344983025133265E+002 0.0000000000000000E+000 0.3976977810542197E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -282 0.1349598387413203E+002 0.0000000000000000E+000 0.2905985782358159E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -281 0.1354277699867723E+002 0.0000000000000000E+000 0.2203894291764412E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -280 0.1358902822245654E+002 0.0000000000000000E+000 0.1579013950962454E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -279 0.1363475365562623E+002 0.0000000000000000E+000 0.1232999655089419E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -278 0.1367967669301039E+002 0.0000000000000000E+000 0.1062971310364175E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -277 0.1372555659942436E+002 0.0000000000000000E+000 0.1114078327721241E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -276 0.1377157481382638E+002 0.0000000000000000E+000 0.1175759258660648E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -275 0.1381837000318044E+002 0.0000000000000000E+000 0.1618513742257797E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -274 0.1386780271621265E+002 0.0000000000000000E+000 0.1992417527271528E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -273 0.1391596691681171E+002 0.0000000000000000E+000 0.2549674829875056E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -272 0.1396725412503404E+002 0.0000000000000000E+000 0.3366905882027155E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -271 0.1401785662385065E+002 0.0000000000000000E+000 0.4022219950524587E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -270 0.1406992180080953E+002 0.0000000000000000E+000 0.5208653721096014E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -269 0.1412192851114569E+002 0.0000000000000000E+000 0.5779366916420291E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -268 0.1417376294661537E+002 0.0000000000000000E+000 0.7172414067536612E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -267 0.1422826765525541E+002 0.0000000000000000E+000 0.8205119046883416E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -266 0.1428152576037873E+002 0.0000000000000000E+000 0.8839042598156082E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -265 0.1433547841695120E+002 0.0000000000000000E+000 0.1047540682869508E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -264 0.1439017585245966E+002 0.0000000000000000E+000 0.1158131197276701E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -263 0.1444509083178941E+002 0.0000000000000000E+000 0.1223200424493971E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -262 0.1450047190588178E+002 0.0000000000000000E+000 0.1338544403168038E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -261 0.1455706757567246E+002 0.0000000000000000E+000 0.1455422020438962E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -260 0.1461293826774758E+002 0.0000000000000000E+000 0.1527473231865209E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -259 0.1466979370656163E+002 0.0000000000000000E+000 0.1674272524659936E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -258 0.1472714930740322E+002 0.0000000000000000E+000 0.1632779025107677E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -257 0.1478406639926454E+002 0.0000000000000000E+000 0.1812414895887492E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -256 0.1484386819168850E+002 0.0000000000000000E+000 0.1867313670377456E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -255 0.1490098868386128E+002 0.0000000000000000E+000 0.1788711297591384E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -254 0.1496016187172854E+002 0.0000000000000000E+000 0.2044869666163970E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -253 0.1502045345309348E+002 0.0000000000000000E+000 0.1902324433720853E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -252 0.1507874968573158E+002 0.0000000000000000E+000 0.1975464311373990E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -251 0.1513949799038438E+002 0.0000000000000000E+000 0.2054613764328833E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -250 0.1519828263462466E+002 0.0000000000000000E+000 0.2002137637278007E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -249 0.1525621103346833E+002 0.0000000000000000E+000 0.2226924926544749E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -248 0.1530857213715710E+002 0.0000000000000000E+000 0.2482857010313358E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -247 0.1535302217953981E+002 0.0000000000000000E+000 0.2507158058566036E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -246 0.1540332297382962E+002 0.0000000000000000E+000 0.2187150031481648E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -245 0.1546181519699536E+002 0.0000000000000000E+000 0.1978806294435105E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -244 0.1552348077237442E+002 0.0000000000000000E+000 0.1800258122173560E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -243 0.1558625199372635E+002 0.0000000000000000E+000 0.1717090741098889E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -242 0.1565027785803439E+002 0.0000000000000000E+000 0.1733714769568442E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -241 0.1571527256262888E+002 0.0000000000000000E+000 0.1579583952512335E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -240 0.1577982071946269E+002 0.0000000000000000E+000 0.1504596059025553E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -239 0.1584735823938166E+002 0.0000000000000000E+000 0.1427644128822107E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -238 0.1591318844063728E+002 0.0000000000000000E+000 0.1281945602838702E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -237 0.1598031778215129E+002 0.0000000000000000E+000 0.1319331029659555E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -236 0.1604901264100746E+002 0.0000000000000000E+000 0.1103062976057693E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -235 0.1611584204925517E+002 0.0000000000000000E+000 0.1017696313225062E-001 0.2138764005532469E+003 0.2409948417554577E+000
1 -234 0.1618650208675939E+002 0.0000000000000000E+000 0.9735017309563939E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -233 0.1625554837540852E+002 0.0000000000000000E+000 0.8388857361406050E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -232 0.1632559817491050E+002 0.0000000000000000E+000 0.7684626715272102E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -231 0.1639618769307610E+002 0.0000000000000000E+000 0.6385354000441633E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -230 0.1646709285925568E+002 0.0000000000000000E+000 0.5943432888326005E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -229 0.1654027446915771E+002 0.0000000000000000E+000 0.5026292095057887E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -228 0.1661135299326800E+002 0.0000000000000000E+000 0.4044159925436232E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -227 0.1668377423127606E+002 0.0000000000000000E+000 0.3483239844422429E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -226 0.1675762501813955E+002 0.0000000000000000E+000 0.2909545920949543E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -225 0.1683207909694682E+002 0.0000000000000000E+000 0.2214700296091440E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -224 0.1690429225890382E+002 0.0000000000000000E+000 0.1677772092354442E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -223 0.1697913119916749E+002 0.0000000000000000E+000 0.1354609795232065E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -222 0.1705461944247425E+002 0.0000000000000000E+000 0.8943507282983964E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -221 0.1712771451212824E+002 0.0000000000000000E+000 0.6391796382111414E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -220 0.1720179479751481E+002 0.0000000000000000E+000 0.4116676516200017E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -219 0.1727267469981545E+002 0.0000000000000000E+000 0.2586323095403868E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -218 0.1734405858017071E+002 0.0000000000000000E+000 0.1640508920687144E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -217 0.1740946411285692E+002 0.0000000000000000E+000 0.1234964863158923E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -216 0.1747556924856177E+002 0.0000000000000000E+000 0.1371818822749238E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -215 0.1754132597705837E+002 0.0000000000000000E+000 0.1707986417363937E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -214 0.1761115069930905E+002 0.0000000000000000E+000 0.2773438521104626E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -213 0.1768358960158892E+002 0.0000000000000000E+000 0.4032858433736992E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -212 0.1775550959106712E+002 0.0000000000000000E+000 0.5684925032935628E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -211 0.1783126310106250E+002 0.0000000000000000E+000 0.7488281645379737E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -210 0.1790900565664254E+002 0.0000000000000000E+000 0.9830226558728778E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -209 0.1798905374265921E+002 0.0000000000000000E+000 0.1247458629923028E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -208 0.1806822408215152E+002 0.0000000000000000E+000 0.1428745719199389E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -207 0.1815040907594439E+002 0.0000000000000000E+000 0.1768991277045846E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -206 0.1823427092437350E+002 0.0000000000000000E+000 0.1980183800649891E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -205 0.1831608844159482E+002 0.0000000000000000E+000 0.2318970811358614E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -204 0.1840084835323940E+002 0.0000000000000000E+000 0.2516698799430657E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -203 0.1848482397820783E+002 0.0000000000000000E+000 0.2644227993098377E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -202 0.1857128347614927E+002 0.0000000000000000E+000 0.3092493339220353E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -201 0.1865779530409689E+002 0.0000000000000000E+000 0.3022333079151656E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -200 0.1874415010353835E+002 0.0000000000000000E+000 0.3207483674108338E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -199 0.1883585236700531E+002 0.0000000000000000E+000 0.3331998537726194E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -198 0.1892507795185491E+002 0.0000000000000000E+000 0.3269047602926676E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -197 0.1901686803081001E+002 0.0000000000000000E+000 0.3427053671695211E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -196 0.1911095330929398E+002 0.0000000000000000E+000 0.3189609350013398E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -195 0.1920626354398723E+002 0.0000000000000000E+000 0.3253082037365268E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -194 0.1930313631594547E+002 0.0000000000000000E+000 0.3052070056951632E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -193 0.1939837146258252E+002 0.0000000000000000E+000 0.2962627440877808E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -192 0.1949922702323961E+002 0.0000000000000000E+000 0.2774740299699339E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -191 0.1959898704407686E+002 0.0000000000000000E+000 0.2454818259419666E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -190 0.1970050629870297E+002 0.0000000000000000E+000 0.2374445299460141E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -189 0.1980182407978575E+002 0.0000000000000000E+000 0.2049124291394741E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -188 0.1990534171745891E+002 0.0000000000000000E+000 0.1853173401112040E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -187 0.2001185119018169E+002 0.0000000000000000E+000 0.1476219823865288E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -186 0.2011585911682414E+002 0.0000000000000000E+000 0.1320292164416708E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -185 0.2022437600149639E+002 0.0000000000000000E+000 0.1084161755697578E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -184 0.2032983802045989E+002 0.0000000000000000E+000 0.7993495931480280E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -183 0.2044005545160609E+002 0.0000000000000000E+000 0.6310634418515095E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -182 0.2054988469559698E+002 0.0000000000000000E+000 0.4336612888368718E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -181 0.2065900593688694E+002 0.0000000000000000E+000 0.3093892512053741E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -180 0.2076755649480905E+002 0.0000000000000000E+000 0.1841885846714147E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -179 0.2087459824787738E+002 0.0000000000000000E+000 0.1145550265694981E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -178 0.2097981973292626E+002 0.0000000000000000E+000 0.6645596785554015E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -177 0.2107843883970493E+002 0.0000000000000000E+000 0.5744398033717620E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -176 0.2117907714266658E+002 0.0000000000000000E+000 0.6591707274568018E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -175 0.2128421179692051E+002 0.0000000000000000E+000 0.1078554381643786E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -174 0.2139873868908994E+002 0.0000000000000000E+000 0.1676043441876241E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -173 0.2151643990959531E+002 0.0000000000000000E+000 0.2656640538714559E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -172 0.2163791219668865E+002 0.0000000000000000E+000 0.3616906843650331E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -171 0.2175745626281292E+002 0.0000000000000000E+000 0.5169139669206885E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -170 0.2188621303767576E+002 0.0000000000000000E+000 0.6394569739574997E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -169 0.2201008099600157E+002 0.0000000000000000E+000 0.7200375882944964E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -168 0.2214110203212313E+002 0.0000000000000000E+000 0.9442260305218821E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -167 0.2227073816588644E+002 0.0000000000000000E+000 0.9895955022336917E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -166 0.2240146833939516E+002 0.0000000000000000E+000 0.1142235298503457E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -165 0.2253648062115329E+002 0.0000000000000000E+000 0.1127510305345983E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -164 0.2267040549404511E+002 0.0000000000000000E+000 0.1266127273279326E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -163 0.2281069435657684E+002 0.0000000000000000E+000 0.1249789109578165E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -162 0.2294483510199273E+002 0.0000000000000000E+000 0.1218297595351297E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -161 0.2308814143391769E+002 0.0000000000000000E+000 0.1205389468846965E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -160 0.2322833070105719E+002 0.0000000000000000E+000 0.1076260929800628E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -159 0.2337455175866149E+002 0.0000000000000000E+000 0.1071558444752645E-002 0.2138764005532469E+003 0.2409948417554577E+000
1 -158 0.2351835872839757E+002 0.0000000000000000E+000 0.8458677006232659E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -157 0.2366342322588315E+002 0.0000000000000000E+000 0.8172292444667980E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -156 0.2381725399996874E+002 0.0000000000000000E+000 0.5965095015202089E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -155 0.2396250454163307E+002 0.0000000000000000E+000 0.4946631799330436E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -154 0.2411998168562292E+002 0.0000000000000000E+000 0.3685648003234194E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -153 0.2426724964177781E+002 0.0000000000000000E+000 0.2364501594969726E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -152 0.2442468799729902E+002 0.0000000000000000E+000 0.1567601163796756E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -151 0.2457463393877306E+002 0.0000000000000000E+000 0.7784346616549883E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -150 0.2472341947095929E+002 0.0000000000000000E+000 0.4272833760584302E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -149 0.2485411717873309E+002 0.0000000000000000E+000 0.2388220816062405E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -148 0.2498821500564327E+002 0.0000000000000000E+000 0.2983889708622242E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -147 0.2512618136344866E+002 0.0000000000000000E+000 0.5012893041582910E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -146 0.2529094237639740E+002 0.0000000000000000E+000 0.1082083712689297E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -145 0.2545035823924795E+002 0.0000000000000000E+000 0.1568976708792406E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -144 0.2562226392852583E+002 0.0000000000000000E+000 0.2826465071490335E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -143 0.2579504561685816E+002 0.0000000000000000E+000 0.3151110123229484E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -142 0.2596780702028743E+002 0.0000000000000000E+000 0.4737852105990698E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -141 0.2615523220094523E+002 0.0000000000000000E+000 0.4916187702750754E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -140 0.2632632238906169E+002 0.0000000000000000E+000 0.5953130816805695E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -139 0.2652132044461053E+002 0.0000000000000000E+000 0.6611878778602524E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -138 0.2669912813725995E+002 0.0000000000000000E+000 0.5999860111982030E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -137 0.2689545729789593E+002 0.0000000000000000E+000 0.6979354143672368E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -136 0.2708534447751221E+002 0.0000000000000000E+000 0.5485613950416665E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -135 0.2727931124742140E+002 0.0000000000000000E+000 0.5941026268419671E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -134 0.2748155774822842E+002 0.0000000000000000E+000 0.4093905927339977E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -133 0.2767269390906219E+002 0.0000000000000000E+000 0.3898085600764998E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -132 0.2788780749986774E+002 0.0000000000000000E+000 0.2539737075024220E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -131 0.2807402292278838E+002 0.0000000000000000E+000 0.1749601964501873E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -130 0.2827533778449683E+002 0.0000000000000000E+000 0.1592607892292288E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -129 0.2834246810439716E+002 0.0000000000000000E+000 0.2113694165524780E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -128 0.2848862699624810E+002 0.0000000000000000E+000 0.4736578624655340E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -127 0.2868428916556834E+002 0.0000000000000000E+000 0.2564075565606301E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -126 0.2883848536230882E+002 0.0000000000000000E+000 0.2115402127765785E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -125 0.2903575920946887E+002 0.0000000000000000E+000 0.3950800174559372E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -124 0.2922164588155148E+002 0.0000000000000000E+000 0.7400170667082327E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -123 0.2946008744988242E+002 0.0000000000000000E+000 0.1522164312034884E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -122 0.2967111196605526E+002 0.0000000000000000E+000 0.1914545576648665E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -121 0.2991379486093227E+002 0.0000000000000000E+000 0.3471781357837647E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -120 0.3014380239811487E+002 0.0000000000000000E+000 0.3213048672034431E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -119 0.3038136390249871E+002 0.0000000000000000E+000 0.4899598951075211E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -118 0.3063867299696982E+002 0.0000000000000000E+000 0.3961926217690662E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -117 0.3087185825991489E+002 0.0000000000000000E+000 0.4936459011770886E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -116 0.3114752466167586E+002 0.0000000000000000E+000 0.4108210510277197E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -115 0.3137735528044857E+002 0.0000000000000000E+000 0.3396055174528571E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -114 0.3166098912517022E+002 0.0000000000000000E+000 0.3129277868257499E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -113 0.3191079490226252E+002 0.0000000000000000E+000 0.1658480891035238E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -112 0.3218450736440232E+002 0.0000000000000000E+000 0.1389122909809094E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -111 0.3244655956125402E+002 0.0000000000000000E+000 0.4722945611562252E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -110 0.3268981508470291E+002 0.0000000000000000E+000 0.2727993232709299E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -109 0.3290182718684483E+002 0.0000000000000000E+000 0.1269187487197612E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -108 0.3313577902184445E+002 0.0000000000000000E+000 0.2883963076654751E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -107 0.3337690087450436E+002 0.0000000000000000E+000 0.4936229250972489E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -106 0.3367805883594815E+002 0.0000000000000000E+000 0.1386893220466034E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -105 0.3396448975241481E+002 0.0000000000000000E+000 0.1428802197164840E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -104 0.3426579158270432E+002 0.0000000000000000E+000 0.2917539642276087E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -103 0.3459301403344207E+002 0.0000000000000000E+000 0.2346095483246814E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -102 0.3487428720388488E+002 0.0000000000000000E+000 0.3135271727027546E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -101 0.3523155165498942E+002 0.0000000000000000E+000 0.2564341547044310E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -100 0.3552097175037599E+002 0.0000000000000000E+000 0.1962235898964196E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -99 0.3588730101088341E+002 0.0000000000000000E+000 0.1641664196608908E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -98 0.3618763813472739E+002 0.0000000000000000E+000 0.7061552131549706E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -97 0.3652135614162039E+002 0.0000000000000000E+000 0.4328004490244962E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -96 0.3679496024874178E+002 0.0000000000000000E+000 0.1242097090589950E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -95 0.3707211506703560E+002 0.0000000000000000E+000 0.1647411125445890E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -94 0.3733419011935639E+002 0.0000000000000000E+000 0.3316079433757980E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -93 0.3771377390348491E+002 0.0000000000000000E+000 0.1043043042273232E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -92 0.3802982956840764E+002 0.0000000000000000E+000 0.1182974926157574E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -91 0.3841268930798253E+002 0.0000000000000000E+000 0.2574021024129380E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -90 0.3879724551344389E+002 0.0000000000000000E+000 0.1867396016726409E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -89 0.3914219900748182E+002 0.0000000000000000E+000 0.2896472560032360E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -88 0.3955666963093740E+002 0.0000000000000000E+000 0.1934332416116209E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -87 0.3988948093625195E+002 0.0000000000000000E+000 0.1293873366464810E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -86 0.4032725457956968E+002 0.0000000000000000E+000 0.8937604749235290E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -85 0.4067126449385742E+002 0.0000000000000000E+000 0.2564274978241059E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -84 0.4101751754805992E+002 0.0000000000000000E+000 0.1428001719634159E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -83 0.4127419500846604E+002 0.0000000000000000E+000 0.1380491987560599E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -82 0.4170672500371013E+002 0.0000000000000000E+000 0.5088481295800327E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -81 0.4207771437213027E+002 0.0000000000000000E+000 0.8095468439697513E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -80 0.4257111417881534E+002 0.0000000000000000E+000 0.1930865289325591E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -79 0.4301167465968543E+002 0.0000000000000000E+000 0.1380578881175937E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -78 0.4349588960599065E+002 0.0000000000000000E+000 0.2139796644164670E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -77 0.4400815261757698E+002 0.0000000000000000E+000 0.9317695315407770E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -76 0.4444389281701847E+002 0.0000000000000000E+000 0.7563948252888947E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -75 0.4495978650384745E+002 0.0000000000000000E+000 0.1556531003007308E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -74 0.4529449160294021E+002 0.0000000000000000E+000 0.1182430573131409E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -73 0.4568159474346020E+002 0.0000000000000000E+000 0.1835075573852385E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -72 0.4621548363123025E+002 0.0000000000000000E+000 0.9195123298095074E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -71 0.4674381047627776E+002 0.0000000000000000E+000 0.9163424228366840E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -70 0.4728370319562774E+002 0.0000000000000000E+000 0.1817955020493219E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -69 0.4793445055722371E+002 0.0000000000000000E+000 0.9171329815030431E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -68 0.4846920676074593E+002 0.0000000000000000E+000 0.7672902269321629E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -67 0.4909895981576726E+002 0.0000000000000000E+000 0.1647875213119624E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -66 0.4947328781923466E+002 0.0000000000000000E+000 0.1119821962478384E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -65 0.4990813421574985E+002 0.0000000000000000E+000 0.1690965191169767E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -64 0.5053645292107890E+002 0.0000000000000000E+000 0.8492918132590230E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -63 0.5120203534677099E+002 0.0000000000000000E+000 0.8119709950349225E-004 0.2138764005532469E+003 0.2409948417554577E+000
1 -62 0.5188050649078963E+002 0.0000000000000000E+000 0.1486138094592853E-003 0.2138764005532469E+003 0.2409948417554577E+000
1 -61 0.5260098989302991E+002 0.0000000000000000E+000 0.5214406729267455E-004 0.2138764005532469E+003 0.2409948417554577E+000
Attachments
deltap.pdf
(187.52 KiB) Downloaded 102 times
profile.txt
(5.63 MiB) Downloaded 128 times

User avatar
rhtownsend
Site Admin
Posts: 397
Joined: Sun Mar 31, 2013 4:22 pm

Re: Lower asymptotic period spacing than expected

Post by rhtownsend » Mon May 31, 2021 7:34 pm

francisco wrote:
Mon May 31, 2021 6:47 pm
Hi Rich,
I returned to this problem because I find this behaviour in other calculations.
Again I found Deltag values quite lower than expected. So I write the N on the MESA
file and proceed to compute "by hand" the intergral to obtain the value of the asimptotic period spacing with the
Tassoul formula.
Doing this, I found a more reasonable value for the deltag. In the figure attached I show the forward period spacing, the asymptotic "by hand" (blue) and the one from the Gyre (green). So now I suspect that im not reading well the MESA output with Gyre.
I dont think so but, could be a units problem? some imcompatibility with the MESA version (r12778)?
Thanks!
Hi Francisco --

From your summary file, it looks like you're interpreting the Delta_g output item as the g-mode period separation, in seconds. But that's not correct; Delta_g is the *inverse* of the period separation, expressed in dimensionless units. To calculate the period separation, you need to use the following formula:

Delta(P) = 1/Delta_g * sqrt(R^3/GM)

cheers,

Rich

Post Reply