# Python routines to handle MESA-Web files # # Version: 2019-Oct-8 # Imports import numpy as np import re # Utility routines def num (s): try: return int(s) except ValueError: sr = re.sub(r'D', 'E', s) sr = re.sub(r'(\d)([+-])', r'\1E\2', sr) try: return float(sr) except ValueError: return s # Read a MESA-Web history file def read_history (filename) : """Read data from a MESA-Web history file Parameters ---------- filename : string giving name of history file Returns ------- hist_data: dict containing header and history data (see below for details) Header Data ----------- The following keys/value pairs in the data dict contain header data -- i.e., scalars describing time-independent properties of the star. Where applicable, units are given in square brackets []. version_number -- version number of MESA initial_mass -- initial mass [Msun] initial_z -- initial metal mass fraction burn_min1 -- 1st limit for reported burning [erg/g/s] burn_min2 -- 2nd limit for reported burning [erg/g/s] History Data ------------ The following keys/value pairs in the data dict contain history data -- i.e., arrays describing global properties of the star over a sequence of time-steps. Where applicable, units are given in square brackets []. model_number -- model number star_age -- stellar age [years] star_mass -- stellar mass [Msun] log_L -- log10(stellar luminosity [Lsun]) log_R -- log10(stellar radius [Rsun]) log_Teff -- log10(effective temperature [K]) log_center_T -- log10(center temperature [K]) log_center_Rho -- log10(center density [g/cm^3]) log_center_P -- log10(center pressure [dyn/cm^2]) center_h1 -- center 1H mass fraction center_he3 -- center 3He mass fraction center_he4 -- center 4He mass fraction center_c12 -- center 12C mass fraction center_n14 -- center 14N mass fraction center_o16 -- center 16O mass fraction center_ne20 -- center 20Ne mass fraction center_mg24 -- center 24Mg mass fraction center_si28 -- center 28Si mass fraction center_s32 -- center 32S mass fraction center_ar36 -- center 36Ar mass fraction center_ca40 -- center 40Ca mass fraction center_ti44 -- center 44Ti mass fraction center_cr48 -- center 48Cr mass fraction center_fe52 -- center 52Fe mass fraction center_fe54 -- center 54Fe mass fraction center_fe56 -- center 56Fe mass fraction center_ni56 -- center 56Ni mass fraction center_degeneracy -- center electron chemical potential [kB*T] center_ye -- center average charge per baryon [e] center_entropy -- center entropy [kB] compactness_parameter -- (m/Msun)/(R(m)/1000km) for m = 2.5 Msun dynamic_timescale -- dynamical timescale [s] kh_timescale -- Kelvin-Helmholtz timescale [s] nuc_timescale -- nuclear timescale [s] pp -- log10(total pp luminosity [Lsun]) cno -- log10(total CNO luminosity [Lsun]) tri_alfa -- log10(total triple-alpha luminosity [Lsun]) log_LH -- log10(total H-burning luminosity, excluding neutrinos [Lsun]) log_LHe -- log10(total He-burning luminosity, excluding neutrinos [Lsun]) log_LZ -- log10(total metal-burning luminosity, excluding neutrinos [Lsun]) log_Lneu -- log10(total neutrino luminosity [Lsun]) he_core_mass -- mass of helium core [Msun] c_core_mass -- mass of carbon core [Msun] o_core_mass -- mass of oxygen core [Msun] si_core_mass -- mass of silicon core [Msun] fe_core_mass -- mass of iron core [Msun] he_core_radius -- radius of helium core [Rsun] c_core_radius -- radius of carbon core [Rsun] o_core_radius -- radius of oxygen core [Rsun] si_core_radius -- radius of silicon core [Rsun] fe_core_radius -- radius of iron core [Rsun] max_abs_v_velocity -- maximum absolute velocity surf_avg_omega_div_omega_crit -- surface average rotation angular frequency [Omega_crit] log_total_angular_momentum -- log10(total angular momentum [cm^2 g/s] surf_avg_omega -- surface average rotation angular frequency [rad/s] surf_avg_v_rot -- surface average rotation velocity [km/s] star_mdot -- mass-loss rate [Msun/year] """ return read_data(filename, rev=False) # Read a MESA-Web profile file def read_profile (filename) : """Read data from a MESA-Web profile file Parameters ---------- filename : string giving name of profile file Returns ------- prof_data: dict containing header and profile data (see below for details) Header Data ----------- The following keys/value pairs in the data dict contain header data -- i.e., scalars describing position-independent properties of the star. Where applicable, units are given in square brackets []. star_mdot -- mass-loss rate [Msun/year] model_number -- model number num_zones -- number of zones initial_mass -- initial mass [Msun] initial_z -- initial metal mass fraction star_age -- stellar age [years] time_step -- current time-step [s] Teff -- effective temperature [K] photosphere_L -- photospheric luminosity [Lsun] photosphere_r -- photospheric radius [Rsun] center_eta -- center electron chemical potential [kB*T] center_h1 -- center 1H mass fraction center_he3 -- center 3He mass fraction center_he4 -- center 4He mass fraction center_c12 -- center 12C mass fraction center_n14 -- center 14N mass fraction center_o16 -- center 16O mass fraction center_ne20 -- center 20Ne mass fraction star_age -- stellar age [years] star_mass -- stellar mass [Msun] star_mdot -- mass-loss rate [Msun/year] star_mass_h1 -- total mass in 1H [Msun] star_mass_he3 -- total mass in 3He [Msun] star_mass_he4 -- total mass in 4He [Msun] star_mass_c12 -- total mass in 12C [Msun] star_mass_n14 -- total mass in 14N [Msun] star_mass_o16 -- total mass in 16O [Msun] star_mass_ne20 -- total mass in 20Ne [Msun] he_core_mass -- mass of helium core [Msun] c_core_mass -- mass of carbon core [Msun] o_core_mass -- mass of oxygen core [Msun] si_core_mass -- mass of silicon core [Msun] fe_core_mass -- mass of iron core [Msun] neutron_rich_core_mass -- mass of neutron-rich core [Msun] tau10_mass -- mass coordinate of optical depth 10 [Msun] tau10_radius -- radius coordinate of optical depth 10 [Rsun] tau100_mass -- mass coordinate of optical depth 100 [Msun] tau100_radius -- radius coordinate of optical depth 100 [Rsun] dynamic_timescale -- dynamical timescale [s] kh_timescale -- Kelvin-Helmholtz timescale [s] nuc_timescale -- nuclear timescale [s] log_LH -- log10(total H-burning luminosity, excluding neutrinos [Lsun]) log_LHe -- log10(total He-burning luminosity, excluding neutrinos [Lsun]) log_LZ -- log10(total metal-burning luminosity, excluding neutrinos [Lsun]) power_nuc_burn -- total nuclear burning luminosity, excluding photodisintegrations [Lsun] power_h_burn -- total H-burning luminosity, excluding neutrinos [Lsun] power_he_burn -- total He-burning luminosity, excluding neutrinos [Lsun] power_neu -- total neutrino luminosity [Lsun] burn_min1 -- 1st limit for reported burning [erg/g/s] burn_min2 -- 2nd limit for reported burning [erg/g/s] Profile Data ------------ The following keys/value pairs in the data dict contain profile data -- i.e., describing local properties of the star over a sequence of spatial zones. Where applicable, units are given in square brackets []. mass -- mass coordinate [Msun] radius -- radius coordinate [Rsun] luminosity -- luminosity [Lsun] pressure -- pressure [dyn/cm^2] logRho -- log10(density [g/cm^3]) logT -- log10(temperature [K]) energy -- log10(specific internal energy [erg/g]) entropy -- log10(specific entropy [kB*N_avo]) cp -- specific heat at constant pressure [erg/K/g] gamma1 -- first adiabatic exponent grada -- adiabatic temperature gradient mu -- mean molecular weight free_e -- free electrons per nucleon ye -- average charge per baryon [e] pgas -- gas pressure [dyn/cm^2] prad -- radiation pressure [dyn/cm^2] gradr -- radiative temperature gradient gradT -- physical temperature gradient velocity -- velocity [km/s] conv_vel -- convective velocity [km/s] opacity -- opacity [cm^2/g] eps_nuc -- nuclear energy release rate, excluding neutrinos [erg/s/g] pp -- pp energy release rate [erg/s/g] cno -- CNO energy release rate [erg/s/g] tri_alfa -- triple-alpha energy release rate [erg/s/g] eps_nuc_neu_total -- energy loss rate as nuclear neutrinos [erg/s/g] non_nuc_neu -- energy loss rate as non-nuclear neutrinos [erg/s/g] eps_grav -- thermal energy release rate [erg/s/g] h1 -- 1H mass fraction he3 -- 3He mass fraction he4 -- 4He mass fraction c12 -- 12C mass fraction n14 -- 14N mass fraction o16 -- 16O mass fraction ne20 -- 20Ne mass fraction mg24 -- 24Mg mass fraction si28 -- 28Si mass fraction s32 -- 32S mass fraction ar36 -- 36Ar mass fraction ca40 -- 40Ca mass fraction ti44 -- 44Ti mass fraction cr48 -- 48Cr mass fraction fe52 -- 52Fe mass fraction fe54 -- 54Fe mass fraction fe56 -- 56Fe mass fraction ni56 -- 56Ni mass fraction eta -- electron chemical potential [kB*T] log_omega -- log10(rotation angular velocity [rad/s]) v_rot -- rotation velocity [km/s] j_rot -- specific angular momentum [cm^2/s] dynamo_log_B_r -- log10(dynamo-generated radial field strength [gauss]) dynamo_log_B_phi -- log10(dynamo-generated azimuthal field strength [gauss]) log_D_conv -- log10(convective diffusivity [cm^2/s]) log_D_semi -- log10(semiconvective diffusivity [cm^2/s]) log_D_ovr -- log10(overshoot diffusivity [cm^2/s]) log_D_thrm -- log10(thermohaline diffusivity [cm^2/s]) """ return read_data(filename, rev=True) # Read MESA data def read_data (filename, rev) : # Open the file file = open(filename, 'r') # Read header data file.readline() header_names = file.readline().split() header_values = [num(value) for value in file.readline().split()] data = dict(zip(header_names, header_values)) # Read array data file.readline() file.readline() column_names = file.readline().split() lines_values = [] for column_name in column_names: if column_name in data: raise Exception('Duplicate field name: {s}'.format(column_name)) while True : line = file.readline() if not line : break lines_values.append([num(value) for value in line.split()]) file.close() # Set up data arrays if (rev) : for i in range(0,len(column_names)) : data[column_names[i]] = np.array([line_values[i] for line_values in reversed(lines_values)]) else : for i in range(0,len(column_names)) : data[column_names[i]] = np.array([line_values[i] for line_values in lines_values]) # Return the data return data