def map_dbt_lightcone(): for i in range(0, len(redshifts)): # print 'filename: setup_dirs.resultsdir()+'map_dbt_lightcone_'+str('%.3f' % redshifts[i])+'.dat' filename = setup_dirs.resultsdir() + "map_dbt_lightcone_" + str("%.3f" % redshifts[i]) + ".dat" temp_filename = setup_dirs.path() + "Temper3D_" + str("%.3f" % redshifts[i]) + ".bin" xfrac_filename = setup_dirs.path() + "xfrac3d_" + str("%.3f" % redshifts[i]) + ".bin" if i % 2 == 0: density_filename = ( "/research/prace/244Mpc_RT/244Mpc_f2_8.2pS_250/coarser_densities/" + str("%.3f" % redshifts[i]) + "n_all.dat" ) else: density_filename = ( "/research/prace/244Mpc_RT/244Mpc_f2_8.2pS_250/coarser_densities/" + str("%.3f" % redshifts[i - 1]) + "n_all.dat" ) tfile = c2t.TemperFile(temp_filename) xfile = c2t.XfracFile(xfrac_filename).xi if ss == " ": dfile = c2t.DensityFile(density_filename).cgs_density else: dfile = np.ones(ss ** 3).reshape(ss, ss, ss) * 1.981e-10 * (1 + redshifts[i]) ** 3 dT_box = c2t.calc_dt_full_lightcone(xfile, tfile, dfile, redshifts[i]) IO.writemap(dT_box, filename) print "Written map to " + filename
def map_temp(): for i in range(start, len(redshifts)): filename = setup_dirs.resultsdir() + "map_temper_" + str("%.3f" % redshifts[i]) + ".dat" temp_filename = setup_dirs.path() + "Temper3D_" + str("%.3f" % redshifts[i]) + ".bin" tfile = c2t.TemperFile(temp_filename) IO.writemap(tfile.temper, filename) print "Writen map to " + filename
def map_xfrac(id): # for i in range(len(redshifts)-6,len(redshifts)-4): print len(redshifts) for i in range(start, len(redshifts)): filename = setup_dirs.resultsdir() + "map_xfrac" + id + "_" + str("%.3f" % redshifts[i]) + ".dat" xfrac_filename = setup_dirs.path() + "xfrac3d" + id + "_" + str("%.3f" % redshifts[i]) + ".bin" xfile = c2t.XfracFile(xfrac_filename) IO.writemap(xfile.xi, filename) print "Written map to " + filename