Initial commit

1 files changed, 258 insertions, 0 deletions

diff --git a/noao/astcat/src/attools/atoutpars.x b/noao/astcat/src/attools/atoutpars.x
new file mode 100644
index 00000000..6c2bf213
--- /dev/null
+++ b/noao/astcat/src/attools/atoutpars.x
@@ -0,0 +1,258 @@
+include "../../lib/astrom.h"
+include "../../lib/acatalog.h"
+include "../../lib/aimpars.h"
+
+
+# AT_PRCPSET -- Write the current parameter values out to the region
+# parameters set. 
+
+procedure at_prcpset (psetname, at)
+
+char    psetname[ARB]           #I the parameter set name
+pointer at                      #I the pointer to the main astrom structure
+
+pointer	sp, str, pp
+int	ival
+double	at_statd()
+pointer	clopset()
+int	at_wrdstr(), at_stati()
+
+begin
+        call smark (sp)
+	call salloc (str, SZ_FNAME, TY_CHAR)
+
+	pp = clopset ("psetname")
+
+	# Update the field center and field width parameters.
+	call clppsetd (pp, "rcra", at_statd (at, RCRA))
+	call clppsetd (pp, "rcdec", at_statd (at, RCDEC))
+	call clppsetd (pp, "rrawidth", at_statd (at, RCRAWIDTH))
+	call clppsetd (pp, "rdecwidth", at_statd (at, RCDECWIDTH))
+
+	# Update the units parameters.
+	ival = at_stati (at, RCRAUNITS)
+	if (ival <= 0)
+	    Memc[str] = EOS
+	else if (at_wrdstr (ival,  Memc[str], SZ_FNAME, AT_RA_UNITS) <= 0)
+	    Memc[str] = EOS
+	call clppset (pp, "rcraunits", Memc[str])
+
+	ival = at_stati (at, RCDECUNITS)
+	if (ival <= 0)
+	    Memc[str] = EOS
+	else if (at_wrdstr (ival,  Memc[str], SZ_FNAME, AT_DEC_UNITS) <= 0)
+	    Memc[str] = EOS
+	call clppset (pp, "rcdecunits", Memc[str])
+
+	# Update the celestial coordinate system.
+	call at_stats (at, RCSYSTEM, Memc[str], SZ_FNAME)
+	call clppset (pp, "rcsystem", Memc[str])
+
+	call clcpset (pp)
+
+	call sfree (sp)
+end
+
+
+# AT_PFSPSET -- Write the current parameter values out to the filtering 
+# parameters parameter set. 
+
+procedure at_pfspset (psetname, at)
+
+char    psetname[ARB]           #I the parameter set name
+pointer at                      #I the pointer to the main astrom structure
+
+pointer	sp, str, pp
+int	ival
+pointer	clopset()
+int	at_stati(), at_wrdstr()
+bool	itob()
+
+begin
+        call smark (sp)
+	call salloc (str, SZ_LINE, TY_CHAR)
+
+	pp = clopset ("psetname")
+
+	call at_stats (pp, FIELDS, Memc[str], SZ_LINE)
+	call clppset (pp, "fields", Memc[str])
+	call at_stats (pp, FEXPR, Memc[str], SZ_LINE)
+	call clppset (pp, "fexpr", Memc[str])
+	call at_stats (pp, FNAMES, Memc[str], SZ_LINE)
+	call clppset (pp, "fnames", Memc[str])
+	call at_stats (pp, FNTYPES, Memc[str], SZ_LINE)
+	call clppset (pp, "fntypes", Memc[str])
+	call at_stats (pp, FNUNITS, Memc[str], SZ_LINE)
+	call clppset (pp, "fnunits", Memc[str])
+	call at_stats (pp, FNFORMATS, Memc[str], SZ_LINE)
+	call clppset (pp, "fnformats", Memc[str])
+
+	call at_stats (pp, FSORT, Memc[str], SZ_LINE)
+	call clppset (pp, "fsort", Memc[str])
+	call clppsetb (pp, "freverse", itob(at_stati(at,FREVERSE)))
+
+	#call clppsetb (pp, "freplace", itob(at_stati(at,FREPLACE)))
+	call at_stats (pp, FOSYSTEM, Memc[str], SZ_LINE)
+	call clppset (pp, "fosystem", Memc[str])
+	call at_stats (pp, FIRA, Memc[str], SZ_LINE)
+	call clppset (pp, "fira", Memc[str])
+	call at_stats (pp, FIDEC, Memc[str], SZ_LINE)
+	call clppset (pp, "fidec", Memc[str])
+	ival = at_wrdstr (at_stati(at,FORAUNITS),  Memc[str], SZ_FNAME,
+	    AT_RA_UNITS)
+	if (ival <= 0)
+	    call clppset (pp, "foraunits", "") 
+	else
+	    call clppset (pp, "foraunits", Memc[str])
+	ival = at_wrdstr (at_stati(at,FODECUNITS),  Memc[str], SZ_FNAME,
+	    AT_DEC_UNITS)
+	if (ival <= 0)
+	    call clppset (pp, "fodecunits", "") 
+	else
+	    call clppset (pp, "fodecunits", Memc[str])
+	call at_stats (pp, FORAFORMAT, Memc[str], SZ_LINE)
+	call clppset (pp, "foraformats", Memc[str])
+	call at_stats (pp, FODECFORMAT, Memc[str], SZ_LINE)
+	call clppset (pp, "fodecformats", Memc[str])
+
+	call at_stats (pp, FIXP, Memc[str], SZ_LINE)
+	call clppset (pp, "fixp", Memc[str])
+	call at_stats (pp, FIYP, Memc[str], SZ_LINE)
+	call clppset (pp, "fiyp", Memc[str])
+	call at_stats (pp, FIXC, Memc[str], SZ_LINE)
+	call clppset (pp, "fixc", Memc[str])
+	call at_stats (pp, FIYC, Memc[str], SZ_LINE)
+	call clppset (pp, "fiyc", Memc[str])
+	call at_stats (pp, FOXFORMAT, Memc[str], SZ_LINE)
+	call clppset (pp, "foxformat", Memc[str])
+	call at_stats (pp, FOYFORMAT, Memc[str], SZ_LINE)
+	call clppset (pp, "foyformat", Memc[str])
+
+	call clcpset (pp)
+
+	call sfree (sp)
+end
+
+
+# AT_PWCPSET -- Write the current parameter values out to the default WCS
+# parameters parameter set. 
+
+procedure at_pwcpset (psetname, at)
+
+char    psetname[ARB]           #I the parameter set name
+pointer at                      #I the pointer to the main astrom structure
+
+pointer	pp, st, sym
+pointer	clopset(), at_statp(), stfind()
+
+
+begin
+	pp = clopset ("psetname")
+	st = at_statp (at, WCST)
+
+	sym = stfind (st, "wxref")
+	call clppset (pp, "wxref", AT_WCSTKVAL(sym))
+
+	sym = stfind (st, "wyref")
+	call clppset (pp, "wyref", AT_WCSTKVAL(sym))
+
+	sym = stfind (st, "wxmag")
+	call clppset (pp, "wxmag", AT_WCSTKVAL(sym))
+
+	sym = stfind (st, "wymag")
+	call clppset (pp, "wymag", AT_WCSTKVAL(sym))
+
+	sym = stfind (st, "wxrot")
+	call clppset (pp, "wxrot", AT_WCSTKVAL(sym))
+
+	sym = stfind (st, "wyrot")
+	call clppset (pp, "wyrot", AT_WCSTKVAL(sym))
+
+	sym = stfind (st, "wraref")
+	call clppset (pp, "wraref", AT_WCSTKVAL(sym))
+
+	sym = stfind (st, "wdecref")
+	call clppset (pp, "wdecref", AT_WCSTKVAL(sym))
+
+	sym = stfind (st, "wraunits")
+	call clppset (pp, "wraunits", AT_WCSTKVAL(sym))
+
+	sym = stfind (st, "wdecunits")
+	call clppset (pp, "wdecunits", AT_WCSTKVAL(sym))
+
+	sym = stfind (st, "wproj")
+	call clppset (pp, "wproj", AT_WCSTKVAL(sym))
+
+	sym = stfind (st, "wsystem")
+	call clppset (pp, "wsystem", AT_WCSTKVAL(sym))
+
+	call clcpset (pp)
+end
+
+
+# AT_PIMPSET -- Write the current parameter values out to the default WCS
+# parameters parameter set. 
+
+procedure at_pimpset (psetname, at)
+
+char    psetname[ARB]           #I the parameter set name
+pointer at                      #I the pointer to the main astrom structure
+
+pointer	pp, st, sym
+pointer	clopset(), at_statp(), stfind()
+
+
+begin
+	pp = clopset ("psetname")
+	st = at_statp (at, IMST)
+
+	sym = stfind (st, "esitelng")
+	call clppset (pp, "esitelng", AT_IMSTKVAL(sym))
+
+	sym = stfind (st, "esitelat")
+	call clppset (pp, "esitelat", AT_IMSTKVAL(sym))
+
+	sym = stfind (st, "emjdobs")
+	call clppset (pp, "emjdobs", AT_IMSTKVAL(sym))
+
+	#sym = stfind (st, "ut")
+	#call clppset (pp, "ut", AT_IMSTKVAL(sym))
+
+	sym = stfind (st, "esitealt")
+	call clppset (pp, "esitealt", AT_IMSTKVAL(sym))
+
+	sym = stfind (st, "esitetz")
+	call clppset (pp, "esitetz", AT_IMSTKVAL(sym))
+
+	#sym = stfind (st, "exposure")
+	#call clppset (pp, "exposure", AT_IMSTKVAL(sym))
+
+	sym = stfind (st, "edatamin")
+	call clppset (pp, "edatamin", AT_IMSTKVAL(sym))
+
+	sym = stfind (st, "edatamax")
+	call clppset (pp, "edatamax", AT_IMSTKVAL(sym))
+
+	sym = stfind (st, "egain")
+	call clppset (pp, "egain", AT_IMSTKVAL(sym))
+
+	sym = stfind (st, "erdnoise")
+	call clppset (pp, "erdnoise", AT_IMSTKVAL(sym))
+
+	sym = stfind (st, "ewavlen")
+	call clppset (pp, "ewavlen", AT_IMSTKVAL(sym))
+
+	sym = stfind (st, "etemp")
+	call clppset (pp, "etemp", AT_IMSTKVAL(sym))
+
+	sym = stfind (st, "epress")
+	call clppset (pp, "epress", AT_IMSTKVAL(sym))
+
+	sym = stfind (st, "observat")
+	call clppset (pp, "observat", AT_IMSTKVAL(sym))
+
+	#sym = stfind (st, "dateobs")
+	#call clppset (pp, "dateobs", AT_IMSTKVAL(sym))
+
+	call clcpset (pp)
+end