ps1d_eorwindow.py 文件源码

def save(self, outfile):
        if self.ps1d_normalized:
            ps1d_desc = "normalized power [K^2]"
        else:
            ps1d_desc = "power [K^2 Mpc^3]"
        header = [
            "EoR window definition:",
            "+ FoV: %f [deg]" % self.ps2d.fov,
            "+ e_ConvWidth: %f" % self.ps2d.e,
            "+ k_perp_min: %f [Mpc^-1]" % self.ps2d.k_perp_min,
            "+ k_perp_max: %f [Mpc^-1]" % self.ps2d.k_perp_max,
            "+ k_los_min: %f [Mpc^-1]" % self.ps2d.k_los_min,
            "+ k_los_max: %f [Mpc^-1]" % self.ps2d.k_los_max,
            "",
            "Columns:",
            "1. k: wavenumber [Mpc^-1]",
            "2. ps1d: %s" % ps1d_desc,
            "ps1d_err: power errors",
            "",
            "k   ps1d   ps1d_err",
        ]
        np.savetxt(outfile, self.ps1d, header="\n".join(header))
        print("Saved 1D power spectrum to file: %s" % outfile)