From fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 Mon Sep 17 00:00:00 2001
From: Joseph Hunkeler <jhunkeler@gmail.com>
Date: Wed, 8 Jul 2015 20:46:52 -0400
Subject: Initial commit

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

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

diff --git a/noao/nproto/ir/iriinit.x b/noao/nproto/ir/iriinit.x
new file mode 100644
index 00000000..a97ade8e
--- /dev/null
+++ b/noao/nproto/ir/iriinit.x
@@ -0,0 +1,28 @@
+# IR_VECINIT -- Procedure to initialize the intensity matching algorithm.
+# If the ranges are undefined and no matching is to take place the
+# ishifts are set to INDEFR and the routine returns. Otherwise the shifts
+# are all initialized to zero and shifts for the missing subrasters are
+# set to INDEFR.
+
+procedure ir_vecinit (deltai, nsubrasters, ranges)
+
+real	deltai[ARB]			# intensity shifts
+int	nsubrasters			# number of subrasters
+int	ranges[ARB]			# ranges of missing subrasters
+
+int	num
+int	get_next_number()
+
+begin
+	# Initialize the shifts to INDEFR.
+	call amovkr (INDEFR, deltai, nsubrasters)
+	if (ranges[1] == NULL)
+	    return
+
+	num = 0
+	while (get_next_number (ranges, num) != EOF) {
+	    if (num  > nsubrasters)
+		break
+	    deltai[num] = 0.0
+	}
+end
-- 
cgit