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: