Initial commit

1 files changed, 29 insertions, 0 deletions

diff --git a/noao/imred/vtel/imratio.x b/noao/imred/vtel/imratio.x
new file mode 100644
index 00000000..5586d204
--- /dev/null
+++ b/noao/imred/vtel/imratio.x
@@ -0,0 +1,29 @@
+# IMRATIO -- Divide two images and return the result in a third image.
+
+procedure imratio (numerator, denominator, ratio, xdim, ydim)
+
+real	numerator[xdim, ydim]		# input numerator
+real	denominator[xdim, ydim]		# input denominator
+real	ratio[xdim, ydim]		# output ratio image
+int	xdim, ydim			# dimensions of the image
+
+int	i
+real	ezero()
+extern	ezero()
+
+begin
+	do i = 1, ydim {
+	    call arltr (denominator[1,i], xdim, 1E-10, 0.0)
+	    call advzr (numerator[1,i], denominator[1,i], ratio[1,i], xdim,
+		ezero)
+	}
+end
+
+
+real procedure ezero (input)
+
+real	input
+
+begin
+	return (0.0)
+end