diff options
| -rwxr-xr-x | astro/21cm/tile_slice.py | 19 |
1 files changed, 10 insertions, 9 deletions
diff --git a/astro/21cm/tile_slice.py b/astro/21cm/tile_slice.py index 4f217e8..cd66896 100755 --- a/astro/21cm/tile_slice.py +++ b/astro/21cm/tile_slice.py @@ -33,20 +33,21 @@ def main(): parser.add_argument("-C", "--clobber", dest="clobber", action="store_true", help="overwrite existing files") - parser.add_argument("-F", "--fov", dest="fov", default=5.0, type=float, + parser.add_argument("-F", "--fov", dest="fov", default=10.0, type=float, help="required FoV [deg] of the output slice " + - "(default: 5.0 deg)") - parser.add_argument("-N", "--n-side", dest="Nside", default=600, type=int, + "(default: 10.0 [deg])") + parser.add_argument("-N", "--n-side", dest="Nside", default=1800, type=int, help="required image size of output slice " + - "(default: 600)") + "(default: 1800)") parser.add_argument("-i", "--infile", dest="infile", required=True, help="input slice") parser.add_argument("-o", "--outfile", dest="outfile", default=outfile_default, help="output tiled slice (default: %s)" % outfile_default) - parser.add_argument("-p", "--prefix", dest="prefix", required=True, - help="prefix for the output tiled slice") + parser.add_argument("-p", "--prefix", dest="prefix", default="deltaTb", + help="prefix for the output tiled slice " + + "(default: 'deltaTb')") exgrp = parser.add_mutually_exclusive_group(required=True) exgrp.add_argument("-z", "--redshift-c", dest="zc", type=float, help="central redshift of the selected cube") @@ -80,9 +81,9 @@ def main(): img_out = scipy.ndimage.zoom(img2, zoom=args.Nside/Nside2, order=1) # Record information to header header["Z_C"] = (zc, "Central redshift") - header["FREQ_C"] = (fc, "Frequency [MHz] w.r.t. to the redshift") - header["FOV"] = (args.fov, "FoV [deg] of this slice") - header["PixSize"] = (60.0*args.fov/args.Nside, "Pixel size [arcmin]") + header["FREQ_C"] = (fc, "[MHz] Frequency w.r.t. to central redshift") + header["FoV"] = (args.fov, "[deg] FoV of this slice") + header["PixSize"] = (3600.0*args.fov/args.Nside, "[arcsec] Pixel size") header.add_history(" ".join(sys.argv)) hdu = fits.PrimaryHDU(data=img_out, header=header) outfile = args.outfile.format(prefix=args.prefix, freq=freq, |