aboutsummaryrefslogtreecommitdiff
path: root/noao/digiphot/daophot/daoedit/dpeconfirm.x
diff options
context:
space:
mode:
Diffstat (limited to 'noao/digiphot/daophot/daoedit/dpeconfirm.x')
-rw-r--r--noao/digiphot/daophot/daoedit/dpeconfirm.x433
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
generated by cgit (git 2.34.1) at 2025-09-20 07:22:11 -0400