aboutsummaryrefslogtreecommitdiff
path: root/noao/digiphot/daophot/daoedit/dpecolon.x
diff options
context:
space:
mode:
Diffstat (limited to 'noao/digiphot/daophot/daoedit/dpecolon.x')
-rw-r--r--noao/digiphot/daophot/daoedit/dpecolon.x441
1 files changed, 441 insertions, 0 deletions
diff --git a/noao/digiphot/daophot/daoedit/dpecolon.x b/noao/digiphot/daophot/daoedit/dpecolon.x
new file mode 100644
index 00000000..36bccdf9
--- /dev/null
+++ b/noao/digiphot/daophot/daoedit/dpecolon.x
@@ -0,0 +1,441 @@
+include <error.h>
+include "daoedit.h"
+
+# DP_ECOLON -- Show/set the DAOPHOT algorithm parameters.
+
+procedure dp_ecolon (cmdstr, gd, redraw)
+
+char cmdstr[ARB] # input colon command
+int gd # pointer to the graphics stream
+int redraw # redraw the radial profile plot
+
+bool bval
+int cmd, pset, ival
+pointer sp, incmd, param, value
+real rval
+bool clgetb()
+int strdic(), nscan(), clgeti()
+real clgetr()
+string datapars "datapars"
+string centerpars "centerpars"
+string fitskypars "fitskypars"
+string photpars "photpars"
+string daopars "daopars"
+string findpars "findpars"
+
+begin
+ call smark (sp)
+ call salloc (incmd, SZ_LINE, TY_CHAR)
+ call salloc (param, SZ_LINE, TY_CHAR)
+ call salloc (value, SZ_LINE, TY_CHAR)
+
+ # Get the command.
+ call sscan (cmdstr)
+ call gargwrd (Memc[incmd], SZ_LINE)
+ if (Memc[incmd] == EOS) {
+ call sfree (sp)
+ return
+ }
+
+ cmd = strdic (Memc[incmd], Memc[incmd], SZ_LINE, CMD_LIST)
+ switch (cmd) {
+ case CMD_LPAR:
+ call gargwrd (Memc[incmd], SZ_LINE)
+ pset = strdic (Memc[incmd], Memc[incmd], SZ_LINE, PSET_LIST)
+ if (pset != 0)
+ call gdeactivate (gd, 0)
+ switch (pset) {
+ case PSET_DATAPARS:
+ call clcmdw ("lparam datapars")
+ case PSET_CENTERPARS:
+ call clcmdw ("lparam centerpars")
+ case PSET_FITSKYPARS:
+ call clcmdw ("lparam fitskypars")
+ case PSET_PHOTPARS:
+ call clcmdw ("lparam photpars")
+ case PSET_DAOPARS:
+ call clcmdw ("lparam daopars")
+ case PSET_FINDPARS:
+ call clcmdw ("lparam findpars")
+ default:
+ call printf ("Unknown parameter set\7\n")
+ }
+
+ case CMD_EPAR:
+ call gargwrd (Memc[incmd], SZ_LINE)
+ pset = strdic (Memc[incmd], Memc[incmd], SZ_LINE, PSET_LIST)
+ if (pset != 0)
+ call gdeactivate (gd, 0)
+ switch (pset) {
+ case PSET_DATAPARS:
+ call clcmdw ("eparam datapars")
+ case PSET_CENTERPARS:
+ call clcmdw ("eparam centerpars")
+ case PSET_FITSKYPARS:
+ call clcmdw ("eparam fitskypars")
+ case PSET_PHOTPARS:
+ call clcmdw ("eparam photpars")
+ case PSET_DAOPARS:
+ call clcmdw ("eparam daopars")
+ case PSET_FINDPARS:
+ call clcmdw ("eparam findpars")
+ default:
+ call printf ("Unknown parameter set\7\n")
+ }
+ if (pset != 0)
+ redraw = YES
+
+ case CMD_UNLEARN:
+ call gargwrd (Memc[incmd], SZ_LINE)
+ pset = strdic (Memc[incmd], Memc[incmd], SZ_LINE, PSET_LIST)
+ switch (pset) {
+ case PSET_DATAPARS:
+ call clcmdw ("unlearn datapars")
+ case PSET_CENTERPARS:
+ call clcmdw ("unlearn centerpars")
+ case PSET_FITSKYPARS:
+ call clcmdw ("unlearn fitskypars")
+ case PSET_PHOTPARS:
+ call clcmdw ("unlearn photpars")
+ case PSET_DAOPARS:
+ call clcmdw ("unlearn daopars")
+ case PSET_FINDPARS:
+ call clcmdw ("unlearn findpars")
+ default:
+ call printf ("Unknown parameter set\7\n")
+ }
+ if (pset != 0)
+ redraw = YES
+
+ case CMD_SCALE:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (datapars)
+ call pargstr (Memc[incmd])
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g\n")
+ call pargstr (Memc[incmd])
+ call pargr (clgetr (Memc[param]))
+ } else {
+ call clputr (Memc[param], rval)
+ redraw = YES
+ }
+
+ case CMD_FWHMPSF, CMD_SIGMA, CMD_DATAMIN, CMD_DATAMAX,
+ CMD_READNOISE, CMD_EPADU, CMD_ITIME, CMD_XAIRMASS:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (datapars)
+ call pargstr (Memc[incmd])
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g\n")
+ call pargstr (Memc[incmd])
+ call pargr (clgetr (Memc[param]))
+ } else
+ call clputr (Memc[param], rval)
+
+ case CMD_EMISSION:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (datapars)
+ call pargstr (Memc[incmd])
+ call gargb (bval)
+ if (nscan() == 1) {
+ call printf ("%s = %b\n")
+ call pargstr (Memc[incmd])
+ call pargb (clgetb (Memc[param]))
+ } else
+ call clputb (Memc[param], bval)
+
+ case CMD_NOISE:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (datapars)
+ call pargstr (Memc[incmd])
+ call clgstr (Memc[param], Memc[value], SZ_LINE)
+ call printf ("%s = %s\n")
+ call pargstr (Memc[incmd])
+ call pargstr (Memc[value])
+
+ case CMD_CCDREAD, CMD_GAIN, CMD_EXPOSURE, CMD_AIRMASS, CMD_FILTER,
+ CMD_OBSTIME, CMD_IFILTER, CMD_OTIME:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (datapars)
+ call pargstr (Memc[incmd])
+ call gargwrd (Memc[value], SZ_LINE)
+ if (nscan() == 1) {
+ call clgstr (Memc[param], Memc[value], SZ_LINE)
+ call printf ("%s = \"%s\"\n")
+ call pargstr (Memc[incmd])
+ call pargstr (Memc[value])
+ } else
+ call clpstr (Memc[param], Memc[value])
+
+ case CMD_CALGORITHM:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (centerpars)
+ call pargstr (Memc[incmd])
+ call gargwrd (Memc[value], SZ_LINE)
+ if (nscan() == 1) {
+ call clgstr (Memc[param], Memc[value], SZ_LINE)
+ call printf ("%s = \"%s\"\n")
+ call pargstr (Memc[incmd])
+ call pargstr (Memc[value])
+ } else {
+ call clpstr (Memc[param], Memc[value])
+ }
+
+ case CMD_CBOX:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (centerpars)
+ call pargstr (Memc[incmd])
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g\n")
+ call pargstr (Memc[incmd])
+ call pargr (clgetr (Memc[param]))
+ } else {
+ call clputr (Memc[param], rval)
+ redraw = YES
+ }
+
+ case CMD_CTHRESHOLD, CMD_MAXSHIFT, CMD_MINSNRATIO, CMD_RCLEAN,
+ CMD_RCLIP, CMD_KCLEAN:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (centerpars)
+ call pargstr (Memc[incmd])
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g\n")
+ call pargstr (Memc[incmd])
+ call pargr (clgetr (Memc[param]))
+ } else
+ call clputr (Memc[param], rval)
+
+ case CMD_CMAXITER:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (centerpars)
+ call pargstr (Memc[incmd])
+ call gargi (ival)
+ if (nscan() == 1) {
+ call printf ("%s = %d\n")
+ call pargstr (Memc[incmd])
+ call pargi (clgeti (Memc[param]))
+ } else
+ call clputi (Memc[param], ival)
+
+ case CMD_CLEAN, CMD_MKCENTER:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (centerpars)
+ call pargstr (Memc[incmd])
+ call gargb (bval)
+ if (nscan() == 1) {
+ call printf ("%s = %b\n")
+ call pargstr (Memc[incmd])
+ call pargb (clgetb (Memc[param]))
+ } else
+ call clputb (Memc[param], bval)
+
+ case CMD_SALGORITHM:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (fitskypars)
+ call pargstr (Memc[incmd])
+ call gargwrd (Memc[value], SZ_LINE)
+ if (nscan() == 1) {
+ call clgstr (Memc[param], Memc[value], SZ_LINE)
+ call printf ("%s = \"%s\"\n")
+ call pargstr (Memc[incmd])
+ call pargstr (Memc[value])
+ } else {
+ call clpstr (Memc[param], Memc[value])
+ call clgstr (Memc[param], Memc[value], SZ_LINE)
+ }
+
+ case CMD_ANNULUS, CMD_DANNULUS:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (fitskypars)
+ call pargstr (Memc[incmd])
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g\n")
+ call pargstr (Memc[incmd])
+ call pargr (clgetr (Memc[param]))
+ } else {
+ call clputr (Memc[param], rval)
+ redraw = YES
+ }
+
+ case CMD_SLOCLIP, CMD_SHICLIP, CMD_SLOREJECT, CMD_SHIREJECT,
+ CMD_KHIST, CMD_BINSIZE, CMD_SKYVALUE, CMD_RGROW:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (fitskypars)
+ call pargstr (Memc[incmd])
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g\n")
+ call pargstr (Memc[incmd])
+ call pargr (clgetr (Memc[param]))
+ } else
+ call clputr (Memc[param], rval)
+
+ case CMD_SMAXITER, CMD_SNREJECT:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (fitskypars)
+ call pargstr (Memc[incmd])
+ call gargi (ival)
+ if (nscan() == 1) {
+ call printf ("%s = %d\n")
+ call pargstr (Memc[incmd])
+ call pargi (clgeti (Memc[param]))
+ } else
+ call clputi (Memc[param], ival)
+
+ case CMD_SMOOTH, CMD_MKSKY:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (fitskypars)
+ call pargstr (Memc[incmd])
+ call gargb (bval)
+ if (nscan() == 1) {
+ call printf ("%s = %b\n")
+ call pargstr (Memc[incmd])
+ call pargb (clgetb (Memc[param]))
+ } else
+ call clputb (Memc[param], bval)
+
+ case CMD_WEIGHTING:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (photpars)
+ call pargstr (Memc[incmd])
+ call clgstr (Memc[param], Memc[value], SZ_LINE)
+ call printf ("%s = %s\n")
+ call pargstr (Memc[incmd])
+ call pargstr (Memc[value])
+
+ case CMD_APERTURES:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (photpars)
+ call pargstr (Memc[incmd])
+ call gargwrd (Memc[value], SZ_LINE)
+ if (nscan() == 1) {
+ call clgstr (Memc[param], Memc[value], SZ_LINE)
+ call printf ("%s = \"%s\"\n")
+ call pargstr (Memc[incmd])
+ call pargstr (Memc[value])
+ } else {
+ call clpstr (Memc[param], Memc[value])
+ redraw = YES
+ }
+
+ case CMD_ZMAG:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (photpars)
+ call pargstr (Memc[incmd])
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g\n")
+ call pargstr (Memc[incmd])
+ call pargr (clgetr (Memc[param]))
+ } else {
+ call clputr (Memc[param], rval)
+ }
+
+ case CMD_MKAPERT:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (photpars)
+ call pargstr (Memc[incmd])
+ call gargb (bval)
+ if (nscan() == 1) {
+ call printf ("%s = %b\n")
+ call pargstr (Memc[incmd])
+ call pargb (clgetb (Memc[param]))
+ } else
+ call clputb (Memc[param], bval)
+
+ case CMD_FUNCTION:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (daopars)
+ call pargstr (Memc[incmd])
+ call clgstr (Memc[param], Memc[value], SZ_LINE)
+ call printf ("%s = %s\n")
+ call pargstr (Memc[incmd])
+ call pargstr (Memc[value])
+
+ case CMD_RECENTER, CMD_FITSKY, CMD_GROUPSKY, CMD_SATURATED:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (daopars)
+ call pargstr (Memc[incmd])
+ call gargb (bval)
+ if (nscan() == 1) {
+ call printf ("%s = %b\n")
+ call pargstr (Memc[incmd])
+ call pargb (clgetb (Memc[param]))
+ } else
+ call clputb (Memc[param], bval)
+
+ case CMD_PSFRAD, CMD_FITRAD:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (daopars)
+ call pargstr (Memc[incmd])
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g\n")
+ call pargstr (Memc[incmd])
+ call pargr (clgetr (Memc[param]))
+ } else
+ call clputr (Memc[param], rval)
+
+ case CMD_MATCHRAD, CMD_SANNULUS, CMD_WSANNULUS, CMD_FLATERR,
+ CMD_PROFERR, CMD_CRITOVERLAP, CMD_CLIPRANGE:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (daopars)
+ call pargstr (Memc[incmd])
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g\n")
+ call pargstr (Memc[incmd])
+ call pargr (clgetr (Memc[param]))
+ } else
+ call clputr (Memc[param], rval)
+
+ case CMD_VARORDER, CMD_NCLEAN, CMD_MAXITER, CMD_MAXGROUP,
+ CMD_MAXNSTAR, CMD_CLIPEXP:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (daopars)
+ call pargstr (Memc[incmd])
+ call gargi (ival)
+ if (nscan() == 1) {
+ call printf ("%s = %d\n")
+ call pargstr (Memc[incmd])
+ call pargi (clgeti (Memc[param]))
+ } else
+ call clputi (Memc[param], ival)
+
+ case CMD_THRESHOLD, CMD_NSIGMA, CMD_RATIO, CMD_THETA, CMD_SHARPLO,
+ CMD_SHARPHI, CMD_ROUNDLO, CMD_ROUNDHI:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (findpars)
+ call pargstr (Memc[incmd])
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g\n")
+ call pargstr (Memc[incmd])
+ call pargr (clgetr (Memc[param]))
+ } else
+ call clputr (Memc[param], rval)
+
+ case CMD_MKDETECTIONS:
+ call sprintf (Memc[param], SZ_LINE, "%s.%s")
+ call pargstr (findpars)
+ call pargstr (Memc[incmd])
+ call gargb (bval)
+ if (nscan() == 1) {
+ call printf ("%s = %b\n")
+ call pargstr (Memc[incmd])
+ call pargb (clgetb (Memc[param]))
+ } else
+ call clputb (Memc[param], bval)
+
+ default:
+ call printf ("Unknown or ambiguous colon command\7\n")
+ }
+
+ call sfree (sp)
+end
generated by cgit (git 2.34.1) at 2025-09-20 09:39:56 -0400