Initial commit

1 files changed, 375 insertions, 0 deletions

diff --git a/noao/astcat/src/attools/atshow.x b/noao/astcat/src/attools/atshow.x
new file mode 100644
index 00000000..e3a9b26b
--- /dev/null
+++ b/noao/astcat/src/attools/atshow.x
@@ -0,0 +1,375 @@
+include "../../lib/astrom.h"
+include "../../lib/acatalog.h"
+include "../../lib/aimpars.h"
+
+# AT_IMSHOW -- Print the current default WCS parameters.
+
+procedure at_imshow (at)
+
+pointer	at			#I the astrometry package descriptor
+
+pointer	sp, str1, str2
+double	at_statd()
+real	at_statr()
+pointer	at_statp()
+
+begin
+	if (at_statp(at, PIMPARS) == NULL)
+	    return
+
+	call smark (sp)
+	call salloc (str1, SZ_FNAME, TY_CHAR)
+	call salloc (str2, SZ_FNAME, TY_CHAR)
+
+	call printf ("\nDefault Image Data Parameters\n")
+	call at_stats (at, OBSERVAT, Memc[str1], SZ_FNAME)
+	call printf ("    The observatory id: %s  MJD: %0.5f\n")
+	    call pargstr (Memc[str1])
+	    call pargd (at_statd (at, EMJDOBS))
+	call printf (
+	"    The site longitude and latitude: %0.2h %0.2h (degrees degrees)\n")
+	    call pargd (at_statd(at, ESITELNG))
+	    call pargd (at_statd(at, ESITELAT))
+	call printf (
+	    "    The site altitude and time zone: %0.1f %0.1f (m hours)\n")
+	    call pargr (at_statr(at, ESITEALT))
+	    call pargr (at_statr(at, ESITETZ))
+	call printf (
+	    "    Effective wavelength: %0.2f (microns)\n")
+	    call pargr (at_statr (at, EWAVLEN))
+	call printf (
+	    "    Effective tempaerature and pressure: %0.2f %0.2f (K mbars)\n")
+	    call pargr (at_statr (at, ETEMP))
+	    call pargr (at_statr (at, EPRESS))
+	call printf (
+	    "    Effective gain and readout noise: %0.2f %0.2f (e-/ADU  e-)\n")
+	    call pargr (at_statr (at, EGAIN))
+	    call pargr (at_statr (at, ERDNOISE))
+	call printf (
+	    "    Low and high good data limits: %0.2f %0.2f (ADU ADU)\n")
+	    call pargr (at_statr (at, EDATAMIN))
+	    call pargr (at_statr (at, EDATAMAX))
+
+	call sfree (sp)
+end
+
+
+# AT_WCSHOW -- Print the current default WCS parameters.
+
+procedure at_wcshow (at)
+
+pointer	at			#I the astrometry package descriptor
+
+pointer	sp, str1, str2
+int	ival
+double	at_statd()
+pointer	at_statp()
+int	at_stati(), at_wrdstr()
+
+begin
+	if (at_statp(at, PWCS) == NULL)
+	    return
+
+	call smark (sp)
+	call salloc (str1, SZ_FNAME, TY_CHAR)
+	call salloc (str2, SZ_FNAME, TY_CHAR)
+
+	call printf ("\nDefault WCS Parameters\n")
+	call at_stats (at, WPROJ, Memc[str1], SZ_FNAME)
+	call printf ("    Sky projection geometry: %s\n")
+	    call pargstr (Memc[str1])
+	ival = at_stati (at, WRAUNITS)
+	if (ival <= 0)
+	    #Memc[str1] = EOS
+	    call strcpy ("default", Memc[str1], SZ_FNAME)
+	else if (at_wrdstr (ival, Memc[str1], SZ_FNAME, AT_RA_UNITS) <= 0)
+	    #Memc[str1] = EOS
+	    call strcpy ("default", Memc[str1], SZ_FNAME)
+	ival = at_stati (at, WDECUNITS)
+	if (ival <= 0)
+	    #Memc[str2] = EOS
+	    call strcpy ("default", Memc[str2], SZ_FNAME)
+	else if (at_wrdstr (ival, Memc[str2], SZ_FNAME, AT_DEC_UNITS) <= 0)
+	    #Memc[str2] = EOS
+	    call strcpy ("default", Memc[str2], SZ_FNAME)
+	call printf ("    Reference point: %0.3h %0.2h (%s %s)\n")
+	    call pargd (at_statd(at, WRAREF))
+	    call pargd (at_statd(at, WDECREF))
+	    call pargstr (Memc[str1])
+	    call pargstr (Memc[str2])
+	call printf ("    Reference point: %0.3f %0.3f (pixels pixels)\n")
+	    call pargd (at_statd(at, WXREF))
+	    call pargd (at_statd(at, WYREF))
+	call printf (
+	    "    X and Y scale: %0.3f %0.3f (arcsec/pixel arsec/pixel)\n")
+	    call pargd (at_statd(at, WXMAG))
+	    call pargd (at_statd(at, WYMAG))
+	call printf (
+	    "    X and Y axis rotation: %0.3f %0.3f (degrees degrees)\n")
+	    call pargd (at_statd(at, WXROT))
+	    call pargd (at_statd(at, WYROT))
+	call at_stats (at, WSYSTEM, Memc[str1], SZ_FNAME)
+	call printf ("    System: %s\n")
+	    call pargstr (Memc[str1])
+
+	call sfree (sp)
+end
+
+
+# AT_IOSHOW -- Print the current i/o parameters.
+
+procedure at_ioshow (at)
+
+pointer	at			#I the astrometry package descriptor
+
+pointer	sp, str1, str2
+pointer	at_statp()
+
+begin
+	if (at_statp(at, PIO) == NULL)
+	    return
+
+	call smark (sp)
+	call salloc (str1, SZ_FNAME, TY_CHAR)
+	call salloc (str2, SZ_FNAME, TY_CHAR)
+
+	call printf ("\nI/O Parameters\n")
+
+	call at_stats (at, CATDB, Memc[str1], SZ_FNAME)
+	call printf ("    catdb = %s\n")
+	    call pargstr (Memc[str1])
+	call at_stats (at, CATALOGS, Memc[str1], SZ_FNAME)
+	call at_stats (at, CATNAME, Memc[str2], SZ_FNAME)
+	call printf ("    catalog = %s  catname = %s\n")
+	    call pargstr (Memc[str1])
+	    call pargstr (Memc[str2])
+
+	call at_stats (at, IMDB, Memc[str1], SZ_FNAME)
+	call printf ("    imdb = %s\n")
+	    call pargstr (Memc[str1])
+	call at_stats (at, SURVEYS, Memc[str1], SZ_FNAME)
+	call at_stats (at, SVNAME, Memc[str2], SZ_FNAME)
+	call printf ("    survey = %s  svname = %s\n")
+	    call pargstr (Memc[str1])
+	    call pargstr (Memc[str2])
+
+	call at_stats (at, IMAGES, Memc[str1], SZ_FNAME)
+	call at_stats (at, IMNAME, Memc[str2], SZ_FNAME)
+	call printf ("    images = %s  imname = %s\n")
+	    call pargstr (Memc[str1])
+	    call pargstr (Memc[str2])
+
+	call at_stats (at, INPUT, Memc[str1], SZ_FNAME)
+	call at_stats (at, INFNAME, Memc[str2], SZ_FNAME)
+	call printf ("    input = %s  infname = %s\n")
+	    call pargstr (Memc[str1])
+	    call pargstr (Memc[str2])
+
+	call at_stats (at, OUTPUT, Memc[str1], SZ_FNAME)
+	call at_stats (at, OUTFNAME, Memc[str2], SZ_FNAME)
+	call printf ("    output = %s  outfname = %s\n")
+	    call pargstr (Memc[str1])
+	    call pargstr (Memc[str2])
+
+	call sfree (sp)
+end
+
+
+# AT_RCSHOW -- Print the current field center parameters.
+
+procedure at_rcshow (at)
+
+pointer	at			#I the astrometry package descriptor
+
+pointer	sp, str1, str2
+int	ival
+double	at_statd()
+pointer	at_statp()
+int	at_wrdstr(), at_stati()
+
+begin
+	if (at_statp(at, PRCENTER) == NULL)
+	    return
+
+	call smark (sp)
+	call salloc (str1, SZ_FNAME, TY_CHAR)
+	call salloc (str2, SZ_FNAME, TY_CHAR)
+
+	call printf ("\nField Center Parameters\n")
+	call at_stats (at, RCSOURCE, Memc[str1], SZ_FNAME)
+	call at_stats (at, RCSYSTEM, Memc[str2], SZ_FNAME)
+	call printf ("    rcsource = %s  rcsystem = %s\n")
+	    call pargstr (Memc[str1])
+	    call pargstr (Memc[str2])
+
+	call printf ("    rcra = %h  rcdec = %h\n")
+	    call pargd (at_statd (at, RCRA))
+	    call pargd (at_statd (at, RCDEC))
+	call printf ("    rrawidth = %0.2f  rdecwidth = %0.2f\n")
+	    call pargd (at_statd (at, RCRAWIDTH))
+	    call pargd (at_statd (at, RCDECWIDTH))
+
+	ival = at_stati (at, RCRAUNITS)
+	if (ival <= 0)
+	    #Memc[str1] = EOS
+	    call strcpy ("default", Memc[str1], SZ_FNAME)
+	else if (at_wrdstr (ival, Memc[str1], SZ_FNAME, AT_RA_UNITS) <= 0)
+	    #Memc[str1] = EOS
+	    call strcpy ("default", Memc[str1], SZ_FNAME)
+	ival = at_stati (at, RCDECUNITS)
+	if (ival <= 0)
+	    #Memc[str2] = EOS
+	    call strcpy ("default", Memc[str2], SZ_FNAME)
+	else if (at_wrdstr (ival, Memc[str2], SZ_FNAME, AT_DEC_UNITS) <= 0)
+	    #Memc[str2] = EOS
+	    call strcpy ("default", Memc[str2], SZ_FNAME)
+	call printf ("    rcraunits = %s  rcdecunits = %s\n")
+	    call pargstr (Memc[str1])
+	    call pargstr (Memc[str2])
+
+	call printf ("\n")
+
+	call sfree (sp)
+end
+
+
+# AT_FSSHOW -- Print the current filtering parameters.
+
+procedure at_fsshow (at)
+
+pointer	at			#I the astrometry package descriptor
+
+pointer	sp, str1, str2
+int	ival
+pointer	at_statp()
+int	at_stati(), at_wrdstr()
+bool	itob()
+
+begin
+	if (at_statp(at, PFILTER) == NULL)
+	    return
+
+	call smark (sp)
+	call salloc (str1, SZ_FNAME, TY_CHAR)
+	call salloc (str2, SZ_FNAME, TY_CHAR)
+
+	call printf ("\nRecord Filtering Parameters\n")
+	call at_stats (at, FIELDS, Memc[str1], SZ_FNAME)
+	call at_stats (at, FEXPR, Memc[str2], SZ_FNAME)
+	call printf ("    fields = %s  fexpr = %s\n")
+	    call pargstr (Memc[str1])
+	    call pargstr (Memc[str2])
+	call at_stats (at, FNAMES, Memc[str1], SZ_FNAME)
+	call at_stats (at, FNTYPES, Memc[str2], SZ_FNAME)
+	call printf ("    fnames = %s fntypes = %s\n")
+	    call pargstr (Memc[str1])
+	    call pargstr (Memc[str2])
+	call at_stats (at, FNUNITS, Memc[str1], SZ_FNAME)
+	call at_stats (at, FNFORMATS, Memc[str2], SZ_FNAME)
+	call printf ("    fnunits = %s fnformats = %s\n")
+	    call pargstr (Memc[str1])
+	    call pargstr (Memc[str2])
+
+	call at_stats (at, FSORT, Memc[str1], SZ_FNAME)
+	call printf ("    fsort = %s  freverse = %b\n")
+	    call pargstr (Memc[str1])
+	    call pargb (itob(at_stati(at, FREVERSE)))
+
+	call at_stats (at, FOSYSTEM, Memc[str1], SZ_FNAME)
+	call printf ("    freplace = %b  fosystem = %s\n")
+	    call pargb (itob(at_stati(at, FREPLACE)))
+	    call pargstr (Memc[str1])
+	call at_stats (at, FIRA, Memc[str1], SZ_FNAME)
+	call at_stats (at, FIDEC, Memc[str2], SZ_FNAME)
+	call printf ("    fira = %s  fidec = %s\n")
+	    call pargstr (Memc[str1])
+	    call pargstr (Memc[str2])
+
+	ival = at_stati (at, FORAUNITS)
+	if (ival <= 0)
+	    Memc[str1] = EOS
+	else if (at_wrdstr (ival, Memc[str1], SZ_FNAME, AT_RA_UNITS) <= 0)
+	    Memc[str1] = EOS
+	ival = at_stati (at, FODECUNITS)
+	if (ival <= 0)
+	    Memc[str2] = EOS
+	else if (at_wrdstr (ival, Memc[str2], SZ_FNAME, AT_DEC_UNITS) <= 0)
+	    Memc[str2] = EOS
+	call printf ("    foraunits = %s  fodecunits = %s\n")
+	    call pargstr (Memc[str1])
+	    call pargstr (Memc[str2])
+	call at_stats (at, FORAFORMAT, Memc[str1], SZ_FNAME)
+	call at_stats (at, FODECFORMAT, Memc[str2], SZ_FNAME)
+	call printf ("    foraformat = %s  fodecformat = %s\n")
+	    call pargstr (Memc[str1])
+	    call pargstr (Memc[str2])
+
+	call at_stats (at, FIXP, Memc[str1], SZ_FNAME)
+	call at_stats (at, FIYP, Memc[str2], SZ_FNAME)
+	call printf ("    fixp = %s  fiyp = %s\n")
+	    call pargstr (Memc[str1])
+	    call pargstr (Memc[str2])
+	call at_stats (at, FIXC, Memc[str1], SZ_FNAME)
+	call at_stats (at, FIYC, Memc[str2], SZ_FNAME)
+	call printf ("    fixc = %s  fiyc = %s\n")
+	    call pargstr (Memc[str1])
+	    call pargstr (Memc[str2])
+	call at_stats (at, FOXFORMAT, Memc[str1], SZ_FNAME)
+	call at_stats (at, FOYFORMAT, Memc[str2], SZ_FNAME)
+	call printf ("    foxformat = %s  foyformat = %s\n")
+	    call pargstr (Memc[str1])
+	    call pargstr (Memc[str2])
+
+	call printf ("\n")
+
+	call sfree (sp)
+end
+
+
+# AT_STSHOW -- Dump the contents of the field center symbol table
+
+procedure at_stshow (at)
+
+pointer	at			#I the astrometry package descriptor
+
+pointer	st, symlist, symbol
+int	i, nfields
+pointer	at_statp(), sthead(), stnext()
+int	stnsymbols()
+
+begin
+	if (at_statp(at, PRCENTER) == NULL)
+	    return
+
+	st = at_statp(at, RCST)
+	if (st == NULL)
+	    return
+
+	nfields = stnsymbols (st, 0)
+	if (nfields <= 0)
+	    return
+
+	call malloc (symlist, nfields, TY_INT)
+	symbol = sthead (st)
+	do i = 1, nfields {
+	    Memi[symlist+i-1] = symbol
+	    symbol = stnext (st, symbol)
+	}
+
+	do i = nfields, 1, -1 {
+	    symbol = Memi[symlist+i-1]
+	    call printf ("%s %s\n")
+		call pargstr (AT_RCSTSOURCE(symbol))
+		call pargstr (AT_RCSTNAME(symbol))
+	    call printf ("    %h %h  %8.3f %8.3f  %d %d  %s\n")
+		call pargd (AT_RCSTRA(symbol))
+		call pargd (AT_RCSTDEC(symbol))
+		call pargd (AT_RCSTRAWIDTH(symbol))
+		call pargd (AT_RCSTDECWIDTH(symbol))
+		call pargi (AT_RCSTRAUNITS(symbol))
+		call pargi (AT_RCSTDECUNITS(symbol))
+		call pargstr (AT_RCSTSYSTEM(symbol))
+	}
+
+	call mfree (symlist, TY_INT)
+end