Initial commit

1 files changed, 305 insertions, 0 deletions

diff --git a/noao/astcat/src/attools/atdefpars.x b/noao/astcat/src/attools/atdefpars.x
new file mode 100644
index 00000000..0b8b2cb7
--- /dev/null
+++ b/noao/astcat/src/attools/atdefpars.x
@@ -0,0 +1,305 @@
+include "../../lib/astrom.h"
+include "../../lib/acatalog.h"
+include "../../lib/aimpars.h"
+
+
+# AT_DRCPSET -- Reset the region definition parameters to their default values.
+
+procedure at_drcpset (at)
+
+pointer at                      #I the pointer to the main astrom structure
+
+begin
+	# Set the default field center and width.
+	call at_setd (at, RCRA, 0.0d0)
+	call at_setd (at, RCDEC, 0.0d0)
+	call at_setd (at, RCRAWIDTH, 10.0d0)
+	call at_setd (at, RCDECWIDTH, 10.0d0)
+
+	# Set the default field center units.
+	call at_seti (at, RCRAUNITS, 0)
+	call at_seti (at, RCDECUNITS, 0)
+
+	# Set the default field center coordinate system.
+	call at_sets (at, RCSYSTEM, "")
+end
+
+
+# AT_DFSPSET -- Reset the filtering parameters to their default values.
+
+procedure at_dfspset (at)
+
+pointer at                      #I the pointer to the main astrom structure
+
+begin
+	call at_sets (at, FIELDS, "f[*]")
+	call at_sets (at, FEXPR, "yes")
+	call at_sets (at, FNAMES, "")
+	call at_sets (at, FNTYPES, "")
+	call at_sets (at, FNUNITS, "")
+	call at_sets (at, FNFORMATS, "")
+	call at_sets (at, FSORT, "")
+	call at_seti (at, FREVERSE, NO)
+	call at_sets (at, FOSYSTEM, "")
+	call at_sets (at, FIRA, "ra")
+	call at_sets (at, FIDEC, "dec")
+	call at_seti (at, FORAUNITS, 0)
+	call at_seti (at, FODECUNITS, 0)
+	call at_sets (at, FORAFORMAT, "")
+	call at_sets (at, FODECFORMAT, "")
+	call at_sets (at, FIXP, "xp")
+	call at_sets (at, FIYP, "yp")
+	call at_sets (at, FIXC, "xc")
+	call at_sets (at, FIYC, "yc")
+	call at_sets (at, FOXFORMAT, "%10.3f")
+	call at_sets (at, FOYFORMAT, "%10.3f")
+end
+
+
+# AT_DWCPSET -- Reset the wcs parameters to their default values.
+
+procedure at_dwcpset (at)
+
+pointer at                      #I the pointer to the main astrom structure
+
+double	dval
+pointer	st, sym
+int	ip
+pointer at_statp(), stfind()
+int	ctod()
+
+begin
+	st = at_statp (at, WCST)
+	if (st == NULL)
+	    return
+
+	sym = stfind (st, "wxref")
+	if (sym != NULL) {
+	    call strcpy ("INDEF", AT_WCSTKVAL(sym), LEN_WCST_STRUCT)
+	    ip = 1
+	    if (ctod ("INDEF", ip, dval) <= 0)
+	        dval = INDEFD
+	    call at_setd (at, WXREF, dval)
+	}
+
+	sym = stfind (st, "wyref")
+	if (sym != NULL) {
+	    call strcpy ("INDEF", AT_WCSTKVAL(sym), LEN_WCST_STRUCT)
+	    ip = 1
+	    if (ctod ("INDEF", ip, dval) <= 0)
+	        dval = INDEFD
+	    call at_setd (at, WYREF, dval)
+	}
+
+	sym = stfind (st, "wxmag")
+	if (sym != NULL) {
+	    call strcpy ("INDEF", AT_WCSTKVAL(sym), LEN_WCST_STRUCT)
+	    ip = 1
+	    if (ctod ("INDEF", ip, dval) <= 0)
+	        dval = INDEFD
+	    call at_setd (at, WXMAG, dval)
+	}
+
+	sym = stfind (st, "wymag")
+	if (sym != NULL) {
+	    call strcpy ("INDEF", AT_WCSTKVAL(sym), LEN_WCST_STRUCT)
+	    ip = 1
+	    if (ctod ("INDEF", ip, dval) <= 0)
+	        dval = INDEFD
+	    call at_setd (at, WYMAG, dval)
+	}
+
+	sym = stfind (st, "wxrot")
+	if (sym != NULL) {
+	    call strcpy ("180.0", AT_WCSTKVAL(sym), LEN_WCST_STRUCT)
+	    ip = 1
+	    if (ctod ("180.0", ip, dval) <= 0)
+	        dval = INDEFD
+	    call at_setd (at, WXROT, dval)
+	}
+
+	sym = stfind (st, "wyrot")
+	if (sym != NULL) {
+	    call strcpy ("0.0", AT_WCSTKVAL(sym), LEN_WCST_STRUCT)
+	    ip = 1
+	    if (ctod ("0.0", ip, dval) <= 0)
+	        dval = INDEFD
+	    call at_setd (at, WYROT, dval)
+	}
+
+	sym = stfind (st, "wraref")
+	if (sym != NULL) {
+	    call strcpy ("RA", AT_WCSTKVAL(sym), LEN_WCST_STRUCT)
+	    ip = 1
+	    if (ctod ("RA", ip, dval) <= 0)
+	        dval = INDEFD
+	    call at_setd (at, WRAREF, dval)
+	}
+
+	sym = stfind (st, "wdecref")
+	if (sym != NULL) {
+	    call strcpy ("DEC", AT_WCSTKVAL(sym), LEN_WCST_STRUCT)
+	    ip = 1
+	    if (ctod ("DEC", ip, dval) <= 0)
+	        dval = INDEFD
+	    call at_setd (at, WDECREF, dval)
+	}
+
+	sym = stfind (st, "wraunits")
+	if (sym != NULL) {
+	    call strcpy ("", AT_WCSTKVAL(sym), LEN_WCST_STRUCT)
+	    call at_seti (at, WRAUNITS, 0)
+	}
+
+	sym = stfind (st, "wdecunits")
+	if (sym != NULL) {
+	    call strcpy ("", AT_WCSTKVAL(sym), LEN_WCST_STRUCT)
+	    call at_seti (at, WDECUNITS, 0)
+	}
+
+	sym = stfind (st, "wproj")
+	if (sym != NULL) {
+	    call strcpy ("INDEF", AT_WCSTKVAL(sym), LEN_WCST_STRUCT)
+	    call at_sets (at, WPROJ, "tan")
+	}
+
+	sym = stfind (st, "wsystem")
+	if (sym != NULL) {
+	    call strcpy ("EQUINOX", AT_WCSTKVAL(sym), LEN_WCST_STRUCT)
+	    call at_sets (at, WSYSTEM, "EQUINOX")
+	}
+end
+
+
+# AT_DIMPSET -- Read in the default image data parameters.
+
+procedure at_dimpset (at)
+
+pointer at                      #I the pointer to the main astrom structure
+
+double	dval
+real	rval
+pointer	st, sym
+int	ip
+pointer	at_statp(), stfind()
+int	ctod(), ctor()
+
+begin
+	st = at_statp (at, IMST)
+	if (st == NULL)
+	    return
+
+	sym = stfind (st, "esitelng")
+	if (sym != NULL) {
+	    call strcpy ("INDEF", AT_IMSTKVAL(sym), LEN_IMST_STRUCT)
+	    ip = 1
+	    if (ctod ("INDEF", ip, dval) <= 0)
+	        dval = INDEFD
+	    call at_setd (at, ESITELNG, dval)
+	}
+
+	sym = stfind (st, "esitelat")
+	if (sym != NULL) {
+	    call strcpy ("INDEF", AT_IMSTKVAL(sym), LEN_IMST_STRUCT)
+	    ip = 1
+	    if (ctod ("INDEF", ip, dval) <= 0)
+	        dval = INDEFD
+	    call at_setd (at, ESITELAT, dval)
+	}
+
+	sym = stfind (st, "esitealt")
+	if (sym != NULL) {
+	    call strcpy ("INDEF", AT_IMSTKVAL(sym), LEN_IMST_STRUCT)
+	    ip = 1
+	    if (ctor ("INDEF", ip, rval) <= 0)
+	        rval = INDEFR
+	    call at_setr (at, ESITEALT, rval)
+	}
+
+	sym = stfind (st, "esitetz")
+	if (sym != NULL) {
+	    call strcpy ("INDEF", AT_IMSTKVAL(sym), LEN_IMST_STRUCT)
+	    ip = 1
+	    if (ctor ("INDEF", ip, rval) <= 0)
+	        rval = INDEFR
+	    call at_setr (at, ESITETZ, rval)
+	}
+
+	sym = stfind (st, "emjdobs")
+	if (sym != NULL) {
+	    call strcpy ("MJD-OBS", AT_IMSTKVAL(sym), LEN_IMST_STRUCT)
+	    ip = 1
+	    if (ctod ("MJD-OBS", ip, dval) <= 0)
+	        dval = INDEFD
+	    call at_setd (at, EMJDOBS, dval)
+	}
+
+	sym = stfind (st, "edatamin")
+	if (sym != NULL) {
+	    call strcpy ("INDEF", AT_IMSTKVAL(sym), LEN_IMST_STRUCT)
+	    ip = 1
+	    if (ctor ("INDEF", ip, rval) <= 0)
+	        rval = INDEFR
+	    call at_setr (at, EDATAMIN, rval)
+	}
+
+	sym = stfind (st, "edatamax")
+	if (sym != NULL) {
+	    call strcpy ("INDEF", AT_IMSTKVAL(sym), LEN_IMST_STRUCT)
+	    ip = 1
+	    if (ctor ("INDEF", ip, rval) <= 0)
+	        rval = INDEFR
+	    call at_setr (at, EDATAMAX, rval)
+	}
+
+	sym = stfind (st, "egain")
+	if (sym != NULL) {
+	    call strcpy ("GAIN", AT_IMSTKVAL(sym), LEN_IMST_STRUCT)
+	    ip = 1
+	    if (ctor ("GAIN", ip, rval) <= 0)
+	        rval = INDEFR
+	    call at_setr (at, EGAIN, rval)
+	}
+
+	sym = stfind (st, "erdnoise")
+	if (sym != NULL) {
+	    call strcpy ("RDNOISE", AT_IMSTKVAL(sym), LEN_IMST_STRUCT)
+	    ip = 1
+	    if (ctor ("RDNOISE", ip, rval) <= 0)
+	        rval = INDEFR
+	    call at_setr (at, ERDNOISE, rval)
+	}
+
+	sym = stfind (st, "ewavlen")
+	if (sym != NULL) {
+	    call strcpy ("INDEF", AT_IMSTKVAL(sym), LEN_IMST_STRUCT)
+	    ip = 1
+	    if (ctor ("INDEF", ip, rval) <= 0)
+	        rval = INDEFR
+	    call at_setr (at, EWAVLEN, rval)
+	}
+
+	sym = stfind (st, "etemp")
+	if (sym != NULL) {
+	    call strcpy ("INDEF", AT_IMSTKVAL(sym), LEN_IMST_STRUCT)
+	    ip = 1
+	    if (ctor ("INDEF", ip, rval) <= 0)
+	        rval = INDEFR
+	    call at_setr (at, ETEMP, rval)
+	}
+
+	sym = stfind (st, "epress")
+	if (sym != NULL) {
+	    call strcpy ("INDEF", AT_IMSTKVAL(sym), LEN_IMST_STRUCT)
+	    ip = 1
+	    if (ctor ("INDEF", ip, rval) <= 0)
+	        rval = INDEFR
+	    call at_setr (at, EPRESS, rval)
+	}
+
+	sym = stfind (st, "observat")
+	if (sym != NULL) {
+	    call strcpy ("OBSERVAT", AT_WCSTKVAL(sym), LEN_WCST_STRUCT)
+	    call at_sets (at, OBSERVAT, "OBSERVAT")
+	}
+end