From 40e5a5811c6ffce9b0974e93cdd927cbcf60c157 Mon Sep 17 00:00:00 2001
From: Joe Hunkeler <jhunkeler@gmail.com>
Date: Tue, 11 Aug 2015 16:51:37 -0400
Subject: Repatch (from linux) of OSX IRAF

---
 noao/nproto/ir/irmedr.x | 35 +++++++++++++++++++++++++++++++++++
 1 file changed, 35 insertions(+)
 create mode 100644 noao/nproto/ir/irmedr.x

(limited to 'noao/nproto/ir/irmedr.x')

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
-- 
cgit