[mesa-users] Density in profile files

[Bill Paxton]

Great!   Please go ahead with updating the default files with whatever you think should be included.  And send me your new versions.
And add more things as you think appropriate.  Redundancy is fine for this.  For example, I added "density" as an alternative to "rho".

-b


On Mar 12, 2015, at 7:48 AM, Josiah Schwab wrote:

> 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
> ------------------------------------------------------------------------