radial order anomaly like -47 -46 -33 -42 -43 -42

[rhtownsend]

The dips are too sharp aren't they? The previous literature ( Bouabid 2013, Miglio 2008 for theory and Van Reeth 2015 for observation) never show such sharp and deep glitches. And the diagram in profile926.pdf should be similar to other two diagrams.

Cheers,

Gang

Hi Gang --

It looks like the sharp dip in the period spacings, for 925 and 927, are a result of an avoided crossing between modes trapped in different parts of the star. I'm not altogether sure what's happening with 926 -- you might want to inspect eigenfunctions for some of the modes to see if they look reasonable.

But more generally, I'm yet to be convinced that the behavior you're seeing is due to a bug in GYRE, or -- as was the case previously -- due to an inadequacy in the input model. If you can present clear evidence that there's a bug in GYRE, then I'm happy to work on fixing it. However, I'm afraid I don't have the time to fix issues with the input model -- you'll have to work on it yourself, and/or seek help elsewhere.

cheers,

Rich

Hi Rich,

I solved that problem! We should use 'mesh grid' correctly to avoid it. The MESA inlist from C. Aerts in MESA summer school 2016 is a good guidance.

Cheers

Gang

Hi Gang --

I'm glad you've solved the problem. But what do you mean by 'mesh grid'?

cheers,

Rich

[gali8292]

Hi Rich,

Here is the inlist from C. Aerts.

Cheers,

LI Gang

! adjustment of the mesh grid
! This is important to properly resolve the gravity modes near the convective core boundary

max_allowed_nz = 10000

! global mesh resolution factor
mesh_delta_coeff = 0.2

! Additional resolution based on the pressure and temperature profiles
P_function_weight = 30
T_function1_weight = 75

! Additional resolution depending on the chemical abundances and gradients
xa_function_species(1) = 'he4' ! name of nuclide as defined in chem_def
xa_function_weight(1) = 80
xa_function_param(1) = 1d-2

xa_function_species(2) = 'he3' ! name of nuclide as defined in chem_def
xa_function_weight(2) = 20
xa_function_param(2) = 1d-5

mesh_dlogX_dlogP_extra = 0.15 ! resol coeff for chemical gradients
mesh_dlogX_dlogP_full_on = 1d-6 ! additional resol on for gradient larger than this
mesh_dlogX_dlogP_full_off = 1d-12 ! additional resol off for gradient smaller than this

mesh_logX_species(1) = 'he4' ! taking into account abundance of He4


! Additional resolution near the boundaries of the convective regions
xtra_coef_czb_full_on = 0.0
xtra_coef_czb_full_off = 1.0

xtra_coef_a_l_nb_czb = 0.015 ! resol coeff above lower nonburn convective boundary
xtra_dist_a_l_nb_czb = 10 ! distance above lower nonburn convective boundary

xtra_coef_b_l_nb_czb = 0.075 ! resol coeff below lower nonburn convective boundary
xtra_dist_b_l_nb_czb = 3 ! distance below lower nonburn convective boundary

! Properly constrain the location of the convective boundaries
convective_bdy_weight = 1
convective_bdy_dq_limit = 1d-5
convective_bdy_min_dt_yrs = 1d-3