[mesa-users] Eq. 9 in Paxton et al. 2011

[Ehsan Moravveji]

Thanks Bill for your detailed explanations.

>         dP0 = s% surf_bc_offset_factor*s% cgrav(1)*s% mstar*s% dm(1)/(8*pi*r*r*r*r)
>         dT0 = dP0*s% gradT(1)*s% T(1)/s% P(1)
There still remains a puzzle for me to resolve:
I understand the necessity of the interlaced mesh for intensive and extensive quantities for the sake of numerical efficiencies. dP0 equation is Okay.
However, what I see is the blend of two variable types in dT0 equation above; here, gradT and dP0 are defined at the face, whereas P and T are defined as cell averages. But, both are used on the RHS of dT0 equation.
Why is this an allowed operation?
Am I missing a point here?

I will try test_suite/low_z with/without s% surf_bc_offset_factor and let you know.

Happy weekend,
Ehsan.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20140914/1087566a/attachment.html>