[Mesa-users] length of cell

[Ebrahim Hassani]

Dear Mesa-users

Thank you for your responses to my question.

I guessed there might be a parameter in the 'star_data.inc' file that
represents the length of each cell. But as you mentioned there is not.
To solve my problem I followed Dr. Marchant's advice and calculated the
length of the last cell outside the do loop as follow:

      part3_h1 = 0
      do i = s% nz, 2, -1
      part3_h1 = (s% r(i)) * ( (s% r(i-1)) - (s% r(i)) )
      end do
      part3_h1 = part3_h1 + (s% r(1)) * ((s% r(1)) - (s% r(2)))

In fact, I do the loop from i = s%nz until i = 2 (a cell before the last
cell) and then calculated the length of the last cell after the do loop.

This time the star evolved without any error.

Thank you again for your time and assistance.

Sincerely yours,
Ebrahim
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20200705/86d02d6b/attachment.htm>