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/digiphot/apphot/fitsky/apreadsky.x | 40 +++++++++++++++++++++++++++++++++
 1 file changed, 40 insertions(+)
 create mode 100644 noao/digiphot/apphot/fitsky/apreadsky.x

(limited to 'noao/digiphot/apphot/fitsky/apreadsky.x')

diff --git a/noao/digiphot/apphot/fitsky/apreadsky.x b/noao/digiphot/apphot/fitsky/apreadsky.x
new file mode 100644
index 00000000..3db9a183
--- /dev/null
+++ b/noao/digiphot/apphot/fitsky/apreadsky.x
@@ -0,0 +1,40 @@
+# AP_READSKY -- Procedure to read sky values, sigma and skew from a file
+# The x and y positions, sky mode, sigma, and skew values, number of sky
+# pixels and number of rejected sky pixels are assumed to be columns
+# 1 to 7 respectively.
+
+int procedure ap_readsky (fd, x, y, sky_mode, sky_sigma, sky_skew, nsky,
+    nsky_reject)
+
+int	fd		# sky file descriptor
+real	x, y		# center of sky annulus
+real	sky_mode	# sky valye
+real	sky_sigma	# sky sigma
+real	sky_skew	# skew of sky pixels
+int	nsky		# number of sky pixels
+int	nsky_reject	# number of rejected pixesl
+
+int	stat
+int	fscan(), nscan()
+
+begin
+	# Initialize.
+	sky_mode = INDEFR
+	sky_sigma = INDEFR
+	sky_skew = INDEFR
+	nsky = 0
+	nsky_reject = 0
+
+	# Read in and decode a sky file text line.
+	stat = fscan (fd)
+	if (stat == EOF)
+	    return (EOF)
+	call gargr (x)
+	call gargr (y)
+	call gargr (sky_mode)
+	call gargr (sky_sigma)
+	call gargr (sky_skew)
+	call gargi (nsky)
+	call gargi (nsky_reject)
+	return (nscan ())
+end
-- 
cgit