import pylab as plt
from atpy import Table
prefix = "/Users/ncrighton/Projects/MPIA_QSO_LBG/Cloudy/J0004/comp1"
name, = glob(prefix + "/uvb_k04/*tilted*tbl")
fig1 = pl.figure(figsize=(7.1,4.7))
fig1.subplots_adjust(left=0.15, bottom=0.15, top=0.96,right=0.96)
ax = fig1.add_subplot(111)
N0 = -21
ref = Table(name)
ind = ref.energy.searchsorted(0.99)
ref.log10jnu = ref.log10jnu - ref.log10jnu[ind] + N0
ax.plot(np.log10(ref.energy), ref.log10jnu - ref.log10jnu[ind] + N0, '-', color='k')
xvals = np.log10((IP * u.eV).to(u.Ry).value)
y = np.interp(xvals, np.log10(ref.energy), ref.log10jnu)
from barak.absorb import split_trans_name
xdiff = np.diff(xvals)
direction = ['up']
previousdown = False
for i in xrange(len(xvals) - 1):
if xdiff[i] < 0.02 and not previousdown:
import pylab as plt
from atpy import Table
prefix = "/Users/ncrighton/Projects/MPIA_QSO_LBG/Cloudy/J0004/comp1"
name, = glob(prefix + "/uvb_k04/*tilted*tbl")
fig1 = pl.figure(figsize=(7.1, 4.7))
fig1.subplots_adjust(left=0.15, bottom=0.15, top=0.96, right=0.96)
ax = fig1.add_subplot(111)
N0 = -21
ref = Table(name)
ind = ref.energy.searchsorted(0.99)
ref.log10jnu = ref.log10jnu - ref.log10jnu[ind] + N0
ax.plot(np.log10(ref.energy), ref.log10jnu - ref.log10jnu[ind] + N0, "-", color="k")
xvals = np.log10((IP * u.eV).to(u.Ry).value)
y = np.interp(xvals, np.log10(ref.energy), ref.log10jnu)
from barak.absorb import split_trans_name
xdiff = np.diff(xvals)
direction = ["up"]
previousdown = False
for i in xrange(len(xvals) - 1):
if xdiff[i] < 0.02 and not previousdown:
