Initial commit

1 files changed, 26 insertions, 0 deletions

diff --git a/noao/twodspec/apextract/apcenter.x b/noao/twodspec/apextract/apcenter.x
new file mode 100644
index 00000000..88f089d1
--- /dev/null
+++ b/noao/twodspec/apextract/apcenter.x
@@ -0,0 +1,26 @@
+include	<pkg/center1d.h>
+
+# AP_CENTER -- Locate the center of an emission profile.  This is done
+# using the CENTER1D algorithm.  The procedure gets the centering
+# parameters using CL queries.  If the center is not found because of the
+# RADIUS or THRESHOLD centering criteria then INDEF is returned.
+
+real procedure ap_center (x, data, npts)
+
+real	x		# Initial guess
+real	data[npts]	# Data
+int	npts		# Number of data points
+
+real	width		# Centering width
+real	radius		# Centering radius
+real	threshold	# Detection threshold
+
+real	apgetr(), center1d()
+
+begin
+	width = apgetr ("width")
+	radius = apgetr ("radius")
+	threshold = apgetr ("threshold")
+
+	return (center1d (x, data, npts, width, EMISSION, radius, threshold))
+end