Initial commit

1 files changed, 76 insertions, 0 deletions

diff --git a/noao/digiphot/apphot/wphot/apwremag.x b/noao/digiphot/apphot/wphot/apwremag.x
new file mode 100644
index 00000000..5f611ec8
--- /dev/null
+++ b/noao/digiphot/apphot/wphot/apwremag.x
@@ -0,0 +1,76 @@
+include <mach.h>
+include "../lib/apphotdef.h"
+include "../lib/apphot.h"
+include "../lib/noisedef.h"
+include "../lib/photdef.h"
+include "../lib/phot.h"
+
+# AP_WREMAG -- Procedure to recompute the magnitudes inside a set of apertures
+# given that the sums and effective areas have already been computed.
+
+int procedure ap_wremag (ap, im, positive, skyval, skysig, nsky)
+
+pointer	ap		# pointer to the apphot structure
+pointer	im		# the input image descriptor
+int	positive	# emission or absorption features
+real	skyval		# sky value
+real	skysig		# sigma of sky
+int	nsky		# number of sky pixels
+
+int	nap
+pointer	nse, phot
+real	zmag
+
+begin
+	# Initalize.
+	phot = AP_PPHOT(ap)
+	nse = AP_NOISE(ap)
+	call amovkr (INDEFR, Memr[AP_MAGS(phot)], AP_NAPERTS(phot))
+	call amovkr (INDEFR, Memr[AP_MAGERRS(phot)], AP_NAPERTS(phot))
+        if (IS_INDEFR(AP_PXCUR(phot)) || IS_INDEFR(AP_PYCUR(phot))) {
+            AP_OPXCUR(phot) = AP_PXCUR(phot)
+            AP_OPYCUR(phot) = AP_PYCUR(phot)
+        } else {
+            switch (AP_WCSOUT(ap)) {
+            case WCS_WORLD, WCS_PHYSICAL:
+                call ap_ltoo (ap, AP_PXCUR(phot), AP_PYCUR(phot),
+                    AP_OPXCUR(phot), AP_OPYCUR(phot), 1)
+            case WCS_TV:
+                call ap_ltov (im, AP_PXCUR(phot), AP_PYCUR(phot),
+                    AP_OPXCUR(phot), AP_OPYCUR(phot), 1)
+            default:
+                AP_OPXCUR(phot) = AP_PXCUR(phot)
+                AP_OPYCUR(phot) = AP_PYCUR(phot)
+            }
+        }
+
+	if (IS_INDEFR(AP_PXCUR(phot)) || IS_INDEFR(AP_PYCUR(phot)))
+	    return (AP_APERT_NOAPERT)
+	if (IS_INDEFR(skyval))
+	    return (AP_APERT_NOSKYMODE)
+
+	nap = min (AP_NMINAP(phot) - 1, AP_NMAXAP(phot))
+
+	# Compute the magnitudes and errors.
+	if (positive == YES)
+	    call apcopmags (Memd[AP_SUMS(phot)], Memd[AP_AREA(phot)],
+	        Memr[AP_MAGS(phot)], Memr[AP_MAGERRS(phot)], nap,
+	        skyval, skysig, nsky, AP_ZMAG(phot), AP_NOISEFUNCTION(nse),
+		AP_EPADU(nse))
+	else
+	    call apconmags (Memd[AP_SUMS(phot)], Memd[AP_AREA(phot)],
+	        Memr[AP_MAGS(phot)], Memr[AP_MAGERRS(phot)], nap,
+	        skyval, skysig, nsky, AP_ZMAG(phot), AP_NOISEFUNCTION(nse),
+		AP_EPADU(nse), AP_READNOISE(nse))
+
+	# Correct for itime.
+	zmag = 2.5 * log10 (AP_ITIME(ap))
+	call aaddkr (Memr[AP_MAGS(phot)], zmag, Memr[AP_MAGS(phot)], nap)
+
+	if (AP_NMINAP(phot) <= AP_NMAXAP(phot))
+	    return (AP_APERT_BADDATA)
+	else if (AP_NMAXAP(phot) < AP_NAPERTS(phot))
+	    return (AP_APERT_OUTOFBOUNDS)
+	else
+	    return (AP_OK)
+end