search

aarondotter / iso

MESA isochrone code
27 stars 4 forks source link

Making EEP tracks from MESA history file #8

Closed poojanagrawal closed 6 years ago

poojanagrawal commented 6 years ago

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 commented 6 years ago

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 commented 6 years ago

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 commented 6 years ago

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 notifications@github.com 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 https://github.com/aarondotter/iso/issues/8#issuecomment-375144701, or mute the thread https://github.com/notifications/unsubscribe-auth/ADFuL8JemhhDarbjPG0IUXtvH_1hhSsTks5tgviJgaJpZM4SxJ__ .

poojanagrawal commented 6 years ago

I will check it.

Thanks again, Poojan