[mesa-users] Density in profile files

[Josiah Schwab]

Hi Bill,

> use "rho" instead of "density".  The profile_columns.list in defaults
> is incomplete and doesn't include "rho".  But it is
> profile_getval.f90.

This made me wonder how many MESA columns I'm missing out on :)

I ran the following (painful to read) one-liners, which see which
profile (history) columns that are in profile_getval.f90 (history.f90)
are not in profile_columns.list (history_columns.list).

This makes the assumption that the flag for history (profile) column
"name" is h_name (p_name), which is at least mostly correct.

Some of these columns appear to be for debugging and so probably
shouldn't be publicly exposed, but many (e.g., Rich's rho) should be
included in the *.list files in a future release.

Josiah


grep -oP "case \(p_\K[A-Za-z0-9_]+" profile_getval.f90 | parallel "grep -qE '^[ \t]+\!?[ \t]?{}' ../defaults/profile_columns.list || echo {}"

------------------------------------------------------------------------
lum_rad_div_L_Edd_div_phi_Joss
log_Lconv
eps_nuc_minus_non_nuc_neu
log_abs_eps_nuc_div_non_nuc_neu
rmid_div_R
sum_check_dEdt_expected
dEdt_actual
check_dEdt_actual
sum_check_dEdt_actual
kap_frac_Type2
rho
scale_height_start
nonnucneu_reco
d_nonnucneu_dlnd
d_nonnucneu_dlnT
extra_jdot
extra_omegadot
eps_grav_composition_term
area
area_times_v
div_v
d_v_div_r_dm
d_v_div_r_dr
dm_div_4pi_rmid2_rho_dr
rmid3_rho
log_rmid3_rho
eta_visc
dP_dm_visc
eps_visc
qvisc
dvdt_visc
dvdt_grav
dvdt_dPdm
dm_eps_grav
env_eps_grav
mlt_Gamma
mlt_Zeta
adjust_mlt_gradT_fraction
burn_num_iters
burn_num_substeps
num_steps
mtx_solve
mtx_factor
avg_order
gradL_sub_gradr
grada_sub_gradr
eps_nuc_start
sig_raw_div_sig
rp_div_GMrho_face
conv_mixing_type
grada_sub_actual_gradT
log_dr_div_scale_height
omega_start
fp_rot_minus_1
ft_rot_minus_1
am_domega_dlnR
gradr_start
grada_at_face_start
gradL_start
gradr_sub_gradL_start
mlt_mixing_length_start
mlt_mixing_type_start
mlt_D_start
mlt_vc_start
mlt_Gamma_start
mlt_cdc_start
dlnR_dk
dlnP_dm_visc
log_maxabs_dt_d_dxdt_dt
energy_check
vcheck
echeck
dm_eps_grav_plus_eps_visc
sum_echeck_from_surf
dlnH1_dlnP
dlnHe3_dlnP
dlnHe4_dlnP
dlnC12_dlnP
dlnC13_dlnP
dlnN14_dlnP
dlnO16_dlnP
dlnNe20_dlnP
dlnMg24_dlnP
dlnSi28_dlnP
dlnY_dr
dlnRho_dr
brunt_g_r_div_cs2
lamb_Sl1
log_lamb_Sl1
log_lamb_Sl10
log_mdot_cs
log_mdot_v
log_L_div_CpTMdot
------------------------------------------------------------------------



grep -oP "case \(h_\K[A-Za-z0-9_]+" history.f90 | parallel "grep -qE '^[ \t]+\!?[ \t]?{}' ../defaults/history_columns.list || echo {}"

------------------------------------------------------------------------
kh_mdot_limit
min_L
min_dL_dm
min_dL_dm_m
k_for_test_CpT_absMdot_div_L
q_for_test_CpT_absMdot_div_L
logxq_for_test_CpT_absMdot_div_L
split_mixing_choice
burn_max_num_substeps
burn_max_num_iters
cells_using_burn_substeps
total_radiation
total_energy_plus_total_radiation
------------------------------------------------------------------------