<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Dear Noam,<div class=""><br class=""></div><div class="">Regarding how the hydrostatic equlibrium is satisfied in your mesa model, you can follow the following page (exercises for master students at KU Leuven), and perhaps that walks you through the necessary steps to compute dP/dr.</div><div class=""><a href="http://www.ster.kuleuven.be/~ehsan/sse/exercise05/exercise05.html" class="">http://www.ster.kuleuven.be/~ehsan/sse/exercise05/exercise05.html</a></div><div class=""><br class=""></div><div class="">Regarding the meshing in MESA, please dive into the first instrument paper (Paxton et al. 2011). You will see quite some light there.</div><div class=""><br class=""></div><div class="">Best,</div><div class="">Ehsan.</div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On 05 Oct 2016, at 09:43, Noam Chai <<a href="mailto:noamchai@gmail.com" class="">noamchai@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="rtl" class=""><div style="direction:ltr" class="">Hi,</div><div style="direction:ltr" class=""><br class=""></div><div style="direction:ltr" class="">I'm having a trouble calculating the hydrostatic equilibrium for a simple star with the mass of the sun. I use matlab 2014a to calculate the gradient of dpdr, by using two methods: </div><div style="direction:ltr" class="">1. the inner function 'gradient'.</div><div style="direction:ltr" class="">2. with a for loop - (p(i+1)-p(i-1))/(r(i+1)-r(i-1))</div><div style="direction:ltr" class=""><br class=""></div><div style="direction:ltr" class="">The plot from both methods have lots of significant spikes that looks like noise in it. Did anyone struggled with this problem? why does it happend? </div><div style="direction:ltr" class=""><br class=""></div><div style="direction:ltr" class="">Another question: </div><div style="direction:ltr" class="">How does mesa decides the intervals of the radius? it seems that the intervals are not constant. How can I change the length of the vector r and how can I make it an evenly spaced vector?</div><div style="direction:ltr" class=""><br class=""></div><div style="direction:ltr" class="">Thank you all,</div><div style="direction:ltr" class=""><br class=""></div><div style="direction:ltr" class="">Noam</div></div>
------------------------------------------------------------------------------<br class="">Check out the vibrant tech community on one of the world's most <br class="">engaging tech sites, <a href="http://SlashDot.org" class="">SlashDot.org</a>! <a href="http://sdm.link/slashdot_______________________________________________" class="">http://sdm.link/slashdot_______________________________________________</a><br class="">mesa-users mailing list<br class=""><a href="mailto:mesa-users@lists.sourceforge.net" class="">mesa-users@lists.sourceforge.net</a><br class="">https://lists.sourceforge.net/lists/listinfo/mesa-users<br class=""></div></blockquote></div><br class=""></div></body></html>