def header(self):
dDc = self.configs.Dc_cell
Dc_min, Dc_max = self.configs.Dc_limit
header = fits.Header()
header["BUNIT"] = (self.configs.unit, "Data unit")
header["zmin"] = (self.configs.zmin, "HI simulation minimum redshift")
header["zmax"] = (self.configs.zmax, "HI simulation maximum redshift")
header["dz"] = (self.configs.dz, "HI simulation redshift step size")
header["Dc_min"] = (Dc_min, "[cMpc] comoving distance at zmin")
header["Dc_max"] = (Dc_max, "[cMpc] comoving distance at zmax")
header["Dc_step"] = (dDc, "[cMpc] comoving distance between slices")
header["Lside"] = (self.configs.Lside, "[cMpc] Simulation side length")
header["Nside"] = (self.configs.Nside, "Number of cells at each side")
header["DATE"] = (datetime.now(timezone.utc).astimezone().isoformat(),
"File creation date")
header.add_history(" ".join(sys.argv))
header.extend(self.wcs.to_header(), update=True)
return header
评论列表
文章目录