[mesa-users] Kjeldsen correction in astero module
Kuldeep Verma
kuldeepv89 at gmail.com
Sun Feb 15 01:48:13 EST 2015
Hi Santiago,
Check if a_div_r is zero in sample_#.data. If so, value of "correction_a"
would have come out to be positive for the model. Take a look at around
line number 222 of "mesa/star/astero/src/astero_support.f".
Try calculating observed delta_nu in exactly same way as MESA does for the
model (slope of linear least-squares fit to l=0 modes, Equation (8) of
Kjeldsen et al. (2008)).
With Best Wishes,
Kuldeep
On Sat, Feb 14, 2015 at 10:51 PM, Santiago Andres Triana <repepo at gmail.com>
wrote:
> Hi MESA users,
>
> Playing around with the astero module and simplex, I noticed that the
> frequencies, as reported in the sample_#.data files, are not
> corrected. They seem to be the raw, uncorrected frequencies from gyre
> although I have specified the 'kjeldsen' method in
> inlist_astero_search_controls.
>
> The inlist is attached. The relevant snippet is this:
>
> ! surface corrections
> correction_scheme = 'kjeldsen' ! options are
> ! 'kjeldsen' Correction of Kjeldsen et al. (2008)
> ! 'cubic' Cubic correction of Ball & Gizon (2014, eqn 3)
> ! 'combined' Combined correction of Ball & Gizon (2014, eq 4)
> ! '' no corrections
> correction_factor = 1
> ! use this fraction of the correction; set to 0 to skip doing
> corrections.
> l0_n_obs(:) = -1 ! the observed radial orders (ignored if < 0)
> ! the observed radial orders are used in calculating surface
> corrections
> ! if <= 0, use default calculation for radial orders
> correction_b = 4.25d0
>
> save_next_best_at_higher_frequency = .true.
> save_next_best_at_lower_frequency = .true.
>
>
> I'm convinced that I'm doing something silly somewhere but cannot
> figure out where... any hints or comments are very much appreciated.
> Many thanks in advance!
>
> Cheers,
>
> Santiago
>
>
> ------------------------------------------------------------------------------
> Dive into the World of Parallel Programming. The Go Parallel Website,
> sponsored by Intel and developed in partnership with Slashdot Media, is
> your
> hub for all things parallel software development, from weekly thought
> leadership blogs to news, videos, case studies, tutorials and more. Take a
> look and join the conversation now. http://goparallel.sourceforge.net/
> _______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150215/9d4499ba/attachment.html>
More information about the Mesa-users
mailing list