[mesa-users] convergence problems

[Josiah Schwab]

Hi Zoe,

> I am trying to add an expression to "extra_dpdm" in other_momentum
> subroutine which has been added to run_star_extras.f of my work
> directory.

In general, it is a good idea to briefly describe your physical
motivation.  For example, that lets us understand whether you need to
include this term during the pre-main sequence evolution of the star.
As indicated in the error message, MESA is struggling while constructing
the PMS model.

> Here is my inlists http://www.uploadman.com/5482b953b1dbf96e or
> http://dsa.uploadman.com/5482b953b1dbf96e?download_token=d735d3feb5bc080405aafe5c3f265a8327d5c4ea385d0ad7685f051395bf809d

I looked at your run_star_extras.f.

I would venture a guess that if you print the value of s% extra_dPdm
you'll be surprised :)

I'll give you two hints to start.  The line

  integer :: k, Y=6d-6

and the lines like 

  rr= (s% r(k+1) - s% r(k))*(1/4)*(pi)*s% r(k)**2

almost certainly don't do what you think they do.

Why don't you try printing the values of Y and rr?  

Josiah

P.S.  There's no need to include the full contents of your working
directory.  We don't need the 20 MB of executable and restart files.
The inlists and run_star_extras files are only a few KB.  That means you
will be able to directly attach them to your messages, making it even
easier for people to take a look.