[mesa-users] Printing planet profile properties

[Gabriel-Dominique Marleau]

Dear MESA/star users,

With in mind the modeling of giant gas planets (which seems to be garnering a healthy amount of attention in this mailing list), would someone know how to:

1. Print log*.data files of the "on-the-side" models calculated while adding a core with relax_M_center and relax_R_center on (as explained in star/public/run_star_defaults.dek)?

2. Add to the log files the composition in terms of the H, H_2, He mass fractions, i.e. track their dissociation and ionisation fractions within the star? I could not find a relevant column name in mesa/data/star_data/profile_columns.list but I may be missing something simple.

3. Get the following columns to print?
- chiY ! dlnP/dlnY used in calculate of brunt_B term
- dlnY_dlnP ! the actual dlnY/dlnP in the model as used in brunt_B
- dlnmu_dr ! for salt-finger mixing
- semiconvection_criterion ! positive means stable
Including them in profile_columns.list with "use_Ledoux_criterion = .true." and values for alpha_semiconvection and thermo_haline_coeff yields for each a "failed to recognize item" error.

Thank you for any help,

Gabriel-Dominique Marleau
Department of Physics
McGill University