unlogg_data = True
tag = ftag
try:
fit_region_R = cf.get_MESA_profile_edge(MESA_file,
quantity='logR',
masscut=masscut,
strip=False)
fit_region_M = cf.get_MESA_profile_edge(MESA_file,
quantity='mass',
masscut=masscut,
strip=False)
if unlogg_data:
fit_region_R = cf.unlog(fit_region_R) * R_to_solar
fit_region_M = fit_region_M * M_to_solar
plt.plot(fit_region_R, fit_region_M, 'm-', color='darkblue', label=tag)
if unlogg_data:
plt.xlabel('R (cm)')
plt.ylabel('M (g)')
else:
plt.xlabel('logR (R/Rsolar)')
plt.ylabel('M (Msolar))')
plt.title('data from profile_' + ftag + ' at ' +
str((1.0 - masscut) * 100.0) + r'% surface depth')
plt.legend(loc=1)
if savefig:
plt.savefig(tag + '_mass_profile_' + str(masscut) + '.png')
plt.show()
def MESA_E(MESA_file, masscut):
## cgs
fit_region_E=cf.get_MESA_profile_edge(MESA_file, quantity='logE', masscut=masscut, strip=False)
fit_region_E=cf.unlog(fit_region_E)
return fit_region_E
def MESA_rho(MESA_file,masscut):
## returns rho in cgs
fit_region_rho = cf.get_MESA_profile_edge(MESA_file, quantity='logRho', masscut=masscut ,strip=False)
fit_region_rho=cf.unlog(fit_region_rho)
return fit_region_rho
def MESA_P(MESA_file, masscut):
## returns pressure in cgs
fit_region_P=cf.get_MESA_profile_edge(MESA_file, quantity='logP', masscut=masscut, strip=False)
fit_region_P=cf.unlog(fit_region_P)
return fit_region_P
def MESA_r(MESA_file, masscut):
## returns R in cm
fit_region_R = cf.get_MESA_profile_edge(MESA_file, quantity='logR', masscut=masscut,strip=False)
fit_region_R=cf.unlog(fit_region_R)*R_to_solar
return fit_region_R