1 files changed, 77 insertions, 0 deletions
diff --git a/pkg/images/immatch/src/wcsmatch/rgmatchio.x b/pkg/images/immatch/src/wcsmatch/rgmatchio.x
new file mode 100644
index 00000000..1a0de167
--- /dev/null
+++ b/pkg/images/immatch/src/wcsmatch/rgmatchio.x
@@ -0,0 +1,77 @@
+include "wcsxymatch.h"
+
+define	DEF_BUFSIZE	200
+
+# RG_RDXY -- Read in the x and y coordinates from a file.
+
+int procedure rg_rdxy (fd, x, y, wcs, xcolumn, ycolumn, xunits, yunits)
+
+int	fd			#I the input file descriptor
+pointer	x			#U pointer to the x coordinates
+pointer	y			#U pointer to the y coordinates
+int	wcs			#I the world coordinate system
+int	xcolumn			#I column containing the x coordinate
+int	ycolumn			#I column containing the y coordinate
+int	xunits			#I the x coordinate units
+int	yunits			#I the y coordinate units
+
+double	xval, yval
+int	i, ip, bufsize, maxcols, npts
+pointer	sp, str
+int	fscan(), nscan(), ctod()
+
+begin
+	call smark (sp)
+	call salloc (str, SZ_FNAME, TY_CHAR)
+
+	bufsize = DEF_BUFSIZE
+	call malloc (x, bufsize, TY_DOUBLE)
+	call malloc (y, bufsize, TY_DOUBLE)
+	maxcols = max (xcolumn, ycolumn)
+
+	npts = 0
+	while (fscan(fd) != EOF) {
+
+	    xval = INDEFD
+	    yval = INDEFD
+	    do i = 1, maxcols {
+		call gargwrd (Memc[str], SZ_FNAME)
+		if (i != nscan())
+		    break
+		ip = 1
+		if (i == xcolumn) {
+		    if (ctod (Memc[str], ip, xval) <= 0)
+			xval = INDEFD
+		} else if (i == ycolumn) {
+		    if (ctod (Memc[str], ip, yval) <= 0)
+			yval = INDEFD
+		}
+	    }
+	    if (IS_INDEFD(xval) || IS_INDEFD(yval))
+		next
+
+	    Memd[x+npts] = xval
+	    Memd[y+npts] = yval
+	    npts = npts + 1
+	    if (npts >= bufsize) {
+		bufsize = bufsize + DEF_BUFSIZE
+		call realloc (x, bufsize, TY_DOUBLE)
+		call realloc (y, bufsize, TY_DOUBLE)
+	    }
+	}
+
+	# Convert the coordinates if necessary.
+	switch (wcs) {
+	case RG_WORLD:
+	    if (xunits == RG_UHOURS)
+		call amulkd (Memd[x], 15.0d0, Memd[x], npts)
+	    if (yunits == RG_UHOURS)
+		call amulkd (Memd[y], 15.0d0, Memd[y], npts)
+	default:
+	    ;
+	}
+
+	call sfree (sp)
+
+	return (npts)
+end