try:
# Declare the script parameters. Their positions in this list define
# their expected position on the script command line. They can also be
# specified by keyword on the command line. No validation of default
# values or values supplied on the command line is performed until the
# parameter value is first accessed within the script, at which time the
# user is prompted for a value if necessary. The parameters "MSG_FILTER",
# "ILEVEL", "GLEVEL" and "LOGFILE" are added automatically by the ParSys
# constructor.
params = []
params.append(
starutil.ParNDG(
"IN", "Template POL2 time series NDFs",
starutil.get_task_par("DATA_ARRAY",
"GLOBAL",
default=Parameter.UNSET)))
params.append(starutil.Par0S("OUT", "Output simulated POL2 data"))
params.append(
starutil.Par0L("NEWART", "Create new artificial I, Q and U maps?"))
params.append(
starutil.ParChoice("ARTFORM", ("0", "1"),
"Form of artificial I, Q and U to create",
"0",
noprompt=True))
params.append(starutil.ParNDG("ARTI", "Artificial I map", maxsize=1))
params.append(starutil.ParNDG("ARTQ", "Artificial Q map", maxsize=1))
params.append(starutil.ParNDG("ARTU", "Artificial U map", maxsize=1))
params.append(
starutil.ParNDG("INCOM",
try:
# Declare the script parameters. Their positions in this list define
# their expected position on the script command line. They can also be
# specified by keyword on the command line. No validation of default
# values or values supplied on the command line is performed until the
# parameter value is first accessed within the script, at which time the
# user is prompted for a value if necessary. The parameters "MSG_FILTER",
# "ILEVEL", "GLEVEL" and "LOGFILE" are added automatically by the ParSys
# constructor.
params = []
params.append(
starutil.ParNDG(
"IN", "The input NDFs",
starutil.get_task_par("DATA_ARRAY",
"GLOBAL",
default=Parameter.UNSET)))
params.append(
starutil.ParChoice("INSTRUMENT",
["SCUBA-2(450)", "SCUBA-2(850)", "ACSIS", "DAS"],
"The JCMT instrument", "SCUBA-2(850)"))
params.append(
starutil.Par0L("JSA", "Are the input NDFs on the JSA "
"all-sky pixel grid?",
True,
noprompt=True))
params.append(
# Catch any exception so that we can always clean up, even if control-C
# is pressed.
try:
# Declare the script parameters. Their positions in this list define
# their expected position on the script command line. They can also be
# specified by keyword on the command line. No validation of default
# values or values supplied on the command line is performed until the
# parameter value is first accessed within the script, at which time the
# user is prompted for a value if necessary. The parameters "MSG_FILTER",
# "ILEVEL", "GLEVEL" and "LOGFILE" are added automatically by the ParSys
# constructor.
params = []
params.append(starutil.ParNDG("INQ", "The input Q images"))
params.append(starutil.ParNDG("INU", "The input U images"))
params.append(starutil.ParNDG("INI", "The input I images"))
params.append(starutil.Par0S("CAT", "The output FITS vector catalogue",
"out.FIT"))
params.append(starutil.ParNDG("PI", "The output polarised intensity map",
default=None, exists=False, minsize=0, maxsize=1 ))
params.append(starutil.Par0L("RETAIN", "Retain temporary files?", False,
noprompt=True))
params.append(starutil.Par0S("QUI", "An HDS container file in which to "
"store the 2D Q, U and I images",
default=None ))
global retain
ParSys.cleanup()
if retain:
msg_out("Retaining temporary files in {0}".format(NDG.tempdir))
else:
NDG.cleanup()
# Catch any exception so that we can always clean up, even if control-C
# is pressed.
for zz in range(1):
params = []
params.append(
starutil.ParNDG(
"IN",
"The input POL2 time series NDFs (with the dome open and calibrator out)",
starutil.get_task_par("DATA_ARRAY",
"GLOBAL",
default=Parameter.UNSET)))
params.append(
starutil.Par0S("OUT",
"The output .sdf file containing the IPT parameters",
"ipdata"))
params.append(
starutil.ParNDG(
"DomeClosedCalIn",
"The input POL2 time series NDFs (with the dome closed and calibrator in)",
starutil.get_task_par("DATA_ARRAY",
"GLOBAL",
default=Parameter.UNSET)))
params.append(
starutil.ParNDG(
# specified by keyword on the command line. No validation of default
# values or values supplied on the command line is performed until the
# parameter value is first accessed within the script, at which time the
# user is prompted for a value if necessary. The parameters "MSG_FILTER",
# "ILEVEL", "GLEVEL" and "LOGFILE" are added automatically by the ParSys
# constructor.
params = []
params.append(
starutil.Par0S("REGION", "The required region (an NDF or "
"AST Region)"))
params.append(
starutil.ParNDG("OUT",
"The output NDF",
default=None,
exists=False,
minsize=0,
maxsize=1))
params.append(
starutil.ParChoice("INSTRUMENT",
["SCUBA-2(450)", "SCUBA-2(850)", "ACSIS", "DAS"],
"The JCMT instrument", "SCUBA-2(850)"))
params.append(
starutil.Par0L("RETAIN",
"Retain temporary files?",
False,
noprompt=True))
# Initialise the parameters to hold any values supplied on the command
try:
# Declare the script parameters. Their positions in this list define
# their expected position on the script command line. They can also be
# specified by keyword on the command line. No validation of default
# values or values supplied on the command line is performed until the
# parameter value is first accessed within the script, at which time the
# user is prompted for a value if necessary. The parameters "MSG_FILTER",
# "ILEVEL", "GLEVEL" and "LOGFILE" are added automatically by the ParSys
# constructor.
params = []
params.append(
starutil.ParNDG("IN",
"The input NDF",
default=None,
minsize=0,
maxsize=1))
params.append(
starutil.Par0S("OUT", "The basename for the output NDFs", None))
params.append(
starutil.ParChoice("INSTRUMENT",
["SCUBA-2(450)", "SCUBA-2(850)", "ACSIS", "DAS"],
"The JCMT instrument", "SCUBA-2(850)"))
params.append(
starutil.Par0I("TRIM", "How to trim the output NDFs", 2,
noprompt=True))
# Catch any exception so that we can always clean up, even if control-C
# is pressed.
try:
# Declare the script parameters. Their positions in this list define
# their expected position on the script command line. They can also be
# specified by keyword on the command line. No validation of default
# values or values supplied on the command line is performed until the
# parameter value is first accessed within the script, at which time the
# user is prompted for a value if necessary. The parameters "MSG_FILTER",
# "ILEVEL", "GLEVEL" and "LOGFILE" are added automatically by the ParSys
# constructor.
params = []
params.append(starutil.ParNDG("TILES", "The available JSA tiles"))
params.append(
starutil.ParNDG("OUT", "The output image", exists=False, maxsize=1))
params.append(
starutil.Par0S("CENTRE1", "Longitude at centre of required "
"circular sky region",
None,
noprompt=True))
params.append(
starutil.Par0S("CENTRE2", "Latitude at centre of required "
"circular sky region"))
params.append(
# Catch any exception so that we can always clean up, even if control-C
# is pressed.
try:
# Declare the script parameters. Their positions in this list define
# their expected position on the script command line. They can also be
# specified by keyword on the command line. No validation of default
# values or values supplied on the command line is performed until the
# parameter value is first accessed within the script, at which time the
# user is prompted for a value if necessary. The parameters "MSG_FILTER",
# "ILEVEL", "GLEVEL" and "LOGFILE" are added automatically by the ParSys
# constructor.
params = []
params.append(starutil.ParNDG("IN", "The input NDF"))
params.append(starutil.Par0S("REGION", "Output text file", default=None))
params.append(
starutil.Par0L("TILES", "Display JSA tiles?", False, noprompt=True))
params.append(
starutil.ParChoice("INSTRUMENT",
["SCUBA-2(450)", "SCUBA-2(850)", "ACSIS", "DAS"],
"The JCMT instrument", "SCUBA-2(850)"))
params.append(
starutil.Par0L("RETAIN",
"Retain temporary files?",
False,
# Catch any exception so that we can always clean up, even if control-C
# is pressed.
try:
# Declare the script parameters. Their positions in this list define
# their expected position on the script command line. They can also be
# specified by keyword on the command line. No validation of default
# values or values supplied on the command line is performed until the
# parameter value is first accessed within the script, at which time the
# user is prompted for a value if necessary. The parameters "MSG_FILTER",
# "ILEVEL", "GLEVEL" and "LOGFILE" are added automatically by the ParSys
# constructor.
params = []
params.append(starutil.ParNDG("IN1", "The first SCUBA-2 map"))
params.append(starutil.ParNDG("IN2", "The second SCUBA-2 map"))
params.append(
starutil.Par0S("REPORT", "Text file in which to "
"put the report",
"sc2compare.rep",
noprompt=True))
params.append(
starutil.Par0L("RETAIN",
"Retain temporary files?",
False,
noprompt=True))
# Set he default value for GLEVEL parameter, created by the ParSys
# constructor. This means that no logfile will be created by default.
starutil.glevel = starutil.NONE
# Catch any exception so that we can always clean up, even if control-C
# is pressed.
try:
# Declare the script parameters. Their positions in this list define
# their expected position on the script command line. They can also be
# specified by keyword on the command line. If no value is supplied on
# the command line, the user is prompted for a value when the parameter
# value is first accessed within this script. The parameters "MSG_FILTER",
# "ILEVEL", "GLEVEL" and "LOGFILE" are added automatically by the ParSys
# constructor.
params = []
params.append(starutil.ParNDG("IN", "The input time series NDFs",
starutil.get_task_par("DATA_ARRAY",
"GLOBAL",
default=starutil.Parameter.UNSET)))
params.append(starutil.ParNDG("OUT", "The output map", default=None,
exists=False, minsize=0, maxsize=1 ))
params.append(starutil.Par0I("NITER", "No. of iterations to perform",
0, noprompt=True))
params.append(starutil.Par0F("PIXSIZE", "Pixel size (arcsec)", None,
maxval=1000, minval=0.01))
params.append(starutil.Par0S("CONFIG", "Map-maker tuning parameters",
"^$STARLINK_DIR/share/smurf/dimmconfig.lis"))
params.append(starutil.ParNDG("ITERMAP", "Output cube holding itermaps",
# Catch any exception so that we can always clean up, even if control-C
# is pressed.
try:
# Declare the script parameters. Their positions in this list define
# their expected position on the script command line. They can also be
# specified by keyword on the command line. If no value is supplied on
# the command line, the user is prompted for a value when the parameter
# value is first accessed within this script. The parameters "MSG_FILTER",
# "ILEVEL", "GLEVEL" and "LOGFILE" are added automatically by the ParSys
# constructor.
params = []
params.append(
starutil.ParNDG(
"IN", "The input time-series NDFs",
"$STARLINK_DIR/share/smurf/s4a20091214_00015_0002.sdf"))
params.append(
starutil.Par0S("OUT", "The base name for the output FITS files"))
params.append(
starutil.ParNDG("FAKEMAP",
"A fake map to add into the data",
default=None,
noprompt=True,
maxsize=1))
params.append(
starutil.Par0L("RETAIN",
"Retain temporary files?",
# is pressed.
try:
# Declare the script parameters. Their positions in this list define
# their expected position on the script command line. They can also be
# specified by keyword on the command line. If no value is supplied on
# the command line, the user is prompted for a value when the parameter
# value is first accessed within this script. The parameters "MSG_FILTER",
# "ILEVEL", "GLEVEL" and "LOGFILE" are added automatically by the ParSys
# constructor.
params = []
params.append(
starutil.ParNDG(
"IN", "The input time series NDFs",
starutil.get_task_par("DATA_ARRAY",
"GLOBAL",
default=starutil.Parameter.UNSET)))
params.append(
starutil.ParNDG("OUT",
"The output map",
default=None,
exists=False,
minsize=0,
maxsize=1))
params.append(
starutil.Par0I("NITER",
"No. of iterations to perform",
0,
