[mesa-users] Extra mesh around convective boundaries

[Ehsan M. <e.m...@gm...>]

Dear mesa users,
I have difficulty to “resolve” the position of the convective boundaries (especially the core), by adding extra mesh points around them, and I would like to ask for your inputs/experience. Maybe some of you have already solved this.

As far as I understand, in the controls.defaults, the following set of options allows lowering mesh_delta_coeff around the boundaries by any desired value. When I change the following controls by giving quite smaller coefficients (than default) over a large distance (in pressure scale height) around the core convective boundary, the total number of mesh points changes very slightly. 

      xtra_coef_czb_full_on = 1
      xtra_coef_czb_full_off = 1

      xtra_coef_a_l_nb_czb = 0.5
      xtra_dist_a_l_nb_czb = 0.5

      xtra_coef_b_l_nb_czb = 0.5
      xtra_dist_b_l_nb_czb = 0.5

      xtra_coef_a_u_nb_czb = 0.5
      xtra_dist_a_u_nb_czb = 0.5

      xtra_coef_b_u_nb_czb = 0.5
      xtra_dist_b_u_nb_czb = 0.5

      xtra_coef_a_l_hb_czb = .01
      xtra_dist_a_l_hb_czb = 5

      xtra_coef_b_l_hb_czb = .01
      xtra_dist_b_l_hb_czb = 5

      xtra_coef_a_u_hb_czb = .01
      xtra_dist_a_u_hb_czb = 5

      xtra_coef_b_u_hb_czb = .01
      xtra_dist_b_u_hb_czb = 5

To give you some numbers, without these controls, the last profile has ~980 mesh, and including them, it reaches around ~1020. However, I was expecting quite a sizeable increase in the total number of zones, since I am requesting a massive amount of mesh to be added!

My question is:
- Am I giving correct values to these controls?
- How is it possible to resolve the convective boundaries by increasing the mesh around them (up to a desired H_p away)?

I have attached my inlist_project (for version 7678).

Best regards
Ehsan.