#{  GETCAT -- Call a Cone service URL and return the raw result.

procedure getcat (resource, fields)

string	resource		{ prompt = "Resource name"		}
string	fields			{ prompt = "Query images/fields"	}
string	pos   = ""		{ prompt = "POS string"			}
real	size  = 0.25		{ prompt = "SIZE value (degrees)"	}

bool	samp     = no		{ prompt = "Broadcast result table(s)?"	}
bool	plot     = no		{ prompt = "Plot result table(s)?"	}
bool	display  = no		{ prompt = "Display field?"		}
bool	overplot  = no		{ prompt = "Overplot display?"		}
string	output = "STDOUT"	{ prompt = "Output filename"		}
string	format = "raw"		{ prompt = "Output format",
				    min="ascii|csv|votable|fits|xml|raw"} 
int	status = 0		{ prompt = "Service status code"	}

begin
	string  lname, lres, lpos, url, loname, tname, ltype, args, tcat
	bool	ldisp, lover, lplot
	real	ra, dec
	int     nread, len, nres, nfields, nrows, ncols, comma


	# Get params to local variables
	#lres    = resource		
	#lname   = fields
	lpos    = ""
	lname   = ""
	print (resource) | translit ("STDIN"," ","+") | scan (lres)
	print (fields) | translit ("STDIN"," ","+") | scan (lname)
	if (fields == "") {
	    lpos = pos
	}

	ltype   = format
	loname  = output
	ldisp   = display
	lover   = overplot
	lplot   = plot

	tcat = mktemp ("tmp$tcat")

	files (lres, sort-) | count ("STDIN") | scan (nres)
	if (substr (lname, 1, 1) == "@") {
	    if (imaccess (substr (lname, 2, strlen (lname))) == yes) {
	        sections (lname, opt="root", > tcat)
		count (tcat) | scan (nfields)
	    }
	} else {
	    if (lpos == "") {
		nfields = 1
	    } else {
	        files (lname, sort-, > tcat)
	        count (tcat) | scan (nfields)
	    }
	}

	# Simple error checking.
	if (loname == "STDOUT" && ltype == "fits") {
	    error (0, "FITS output not allowed to STDOUT")
	}

	if (nres > 1 || nfields > 1) {
	    vodata (lres, "@"//tcat, size=size, samp=samp, count-, all-,
		type="catalog", output=output, format=format)

	    delete (tcat, verify-, >& "dev$null")
	    return
	}

	#  Resolve resource name.
	if (lres == "")
	    error (0, "No resource specified")
	    
	if (substr (lres, 1, 7) == "http://") {
	    url  = lres

	#} else if (access ("uparm$url") == yes) {
	#    list = "uparm$url" ; nread = fscan (list, url) ; list = ""
	} else {
	    regdb ("resolve", lres, type="C", >& "dev$null")
	    if (regdb.status == 0 && regdb.svctype != "C")  {
	        error (0, "Resource '" // lres // "' is not a catalog service")
	    } else if (regdb.status == 1) {
                url   = regResolver (lres)  # not found, query the Registry
            } else {
                url   = regdb.url
	    }

            if (url == "INDEF") {
                error (0, "Resource '" // lres // "' is not a known service")
            }
	}

	#  Make sure the URL has a trailing '/' or '&'
        len = strlen (url)
        if (substr(url,len-4,len) == "&")
            url = substr (url, 1, len-5)
        if (substr(url,len,len) != '?' && substr(url,len,len) != '&') {
            if (strstr ("?", url) == 0)
                url = url // "?"
            else
                url = url // "&"
        }


	# Determine the query params from the image WCS.
	if (imaccess (lname) == no) {
	    if (lpos == "") {
	        #error (0, "Cannot open image '" // lname // "'")
                sesame (lname, verbose-)
                ra  = sesame.ra
                dec = sesame.dec

	        if (ldisp)
		    dss (lname, use_disp+)
                lname = "cache$" // lname // ".fits"
	        args = "RA=" // ra // "&DEC=" // dec  // "&SR=" // size
	    } else {
		comma = stridx (",", lpos)
		ra = real (substr(lpos,1,comma-1))
		dec = real (substr(lpos,comma+1,strlen(lpos)))
	        args = "RA=" // ra 
	        args = args // "&DEC=" // dec 
	        args = args // "&SR=" // size
	    }

	} else {
	    wcsinfo (lname)
	    args = "RA=" // wcsinfo.ra 
	    args = args // "&DEC=" // wcsinfo.dec 
	    args = args // "&SR=" // wcsinfo.size

	    if (ldisp)
		display (lname, 1, >& "dev$null")
	}

	if (vo.runid != "")
	    args = args // "&RUNID=" // vo.runid

	# Create a temporary output name.
	tname = mktemp ("tmp$raw") // ".xml"

	# Call the raw service via the URLGET generic task.
	urlget (url//args, tname, extn="", use_cache+, cache="cache$")

	if (samp)
	    samp ("loadVOTable", tname, id=lres, >& "dev$null")

        # Get the size of the result table.
        votsize (osfn(tname)) |& scan (nrows, ncols)
        if (nrows == 0) {
            print ("No results found.")
            delete (tcat, verify-, >& "dev$null")
            delete (tname, verify-, >& "dev$null")
            status = 1
            return
        }

	# Do the overlay if we displayed an image.
	if (ldisp && lover) {
            votpos (tname, out=tcat, verb-, number+, header-)
            taboverlay (lname, tcat, lab=1, ra=2, dec=3, mkcolor=207)
	}
	if (lplot)
            votpos (tname, out="", verb-, plot+, number-, overplot=lover)

	# If we didn't want a VOTable, convert the file and rename to
	# desired output name.  A null output name means we don't want
	# to see the table at all.
	if (loname != "")
            votcopy (tname, loname, ltype, header+, verb-)

	delete (tname, verify-)
	delete (tcat, verify-,  >& "dev$null")
end