[Mesa-users] A typo in atm module?

[Yixiao Zhou]

Dear MESA developers and users,


There might be a small problem in the grey atmosphere routine "eval_grey" in mod_atm.f90.


The derivative used in "eval_grey", d_logg_dlnR should equals to -2/ln10 if assuming surface gravity g is propotional to GM/R^2. Similraly, d_logg_dlnM = 1/ln10.


This is indeed the case in other part of atmosphere module (e.g., routine "do_atm_get" in mod_atm.f90). However, in "eval_grey" they are set to be

d_logg_dlnR = -2

d_logg_dlnM = 1


Is this a typo or I misunderstood something?


Thanks,

Yixiao Zhou

PhD student

RSAA

ANU
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20181024/1d2a34be/attachment.html>