Initial commit

1 files changed, 289 insertions, 0 deletions

diff --git a/noao/astcat/src/asttest.cl b/noao/astcat/src/asttest.cl
new file mode 100644
index 00000000..23292172
--- /dev/null
+++ b/noao/astcat/src/asttest.cl
@@ -0,0 +1,289 @@
+# ASTTEST - Self testing procedure for the ASTCAT package.
+
+procedure asttest (rootname)
+
+string	rootname	{prompt="Root name of the output test files"}
+string	astlogfile	{"", prompt="Name of the output log file"}
+
+struct	*catlist
+struct	*svlist
+
+begin
+	# Declare local variables.
+	string	troot, tastlog, tcatlist, tcatalog, tcatfile
+	string	tsvlist, tsurvey, timage, ttemp1
+	int tfirst, tlast, tindex
+
+	# Check that the user truly wants to proceed.
+	ttemp1 = ""
+	print ("")
+	print ("ASTTEST initializes the ASTCAT task parameters")
+	print ("Type 'q' or 'Q' to quit, any other key to proceed")
+	if (scan (ttemp1) != EOF) {
+	    if (ttemp1 == "q" || ttemp1 == "Q") {
+		print ("Terminating the ASTTEST task")
+		bye
+	    }
+	}
+	print ("")
+
+	# Define the plot file name.
+
+	troot = rootname
+	tastlog = astlogfile
+	if (tastlog == "") {
+	    tastlog = troot // ".log"
+	}
+
+	# Check for the existence of the test files.
+
+	if (access (tastlog)) {
+	    error (0, "Error: The log file already exists on disk")
+	}
+
+	# Create some temporary files.
+
+	tcatlist = mktemp (troot)
+	tsvlist = mktemp (troot)
+
+	# Initialize the ASTCAT package.
+
+	print ("Initialize the ASTCAT package", >> tastlog)
+	print ("", >> tastlog)
+	print ("")
+	print ("Initialize the ASTCAT package")
+	print ("")
+
+	astcat.catalogs = "usno2@noao"
+	astcat.catdb = "astcat$lib/catdb.dat"
+	astcat.imsurveys = "dss2@cadc"
+	astcat.imdb = "astcat$lib/imdb.dat"
+
+	# Unlearning tasks and psets. Ran into a problem with doing
+	# unlearn on astcat which I will have to work around.
+
+	# unlearn ("astcat")
+	unlearn ("acatpars")
+	unlearn ("aclist")
+	unlearn ("adumpcat")
+	unlearn ("adumpim")
+	unlearn ("afiltcat")
+	unlearn ("afiltpars")
+	unlearn ("agetcat")
+	unlearn ("agetim")
+	unlearn ("ahedit")
+	unlearn ("aimfind")
+	unlearn ("aimpars")
+	unlearn ("aregpars")
+	unlearn ("aslist")
+	unlearn ("awcspars")
+
+	# Do the tests.
+
+	# Test the ACLIST task.
+
+	print ("", >> tastlog)
+	print ("Testing the ACLIST task", >> tastlog)
+	print ("Testing the ACLIST task")
+	print ("", >> tastlog)
+
+	aclist ("*", verbose+, catdb=")_.catdb", >> tastlog)
+	aclist ("*", verbose-, catdb=")_.catdb", > tcatlist)
+	print ("", >> tastlog)
+
+	# Test the ADUMPCAT task.
+
+	print ("", >> tastlog)
+	print ("Testing the ADUMPCAT task", >> tastlog)
+	print ("Testing the ADUMPCAT task")
+	print ("", >> tastlog)
+
+	catlist = tcatlist
+	while (fscan (catlist, tcatalog) != EOF) {
+	    print ("", >> tastlog)
+	    print (tcatalog, >> tastlog)
+	    print ("", >> tastlog)
+	    print ("    ", tcatalog)
+	    if (tcatalog == "tmass@ipac") {
+		adumpcat (tcatalog, "STDOUT", ra="00:00:00.0", dec="00:00:00",
+		    size="0.17", catdb=")_.catdb", >> tastlog)
+	    } else {
+		adumpcat (tcatalog, "STDOUT", ra="00:00:00.0", dec="00:00:00",
+		    size="10.0", catdb=")_.catdb", >> tastlog)
+	    }
+	    print ("", >> tastlog)
+	}
+
+	# Test the AGETCAT and AFILTCAT tasks.
+
+	print ("", >> tastlog)
+	print ("Testing the AGETCAT and AFILTCAT tasks", >> tastlog)
+	print ("Testing the AGETCAT and AFILTCAT tasks")
+	print ("", >> tastlog)
+
+	catlist = tcatlist 
+	while (fscan (catlist, tcatalog) != EOF) {
+	    print ("", >> tastlog)
+	    print (tcatalog, >> tastlog)
+	    print ("", >> tastlog)
+	    print ("    ", tcatalog)
+	    tfirst = 1
+	    tlast = strlen (tcatalog)
+	    tindex = stridx ("@", tcatalog)
+	    tcatfile = substr (tcatalog, tfirst, tindex - 1) // "." //
+		substr (tcatalog, tindex + 1, tlast) // ".cat"
+	    if (access (tcatfile)) {
+		delete (tcatfile, verify-)
+	    }
+	    agetcat ("pars", tcatfile, rcra="00:00:00.00", rcdec="+00:00:00.0",
+	        rrawidth=20.0, rdecwidth=20.0, catalogs=tcatalog, standard+,
+		filter-, update-, verbose+, catdb=")_.catdb", >> tastlog)
+	    print ("", >> tastlog)
+	    type (tcatfile, map_cc+, device="terminal", >> tastlog)
+	    print ("", >> tastlog)
+	    afiltcat (tcatfile, tcatfile, catalogs="filename@noao", standard+,
+	        filter+, fsort="mag1", update-, verbose+, catdb=")_.catdb",
+		>> tastlog)
+	    print ("", >> tastlog)
+	    type (tcatfile, map_cc+, device="terminal", >> tastlog)
+	    print ("", >> tastlog)
+	    delete (tcatfile, verify-)
+	}
+
+	# Test the ASLIST task.
+
+	print ("", >> tastlog)
+	print ("Testing the ASLIST task", >> tastlog)
+	print ("Testing the ASLIST task")
+	print ("", >> tastlog)
+
+	aslist ("*", verbose+, imdb=")_.imdb", >> tastlog)
+	aslist ("*", verbose-, imdb=")_.imdb", > tsvlist)
+	print ("", >> tastlog)
+
+	# Test the ADUMPIM task.
+
+	print ("", >> tastlog)
+	print ("Testing the ADUMPIM task", >> tastlog)
+	print ("Testing the ADUMPIM task")
+	print ("", >> tastlog)
+
+	svlist = tsvlist
+	while (fscan (svlist, tsurvey) != EOF) {
+	    print ("", >> tastlog)
+	    print (tsurvey, >> tastlog)
+	    print ("", >> tastlog)
+	    print ("    ", tsurvey)
+	    tfirst = 1
+	    tlast = strlen (tsurvey)
+	    tindex = stridx ("@", tsurvey)
+	    timage = substr (tsurvey, tfirst, tindex - 1) // "." //
+		substr (tsurvey, tindex + 1, tlast) // ".fits"
+	    if (imaccess (timage)) {
+	        imdelete (timage, verify-)
+	    }
+	    adumpim (tsurvey, timage, ra="14:28:07.0", dec="+34:55:00",
+		    size="10.0", imdb=")_.imdb", >> tastlog)
+	    printf ("    ")
+	    imheader (timage, longheader-, userfields+)
+	    print ("", >> tastlog)
+	    imheader (timage, longheader+, userfields+, >> tastlog)
+	    print ("", >> tastlog)
+	    imdelete (timage, verify-)
+	}
+
+	# Test the AGETIM and AHEDIT tasks.
+
+	print ("", >> tastlog)
+	print ("Testing the AGETIM and AHEDIT tasks", >> tastlog)
+	print ("Testing the AGETIM and AHEDIT tasks")
+	print ("", >> tastlog)
+
+	svlist = tsvlist
+	while (fscan (svlist, tsurvey) != EOF) {
+	    print ("", >> tastlog)
+	    print (tsurvey, >> tastlog)
+	    print ("", >> tastlog)
+	    tfirst = 1
+	    tlast = strlen (tsurvey)
+	    tindex = stridx ("@", tsurvey)
+	    timage = substr (tsurvey, tfirst, tindex - 1) // "." //
+		substr (tsurvey, tindex + 1, tlast) // ".fits"
+	    if (imaccess (timage)) {
+	        imdelete (timage, verify-)
+	    }
+	    agetim ("pars", timage, rcra="14:28:07.00", rcdec="+34:55:00.0",
+	        rrawidth=10.0, rdecwidth=10.0, imsurveys=tsurvey, wcsedit-,
+		hdredit-, update-, verbose+, imdb=")_.imdb", >> tastlog)
+	    printf ("    ")
+	    imheader (timage, longheader-, userfields+)
+	    print ("", >> tastlog)
+	    imheader (timage, longheader+, userfields+, >> tastlog)
+	    print ("", >> tastlog)
+	    ahedit (timage, tsurvey, hupdate+, wcsedit+, wcs="none", hdredit+,
+	        update-, verbose+, imdb=")_.imdb", >> tastlog)
+	    print ("", >> tastlog)
+	    imheader (timage, longheader+, userfields+, >> tastlog)
+	    #imdelete (timage, verify-)
+	}
+
+	# Test the AIMFIND task.
+
+	print ("", >> tastlog)
+	print ("Testing the AIMFIND task", >> tastlog)
+	print ("Testing the AIMFIND task")
+	print ("", >> tastlog)
+
+	# Test the aimfind task using the USNO2 survey
+
+	svlist = tsvlist
+	while (fscan (svlist, tsurvey) != EOF) {
+	    tfirst = 1
+	    tlast = strlen (tsurvey)
+	    tindex = stridx ("@", tsurvey)
+	    timage = substr (tsurvey, tfirst, tindex - 1) // "." //
+		substr (tsurvey, tindex + 1, tlast) // ".fits"
+	    tcatfile = substr (tsurvey, tfirst, tindex - 1) // "." //
+		substr (tsurvey, tindex + 1, tlast) // ".cat"
+	    if (access (tcatfile)) {
+		delete (tcatfile, verify-)
+	    }
+	    aimfind (timage, tcatfile, imfile="", catalogs="usno2@noao",
+	        standard+, filter-, append-, update-, verbose+,
+		catdb=")_.catdb", >> tastlog)
+	}
+
+	# Reinitialize the astcat package.
+
+	# unlearn ("astcat")
+	unlearn ("acatpars")
+	unlearn ("aclist")
+	unlearn ("adumpcat")
+	unlearn ("adumpim")
+	unlearn ("afiltcat")
+	unlearn ("afiltpars")
+	unlearn ("agetcat")
+	unlearn ("agetim")
+	unlearn ("ahedit")
+	unlearn ("aimfind")
+	unlearn ("aimpars")
+	unlearn ("aregpars")
+	unlearn ("aslist")
+	unlearn ("awcspars")
+
+	# Delete some temporary files.
+
+	delete (tsvlist, verify-)
+	delete (tcatlist, verify-)
+	svlist = ""
+	catlist = ""
+
+        print ("", >> tastlog)
+        print ("ASTCAT package tests completed", >> tastlog)
+        print ("", >> tastlog)
+        print ("")
+        print ("ASTCAT package tests completed")
+        print ("")
+
+	bye
+end