Документ взят из кэша поисковой машины. Адрес оригинального документа : http://crydee.sai.msu.ru/mesa/profile_columns_modified.list Дата изменения: Mon Apr 9 16:37:05 2012 Дата индексирования: Mon Oct 1 21:42:38 2012 Кодировка:

! 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