Initial commit

1 files changed, 433 insertions, 0 deletions

diff --git a/noao/digiphot/daophot/daoedit/dpeconfirm.x b/noao/digiphot/daophot/daoedit/dpeconfirm.x
new file mode 100644
index 00000000..31aa359f
--- /dev/null
+++ b/noao/digiphot/daophot/daoedit/dpeconfirm.x
@@ -0,0 +1,433 @@
+include "daoedit.h"
+
+# DP_CBANNER -- Pritnt the confirm banner.
+
+procedure dp_cbanner ()
+
+begin
+	call printf ("\nVERIFY THE NEW VALUES\n")
+end
+
+
+# DP_CFWHMPSF -- Confirm the new value of fwhmpsf.
+
+procedure dp_cfwhmpsf()
+
+real	fwhmpsf, scale, rval
+int	scan(), nscan()
+real	clgetr()
+
+begin
+	# Get the current values.
+	scale = 1.0 / clgetr ("datapars.scale")
+	fwhmpsf = clgetr ("datapars.fwhmpsf")
+
+	# Confirm the fwhmpsf.
+	call printf ("FWHM of features (%g scale units) (CR or value): ")
+	    call pargr (fwhmpsf)
+	call flush (STDOUT)
+	if (scan() != EOF) {
+	    call gargr (rval)
+	    if (nscan() == 1)
+	        fwhmpsf = rval
+	}
+	call printf ("\tNew FWHM of features: %g scale units  %g pixels\n")
+	    call pargr (fwhmpsf)
+	    call pargr (scale * fwhmpsf)
+
+	# Store the new fwhmpsf.
+	call clputr ("datapars.fwhmpsf", fwhmpsf)
+end
+
+
+# DP_CSIGMA -- Confirm the sigma parameters.
+
+procedure dp_csigma()
+
+real	sigma, rval
+int	scan(), nscan()
+real	clgetr()
+
+begin
+	# Get the current value.
+	sigma = clgetr ("datapars.sigma")
+
+	# Confirm the sky sigma.
+	call printf (
+	    "Standard deviation of background (%g counts) (CR or value): ")
+	    call pargr (sigma)
+	call flush (STDOUT)
+	if (scan() != EOF) {
+	    call gargr (rval)
+	    if (nscan () == 1)
+	        sigma = rval
+	}
+	call printf ("\tNew standard deviation of background: %g counts\n")
+	    call pargr (sigma)
+
+	# Store the new sky sigma.
+	call clputr ("datapars.sigma", sigma)
+end
+
+
+# DP_CDMIN -- Confirm the good data minimum value.
+
+procedure dp_cdmin ()
+
+real	datamin, rval
+int	scan(), nscan()
+real	clgetr()
+
+begin
+	# Get the current value.
+	datamin = clgetr ("datapars.datamin")
+
+	# Confirm the new minimum good data value.
+	call printf (
+	"Minimum good data value (%g counts) (CR or value): ")
+	    call pargr (datamin)
+	call flush (STDOUT)
+	if (scan() != EOF) {
+	    call gargr (rval)
+	    if (nscan () == 1)
+	        datamin = rval
+	}
+	call printf ("\tNew good data minimum: %g counts\n")
+	    call pargr (datamin)
+
+	# Store the new good data minimum.
+	call clputr ("datapars.datamin", datamin)
+end
+
+
+# DP_CDMAX -- Confirm the good data maximum value.
+
+procedure dp_cdmax ()
+
+real	datamax, rval
+int	scan(), nscan()
+real	clgetr()
+
+begin
+	# Get the current value.
+	datamax = clgetr ("datapars.datamax")
+
+	# Confirm the new maximum good data value.
+	call printf (
+	"Maximum good data value (%g counts) (CR or value): ")
+	    call pargr (datamax)
+	call flush (STDOUT)
+	if (scan() != EOF) {
+	    call gargr (rval)
+	    if (nscan () == 1)
+	        datamax = rval
+	}
+	call printf ("\tNew good data maximum: %g counts\n")
+	    call pargr (datamax)
+
+	# Store the new maximum good data value.
+	call clputr ("datapars.datamax", datamax)
+end
+
+
+# DP_CCBOX -- Confirm the cbox parameter.
+
+procedure dp_ccbox()
+
+real	scale, cbox, rval
+int	scan(), nscan()
+real	clgetr()
+
+begin
+	# Get the current values.
+	scale = 1.0 / clgetr ("datapars.scale")
+	cbox = clgetr ("centerpars.cbox")
+
+	# Confirm the centering box value.
+	call printf ("Centering box width (%g scale units) (CR or value): ")
+	    call pargr (cbox)
+	call flush (STDOUT)
+	if (scan() != EOF) {
+	    call gargr (rval)
+	    if (nscan () == 1)
+	        cbox = rval
+	}
+	call printf ("\tNew centering box width: %g scale units  %g pixels\n")
+	    call pargr (cbox)
+	    call pargr (scale * cbox)
+
+	# Store the new centering box.
+	call clputr ("centerpars.cbox", cbox)
+end
+
+
+# DP_CRCLEAN -- Confirm the cleaning radius.
+
+procedure dp_crclean()
+
+real	scale, rclean, rval
+int	scan(), nscan()
+real	clgetr()
+
+begin
+	# Get the current values.
+	scale = 1.0 / clgetr ("datapars.scale")
+	rclean = clgetr ("centerpars.rclean")
+
+	# Confirm the cleaning radius..
+	call printf ("Cleaning radius (%g scale units) (CR or value): ")
+	    call pargr (rclean)
+	call flush (STDOUT)
+	if (scan() != EOF) {
+	    call gargr (rval)
+	    if (nscan () == 1)
+	        rclean = rval
+	}
+	call printf ("\tNew cleaning radius: %g scale units  %g pixels\n")
+	    call pargr (rclean)
+	    call pargr (scale * rclean)
+
+	# Store the new cleaning radius.
+	call clputr ("centerpars.rclean", rclean)
+end
+
+
+# DP_CRCLIP -- Confirm the clipping radius.
+
+procedure dp_crclip()
+
+real	scale, rclip, rval
+int	scan(), nscan()
+real	clgetr()
+
+begin
+	# Get the current values.
+	scale = 1.0 / clgetr ("datapars.scale")
+	rclip = clgetr ("centerpars.rclip")
+
+	# Confirm the cleaning radius..
+	call printf ("Clipping radius (%g scale units) (CR or value): ")
+	    call pargr (rclip)
+	call flush (STDOUT)
+	if (scan() != EOF) {
+	    call gargr (rval)
+	    if (nscan () == 1)
+	        rclip = rval
+	}
+	call printf ("\tNew clipping radius: %g scale units  %g pixels\n")
+	    call pargr (rclip)
+	    call pargr (scale * rclip)
+
+	# Store the new clipping radius.
+	call clputr ("centerpars.rclip", rclip)
+end
+
+
+# DP_CANNULUS -- Confirm the inner radius of the sky annulus.
+
+procedure dp_cannulus()
+
+real	scale, annulus, rval
+int	scan(), nscan()
+real	clgetr()
+
+begin 
+	# Get the current values.
+	scale = 1.0 / clgetr ("datapars.scale")
+	annulus = clgetr ("fitskypars.annulus")
+
+	# Confirm the sky annulus.
+	call printf (
+	    "Inner radius of sky annulus (%g scale units) (CR or value): ")
+            call pargr (annulus)
+	call flush (STDOUT)
+	if (scan () != EOF) {
+	    call gargr (rval)
+	    if (nscan () == 1)
+	        annulus = rval
+	}
+	call printf (
+	    "\tNew inner radius of sky annulus: %g scale units %g pixels\n")
+	    call pargr (annulus)
+	    call pargr (scale * annulus)
+
+	# Store the new sky annulus.
+	call clputr ("fitskypars.annulus", annulus)
+end
+
+
+# DP_CDANNULUS -- Confirm the annulus width parameter.
+
+procedure dp_cdannulus()
+
+real	scale, dannulus, rval
+int	scan(), nscan()
+real	clgetr()
+
+begin
+	# Get the current values.
+	scale = 1.0 / clgetr ("datapars.scale")
+	dannulus = clgetr ("fitskypars.dannulus")
+
+	# Confirm the sky annulus width.
+	call printf (
+	    "Width of the sky annulus (%g scale units) (CR or value): ")
+	    call pargr (dannulus)
+	call flush (STDOUT)
+	if (scan() != EOF) {
+	    call gargr (rval)
+	    if (nscan () == 1)
+	        dannulus = rval
+	}
+	call printf (
+	    "\tNew width of the sky annulus: %g scale units %g pixels\n")
+	    call pargr (dannulus)
+	    call pargr (scale * dannulus)
+
+	# Save the new sky annulus width.
+	call clputr ("fitskypars.dannulus", dannulus)
+end
+
+
+# DP_CRGROW -- Confirm the region growing radius.
+
+procedure dp_crgrow()
+
+real	scale, rgrow, rval
+int	scan(), nscan()
+real	clgetr()
+
+begin
+	# Mark the region growing radius.
+	scale = 1.0 / clgetr ("datapars.scale")
+	rgrow = clgetr ("fitskypars.rgrow")
+
+	# Confirm the new region growing radius.
+	call printf (
+	    "Region growing radius (%g scale units) (CR or value): ")
+	    call pargr (rgrow)
+	call flush (STDOUT)
+	if (scan() != EOF) {
+	    call gargr (rval)
+	    if (nscan () == 1)
+	        rgrow = rval
+	}
+	call printf (
+	    "\tNew region growing radius: %g scale units %g pixels\n")
+	    call pargr (rgrow)
+	    call pargr (scale * rgrow)
+
+	# Save the new region growing radius.
+	call clputr ("fitskypars.rgrow", rgrow)
+end
+
+
+# DP_CAPER -- Confirm the aperture string.
+
+procedure dp_caper()
+
+int	i, naperts
+pointer	sp, apstr, newapstr, aperts
+real	scale
+int	scan(), nscan(), dp_gaperts()
+real	clgetr()
+
+begin
+	call smark (sp)
+	call salloc (apstr, SZ_LINE, TY_CHAR)
+	call salloc (newapstr, SZ_LINE, TY_CHAR)
+	call salloc (aperts, MAX_NAPERTS, TY_REAL)
+
+	# Get the current values.
+	scale = 1.0 / clgetr ("datapars.scale")
+	call clgstr ("photpars.apertures", Memc[apstr], SZ_LINE)
+
+	# Confirm the aperture string.
+	call printf (
+	"File/list of aperture radii (%s scale units) (CR or value): ")
+	    call pargstr (Memc[apstr])
+	call flush (STDOUT)
+
+	# Get the new apertures.
+	if (scan() != EOF) {
+	    call gargwrd (Memc[newapstr], SZ_LINE)
+	    if (nscan () == 1)
+		call strcpy (Memc[newapstr], Memc[apstr], SZ_LINE)
+	}
+
+	# Print the new apertures.
+	naperts = dp_gaperts (Memc[apstr], Memr[aperts], MAX_NAPERTS)
+	do i = 1, naperts {
+	    call printf ("\tAperture radius %d: %g scale units %g pixels\n")
+		call pargi (i)
+		call pargr (Memr[aperts+i-1])
+		call pargr (scale * Memr[aperts+i-1])
+	}
+
+	# Save the new aperture string.
+	call clpstr ("photpars.apertures", Memc[apstr])
+
+	call sfree (sp)
+end
+
+
+# DP_CPSFRAD -- Confirm the psf radius.
+
+procedure dp_cpsfrad()
+
+real	scale, psfrad, rval
+int	scan(), nscan()
+real	clgetr()
+
+begin
+	# Get the current values.
+	scale = 1.0 / clgetr ("datapars.scale")
+	psfrad = clgetr ("daopars.psfrad")
+
+	# Confirm the new PSF radius.
+	call printf ("PSF radius (%g scale units) (CR or value): ")
+	    call pargr (psfrad)
+	call flush (STDOUT)
+	if (scan() != EOF) {
+	    call gargr (rval)
+	    if (nscan () == 1)
+	        psfrad = rval
+	}
+	call printf ("\tNew PSF radius: %g scale units  %g pixels\n")
+	    call pargr (psfrad)
+	    call pargr (scale * psfrad)
+
+	# Store the new PSF radius.
+	call clputr ("daopars.psfrad", psfrad)
+end
+
+
+# DP_CFITRAD -- Confirm the fitting radius.
+
+procedure dp_cfitrad ()
+
+real	scale, fitrad, rval
+int	scan(), nscan()
+real	clgetr()
+
+begin
+	# Get the current values.
+	scale = 1.0 / clgetr ("datapars.scale")
+	fitrad = clgetr ("daopars.fitrad")
+
+	# Confirm the new fitting radius.
+	call printf ("Fitting radius (%g scale units) (CR or value): ")
+	    call pargr (fitrad)
+	call flush (STDOUT)
+	if (scan() != EOF) {
+	    call gargr (rval)
+	    if (nscan () == 1)
+	        fitrad = rval
+	}
+	call printf ("\tNew fitting radius: %g scale units  %g pixels\n")
+	    call pargr (fitrad)
+	    call pargr (scale * fitrad)
+
+	# Store the new fitting radius.
+	call clputr ("daopars.fitrad", fitrad)
+end