def make_iso_table(
project,
oiso=None,
logtmin=6.0,
logtmax=10.13,
dlogt=0.05,
z=[0.0152],
iso_fname=None,
):
"""
The isochrone tables are loaded (downloading if necessary)
Parameters
----------
project: str
project name
oiso: isochrone.Isochrone object
contains the full isochrones information
logtmin: float
log-age min
logtmax: float
log-age max
dlogt: float
log-age step to request
z: float or sequence
list of metalicity values, where default (Z=0.152) is adopted Z_sun
for PARSEC/COLIBRI models
Returns
-------
fname: str
name of saved file
oiso: isochrone.Isochrone object
contains the full isochrones information
"""
if iso_fname is None:
iso_fname = "%s/%s_iso.csv" % (project, project)
if not os.path.isfile(iso_fname):
if oiso is None:
oiso = isochrone.PadovaWeb()
t = oiso._get_t_isochrones(max(5.0, logtmin), min(10.13, logtmax),
dlogt, z)
t.header["NAME"] = "{0} Isochrones".format("_".join(
iso_fname.split("_")[:-1]))
print("{0} Isochrones".format("_".join(iso_fname.split("_")[:-1])))
t.write(iso_fname)
# read in the isochrone data from the file
# not sure why this is needed, but reproduces previous ezpipe method
oiso = ezIsoch(iso_fname)
return (iso_fname, oiso)
# Metallicity : list of floats # Here: Z == Z_initial, NOT Z(t) surface abundance # PARSECv1.2S accepts values 1.e-4 < Z < 0.06 # example z = [0.03, 0.019, 0.008, 0.004] # can they be set as [min, max, step]? z = [0.03, 0.019, 0.008, 0.004] # Isochrone Model Grid # Current Choices: Padova or MIST # PadovaWeb() -- `modeltype` param for iso sets from ezpadova # (choices: parsec12s_r14, parsec12s, 2010, 2008, 2002) # MISTWeb() -- `rotation` param (choices: vvcrit0.0=default, vvcrit0.4) # # Default: PARSEC+COLIBRI oiso = isochrone.PadovaWeb() # Alternative: PARSEC1.2S -- old grid parameters #oiso = isochrone.PadovaWeb(modeltype='parsec12s', filterPMS=True) # Alternative: MIST -- v1, no rotation #oiso = isochrone.MISTWeb() # Stellar Atmospheres library definition osl = stellib.Tlusty() + stellib.Kurucz() ################ ### Dust extinction grid definition extLaw = extinction.Gordon16_RvFALaw() # A(V): dust column in magnitudes # acceptable avs > 0.0
# Here: Z == Z_initial, NOT Z(t) surface abundance # PARSECv1.2S accepts values 1.e-4 < Z < 0.06 # example z = [0.03, 0.019, 0.008, 0.004] # can they be set as [min, max, step]? z = [0.03, 0.019, 0.008, 0.004] # Isochrone Model Grid # Current Choices: Padova or MIST # PadovaWeb() -- `modeltype` param for iso sets from ezpadova # (choices: parsec12s_r14, parsec12s, 2010, 2008, 2002) # MISTWeb() -- `rotation` param (choices: vvcrit0.0=default, vvcrit0.4) # # Default: PARSEC+CALIBRI #oiso = isochrone.PadovaWeb() # Alternative: PARSEC1.2S -- old grid parameters oiso = isochrone.PadovaWeb(modeltype='parsec12s', filterPMS=True) # Alternative: MIST -- v1, no rotation #oiso = isochrone.MISTWeb() # Stellar Atmospheres library definition osl = stellib.Tlusty() + stellib.Kurucz() ################ ### Dust extinction grid definition extLaw = extinction.Gordon16_RvFALaw() # A(V): dust column in magnitudes # acceptable avs > 0.0 # example [min, max, step] = [0.0, 10.055, 1.0] avs = [0.0, 10.055, 1.0]
