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')