From fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 Mon Sep 17 00:00:00 2001
From: Joseph Hunkeler <jhunkeler@gmail.com>
Date: Wed, 8 Jul 2015 20:46:52 -0400
Subject: Initial commit

---
 noao/imred/kpnocoude/do3fiber.cl | 60 ++++++++++++++++++++++++++++++++++++++++
 1 file changed, 60 insertions(+)
 create mode 100644 noao/imred/kpnocoude/do3fiber.cl

(limited to 'noao/imred/kpnocoude/do3fiber.cl')

diff --git a/noao/imred/kpnocoude/do3fiber.cl b/noao/imred/kpnocoude/do3fiber.cl
new file mode 100644
index 00000000..649684cb
--- /dev/null
+++ b/noao/imred/kpnocoude/do3fiber.cl
@@ -0,0 +1,60 @@
+# DO3FIBERS -- Process Coude fiber 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 do3fiber (objects)
+
+string	objects = ""		{prompt="List of object spectra"}
+
+file	apref = ""		{prompt="Aperture reference spectrum"}
+file	flat = ""		{prompt="Flat field spectrum"}
+string	arcs = ""		{prompt="List of arc spectra"}
+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 = 3		{prompt="Number of fibers"}
+real	width = 6.		{prompt="Width of profiles (pixels)"}
+string	crval = "INDEF"		{prompt="Approximate central wavelength"}
+string	cdelt = "INDEF"		{prompt="Approximate dispersion"}
+string	objaps = "2"		{prompt="Object apertures"}
+string	arcaps = "1,3"		{prompt="Arc apertures\n"}
+
+bool	scattered = no		{prompt="Subtract scattered light?"}
+bool	fitflat = yes		{prompt="Fit and ratio flat field spectrum?"}
+bool	recenter = yes		{prompt="Recenter object apertures?"}
+bool	edit = no		{prompt="Edit/review object apertures?"}
+bool	clean = no		{prompt="Detect and replace bad pixels?"}
+bool	dispcor = yes		{prompt="Dispersion correct spectra?"}
+bool	splot = yes		{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.radius = width
+	apscript.clean = clean
+	apscript.order = "increasing"
+	proc.datamax = datamax
+
+	proc (objects, apref, flat, "", arcs, "", "",
+	    arctable, fibers, "", crval, cdelt, objaps, "", arcaps, "",
+	    "", "", scattered, fitflat, recenter, edit, no, no, clean,
+	    dispcor, no, no, no, no, no, splot, redo, update, batch, listonly)
+
+	if (proc.dobatch) {
+	    print ("-- Do remaining spectra as a batch job --")
+	    print ("batch&batch") | cl
+	}
+end
-- 
cgit