Issue with output file

[linovais]

Hi I am currently using MESA to obtain stellar models. To compare my model to the observations I am using GYRE both in non rotating stars and rotating stars. It run fine, but when it is time to extract information from the file it gives back an error.
I am using a program developed by a colegue of mine and also stelarmodels package for python, were the result is the same, it doesn't read. Can someone help me?
Thanks,
Liliana
PS: Error from my colegue program:

gyre_ad [(4.0)]
===============

Compiler : GCC version 4.9.4 20150827 (prerelease)
Compiler options : -I /usr/include/x86_64-linux-gnu -I /home/bnsamba/mesasdk/include -march=broadwell -mmmx -mno-3dnow -msse -msse2 -msse3 -mssse3 -mno-sse4a -mcx16 -msahf -mmovbe -maes -mno-sha -mpclmul -mpopcnt -mabm -mno-lwp -mfma -mno-fma4 -mno-xop -mbmi -mno-tbm -mavx2 -msse4.2 -msse4.1 -mlzcnt -mno-rtm -mno-hle -mrdrnd -mf16c -mfsgsbase -mrdseed -mprfchw -madx -mfxsr -mxsave -mxsaveopt -mno-avx512f -mno-avx512er -mno-avx512cd -mno-avx512pf -mno-prefetchwt1 --param l1-cache-size=32 --param l1-cache-line-size=64 --param l2-cache-size=4096 -mtune=generic -mno-avx -mno-bmi2 -O2 -ffpe-summary=none -finit-real=snan -ffpe-trap=invalid,overflow,zero -fbacktrace -fopenmp
OpenMP Threads : 2
Input filename : HD_10700.data.FGONG.in
GYRE_DIR :

Initialization
==============

Reading from FGONG file HD_10700.data.FGONG
Initial points : 834
File version : 1300
At line 146 of file gyre_fgong_file.f90
Fortran runtime error: Bad value during floating point read


Error from stelar models:

Traceback (most recent call last):
File "fgong_reader.py", line 12, in <module>
starg,starl = fileio.read_fgong('HD_10700_rot.data.FGONG')
File "/home/linovais/anaconda/lib/python2.7/site-packages/stellarmodels/fileio.py", line 385, in read_fgong
' wrong IVERS ({})'.format(IVERS))
ValueError: cannot interpret FGONG file: wrong IVERS (1300)

[rhtownsend]

Hi I am currently using MESA to obtain stellar models. To compare my model to the observations I am using GYRE both in non rotating stars and rotating stars. It run fine, but when it is time to extract information from the file it gives back an error.
I am using a program developed by a colegue of mine and also stelarmodels package for python, were the result is the same, it doesn't read. Can someone help me?
Thanks,
Liliana
PS: Error from my colegue program:

gyre_ad [(4.0)]
===============

Compiler : GCC version 4.9.4 20150827 (prerelease)
Compiler options : -I /usr/include/x86_64-linux-gnu -I /home/bnsamba/mesasdk/include -march=broadwell -mmmx -mno-3dnow -msse -msse2 -msse3 -mssse3 -mno-sse4a -mcx16 -msahf -mmovbe -maes -mno-sha -mpclmul -mpopcnt -mabm -mno-lwp -mfma -mno-fma4 -mno-xop -mbmi -mno-tbm -mavx2 -msse4.2 -msse4.1 -mlzcnt -mno-rtm -mno-hle -mrdrnd -mf16c -mfsgsbase -mrdseed -mprfchw -madx -mfxsr -mxsave -mxsaveopt -mno-avx512f -mno-avx512er -mno-avx512cd -mno-avx512pf -mno-prefetchwt1 --param l1-cache-size=32 --param l1-cache-line-size=64 --param l2-cache-size=4096 -mtune=generic -mno-avx -mno-bmi2 -O2 -ffpe-summary=none -finit-real=snan -ffpe-trap=invalid,overflow,zero -fbacktrace -fopenmp
OpenMP Threads : 2
Input filename : HD_10700.data.FGONG.in
GYRE_DIR :

Initialization
==============

Reading from FGONG file HD_10700.data.FGONG
Initial points : 834
File version : 1300
At line 146 of file gyre_fgong_file.f90
Fortran runtime error: Bad value during floating point read


Error from stelar models:

Traceback (most recent call last):
File "fgong_reader.py", line 12, in <module>
starg,starl = fileio.read_fgong('HD_10700_rot.data.FGONG')
File "/home/linovais/anaconda/lib/python2.7/site-packages/stellarmodels/fileio.py", line 385, in read_fgong
' wrong IVERS ({})'.format(IVERS))
ValueError: cannot interpret FGONG file: wrong IVERS (1300)

Hi Liliana --

As I mentioned in my post to the MESA forum, the problem is caused by MESA writing out an extended-precision FGONG file, which has ivers > 1000. This can only be properly read by GYRE 4.0 if you add the following line:

data_format = '(1P,5(X,E26.18E3))’

to the &model namelist.

Your colleague's python routine will also need to be modified in order to cope with the extended precision format; to do this, they may want to look at the PDF document I attach here, which describes the new format.