1 files changed, 316 insertions, 0 deletions
diff --git a/noao/digiphot/apphot/aplib/apstat1.x b/noao/digiphot/apphot/aplib/apstat1.x
new file mode 100644
index 00000000..0df220ae
--- /dev/null
+++ b/noao/digiphot/apphot/aplib/apstat1.x
@@ -0,0 +1,316 @@
+include "../lib/apphotdef.h"
+include "../lib/apphot.h"
+include "../lib/centerdef.h"
+include "../lib/center.h"
+include "../lib/fitskydef.h"
+include "../lib/fitsky.h"
+include "../lib/photdef.h"
+include "../lib/phot.h"
+include "../lib/fitpsfdef.h"
+include "../lib/fitpsf.h"
+
+# AP1STATS -- Procedure to fetch an apphot string parameter.
+
+procedure ap1stats (ap, param, str, maxch)
+
+pointer	ap		# pointer to apphot structure
+int	param		# parameter
+char	str[ARB]	# string
+int	maxch		# maximum number of characters
+
+pointer	cen, sky, phot, psf
+
+begin
+	cen = AP_PCENTER(ap)
+	sky = AP_PSKY(ap)
+	phot = AP_PPHOT(ap)
+	psf = AP_PPSF(ap)
+
+	switch (param) {
+	case IMNAME:
+	    call strcpy (AP_IMNAME(ap), str, maxch)
+	case IMROOT:
+	    call strcpy (AP_IMROOT(ap), str, maxch)
+	case CLNAME:
+	    call strcpy (AP_CLNAME(ap), str, maxch)
+	case CLROOT:
+	    call strcpy (AP_CLROOT(ap), str, maxch)
+	case PLOTFILE:
+	    call strcpy (AP_PLOTFILE(ap), str, maxch)
+	case OUTNAME:
+	    call strcpy (AP_OUTNAME(ap), str, maxch)
+	case EXPOSURE:
+	    call strcpy (AP_EXPOSURE(ap), str, maxch)
+	case AIRMASS:
+	    call strcpy (AP_AIRMASS(ap), str, maxch)
+	case FILTER:
+	    call strcpy (AP_FILTER(ap), str, maxch)
+	case FILTERID:
+	    call strcpy (AP_FILTERID(ap), str, maxch)
+	case OBSTIME:
+	    call strcpy (AP_OBSTIME(ap), str, maxch)
+	case OTIME:
+	    call strcpy (AP_OTIME(ap), str, maxch)
+
+	case CSTRING:
+	    call strcpy (AP_CSTRING(cen), str, maxch)
+
+	case SSTRING:
+	    call strcpy (AP_SSTRING(sky), str, maxch)
+
+	case APSTRING:
+	    call strcpy (AP_APSTRING(phot), str, maxch)
+	case APERTS:
+	    call strcpy (AP_APSTRING(phot), str, maxch)
+	case PWSTRING:
+	    call strcpy (AP_PWSTRING(phot), str, maxch)
+
+	case PSFSTRING:
+	    call strcpy (AP_PSFSTRING(psf), str, maxch)
+
+	default:
+	    call error (0, "APSTATS: Unknown apphot string parameter")
+	}
+end
+
+
+# AP1STATI -- Procedure to set an integer apphot parameter.
+
+int procedure ap1stati (ap, param)
+
+pointer	ap		# pointer to apphot structure
+int	param		# parameter
+
+pointer cen, sky, phot, psf
+
+begin
+	cen = AP_PCENTER(ap)
+	sky = AP_PSKY(ap)
+	phot = AP_PPHOT(ap)
+	psf = AP_PPSF(ap)
+
+	switch (param) {
+	case POSITIVE:
+	    return (AP_POSITIVE(ap))
+	case WCSIN:
+	    return (AP_WCSIN(ap))
+	case WCSOUT:
+	    return (AP_WCSOUT(ap))
+	case MW:
+	    return (AP_MW(ap))
+	case CTIN:
+	    return (AP_CTIN(ap))
+	case CTOUT:
+	    return (AP_CTOUT(ap))
+
+	case CENTERFUNCTION:
+	    return (AP_CENTERFUNCTION(cen))
+	case CLEAN:
+	    return (AP_CLEAN(cen))
+	case CMAXITER:
+	    return (AP_CMAXITER(cen))
+
+	case SKYFUNCTION:
+	    return (AP_SKYFUNCTION(sky))
+	case SMAXITER:
+	    return (AP_SMAXITER(sky))
+	case SNREJECT:
+	    return (AP_SNREJECT(sky))
+	case SMOOTH:
+	    return (AP_SMOOTH(sky))
+	case NSKY:
+	    return (AP_NSKY(sky))
+	case NSKY_REJECT:
+	    return (AP_NSKY_REJECT(sky))
+
+	case NAPERTS:
+	    return (AP_NAPERTS(phot))
+	case PWEIGHTS:
+	    return (AP_PWEIGHTS(phot))
+
+	case MAXNPARS:
+	    return (AP_MAXNPARS(psf))
+	case NPARS:
+	    return (AP_PSFNPARS(psf))
+	case PMAXITER:
+	    return (AP_PMAXITER(psf))
+	case PSFUNCTION:
+	    return (AP_PSFUNCTION(psf))
+	case PNREJECT:
+	    return (AP_PNREJECT(psf))
+
+	default:
+	    call error (0, "APSTATI: Unknown apphot integer parameter")
+	}
+end
+
+
+# AP1STATR -- Procedure to set a real apphot parameter.
+
+real procedure ap1statr (ap, param)
+
+pointer	ap		# pointer to apphot structure
+int	param		# parameter
+
+pointer	cen, sky, phot, psf
+
+begin
+	cen = AP_PCENTER(ap)
+	sky = AP_PSKY(ap)
+	phot = AP_PPHOT(ap)
+	psf = AP_PPSF(ap)
+
+	switch (param) {
+
+	case FWHMPSF:
+	    return (AP_FWHMPSF(ap))
+	case SCALE:
+	    return (AP_SCALE(ap))
+	case WX:
+	    return (AP_WX(ap))
+	case WY:
+	    return (AP_WY(ap))
+	case ITIME:
+	    return (AP_ITIME(ap))
+	case CWX:
+	    return (AP_CWX(ap))
+	case CWY:
+	    return (AP_CWY(ap))
+	case DATAMIN:
+	    return (AP_DATAMIN(ap))
+	case DATAMAX:
+	    return (AP_DATAMAX(ap))
+	case XAIRMASS:
+	    return (AP_XAIRMASS(ap))
+
+	case CDATALIMIT:
+	    return (AP_CDATALIMIT(cen))
+	case XSHIFT:
+	    return (AP_XSHIFT(cen))
+	case YSHIFT:
+	    return (AP_YSHIFT(cen))
+	case OXSHIFT:
+	    return (AP_OXSHIFT(cen))
+	case OYSHIFT:
+	    return (AP_OYSHIFT(cen))
+	case CXCUR:
+	    return (AP_CXCUR(cen))
+	case CYCUR:
+	    return (AP_CYCUR(cen))
+	case CAPERT:
+	    return (AP_CAPERT(cen))
+	case CTHRESHOLD:
+	    return (AP_CTHRESHOLD(cen))
+	case MAXSHIFT:
+	    return (AP_MAXSHIFT(cen))
+	case MINSNRATIO:
+	    return (AP_MINSNRATIO(cen))
+	case RCLEAN:
+	    return (AP_RCLEAN(cen))
+	case RCLIP:
+	    return (AP_RCLIP(cen))
+	case SIGMACLEAN:
+	    return (AP_SIGMACLEAN(cen))
+
+	case OXINIT:
+	    return (AP_OXINIT(cen))
+	case OYINIT:
+	    return (AP_OYINIT(cen))
+	case XCENTER:
+	    return (AP_XCENTER(cen))
+	case YCENTER:
+	    return (AP_YCENTER(cen))
+	case OXCENTER:
+	    return (AP_OXCENTER(cen))
+	case OYCENTER:
+	    return (AP_OYCENTER(cen))
+	case XERR:
+	    return (AP_XERR(cen))
+	case YERR:
+	    return (AP_YERR(cen))
+
+	case SXCUR:
+	    return (AP_SXCUR(sky))
+	case SYCUR:
+	    return (AP_SYCUR(sky))
+	case OSXCUR:
+	    return (AP_OSXCUR(sky))
+	case OSYCUR:
+	    return (AP_OSYCUR(sky))
+	case ANNULUS:
+	    return (AP_ANNULUS(sky))
+	case DANNULUS:
+	    return (AP_DANNULUS(sky))
+	case K1:
+	    return (AP_K1(sky))
+	case SLOREJECT:
+	    return (AP_SLOREJECT(sky))
+	case SHIREJECT:
+	    return (AP_SHIREJECT(sky))
+	case SLOCLIP:
+	    return (AP_SLOCLIP(sky))
+	case SHICLIP:
+	    return (AP_SHICLIP(sky))
+	case BINSIZE:
+	    return (AP_BINSIZE(sky))
+	case RGROW:
+	    return (AP_RGROW(sky))
+	case SKY_BACKGROUND:
+	    return (AP_SKYBACKGROUND(sky))
+	case SKY_MODE:
+	    return (AP_SKY_MODE(sky))
+	case SKY_SIGMA:
+	    return (AP_SKY_SIG(sky))
+	case SKY_SKEW:
+	    return (AP_SKY_SKEW(sky))
+
+	case PXCUR:
+	    return (AP_PXCUR(phot))
+	case PYCUR:
+	    return (AP_PYCUR(phot))
+	case OPXCUR:
+	    return (AP_OPXCUR(phot))
+	case OPYCUR:
+	    return (AP_OPYCUR(phot))
+	case ZMAG:
+	    return (AP_ZMAG(phot))
+
+	case PSFAPERT:
+	    return (AP_PSFAPERT(psf))
+	case PK2:
+	    return (AP_PK2(psf))
+	case PFXCUR:
+	    return (AP_PFXCUR(psf))
+	case PFYCUR:
+	    return (AP_PFYCUR(psf))
+	case OPFXCUR:
+	    return (AP_OPFXCUR(psf))
+	case OPFYCUR:
+	    return (AP_OPFYCUR(psf))
+
+	default:
+	    call error (0, "APSTATR: Unknown apphot real parameter")
+	}
+end
+
+
+# AP1STATD -- Procedure to set a double apphot parameter.
+
+double procedure ap1statd (ap, param)
+
+pointer	ap		# pointer to apphot structure
+int	param		# parameter
+
+pointer	cen, sky, phot, psf
+
+begin
+	cen = AP_PCENTER(ap)
+	sky = AP_PSKY(ap)
+	phot = AP_PPHOT(ap)
+	psf = AP_PPSF(ap)
+
+	switch (param) {
+	default:
+	    call error (0, "APSTATD: Unknown apphot double parameter")
+	}
+end