HELP_tap_irsa.csh

#!/bin/tcsh
#----------------------------------------------------------------------
# Purpose:
#   To launch a query to the IRSA TAP service, producing a FITS FILE
#   
#----------------------------------------------------------------------
# Parameters
#    $1 the name of the field  - purely to determine names of outputs
#    $2-$5 RAMIN, RAMAX, DECMIN, DECMAX - all in decimal degrees
#    $6 table name in IRSA (see irsa_tables.* or get latest versions from 
#    curl -o irsa_tables.fits "http://irsa.ipac.caltech.edu/TAP/sync?QUERY=SELECT+*+FROM+TAP_SCHEMA.tables&FORMAT=FITS"
#
#----------------------------------------------------------------------
# How to Improve:
#    The polygon is bounded by great circles, so this may not be exactly
#    what the user expects by specifying min/max limits.  Better might be
#    to use the BOX geometry but then we'd need to do some calculations in
#    the script and I couldn't be bothered.
#
#    Could have an option to select FITS/VOTABLE etc.  Current a VOT format
#    version is commented out
#
#----------------------------------------------------------------------
# Author(s): 
#    Seb Oliver (University of Sussex)
#----------------------------------------------------------------------
# Date: 
#    29th November 2015 (tested with commented header)
#----------------------------------------------------------------------
# Acknowledgement: 
#    This script was produced as part of the Herschel Extragalactic 
#    Legacy Project: HELP
#    The research leading to these results has received funding from the 
#    European Union Seventh Framework Programme FP7/2007-2013/ under 
#    grant agreement no 607254.”
#    "This publication reflects only the author’s view and the European 
#     Union is not responsible for any use that may be made of the 
#     information contained therein.” 
#----------------------------------------------------------------------
# end of header
#----------------------------------------------------------------------

set field=$1
set ramin=$2
set ramax=$3
set decmin=$4
set decmax=$5
set table=$6

set id=$field"_"$table

# FITS FORMAT
 curl  --stderr $id.log -v "http://irsa.ipac.caltech.edu/TAP/async?QUERY=SELECT+*+FROM+"$table"+WHERE+CONTAINS(POINT('J2000',ra,dec),POLYGON('J2000',"$ramin","$decmin","$ramax","$decmin","$ramax","$decmax","$ramin","$decmax"))=1&FORMAT=FITS&PHASE=RUN" 

## VOT FORMAT
# curl  --stderr $id.log -v "http://irsa.ipac.caltech.edu/TAP/async?QUERY=SELECT+*+FROM+"$table"+WHERE+CONTAINS(POINT('J2000',ra,dec),POLYGON('J2000',#"$ramin","$decmin","$ramax","$decmin","$ramax","$decmax","$ramin","$decmax"))=1&PHASE=RUN" 


echo $id >> HELP_tap_irsa.urls
grep Location $id.log   >> HELP_tap_irsa.urls
echo "-----------------------" >> HELP_tap_irsa.urls