1 files changed, 24 insertions, 0 deletions
diff --git a/noao/onedspec/ecidentify/ecpeak.x b/noao/onedspec/ecidentify/ecpeak.x
new file mode 100644
index 00000000..f797fbac
--- /dev/null
+++ b/noao/onedspec/ecidentify/ecpeak.x
@@ -0,0 +1,24 @@
+include	"ecidentify.h"
+
+# EC_PEAK -- Find the peak value above continuum.
+
+double procedure ec_peak (ec, pix)
+
+pointer	ec			# ID pointer
+double	pix			# Pixel position
+double	peak			# Peak value
+
+int	c, l, u
+
+begin
+	if (IS_INDEFD(pix))
+	    return (INDEFD)
+
+	c = nint (pix)
+	l = max (1, nint (pix - EC_FWIDTH(ec)))
+	u = min (EC_NPTS(ec), nint (pix + EC_FWIDTH(ec)))
+	peak = IMDATA(ec,c) - (IMDATA(ec,l) +
+	    IMDATA(ec,u)) / 2.
+
+	return (peak)
+end