From c73b3b3db67dc5bd59d4b2925e81e3da31354d24 Mon Sep 17 00:00:00 2001
From: Aaron LI <aly@aaronly.me>
Date: Mon, 11 Sep 2017 22:36:05 +0800
Subject: wsclean.py: Add more arguments; as well as add short version

---
 astro/oskar/wsclean.py | 49 +++++++++++++++++++++++++++++++++++--------------
 1 file changed, 35 insertions(+), 14 deletions(-)

diff --git a/astro/oskar/wsclean.py b/astro/oskar/wsclean.py
index 574c1dc..cfef30d 100755
--- a/astro/oskar/wsclean.py
+++ b/astro/oskar/wsclean.py
@@ -44,30 +44,41 @@ def main():
                         help="only create dirty images (by setting niter=0)")
     parser.add_argument("--update-model", dest="update_model",
                         action="store_true",
-                        help="update the MODEL_DATA column in MS")
+                        help="write/update the MODEL_DATA column in MS")
     parser.add_argument("--save-weights", dest="save_weights",
                         action="store_true",
-                        help="save the gridded weights in <name>-weights.fits")
+                        help="save gridded weights in <name>-weights.fits")
+    parser.add_argument("--save-uv", dest="save_uv",
+                        action="store_true",
+                        help="save gridded uv in <name>-uv-{real,imag}.fits")
+    parser.add_argument("--uv-range", dest="uv_range", default=":",
+                        help="uv range [lambda] (i.e., baseline lengths) " +
+                        "used for imaging; syntax: '<min>:<max>' " +
+                        "(default: ':', i.e., all uv/baselines)")
     parser.add_argument("-w", "--weight", dest="weight", default="briggs",
                         choices=["uniform", "natural", "briggs"],
                         help="weighting method (default: 'briggs')")
-    parser.add_argument("--briggs", dest="briggs", type=float, default=0.0,
+    parser.add_argument("-B", "--briggs", dest="briggs",
+                        type=float, default=0.0,
                         help="Briggs robustness parameter (default: 0); " +
                         "-1 (uniform) -> 1 (natural)")
-    parser.add_argument("--niter", dest="niter", type=int, default=100000,
+    parser.add_argument("-#", "--niter", dest="niter",
+                        type=int, default=100000,
                         help="maximum number of CLEAN iterations " +
                         "(default: 100,000)")
     parser.add_argument("--gain", dest="gain", type=float, default=0.1,
-                        help="CLEAN gain for each minor iteration")
+                        help="CLEAN gain for each minor iteration " +
+                        "(default: 0.1)")
     parser.add_argument("--mgain", dest="mgain", type=float, default=0.85,
-                        help="CLEAN gain for major iterations")
+                        help="CLEAN gain for major iterations " +
+                        "(default: 0.85)")
     parser.add_argument("-s", "--size", dest="size", type=int,
                         required=True,
                         help="output image size (pixel number on a side)")
     parser.add_argument("-p", "--pixelsize", dest="pixelsize", type=float,
                         required=True,
                         help="output image pixel size [arcsec]")
-    parser.add_argument("--taper-gaus", dest="taper_gaus", type=float,
+    parser.add_argument("-G", "--taper-gaus", dest="taper_gaus", type=float,
                         help="taper the weights with a Gaussian function " +
                         "to reduce the contribution of long baselines. " +
                         "Gaussian beam size in [arcsec].")
@@ -76,11 +87,12 @@ def main():
                         "spectra with [order] polynomial in normal-space")
     #
     exgrp = parser.add_mutually_exclusive_group()
-    exgrp.add_argument("--threshold-auto", dest="threshold_auto",
-                       type=float, default=1.5,
-                       help="estimate noise level and stop at <sigma>*<std>")
-    exgrp.add_argument("--threshold", dest="threshold", type=float,
-                       help="stopping CLEAN threshold [Jy]")
+    exgrp.add_argument("-S", "--threshold-nsigma", dest="threshold_nsigma",
+                       type=float, default=2.0,
+                       help="estimate the noise level <sigma> and stop at " +
+                       "nsigma*<sigma> (default: 2.0 <sigma>)")
+    exgrp.add_argument("-t", "--threshold", dest="threshold", type=float,
+                       help="stopping CLEAN threshold [mJy]")
     #
     parser.add_argument("-N", "--name", dest="name", required=True,
                         help="filename prefix for the output files")
@@ -125,11 +137,20 @@ def main():
 
     if args.save_weights:
         cmdargs += ["-saveweights"]
+    if args.save_uv:
+        cmdargs += ["-saveuv"]
+
+    # uv/baseline range
+    uvmin, uvmax = args.uv_range.strip().split(":")
+    if uvmin:
+        cmdargs += ["-minuv-l", float(uvmin)]
+    if uvmax:
+        cmdargs += ["-maxuv-l", float(uvmax)]
 
     if args.threshold:
-        cmdargs += ["-threshold", args.threshold]
+        cmdargs += ["-threshold", args.threshold*1e-3]  # [mJy] -> [Jy]
     else:
-        cmdargs += ["-auto-threshold", args.threshold_auto]
+        cmdargs += ["-auto-threshold", args.threshold_nsigma]
 
     if args.taper_gaus:
         cmdargs += ["-taper-gaussian", args.taper_gaus]
-- 
cgit v1.2.2