}
parFormat = ["age"]
readMe = '''
The stellar emission templates are generated with EzGal.
'''
templateDict = {
"tckList": tckList,
"kdTree": kdt,
"parList": XList,
"modelInfo": modelInfo,
"parFormat": parFormat,
"readMe": readMe
}
t = Template(tckList=tckList,
kdTree=kdt,
parList=XList,
modelInfo=modelInfo,
parFormat=parFormat,
readMe=readMe)
t = Template(**templateDict)
fp = open("bc03_kdt.tmplt", "w")
pickle.dump(templateDict, fp)
fp.close()
#Test
cl = ['r', 'g', 'b', 'y', 'c']
for loop in range(nAge):
pars = [ageList[loop]]
f = t(wave, pars)
flux = fluxList[loop]
plt.plot(wave, f, color=cl[loop], linestyle="--")
plt.plot(wave, flux, color=cl[loop], linestyle=":")
readMe = '''
This template is from: http://www.pa.uky.edu/clumpy/
The interpolation is tested well!
'''
templateDict = {
"tckList": tckList,
"kdTree": kdt,
"parList": XList,
"modelInfo": modelInfo,
"parFormat": parFormat,
"readMe": readMe
}
print("haha")
t = Template(tckList=tckList,
kdTree=kdt,
parList=XList,
modelInfo=modelInfo,
parFormat=parFormat,
readMe=readMe)
print("haha")
t = Template(**templateDict)
print("haha")
fp = open("clumpy_kdt.tmplt", "w")
#pickle.dump(t, fp)
pickle.dump(templateDict, fp)
fp.close()
if f_test:
fp = open("clumpy_kdt.tmplt", "r")
tpDict = pickle.load(fp)
fp.close()
"Silicate": silDict,
"Graphite": graDict,
"readMe": templateReadMe
}
fp = open("/Users/jinyi/Work/mcmc/Fitter/template/dust_xl_kdt.tmplt", "w")
pickle.dump(dustModel, fp)
fp.close()
#"""
##Test the KDTree and the interpolation
fp = open("/Users/jinyi/Work/mcmc/Fitter/template/dust_xl_kdt.tmplt", "r")
dustModel = pickle.load(fp)
fp.close()
#For the astronomical silicate
silDict = dustModel["Silicate"]
t = Template(**silDict)
colorList = ["r", "b", "m", "g", "y", "c"]
for nsil in silTList:
tmplDir = modelDir + "{0}.dat".format(silNameList[nsil])
modelSil = np.loadtxt(tmplDir)
size_all = modelSil[:, 0]
wave_all = modelSil[:, 1]
kappa_all = modelSil[:, 2]
sizeList = np.unique(size_all)
for nsz in range(len(sizeList)): # [0, 79, 99]:
sz = sizeList[nsz]
fltr = size_all == sz
wave = wave_all[fltr]
kappa = kappa_all[fltr]
par = [nsil, sz-0.01]
kappa_t = t(wave, par)