def model_klf():
# Read in Geneva tracks
genevaFile = '/u/jlu/work/models/geneva/iso/020/c/'
genevaFile += 'iso_c020_0675.UBVRIJHKLM'
model = asciidata.open(genevaFile)
modMass = model[1].tonumpy()
modV = model[6].tonumpy()
modVK = model[11].tonumpy()
modHK = model[15].tonumpy()
modJLp = model[19].tonumpy()
modJK = model[17].tonumpy()
# genevaFile2 = '/u/jlu/work/models/geneva/iso/020/c/'
# genevaFile2 += 'iso_c020_068.UBVRIJHKLM'
# model = asciidata.open(genevaFile)
# modMass = model[1].tonumpy()
# modV = model[6].tonumpy()
# modVK = model[11].tonumpy()
# modHK = model[15].tonumpy()
# modJLp = model[19].tonumpy()
# modJK = model[17].tonumpy()
# Reddening
aV = 27.0
RV = 2.9
# # cardelli() returns A_L
# aJ = aV * extinction.cardelli(1.248, RV)
# aH = aV * extinction.cardelli(1.6330, RV)
# aKp = aV * extinction.cardelli(2.1245, RV)
# aK = aV * extinction.cardelli(2.196, RV)
# aKs = aV * extinction.cardelli(2.146, RV)
aKs = 2.7
aJ = extinction.nishiyama09(1.248, aKs)
aH = extinction.nishiyama09(1.6330, aKs)
aKp = extinction.nishiyama09(2.1245, aKs)
aK = extinction.nishiyama09(2.196, aKs)
aKs = extinction.nishiyama09(2.146, aKs)
modK = modV - modVK
modH = modK + modHK
modJ = modK + modJK
modLp = modJ - modJLp
modKs = modK + 0.002 + 0.026 * (modJK)
modKp = modK + 0.22 * (modHK)
dist = 8400.0
distMod = -5.0 + 5.0 * math.log10(dist)
modK_extinct = modK + aK + distMod
modKp_extinct = modKp + aKp + distMod
modKs_extinct = modKs + aKs + distMod
return modKp_extinct
def model_klf():
# Read in Geneva tracks
genevaFile = '/u/jlu/work/models/geneva/iso/020/c/'
genevaFile += 'iso_c020_0675.UBVRIJHKLM'
model = asciidata.open(genevaFile)
modMass = model[1].tonumpy()
modV = model[6].tonumpy()
modVK = model[11].tonumpy()
modHK = model[15].tonumpy()
modJLp = model[19].tonumpy()
modJK = model[17].tonumpy()
# genevaFile2 = '/u/jlu/work/models/geneva/iso/020/c/'
# genevaFile2 += 'iso_c020_068.UBVRIJHKLM'
# model = asciidata.open(genevaFile)
# modMass = model[1].tonumpy()
# modV = model[6].tonumpy()
# modVK = model[11].tonumpy()
# modHK = model[15].tonumpy()
# modJLp = model[19].tonumpy()
# modJK = model[17].tonumpy()
# Reddening
aV = 27.0
RV = 2.9
# # cardelli() returns A_L
# aJ = aV * extinction.cardelli(1.248, RV)
# aH = aV * extinction.cardelli(1.6330, RV)
# aKp = aV * extinction.cardelli(2.1245, RV)
# aK = aV * extinction.cardelli(2.196, RV)
# aKs = aV * extinction.cardelli(2.146, RV)
aKs = 2.7
aJ = extinction.nishiyama09(1.248, aKs)
aH = extinction.nishiyama09(1.6330, aKs)
aKp = extinction.nishiyama09(2.1245, aKs)
aK = extinction.nishiyama09(2.196, aKs)
aKs = extinction.nishiyama09(2.146, aKs)
modK = modV - modVK
modH = modK + modHK
modJ = modK + modJK
modLp = modJ - modJLp
modKs = modK + 0.002 + 0.026 * (modJK)
modKp = modK + 0.22 * (modHK)
dist = 8400.0
distMod = -5.0 + 5.0 * math.log10(dist)
modK_extinct = modK + aK + distMod
modKp_extinct = modKp + aKp + distMod
modKs_extinct = modKs + aKs + distMod
return modKp_extinct
def __init__(self, k=3):
# Fetch the extinction curve, pre-interpolate across 1-8 microns
wave = np.arange(0.5, 8.0, 0.001)
# This will eventually be scaled by AKs when you
# call reddening(). Right now, calc for AKs=1
Alambda_scaled = extinction.nishiyama09(wave, 1.0, makePlot=False, k=k)
# Convert wavelength to angstrom
wave *= 10**4
pysynphot.reddening.CustomRedLaw.__init__(self, wave=wave,
waveunits='angstrom',
Avscaled=Alambda_scaled,
name='Nishiyama09',
litref='Nishiyama+ 2009')
def __init__(self):
# Fetch the extinction curve, pre-interpolate across 1-8 microns
wave = np.arange(1.0, 8.0, 0.01)
# This will eventually be scaled by AKs when you
# call reddening(). Right now, calc for AKs=1
Alambda_scaled = extinction.nishiyama09(wave, 1.0, makePlot=False)
# Convert wavelength to angstrom
wave *= 10**4
pysynphot.reddening.CustomRedLaw.__init__(self,
wave=wave,
waveunits='angstrom',
Avscaled=Alambda_scaled,
name='Nishiyama09',
litref='Nishiyama+ 2009')
