Initial commit

1 files changed, 35 insertions, 0 deletions

diff --git a/noao/nproto/ir/irmedr.x b/noao/nproto/ir/irmedr.x
new file mode 100644
index 00000000..436a3673
--- /dev/null
+++ b/noao/nproto/ir/irmedr.x
@@ -0,0 +1,35 @@
+# IRMEDR -- Procedure to compute the median of an array in which some
+# elements are undefined.
+
+real procedure irmedr (a, aindex, npts)
+
+real	a[ARB]		# input array
+int	aindex[ARB]	# definition array
+int	npts		# number of points
+
+int	i, n
+pointer	sp, b
+real	med
+real	asokr()
+
+begin
+	call smark (sp)
+	call salloc (b, npts, TY_REAL)
+
+	n = 0
+	do i = 1, npts {
+	    if (aindex[i] > 0) {
+		Memr[b+n] = a[i]
+		n = n + 1
+	    }
+	}
+
+	if (n == 0)
+	    med = INDEFR
+	else
+	    med = asokr (Memr[b], n, (n + 1) / 2)
+
+	call sfree (sp)
+
+	return (med)
+end