def ps():
for i in range(len(redshifts)):
data = IO.readmap("dbt_" + str("%.3f" % redshifts[i]))
ps = power_spectrum_1d(data)
ps = np.asarray(ps)
print np.shape(ps), type(ps)
IO.writedata(ps[0, :], "data/powerSpectra_dbt_" + str("%.3f" % redshifts[i]) + ".dat")
IO.writedata(ps[1, :], "data/powerSpectraFrequencies_dbt_" + str("%.3f" % redshifts[i]) + ".dat")
