Repatch (from linux) of OSX IRAF

1 files changed, 146 insertions, 0 deletions

diff --git a/vendor/voclient/libvoclient/examples/siap4.c b/vendor/voclient/libvoclient/examples/siap4.c
new file mode 100644
index 00000000..5d5b8def
--- /dev/null
+++ b/vendor/voclient/libvoclient/examples/siap4.c
@@ -0,0 +1,146 @@
+/************************************************************************
+ *   Call a SIAP search service and download resulting datasets.
+ *
+ *   Usage:        siap4 ra dec size [format [maximages [serviceURL]]]
+ *
+ *   Or call with no args for the built-in unit test.
+ *
+ *  M. Fitzpatrick, NOAO, June 2006
+ */
+
+
+#include <stdio.h>
+#include <stdlib.h>
+#include "VOClient.h"
+
+double  ra       = 12.0;			/* default values	*/
+double  dec      = 0.0;
+double  size     = 0.5;
+char	*format  = "ALL";
+int	maximages= 5;
+
+char   *service1  = "http://www-gsss.stsci.edu/gscvo/DSS2.jsp";
+char   *service2  = "http://skyview.gsfc.nasa.gov/cgi-bin/vo/sia.pl?";
+
+char   *server   = "6200:localhost";
+
+static void callSiapService (char *url, double ra, double dec, double size,
+    char *format, int maximages);
+
+
+int main (int argc, char *argv[])
+{
+    char *service = service2;
+
+
+    /* Process command line arguments.
+     */
+    if (argc <= 1) {
+	/* Use builtin defaults. */
+    } else if (argc >= 3) {
+	int arg = 1;
+    
+	/* Look for a server specification. */
+	if (strncmp (argv[arg], "-ds",3) == 0)
+	    server = argv[++arg];
+	ra   = atof (argv[arg++]);
+	dec  = atof (argv[arg++]);
+	size = atof (argv[arg++]);
+	if (arg < argc) format    = argv[arg++];
+	if (arg < argc) maximages = atoi (argv[arg++]);
+	if (arg < argc) service   = argv[arg++];
+
+    } else {
+        fprintf (stderr, 
+	    "Usage: siap4 [-ds server] ra dec size [fmt [maximg [siapURL]]]\n");
+        exit(1);
+    }
+
+    /* Now call the Siap Service and summarize the results.
+     */
+    callSiapService (service, ra, dec, size, format, maximages);
+
+    return (0);
+}
+
+
+/*  Simple test routine to call a Siap search service and summarize results.
+ */
+static void
+callSiapService (char *service_url, double ra, double dec, double size,
+    char *format, int maximages)
+{
+    char *voc_opts = NULL;
+    char *acref    = NULL;
+    char  fname[SZ_FNAME];
+    int   i, nattr = 0, nrec = 0;
+
+    DAL	      siap;				/* DAL Connection	 */
+    Query     query;				/* query handle		 */
+    QResponse qr;                               /* query response handle */
+    QRecord   rec;				/* result record handle	 */
+    QRAttribute v;				/* dataset attribute	 */
+	
+
+    /*  Initialize the VOClient code.  Error messages are printed by the
+     *  interface so we just quit if there is a problem.
+     */
+    if (voc_initVOClient (voc_opts) == ERR) 
+        return;
+
+    /*  Get a new connection to the named service.
+     */
+    siap = voc_openSiapConnection (service_url);    /* open a connection    */
+
+    /*  Form a query.  Here we'll use the one search size we're given for
+     *  both the RA,DEC sizes, and specify a null format.
+     */
+    query = voc_getSiapQuery (siap, ra, dec, size, size, format);
+
+    printf ("Executing Query:\n  %s\n\n", 
+        voc_getQueryString (query, SIAP_CONN, 0));
+    qr = voc_executeQuery (query);                  /* execute the query    */
+
+
+    /* Summarize response.
+     */
+    if ((nrec = voc_getRecordCount (qr)) <= 0) {    
+        fprintf (stderr, "no records matched");
+        exit (1);
+
+    } else {
+        rec = voc_getRecord (qr, 0);
+        nattr = (rec != ERR) ? voc_getAttrCount (rec) : 0;
+
+        printf ("# returns %d records containing %d attributes each\n#\n",
+            nrec, nattr);
+    }
+
+
+    /*  Download the first 'maximages' images.
+     */
+    printf ("Downloading images:\n");
+
+    for (i=0; i < nrec && i < maximages; i++) {
+	rec = voc_getRecord (qr, i);            /* get a row in the table    */
+
+	v = voc_getAttribute (rec, "AccessReference");
+
+	if (v <= 0) 
+	    continue;
+
+	bzero (fname, SZ_FNAME);
+	sprintf (fname, "dataset.%04d", i);
+	printf ("Downloading: %s\n", (acref = voc_stringValue (v)) );
+
+	if ( voc_getDataset (rec, acref, fname) == OK )
+	    printf ("Downloaded %s\n", fname);
+	else
+	    printf ("Download failed\n");
+    }
+
+    voc_closeConnection (siap);			/* close the siap connection */
+    voc_closeVOClient (1);		        /* clean up and shutdown     */
+
+    return;
+}