! profile_columns.list -- determines the contents of star model profiles ! you can use a non-standard version by setting profile_columns_file in your inlist ! units are cgs unless otherwise noted. ! reorder the following names as desired to reorder columns. ! comment out the name to omit a column (fewer columns => less IO => faster running). ! remove '!' to restore a column. ! if you have a situation where you want a non-standard set of columns, ! make a copy of this file, edit as desired, and give the new filename in your inlist ! as profile_columns_file. if you are just adding columns, you can 'include' this file, ! and just list the additions in your file. note: to include the standard default ! version, use include '' -- the 0 length string means include the default file. ! if you need to have something added to the list of options, let me know.... ! the first few lines of the profile contain general info about the model. ! for completeness, those items are described at the end of this file. ! note: you can include another list by doing ! include 'filename' ! include '' means include the default standard list file ! the following lines of the profile contain info for 1 zone per row, surface to center. 1. ! numbers start with 1 at the surface 2. logT ! log10(temperature) at center of zone 3. logRho ! log10(density) at center of zone 4. logP ! log10(pressure) at center of zone 5. logR ! log10(radius/Rsun) at outer boundary of zone !log_column_depth ! log10 column depth, exterior mass / area (g cm^-2) !log_radial_depth ! log10 radial distance to surface (cm) 6. luminosity ! luminosity at outer boundary of zone (in Lsun units) !luminosity_rad ! radiative luminosity at outer boundary of zone (in Lsun units) ! -(4 pi r^2)^2 a c / (3 kap) d(T^4)/dm !lum_minus_lum_rad ! at outer boundary of zone (in Lsun units) !grav ! gravitational acceleration (cm sec^2) 7. eps_grav ! -T*ds/dt (negative for expansion) 8. signed_log_eps_grav ! sign(eps_grav)*log10(max(1,abs(eps_grav))) !net_nuclear_energy ! erg/gm/s from nuclear reactions minus all neutrino losses ! The value plotted is net_nuclear_energy = sign(val)*log10(max(1,abs(val))) ! where val = net nuclear energy minus all neutrino losses. 9. net_energy ! net_energy + eps_grav. ! The value plotted is net_energy = sign(val)*log10(max(1,abs(val))) ! where val = net nuclear energy plus eps_grav minus all neutrino losses. !logL ! log10(max(1d-2,L/Lsun)) !log_Ledd ! log10(Leddington/Lsun) -- local Ledd, 4 pi clight G m / kap !log_L_div_Ledd ! log10(max(1d-12,L/Leddington)) 10. signed_log_power ! sign(L)*log10(max(1,abs(L))) 11. velocity ! velocity at outer boundary of zone -- 0 if no velocity variable 12. entropy ! specific entropy divided by (avo*kerg) 13. mixing_type ! mixing types are defined in mesa/mlt/public/mlt_def 14. csound ! sound speed 15. v_div_csound ! velocity divided by sound speed !binding_energy ! v^2/2 - G m / r + E + P / rho (ergs/gm). negative if bound. !binding_energy_integral ! sum from surface inwards of dm*(v^2/2 - G m / r + E + P/rho) (ergs) !v_div_r ! velocity divided by radius !scale_height ! pressure scale height; P / (g rho) (in Rsun units) 16. eta ! electron degeneracy parameter (eta >> 1 for significant degeneracy) 17. mu ! mean molecular weight per gas particle (ions + free electrons) 18. logdq ! log10(dq) 19. dq_ratio ! dq(k-1)/dq(k) 20. q ! fraction of star mass interior to outer boundary of this zone 21. radius ! radius at outer boundary of zone (in Rsun units) !r_div_R ! fraction of total radius 22. temperature ! temperature at center of zone 23. tau ! optical depth 24. logtau ! log10(optical depth) at center of zone 25. pressure ! total pressure at center of zone (pgas + prad) !pgas ! gas pressure at center of zone (electrons and ions) 26. pgas_div_ptotal ! pgas/pressure 27. logPgas ! log10(pgas) !prad ! radiation pressure at center of zone !energy ! internal energy (ergs/g) !logE ! log10(specific internal energy) at center of zone 28. grada ! dlnT_dlnP at constant S !dE_dRho ! at constant T !cv ! specific heat at constant volume !cp ! specific heat at constant total pressure !logS ! log10(specific entropy) !logS_per_baryon ! log10(specific entropy per baryon / kerg) !gamma1 ! dlnP_dlnRho at constant S !gamma3 ! gamma3 - 1 = dlnT_dlnRho at constant S !theta_e ! electron degeneracy factor for graboske screening !gam ! plasma interaction parameter (> 160 or so means starting crystallization) 29. free_e ! free_e is mean number of free electrons per nucleon !logfree_e ! log10(free_e), free_e is mean number of free electrons per nucleon !chiRho ! dlnP_dlnRho at constant T !chiT ! dlnP_dlnT at constant Rho 30. abar ! average atomic weight (g/mole) !zbar ! average charge !z2bar ! average charge^2 31. ye ! average charge per baryon = proton fraction !opacity ! opacity measured at center of zone 32. log_opacity ! log10(opacity) 33. eps_nuc ! ergs/g/sec from nuclear reactions (reaction neutrinos subtracted) !eps_nuc_neu_total ! erg/gm/sec as neutrinos from nuclear reactions 34. non_nuc_neu ! non-nuclear-reaction neutrino losses !nonnucneu_plas ! plasmon neutrinos (for collective reactions like gamma_plasmon => nu_e + nubar_e) !nonnucneu_brem ! bremsstrahlung (for reactions like e- + (z,a) => e- + (z,a) + nu + nubar) !nonnucneu_phot ! photon neutrinos (for reactions like e- + gamma => e- + nu_e + nubar_e) !nonnucneu_pair ! pair production (for reactions like e+ + e- => nu_e + nubar_e) !extra_heat !log_irradiation_heat 35. mlt_mixing_length ! mixing length for mlt (cm) 36. mlt_mixing_type ! value returned by mlt 37. gradT_sub_grada ! gradT-grada at cell boundary !super_ad ! max(0,gradT-grada) at cell boundary 38. log_D_mix ! log10 diffusion coefficient for mixing in units of cm^2/second (Eulerian) !log_D_mix_non_rotation 39. log_conv_vel ! log10 convection velocity (cm/sec) 40. conv_vel_div_csound ! convection velocity divided by sound speed 41. log_mlt_D_mix ! log10 diffusion coefficient for mixing from mlt (cm^2/sec) 42. pressure_scale_height ! in Rsun units !log_D_conv ! D_mix for regions where mix_type = convective_mixing !log_D_semi ! D_mix for regions where mix_type = semiconvective_mixing !log_D_ovr ! D_mix for regions where mix_type = overshoot_mixing !log_D_th ! D_mix for regions where mix_type = thermo_haline_mixing 43. gradT ! mlt value for required temperature gradient dlnT/dlnP 44. gradr ! dlnT/dlnP required for purely radiative transport !grad_temperature ! smoothed dlnT/dlnP at cell boundary !grad_density ! smoothed dlnRho/dlnP at cell boundary !gradmu ! smoothed dlnmu/dlnP at cell boundary (uses mu from eos) !mu_alt ! calculated from actual composition assuming complete ionization. mu_alt = abar/(1 + zbar) !gradmu_alt ! smoothed dlnmu_alt/dlnP at cell boundary (uses mu from xa assuming complete ionization) !dlnmu_X ! (at cell boundary) based on hydrogen composition !gradmu_X ! smoothed dlnmu_X/dlnP at cell boundary !gradL ! gradient for Ledoux criterion for convection !gradL_composition_term ! = composition term in gradL. negative means thermohaline instability. !sch_stable ! 1 if grada > gradr, 0 otherwise !ledoux_stable ! 1 if gradL > gradr, 0 otherwise !stability_type ! values same as defined in mlt_def for mixing types ! 0 if both sch_stable and ledoux_stable, ! 1 if neither sch_stable and ledoux_stable ! 3 if ledoux_stable but not sch_stable ! 4 if sch_stable but not ledoux_stable ! the "for_mixing" values are the ones used to calculate the mixing diffusion coeffs, ! i.e., from the start of the step. ! the names without "for_mixing" are the values at the end of the step. !gradr_for_mixing !gradT_for_mixing !grada_for_mixing !gradL_for_mixing !gradL_comp_pre_hydro !gradmu_pre_hydro !gradmu_alt_pre_hydro !gradmu_X_pre_hydro 45. dlnd_dt ! time derivative of log(density) at fixed mass coordinate (Langranian) 46. dlnT_dt ! time derivative of log(temperature) at fixed mass coordinate (Langranian) 47. signed_dlnd ! sign(dlnd)*log10(max(1,abs(1d6*dlnd))) 48. signed_dlnT ! sign(dlnT)*log10(max(1,abs(1d6*dlnT))) !dv_dt ! time derivative of velocity at fixed mass coordinate (Langranian) !accel_div_grav ! dv_dt/grav -- only if v_flag is true. 0 otherwise. 49. mass ! m/Msun. mass coordinate of outer boundary of cell. 50. mmid ! mass at midpoint of cell (average of mass coords of the cell boundaries) Msun units. !logM ! log10(m/Msun) !xm ! mass exterior to point (Msun units) !dq ! mass of zone as a fraction of total star mass 51. logxq ! log10(1-q) 52. logxm ! log10(xm) !dr ! cell width (cm) !log_dr ! log10 cell width (cm) !dr_div_cs ! cell sound crossing time (sec) !log_dr_div_cs ! log10 cell sound crossing time (sec) !dr_div_cs_yr ! cell sound crossing time (years) !log_dr_div_cs_yr ! log10 cell sound crossing time (years) !grav_gr_factor ! (1/sqrt(1 - 2Gm/(rc^2)) -- only important for neutron stars !log_grav_gr_factor ! log10(grav_gr_factor) !z ! metallicity !log_z ! metallicity add_abundances ! this adds all of the isos that are in the current net ! NOTE: you can list specific isotopes by giving their names (from chem_def) 53 !h1 54 !he3 55 !he4 56 !c12 57 !n14 58 !o16 !add_log_abundances ! this adds log10 of all of the isos that are in the current net ! NOTE: you can list specific isotopes by giving their names (from chem_def) !log h1 !log he3 !log he4 !log c12 !log n14 !log o16 ! ergs/g/sec for reaction categories add_reaction_categories ! this adds all the reaction categories ! NOTE: you can list specific categories by giving their names (from net_def) 61 !pp 62 !cno 63 !tri_alfa ! equation residuals !equP ! momentum conservation !equT ! energy transport !equR ! mass-volume-density relation !equL ! energy conservation !equv ! velocity definition !equJ ! angular momentum !equ_chem h1 !equ_chem he4 ! any other isotopes you want ! rotation !omega ! angular velocity = j_rot/i_rot !log_omega !log_j_rot !log_J_div_M53 ! J is j*1e-15 integrated from center; M53 is m^(5/3) !log_J_inside ! J_inside is j_rot integrated from center !shear ! -dlnomega/dlnR !log_abs_shear ! log10(abs(dlnomega/dlnR)) !domega_dt ! time derivative of omega at fixed mass coordinate (Langranian) !delta_omega ! omega minus omega_pre_hydro !i_rot ! specific moment of interia at cell boundary !j_rot ! specific angular momentum at cell boundary !v_rot ! rotation velocity at cell boundary (km/sec) !fp_rot ! rotation factor for pressure !ft_rot ! rotation factor for temperature !log_am_nu ! angular momentum diffusion coef at cell boundary ! this is the rotational_viscosity, nu, from Heger 2000 eqn for ang. mom. transport. !r_polar ! (Rsun) !log_r_polar ! log10 (Rsun) !r_equitorial ! (Rsun) !log_r_equitorial ! log10 (Rsun) !r_e_div_r_p ! r_equitorial/r_polar !omega_crit ! breakup angular velocity = sqrt(G M / r_equitorial^3) !omega_div_omega_crit !am_log_sig !am_log_D_DSI ! diffusion coeff for dynamical shear instability !am_log_D_SH ! diffusion coeff for Solberg-Hoiland instability !am_log_D_SSI ! diffusion coeff for secular shear instability !am_log_D_ES ! diffusion coeff for Eddington-Sweet circulation !am_log_D_GSF ! diffusion coeff for Goldreich-Schubert-Fricke instability !am_log_D_ST ! Spruit dynamo mixing diffusivity !am_log_nu_ST ! Spruit dynamo effective viscosity !dynamo_log_B_r ! (Gauss) !dynamo_log_B_phi ! (Gauss) !am_coeff_D_SSI ! strength coeff for D_SSI in [0..1] ! misc !gradr_sub_grada ! gradr - grada; > 0 => Schwarzschild unstable for convection !dlnP_dm ! for structure equation !dlnT_dm ! for structure equation !dL_dm ! for structure equation !dlnP_dlnm ! for structure equation !dlnT_dlnm ! for structure equation !delta_r ! r(outer edge) - r(inner edge); radial extent of cell in cm. !delta_v ! v(inner edge) - v(outer edge); rate at which delta_r is shrinking (cm/sec). !dt_dv_div_dr ! dt*delta_v/delta_r; need to have this << 1 for every cell ! element diffusion velocity for species !edv h1 ! element diffusion velocity for h1 !edv he4 ! element diffusion velocity for he4 !edv o16 ! element diffusion velocity for o16 ! ionization state for given species !ionization he4 !ionization c12 !ionization fe52 !cno_div_z ! abundance of c12, n14, and o16 as a fraction of total z !dlog_h1_dlogP ! (log(h1(k)) - log(h1(k-1)))/(log(P(k)) - log(P(k-1))) !dlog_he3_dlogP !dlog_he4_dlogP !dlog_c12_dlogP !dlog_c13_dlogP !dlog_n14_dlogP !dlog_o16_dlogP !dlog_ne20_dlogP !dlog_mg24_dlogP !dlog_si28_dlogP !dlog_pp_dlogP !dlog_cno_dlogP !dlog_3alf_dlogP !dlog_burn_c_dlogP !dlog_burn_n_dlogP !dlog_burn_o_dlogP !dlog_burn_ne_dlogP !dlog_burn_na_dlogP !dlog_burn_mg_dlogP !dlog_cc_dlogP !dlog_co_dlogP !dlog_oo_dlogP !dlog_burn_si_dlogP !dlog_burn_s_dlogP !dlog_burn_ar_dlogP !dlog_burn_ca_dlogP !dlog_burn_ti_dlogP !dlog_burn_cr_dlog !dlog_burn_fe_dlogP !dlog_pnhe4_dlogP !dlog_photo_dlogP !dlog_other_dlogP !brunt_N2 ! brunt-vaisala frequency squared !brunt_N2_dimensionless ! N2 in units of 3GM/R^3 !brunt_N_dimensionless ! N in units of sqrt(3GM/R^3) !brunt_frequency ! cycles per day !brunt_N ! sqrt(abs(brunt_N2)) !log_brunt_N ! log10(brunt_N) !log_brunt_N2 ! log10(brunt_N2) !sign_brunt_N2 ! sign of brunt_N2 (+1 for Ledoux stable; -1 for Ledoux unstable) !lamb_S ! for l=1: S = sqrt(2)*csound/r !dlnRho_dlnR ! smoothed numerical difference for use with brunt !logQ ! logQ = logRho - 2*logT + 12 !log_L_div_CpTMdot ! (unitless) !cs_at_cell_bdy ! sound speed at cell boundary (csound is at cell center) ! the first few lines of the profile contain general info about the model. ! for completeness, those items are described here. ! initial mass and Z ! initial_mass ! initial_z ! general properties of the current state ! model_number ! num_zones ! star_age ! time_step ! properties at the photosphere ! Teff ! photosphere_L ! photosphere_r ! properties at the outermost zone of the model ! log_surface_L ! log_surface_radius ! log_surface_temp ! properties near the center of the model ! log_center_temp ! log_center_density ! log_center_P ! center_eta ! abundances near the center ! center_h1 ! center_he3 ! center_he4 ! center_c12 ! center_n14 ! center_o16 ! center_ne20 ! information about total mass ! star_mass ! star_mdot ! star_mass_h1 ! star_mass_he3 ! star_mass_he4 ! star_mass_c12 ! star_mass_n14 ! star_mass_o16 ! star_mass_ne20 ! locations of abundance transitions ! h1_boundary_mass ! he4_boundary_mass ! c12_boundary_mass ! location of optical depths 10 and 100 ! tau10_mass ! tau10_radius ! tau100_mass ! tau100_radius ! time scales ! dynamic_time ! kh_timescale ! nuc_timescale ! various kinds of total power ! power_nuc_burn ! power_h_burn ! power_he_burn ! power_neu ! a few control parameter values ! h1_boundary_limit ! he4_boundary_limit ! c12_boundary_limit ! burn_min1 ! burn_min2