aboutsummaryrefslogtreecommitdiff
path: root/noao/digiphot/apphot/fitsky/apskycolon.x
diff options
context:
space:
mode:
Diffstat (limited to 'noao/digiphot/apphot/fitsky/apskycolon.x')
-rw-r--r--noao/digiphot/apphot/fitsky/apskycolon.x367
1 files changed, 367 insertions, 0 deletions
diff --git a/noao/digiphot/apphot/fitsky/apskycolon.x b/noao/digiphot/apphot/fitsky/apskycolon.x
new file mode 100644
index 00000000..35b59aa9
--- /dev/null
+++ b/noao/digiphot/apphot/fitsky/apskycolon.x
@@ -0,0 +1,367 @@
+include "../lib/apphot.h"
+include "../lib/noise.h"
+include "../lib/fitsky.h"
+include "../lib/display.h"
+
+
+# APSKYCOLON -- Procedure to process the fitsky colon commands.
+
+procedure apskycolon (ap, im, cl, out, stid, ltid, cmdstr, newimage,
+ newskybuf, newsky)
+
+pointer ap # pointer to the apphot structure
+pointer im # pointer to the iraf image
+int cl # coordinate file descriptor
+int out # output file descriptor
+int stid # output file sequence number
+int ltid # coord list sequence number
+char cmdstr[ARB] # command string
+int newimage # new image ?
+int newskybuf # new sky buffer ?
+int newsky # new sky fit ?
+
+int junk
+pointer sp, incmd, outcmd
+int strdic()
+
+begin
+ # Get the command.
+ call smark (sp)
+ call salloc (incmd, SZ_LINE, TY_CHAR)
+ call salloc (outcmd, SZ_LINE, TY_CHAR)
+ call sscan (cmdstr)
+ call gargwrd (Memc[incmd], SZ_LINE)
+ if (Memc[incmd] == EOS) {
+ call sfree (sp)
+ return
+ }
+
+ # Process the command.
+ if (strdic (Memc[incmd], Memc[outcmd], SZ_LINE, SCMDS) != 0)
+ call apscolon (ap, out, stid, cmdstr, newskybuf, newsky)
+ else if (strdic (Memc[incmd], Memc[outcmd], SZ_LINE, APCMDS) != 0)
+ call ap_apcolon (ap, im, cl, out, stid, ltid, cmdstr,
+ newimage, junk, junk, newskybuf, newsky, junk, junk)
+ else if (strdic (Memc[incmd], Memc[outcmd], SZ_LINE, NCMDS) != 0)
+ call apnscolon (ap, im, out, stid, cmdstr, junk, junk,
+ newskybuf, newsky, junk, junk)
+ else
+ call ap_simcolon (ap, cmdstr)
+
+ call sfree (sp)
+end
+
+
+# APSCOLON -- Procedure to examine and edit the sky fitting parameters.
+
+procedure apscolon (ap, out, stid, cmdstr, newbuf, newfit)
+
+pointer ap # pointer to the apphot structure
+int out # output file descriptor
+int stid # output file number
+char cmdstr # command string
+int newbuf # new sky buffer
+int newfit # new sky fit
+
+bool bval
+int ncmd, ival, stat
+pointer sp, cmd
+real rval
+
+bool itob()
+int nscan(), strdic(), btoi(), apstati()
+real apstatr()
+
+
+begin
+ # Get the command
+ call smark (sp)
+ call salloc (cmd, SZ_LINE, TY_CHAR)
+ call sscan (cmdstr)
+ call gargwrd (Memc[cmd], SZ_LINE)
+ if (Memc[cmd] == EOS) {
+ call sfree (sp)
+ return
+ }
+
+ # Process the command.
+ ncmd = strdic (Memc[cmd], Memc[cmd], SZ_LINE, SCMDS)
+ switch (ncmd) {
+ case SCMD_ANNULUS:
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g %s\n")
+ call pargstr (KY_ANNULUS)
+ call pargr (apstatr (ap, ANNULUS))
+ call pargstr (UN_SSCALEUNIT)
+ } else {
+ call apsetr (ap, ANNULUS, rval)
+ if (stid > 1)
+ call ap_rparam (out, KY_ANNULUS, rval, UN_SSCALEUNIT,
+ "inner radius of sky annulus")
+ newbuf = YES
+ newfit = YES
+ }
+ case SCMD_DANNULUS:
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g %s\n")
+ call pargstr (KY_DANNULUS)
+ call pargr (apstatr (ap, DANNULUS))
+ call pargstr (UN_SSCALEUNIT)
+ } else {
+ call apsetr (ap, DANNULUS, rval)
+ if (stid > 1)
+ call ap_rparam (out, KY_DANNULUS, rval, UN_SSCALEUNIT,
+ "width of the sky annulus")
+ newbuf = YES
+ newfit = YES
+ }
+ case SCMD_SALGORITHM:
+ call gargwrd (Memc[cmd], SZ_LINE)
+ if (Memc[cmd] == EOS) {
+ call apstats (ap, SSTRING, Memc[cmd], SZ_FNAME)
+ call printf ("%s = %s\n")
+ call pargstr (KY_SSTRING)
+ call pargstr (Memc[cmd])
+ } else {
+ stat = strdic (Memc[cmd], Memc[cmd], SZ_LINE, SFUNCS)
+ if (stat > 0) {
+ call apseti (ap, SKYFUNCTION, stat)
+ call apsets (ap, SSTRING, Memc[cmd])
+ if (stid > 1)
+ call ap_sparam (out, KY_SSTRING, Memc[cmd],
+ UN_SALGORITHM, "sky fitting algorithm")
+ newfit = YES
+ }
+ }
+ case SCMD_KHIST:
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g %s\n")
+ call pargstr (KY_K1)
+ call pargr (apstatr (ap, K1))
+ call pargstr (UN_SSIGMA)
+ } else {
+ call apsetr (ap, K1, rval)
+ if (stid > 1)
+ call ap_rparam (out, KY_K1, rval, UN_SSIGMA,
+ "half width of sky histogram")
+ newfit = YES
+ }
+ case SCMD_SLOREJECT:
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g %s\n")
+ call pargstr (KY_SLOREJECT)
+ call pargr (apstatr (ap, SLOREJECT))
+ call pargstr (UN_SSIGMA)
+ } else {
+ call apsetr (ap, SLOREJECT, rval)
+ if (stid > 1)
+ call ap_rparam (out, KY_SLOREJECT, rval, UN_SSIGMA,
+ "lower k-sigma rejection criterion")
+ newfit = YES
+ }
+ case SCMD_SHIREJECT:
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g %s\n")
+ call pargstr (KY_SHIREJECT)
+ call pargr (apstatr (ap, SHIREJECT))
+ call pargstr (UN_SSIGMA)
+ } else {
+ call apsetr (ap, SHIREJECT, rval)
+ if (stid > 1)
+ call ap_rparam (out, KY_SHIREJECT, rval, UN_SSIGMA,
+ "upper k-sigma rejection criterion")
+ newfit = YES
+ }
+ case SCMD_SLOCLIP:
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g %s\n")
+ call pargstr (KY_SLOCLIP)
+ call pargr (apstatr (ap, SLOCLIP))
+ call pargstr (UN_SPERCENT)
+ } else {
+ call apsetr (ap, SLOCLIP, rval)
+ if (stid > 1)
+ call ap_rparam (out, KY_SLOCLIP, rval, UN_SPERCENT,
+ "lower k-sigma rejection criterion")
+ newfit = YES
+ }
+ case SCMD_SHICLIP:
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g %s\n")
+ call pargstr (KY_SHICLIP)
+ call pargr (apstatr (ap, SHICLIP))
+ call pargstr (UN_SPERCENT)
+ } else {
+ call apsetr (ap, SHICLIP, rval)
+ if (stid > 1)
+ call ap_rparam (out, KY_SHICLIP, rval, UN_SPERCENT,
+ "lower k-sigma rejection criterion")
+ newfit = YES
+ }
+ case SCMD_SMAXITER:
+ call gargi (ival)
+ if (nscan () == 1) {
+ call printf ("%s = %d\n")
+ call pargstr (KY_SMAXITER)
+ call pargi (apstati (ap, SMAXITER))
+ } else {
+ call apseti (ap, SMAXITER, ival)
+ if (stid > 1)
+ call ap_iparam (out, KY_SMAXITER, ival, UN_SNUMBER,
+ "maximum number of iterations")
+ newfit = YES
+ }
+ case SCMD_BINSIZE:
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g %s\n")
+ call pargstr (KY_BINSIZE)
+ call pargr (apstatr (ap, BINSIZE))
+ call pargstr (UN_SSIGMA)
+ } else {
+ call apsetr (ap, BINSIZE, rval)
+ if (stid > 1)
+ call ap_rparam (out, KY_BINSIZE, rval, UN_SSIGMA,
+ "width of the sky histogram bin")
+ newfit = YES
+ }
+ case SCMD_SMOOTH:
+ call gargb (bval)
+ if (nscan () == 1) {
+ call printf ("%s = %b\n")
+ call pargstr (KY_SMOOTH)
+ call pargb (itob (apstati (ap, SMOOTH)))
+ } else {
+ call apseti (ap, SMOOTH, btoi (bval))
+ if (stid > 1)
+ call ap_bparam (out, KY_SMOOTH, bval, UN_SSWITCH,
+ "Lucy smooth the histogram")
+ newfit = YES
+ }
+ case SCMD_RGROW:
+ call gargr (rval)
+ if (nscan() == 1) {
+ call printf ("%s = %g %s\n")
+ call pargstr (KY_RGROW)
+ call pargr (apstatr (ap, RGROW))
+ call pargstr (UN_SSCALEUNIT)
+ } else {
+ call apsetr (ap, RGROW, rval)
+ if (stid > 1)
+ call ap_rparam (out, KY_RGROW, rval, UN_SSCALEUNIT,
+ "region growing radius")
+ newfit = YES
+ }
+ case SCMD_SNREJECT:
+ call gargi (ival)
+ if (nscan () == 1) {
+ call printf ("%s = %d\n")
+ call pargstr (KY_SNREJECT)
+ call pargi (apstati (ap, SNREJECT))
+ } else {
+ call apseti (ap, SNREJECT, ival)
+ if (stid > 1)
+ call ap_iparam (out, KY_SNREJECT, ival, UN_SNUMBER,
+ "maximum number of rejection cycles")
+ newfit = YES
+ }
+ case SCMD_SKYVALUE:
+ call gargr (rval)
+ if (nscan () == 1) {
+ call printf ("%s = %g\n")
+ call pargstr (KY_SKY_BACKGROUND)
+ call pargr (apstatr (ap, SKY_BACKGROUND))
+ call pargstr (UN_SCOUNTS)
+ } else {
+ call apsetr (ap, SKY_BACKGROUND, rval)
+ if (stid > 1)
+ call ap_rparam (out, KY_SKY_BACKGROUND, rval,
+ UN_SCOUNTS, "user supplied sky value")
+ newfit = YES
+ }
+ case SCMD_MKSKY:
+ call gargb (bval)
+ if (nscan () == 1) {
+ call printf ("%s = %b\n")
+ call pargstr (KY_MKSKY)
+ call pargb (itob (apstati (ap, MKSKY)))
+ } else {
+ call apseti (ap, MKSKY, btoi (bval))
+ }
+ default:
+ # do nothing gracefully
+ call printf ("Unrecognized command\7\n")
+ }
+
+ call sfree (sp)
+end
+
+
+# AP_SIMCOLON -- Procedure to process fitsky commands which alter parameters
+# other than the sky fitting parameters themselves.
+
+procedure ap_simcolon (ap, cmdstr)
+
+pointer ap # pointer to the apphot structure
+char cmdstr[ARB] # command string
+
+bool bval
+int ncmd
+pointer sp, cmd
+bool itob()
+int strdic(), nscan(), apstati(), btoi()
+
+begin
+ call smark (sp)
+ call salloc (cmd, SZ_LINE, TY_CHAR)
+
+ # Get the command.
+ call sscan (cmdstr)
+ call gargwrd (Memc[cmd], SZ_LINE)
+ if (Memc[cmd] == EOS) {
+ call sfree (sp)
+ return
+ }
+
+ # Process the commands.
+ ncmd = strdic (Memc[cmd], Memc[cmd], SZ_LINE, MISC)
+ switch (ncmd) {
+ case ACMD_SHOW:
+ call gargwrd (Memc[cmd], SZ_LINE)
+ ncmd = strdic (Memc[cmd], Memc[cmd], SZ_LINE, SSHOWARGS)
+ switch (ncmd) {
+ case SCMD_DATA:
+ call printf ("\n")
+ call ap_nshow (ap)
+ call printf ("\n")
+ case SCMD_SKY:
+ call printf ("\n")
+ call ap_spshow (ap)
+ call printf ("\n")
+ default:
+ call printf ("\n")
+ call ap_sshow (ap)
+ call printf ("\n")
+ }
+ case ACMD_RADPLOTS:
+ call gargb (bval)
+ if (nscan () == 1) {
+ call printf ("%s = %b\n")
+ call pargstr (KY_RADPLOTS)
+ call pargb (itob (apstati (ap, RADPLOTS)))
+ } else {
+ call apseti (ap, RADPLOTS, btoi (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 11:23:47 -0400