1 files changed, 188 insertions, 0 deletions
diff --git a/noao/digiphot/apphot/aplib/apverify2.x b/noao/digiphot/apphot/aplib/apverify2.x
new file mode 100644
index 00000000..b5f02f46
--- /dev/null
+++ b/noao/digiphot/apphot/aplib/apverify2.x
@@ -0,0 +1,188 @@
+include "../lib/apphot.h"
+include "../lib/noise.h"
+include "../lib/find.h"
+include "../lib/fitpsf.h"
+include "../lib/radprof.h"
+
+
+# AP_VTHRESHOLD -- Verify the full detection threshold.
+
+real procedure ap_vthreshold (ap)
+
+pointer	ap		# pointer to the apphot structure
+
+real	skysigma, threshold
+int	scan(), nscan()
+real	apstatr()
+
+begin
+	skysigma = apstatr (ap, SKYSIGMA)
+
+	# Confirm the threshold parameter.
+	call printf (
+	"Detection threshold in sigma (%g) (CR or value): ")
+	    call pargr (apstatr (ap, THRESHOLD))
+	call flush (STDOUT)
+	if (scan() == EOF)
+	    threshold = apstatr (ap, THRESHOLD)
+	else {
+	    call gargr (threshold)
+	    if (nscan () != 1)
+	        threshold = apstatr (ap, THRESHOLD)
+	}
+
+	call printf ("\tNew detection threshold: %g sigma %g counts\n")
+	    call pargr (threshold)
+	call apsetr (ap, THRESHOLD, threshold)
+	if (IS_INDEFR(skysigma))
+	    call pargr (INDEFR)
+	else
+	    call pargr (threshold * skysigma)
+
+	return (threshold)
+end
+
+# AP_VPFSTRING -- Verify the psf fitting function.
+
+procedure ap_vpfstring (ap, str, maxch)
+
+pointer	ap		# pointer to the apphot structure
+char	str[ARB]	# output string
+int	maxch		# maximum number of characters
+
+int	pfunc
+int	scan(), nscan(), strdic()
+
+begin
+	# Print the old string value.
+	call apstats (ap, PSFSTRING, str, maxch)
+	call printf ("Fitting function (%s) (CR or value): ")
+	    call pargstr (str)
+	call flush (STDOUT)
+
+	# Get the new value.
+	if (scan() != EOF) {
+	    call gargwrd (str, maxch)
+	    pfunc = strdic (str, str, maxch, PSFFUNCS)
+	    if (nscan () == 1 && pfunc > 0) {
+		call apseti (ap, PSFUNCTION, pfunc)
+		call apsets (ap, PSFSTRING, str)
+	    }
+	}
+
+	# Print the new value.
+	call apstats (ap, PSFSTRING, str, SZ_LINE)
+	call printf ("\tNew fitting function: %s\n")
+	    call pargstr (str)
+end
+
+
+# AP_VPSFAPERT -- Verify the point spread function fitting aperture.
+
+real procedure ap_vpsfapert (ap)
+
+pointer	ap		# pointer to the apphot strucuture
+
+real	scale, psfapert
+int	scan(), nscan()
+real	apstatr()
+
+begin
+	scale = apstatr (ap, SCALE)
+
+	# Print the old value.
+	call printf ("Fitting box width in scale units (%g) (CR or value): ")
+	    call pargr (2.0 * apstatr (ap, PSFAPERT))
+	call flush (STDOUT)
+
+	# Get the new value.
+	if (scan() == EOF)
+	    psfapert = 2.0 * apstatr (ap, PSFAPERT)
+	else {
+	    call gargr (psfapert)
+	    if (nscan () != 1)
+	        psfapert = 2.0 * apstatr (ap, PSFAPERT)
+	}
+
+	# Print the new  value.
+	call printf ("\tNew fitting box width: %g scale units  %g pixels\n")
+	    call pargr (psfapert)
+	    call pargr (scale 	* psfapert)
+	call apsetr (ap, PSFAPERT, psfapert / 2.0)
+
+	return (psfapert / 2.0)
+end
+
+
+# AP_VRPRADIUS -- Verify the radial profile fitting radius.
+
+real procedure ap_vrpradius (ap)
+
+pointer	ap		# pointer to the apphot strucuture
+
+real	scale, radius
+int	scan(), nscan()
+real	apstatr()
+
+begin
+	scale = apstatr (ap, SCALE)
+
+	# Print the old value.
+	call printf ("Fitting radius in scale units (%g) (CR or value): ")
+	    call pargr (apstatr (ap, RPRADIUS))
+	call flush (STDOUT)
+
+	# Get the new value.
+	if (scan() == EOF)
+	    radius = apstatr (ap, RPRADIUS)
+	else {
+	    call gargr (radius)
+	    if (nscan () != 1)
+	        radius = apstatr (ap, RPRADIUS)
+	}
+
+	# Print the new  value.
+	call printf ("\tNew fitting radius: %g scale units  %g pixels\n")
+	    call pargr (radius)
+	    call pargr (scale 	* radius)
+	call apsetr (ap, RPRADIUS, radius)
+
+	return (radius)
+end
+
+
+# AP_VSTEP -- Verify the profile step size.
+
+real procedure ap_vstep (ap)
+
+pointer	ap		# pointer to the apphot strucuture
+
+real	scale, step
+int	scan(), nscan()
+real	apstatr()
+
+begin
+	scale = apstatr (ap, SCALE)
+
+	# Print the old value.
+	call printf ("Step size in scale units (%g) (CR or value): ")
+	    call pargr (apstatr (ap, RPSTEP))
+	call flush (STDOUT)
+
+	# Get the new value.
+	if (scan() == EOF)
+	    step = apstatr (ap, RPSTEP)
+	else {
+	    call gargr (step)
+	    if (nscan () != 1)
+	        step = apstatr (ap, RPSTEP)
+	}
+
+	# Print the new  value.
+	call printf ("\tNew step size: %g scale units  %g pixels\n")
+	    call pargr (step)
+	    call pargr (scale * step)
+	call apsetr (ap, RPSTEP, step)
+
+	return (step)
+end