Initial commit

1 files changed, 123 insertions, 0 deletions

diff --git a/noao/digiphot/apphot/fitpsf/apppsf.x b/noao/digiphot/apphot/fitpsf/apppsf.x
new file mode 100644
index 00000000..c3604c87
--- /dev/null
+++ b/noao/digiphot/apphot/fitpsf/apppsf.x
@@ -0,0 +1,123 @@
+include "../lib/apphotdef.h"
+include "../lib/fitpsfdef.h"
+include "../lib/apphot.h"
+include "../lib/fitpsf.h"
+
+# AP_PPSF -- Procedure to write the results of fitpsf to the output file.
+
+procedure ap_ppsf (ap, fd, id, lid, ier)
+
+pointer	ap	# pointer to apphot structure
+int	fd	# output file descriptor
+int	id	# sequence number of star
+int	lid	# list id of star
+int	ier	# comment string
+
+real	apstatr()
+
+begin
+	# Initialize.
+	if (fd == NULL)
+	    return
+
+	# Print the stars id.
+	call ap_wid (ap, fd, apstatr (ap, OPFXCUR), apstatr (ap, OPFYCUR), id,
+	    lid, '\\')
+
+	# Print the parameters.
+	call ap_wfres (ap, fd, ier)
+end
+
+
+# AP_QPPSF -- Procedure to print the results of fitpsf on the standard output
+# in short form.
+
+procedure ap_qppsf (ap, ier)
+
+pointer	ap	# pointer to apphot structure
+int	ier	# comment string
+
+pointer	sp, imname, psf
+
+begin
+	# Initialize.
+	call smark (sp)
+	call salloc (imname, SZ_FNAME, TY_CHAR)
+	psf = AP_PPSF(ap)
+
+	# Print the parameters on the standard output.
+	call apstats (ap, IMROOT, Memc[imname], SZ_FNAME)
+	switch (AP_PSFUNCTION(psf)) {
+	case AP_RADGAUSS:
+	    call printf (
+	    "%s  %8.2f %8.2f  %6.3f %8g %8g  %s\n") 
+		call pargstr (Memc[imname])
+		call pargr (Memr[AP_PPARS(psf)+1])
+		call pargr (Memr[AP_PPARS(psf)+2])
+		call pargr (Memr[AP_PPARS(psf)+3])
+		call pargr (Memr[AP_PPARS(psf)])
+		call pargr (Memr[AP_PPARS(psf)+4])
+		if (ier != AP_OK)
+		    call pargstr ("err")
+		else
+		    call pargstr ("ok")
+	case AP_ELLGAUSS:
+	    call printf ("%s  %8.2f %8.2f  %6.3f %6.3f  ")
+		call pargstr (Memc[imname])
+		call pargr (Memr[AP_PPARS(psf)+1])
+		call pargr (Memr[AP_PPARS(psf)+2])
+		call pargr (Memr[AP_PPARS(psf)+3])
+		call pargr (Memr[AP_PPARS(psf)+4])
+	    call printf ("%6.1f %8g %8g  %s\n")
+		call pargr (Memr[AP_PPARS(psf)+5])
+		call pargr (Memr[AP_PPARS(psf)])
+		call pargr (Memr[AP_PPARS(psf)+6])
+		if (ier != AP_OK)
+		    call pargstr ("err")
+		else
+		    call pargstr ("ok")
+	case AP_MOMENTS:
+	    call printf ("%s  %8.2f %8.2f  %6.3f %6.3f ")
+		call pargstr (Memc[imname])
+		call pargr (Memr[AP_PPARS(psf)+1])
+		call pargr (Memr[AP_PPARS(psf)+2])
+		call pargr (Memr[AP_PPARS(psf)+3])
+		call pargr (Memr[AP_PPARS(psf)+4])
+	    call printf ("%6.1f  %8g %8g  %s\n")
+		call pargr (Memr[AP_PPARS(psf)+5])
+		call pargr (Memr[AP_PPARS(psf)])
+		call pargr (Memr[AP_PPARS(psf)+6])
+		if (ier != AP_OK)
+		    call pargstr ("err")
+		else
+		    call pargstr ("ok")
+	}
+
+	call sfree (sp)
+end
+
+
+# AP_PFHDR -- Procedure to write the banner for the fitpsf output file.
+
+procedure ap_pfhdr (ap, fd)
+
+pointer	ap		# pointer to the apphot strucuture
+int	fd		# output file descriptor
+
+int	apstati()
+
+begin
+	if (fd == NULL)
+	    return
+
+	switch (apstati (ap, PSFUNCTION)) {
+	case AP_RADGAUSS:
+	    call radhdr (ap, fd)
+	case AP_ELLGAUSS:
+	    call elhdr (ap, fd)
+	case AP_MOMENTS:
+	    call momhdr (ap, fd)
+	default:
+	    ;
+	}
+end