diff options
| author | Joe Hunkeler <jhunkeler@gmail.com> | 2015-08-11 16:51:37 -0400 |
|---|---|---|
| committer | Joe Hunkeler <jhunkeler@gmail.com> | 2015-08-11 16:51:37 -0400 |
| commit | 40e5a5811c6ffce9b0974e93cdd927cbcf60c157 (patch) | |
| tree | 4464880c571602d54f6ae114729bf62a89518057 /noao/onedspec/ecidentify/ecpeak.x | |
| download | iraf-osx-40e5a5811c6ffce9b0974e93cdd927cbcf60c157.tar.gz | |
Repatch (from linux) of OSX IRAF
Diffstat (limited to 'noao/onedspec/ecidentify/ecpeak.x')
| -rw-r--r-- | noao/onedspec/ecidentify/ecpeak.x | 24 |
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 |