[mesa-users] Rotation : relax_initial_omega_div_omega_crit + preserving omega in several inlists

[Francis Timmes]

hi hector,

i'm happy to see you're poking into the source code.
i'm going to encourage you to delve deeper!

your first task is figure out what three variables are being written out 
that accompany the "relax to omega: ..." print statement. then figure out 
how the third variable, which might be named new_omega, gets a value of ~1.7e-5. 
this will help align/clarify your expectation of ~0.83 * 0.1.

report what you find!

fxt




> On Jun 16, 2016, at 8:48 PM, Héctor MR <hector.mr at pitt.edu> wrote:
> 
>   Dear all,
> 
>   I have a question regarding the &star_job command relax_initial_omega_div_omega_crit. [Disclosure: this is the first time I am implementing rotation in MESA]
> 
>   I am analyzing rotation on an accreting, 1Msun-WD with radius R=8.25E-3*Rsun, so that the break-up frequency sqrt(G*M/R**3) is ~0.83; however, when I use new_omega_div_omega_crit = 0.1, I get
> 
> 
>   num_steps_to_relax_rotation         100
>         relax to omega: wanted-current, current, wanted           1    1.7485786361779884D-05    0.0000000000000000D+00    1.7485786361779884D-05
>         relax to omega: wanted-current, current, wanted           2    1.7330581050599486D-05    1.7485788969702041D-07    1.7505438940296507D-05
>         relax to omega: wanted-current, current, wanted           3    1.7177934014914814D-05    3.5010756599501660D-07    1.7528041580909830D-05
>         relax to omega: wanted-current, current, wanted           4    1.7016693789573586D-05    5.2583873499324927D-07    1.7542532524566834D-05
>         relax to omega: wanted-current, current, wanted           5    1.6850228116883350D-05    7.0169673253464221D-07    1.7551924849417992D-05
>         relax to omega: wanted-current, current, wanted           6    1.6680900334687578D-05    8.7758847614349309D-07    1.7558488810831069D-05
>         relax to omega: wanted-current, current, wanted           7    1.6513051889555891D-05    1.0534860275206089D-06    1.7566537917076500D-05
> 
> 
>   Rather than, let's say, 0.083, as I was expecting. I have done some "research" with "grep -r" in $MESA_DIR/star/private:
> 
>   /relax.f90: integer function relax_omega_check_model
> 
>   /do_one_utils: integer function do_check_limits
> 
>   /star_utils.f90: real(dp) function omega_crit(s, k), subroutine set_surf_avg_rotation_info(s)
> 
>   /hydro_vars.f90: subroutine set_rmid_and_Amid        !definition of rmid
> 
> 
>   And have not found anything weird. Is there anyone who can tell me what I am either doing wrong, obviating or misuderstanding? This sums everything up:
> 
>       new_rotation_flag = .true.
>       change_initial_rotation_flag = .true.
> 
>       !set_initial_omega_div_omega_crit = .true.
>       new_omega_div_omega_crit = 0.1
> 
>       relax_initial_omega_div_omega_crit = .true.
>       num_steps_to_relax_rotation = 100
> 
> 
>   Finally, my inlists are split in two pieces. How can I preserve the last value of omega_div_omega_crit in the first one for the initial angular velocity in the second one? Via x_ctrl(1)? And if so, how?
> 
> 
>   Thank you so much!
> 
> 
>   --
>   Héctor
> ------------------------------------------------------------------------------
> What NetFlow Analyzer can do for you? Monitors network bandwidth and traffic
> patterns at an interface-level. Reveals which users, apps, and protocols are 
> consuming the most bandwidth. Provides multi-vendor support for NetFlow, 
> J-Flow, sFlow and other flows. Make informed decisions using capacity planning
> reports. http://sdm.link/zohomanageengine_______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users