Line emissivity from ATOMDB by python

#!/usr/bin/env python


import pyatomdb
import sys

# He-like Fe
Z = 26
z1 = 25
upind = 7
loind = 1
linefile=""

emis = pyatomdb.atomdb.get_line_emissivity(Z,z1, upind, loind)

for t, e, dens in zip(emis['kT'], emis['epsilon'], emis['dens']):
    print t,dens, e