[mesa-users] adding rotation

[Catherine L. <clo...@mt...>]

Hi everyone,


I'm trying to run some rotating models (1.2-2.2 solar masses), and I've run into a problem.  For now, I'm trying to include rotation in the simplest way possible,  so I've got an inlist that has the following:

&star_job

   show_log_description_at_start = .false.

   eos_file_prefix = 'mesa'
   kappa_file_prefix = 'gs98'

   change_v_flag = .true.
   new_v_flag = .true.

   new_rotation_flag = .true.
   change_rotation_flag = .true.

   relax_initial_surface_rotation_v = .true.
   set_surface_rotation_v = .true.
   new_surface_rotation_v = 100


/ ! end of star_job namelist


When I run a model, it starts out with whatever rotation rate I set on the last line, but around 30 models in, I see:
 turn on rotation
                               new_surface_rotation_v =    4.5000000000000000D+01    3.7745766462731172D-05

and the rotation rate of my star drops.  I can confirm that this is what happens when I plot omega or v as a function of the age of the star. The rotation velocity starts out where I set it, and then drops to close to 45 km/s and evolves from there.

 I get a surface rotation of 45 km/s for my model even if I comment all of the rotation commands out of my inlist, so it seems that my version of mesa is reading a velocity from somewhere else, or there's a default that isn't getting changed somehow.   I've tried deleting all my inlists and making new ones, as well as recompiling the code, but the problem doesn't go away.  Is there another rotation parameter I'm not setting, or another place I should be changing the rotation?  I'm running version 8845 on a Mac OSX 10.9.5.


Thanks,


Catherine


Catherine Lovekin
Assistant Professor
Department of Physics
Dunn 222
Mount Allison University
clo...@mt...
phone: (506) 364-2582