Initial commit

1 files changed, 75 insertions, 0 deletions

diff --git a/noao/imred/hydra/dohydra.cl b/noao/imred/hydra/dohydra.cl
new file mode 100644
index 00000000..74264633
--- /dev/null
+++ b/noao/imred/hydra/dohydra.cl
@@ -0,0 +1,75 @@
+# DOHYDRA -- Process HYDRA spectra from 2D to wavelength calibrated 1D.
+#
+# The task PROC does all of the interactive work and BATCH does the
+# background work.  This procedure is organized this way to minimize the
+# dictionary space when the background task is submitted.
+
+procedure dohydra (objects)
+
+string	objects = ""		{prompt="List of object spectra"}
+
+file	apref = ""		{prompt="Aperture reference spectrum"}
+file	flat = ""		{prompt="Flat field spectrum"}
+file	throughput = ""		{prompt="Throughput file or image (optional)"}
+string	arcs1 = ""		{prompt="List of arc spectra"}
+string	arcs2 = ""		{prompt="List of shift arc spectra"}
+file	arcreplace = ""		{prompt="Special aperture replacements"}
+file	arctable = ""		{prompt="Arc assignment table (optional)\n"}
+
+string	readnoise = "RDNOISE"	{prompt="Read out noise sigma (photons)"}
+string	gain = "GAIN"		{prompt="Photon gain (photons/data number)"}
+real	datamax = INDEF		{prompt="Max data value / cosmic ray threshold"}
+int	fibers = 97		{prompt="Number of fibers"}
+real	width = 12.		{prompt="Width of profiles (pixels)"}
+real	minsep = 8.	{prompt="Minimum separation between fibers (pixels)"}
+real	maxsep = 15.	{prompt="Maximum separation between fibers (pixels)"}
+file	apidtable = ""		{prompt="Aperture identifications"}
+string	crval = "INDEF"		{prompt="Approximate central wavelength"}
+string	cdelt = "INDEF"		{prompt="Approximate dispersion"}
+string	objaps = ""		{prompt="Object apertures"}
+string	skyaps = ""		{prompt="Sky apertures"}
+string	arcaps = ""		{prompt="Arc apertures"}
+string	objbeams = "0,1"	{prompt="Object beam numbers"}
+string	skybeams = "0"		{prompt="Sky beam numbers"}
+string	arcbeams = ""		{prompt="Arc beam numbers\n"}
+
+bool	scattered = no		{prompt="Subtract scattered light?"}
+bool	fitflat = yes		{prompt="Fit and ratio flat field spectrum?"}
+bool	clean = yes		{prompt="Detect and replace bad pixels?"}
+bool	dispcor = yes		{prompt="Dispersion correct spectra?"}
+bool	savearcs = yes		{prompt="Save simultaneous arc apertures?"}
+bool	skyalign = no		{prompt="Align sky lines?"}
+bool	skysubtract = yes	{prompt="Subtract sky?"}
+bool	skyedit = yes		{prompt="Edit the sky spectra?"}
+bool	saveskys = yes		{prompt="Save sky spectra?"}
+bool	splot = no		{prompt="Plot the final spectrum?"}
+bool	redo = no		{prompt="Redo operations if previously done?"}
+bool	update = yes		{prompt="Update spectra if cal data changes?"}
+bool	batch = no		{prompt="Extract objects in batch?"}
+bool	listonly = no		{prompt="List steps but don't process?\n"}
+
+pset	params = ""		{prompt="Algorithm parameters"}
+
+begin
+	apscript.readnoise = readnoise
+	apscript.gain = gain
+	apscript.nfind = fibers
+	apscript.width = width
+	apscript.t_width = width
+	apscript.minsep = minsep
+	apscript.maxsep = maxsep
+	apscript.radius = minsep
+	apscript.clean = clean
+	proc.datamax = datamax
+
+	proc (objects, apref, flat, throughput, arcs1, arcs2, arcreplace,
+	    arctable, fibers, apidtable, crval, cdelt, objaps, skyaps,
+	    arcaps, objbeams, skybeams, arcbeams, scattered, fitflat, no,
+	    no, no, no, clean, dispcor, savearcs, skyalign, skysubtract,
+	    skyedit, saveskys, splot, redo, update, batch, listonly)
+
+	if (proc.dobatch) {
+	    print ("-- Do remaining spectra as a batch job --")
+	    print ("batch&batch") | cl
+	}
+end