[mesa-users] other_momentum subroutine

[sheyda najafi]

Hi MESA users

I am trying to add an extra term (- dx/dm) to the hydrostatic equilibrium
equation as follows:

dp/dm = - G m/4 pi r^4 - dx/dm

 or

(( s%  p(k+1)-    s% p(k))/s% dm_bar(k)) = - (s% standard_cgrav)* s% m(k) /
(4*pi* (s% r(k)**4 )) -((   xn(k+1)-    xn(k)) / s% dm_bar(k))

where x is a new quantity. I am not sure about setting

s% extra_dpdm(k) = - ((   xn(k+1)-    xn(k)) / s% dm_bar(k))

Is it correct? or it must be divided by the interpolated pressure to be
consistent with the hydrostatic equilibrium equation mesa solve.

Another question is :  on which side of the hydrostatic equilibrium
equation, the term (s% extra_dpdm) is added into the mesa code ?
 In other words, what is the sign of (dx/dm) in the (s% extra_dpdm) in my
case?


Thank you in advance

Sheyda Najafi
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20170119/b9388a7e/attachment.html>