def test_morton03(): m03 = parse.parse_morton03(orig=True) # np.testing.assert_allclose(m03['wrest'][5], 930.7482, rtol=1e-7) assert m03['wrest'].unit == u.Angstrom # m03 = parse.mktab_morton03() #m03 = parse.mktab_morton03(do_this=True, outfil='tmp.fits') m03 = parse.mktab_morton03(do_this=True, fits=False, outfil='tmp.vo')
import imp jbvp_path = imp.find_module('joebvp')[1] vernerlist=np.genfromtxt(jbvp_path+'/atomicdata/verner6.txt',dtype=None,delimiter=[10,8,3,4,3,2,9,6]) vernlam=jbg.arrfromcol(vernerlist,0) vernion=jbg.arrfromcol(vernerlist,1) vernzatom=jbg.arrfromcol(vernerlist,2) vernnume=jbg.arrfromcol(vernerlist,3) verngl=jbg.arrfromcol(vernerlist,4) verngu=jbg.arrfromcol(vernerlist,5) vernosc=jbg.arrfromcol(vernerlist,6) vernp=jbg.arrfromcol(vernerlist,7) # A start to using the linetools atomic data framework adata=lilp.parse_morton03() vdata=lilp.parse_verner96() for i in range(len(vernion)): vernion[i]=vernion[i].strip() def setatomicdata(lines,precise=True): #if isinstance(lines,float): # lines = [lines] lam=np.zeros(len(lines)) ; fosc=np.zeros(len(lines)) ; gam=np.zeros(len(lines)) for i,ll in enumerate(lines): try: al=AbsLine(ll*u.AA,closest=True) lam[i]=al.data['wrest'].value ; fosc[i]=al.data['f'] ; gam[i]=al.data['gamma'].value except: idx=jbg.closest(adata['wrest'],ll)
jbvp_path = imp.find_module('joebvp')[1] vernerlist = np.genfromtxt(jbvp_path + '/atomicdata/verner6.txt', dtype=None, delimiter=[10, 8, 3, 4, 3, 2, 9, 6]) vernlam = jbg.arrfromcol(vernerlist, 0) vernion = jbg.arrfromcol(vernerlist, 1) vernzatom = jbg.arrfromcol(vernerlist, 2) vernnume = jbg.arrfromcol(vernerlist, 3) verngl = jbg.arrfromcol(vernerlist, 4) verngu = jbg.arrfromcol(vernerlist, 5) vernosc = jbg.arrfromcol(vernerlist, 6) vernp = jbg.arrfromcol(vernerlist, 7) # A start to using the linetools atomic data framework adata = lilp.parse_morton03() vdata = lilp.parse_verner96() for i in range(len(vernion)): vernion[i] = vernion[i].strip() def setatomicdata(lines, precise=True, linelist=ilist): lam = np.zeros(len(lines)) fosc = np.zeros(len(lines)) gam = np.zeros(len(lines)) for i, ll in enumerate(lines): try: al = AbsLine(ll * u.AA, closest=True, linelist=linelist) lam[i] = al.data['wrest'].value fosc[i] = al.data['f']
def test_morton03(): m03 = parse.parse_morton03(orig=True) # np.testing.assert_allclose(m03['wrest'][5], 930.7482, rtol=1e-7) assert m03['wrest'].unit == u.Angstrom