Making EEP tracks from MESA history file

[poojanagrawal]

Hi,

I am trying to use iso code for converting MESA history files into EEP based track (.eep files) using the code make_eep.

The code runs and produces output; however, it complains that the history files don’t have certain columns e.g. star mass and writes 0.000000 for all entries for those columns in .eep file.

./make_eep input.example
found a tab in line: 	

do not have star_mass
do not have star_mdot
 ..
do not have e_thermal
 history.data          10       10108
       2     180       0       0       0       0       0       0       0       0

I even tried running MESA code with same ‘my_history_columns_full.list’ as is used in iso code, and then using it to make EEP tracks but the problem persists.

./make_eep input.example
  found a tab in line: 	
  found a tab in line: 	
 do not have star_age
 do not have c_core_mass
 …
 do not have max_L_rad_div_Ledd
 history.data          10       10108
       0       0       0       0       0       0       0       0       0       0
  PROBLEM WITH TRACK: NO EEPS DEFINED

Can you advise me on how to resolve this issue?

Thanks,
Poojan

git.zip

[aarondotter]

The format of the history file needs to be specified in the input.format file.

Try changing the 3rd line, which is currently 32, to 40. Remove any .bin files (e.g., history.data.bin) that were previously created and try again.

[poojanagrawal]

Thank you for the help.
Changing the number in 3rd line of input.format to 40 solved the problem.
The code is working fine now.
Can you let me know how did you got this number?

Many Thanks,
Poojan

[aarondotter]

It comes from the format for double precision in MESA. You can find it in mesa/star/defaults/controls.defaults

…

On Wed, 21 Mar 2018 at 9:04 pm poojanagrawal ***@***.***> wrote: Thank you for the help. Changing the number in 3rd line of input.format to 40 solved the problem. The code is working fine now. Can you let me know how did you got this number? Many Thanks, Poojan — You are receiving this because you commented. Reply to this email directly, view it on GitHub <#8 (comment)>, or mute the thread <https://github.com/notifications/unsubscribe-auth/ADFuL8JemhhDarbjPG0IUXtvH_1hhSsTks5tgviJgaJpZM4SxJ__> .

[poojanagrawal]

I will check it.

Thanks again,
Poojan