Initial commit

8 files changed, 833 insertions, 0 deletions

diff --git a/noao/imred/quadred/doc/package.hlp b/noao/imred/quadred/doc/package.hlp
new file mode 100644
index 00000000..ffea9446
--- /dev/null
+++ b/noao/imred/quadred/doc/package.hlp
@@ -0,0 +1,142 @@
+.help package Sep93 quadred
+.ih
+NAME
+quadred -- CCD reductions of images in multi-amp readout format
+.ih
+SYNOPSIS
+This package is a varient of \fBccdred\fR that operates on a
+multi-amplifier data format in which the various amplifier readouts are
+recorded in sections of a regular two-dimensional image.  The CTIO Arcon
+dual or quad readout data is an example of this format.  See help on
+\fBquadformat\fR for details.  Most tasks are the same as in the
+\fBccdred\fR package.  The difference is the version of \fBccdproc\fR in
+this package also works on the multi-amp format.  An alternative to using
+this version of \fBccdproc\fR is \fBquadproc\fR and the alternate
+calibration combining task based on this task.
+.ih
+USAGE
+quadred
+.ih
+PARAMETERS
+The following are "package" parameters.  This means that they apply to
+many of the tasks in this package.
+
+.ls pixeltype = "real real"
+Output pixel datatype and calculation datatype.  When images are processed
+or created the output pixel datatype is determined by this parameter.
+The allowed types are "short" for short integer, and "real" for real
+floating point.  The calculation datatypes are also short and real with a
+default of real if none is specified.
+.le
+.ls verbose = no
+Print log information to the standard output?
+.le
+.ls logfile = "logfile"
+Text log file.  If no filename is specified then no log file is kept.
+.le
+.ls plotfile = ""
+Log metacode plot file for the overscan bias vector fits.  If
+no filename is specified then no metacode plot file is kept.
+.le
+.ls backup = ""
+Backup prefix for backup images.  If no prefix is specified then no backup
+images are kept when processing.  If specified then the backup image
+has the specified prefix.
+.le
+.ls instrument = ""
+CCD instrument translation file.  This is usually set with \fBsetinstrument\fR.
+.le
+.ls ssfile = "subsets"
+Subset translation file used to define the subset identifier.  See
+\fBsubsets\fR for more.
+.le
+.ls graphics = "stdgraph"
+Interactive graphics output device when fitting the overscan bias vector.
+.le
+.ls cursor = ""
+Graphics cursor input.  The default is the standard graphics cursor.
+.le
+.ls version = "Version 1.0 - August 22, 2001"
+Package version.
+.le
+.ih
+DESCRIPTION
+The \fBquadred\fR package contains all basic tasks necessary for the
+reduction of CCD data in single image format.  This includes both single
+amplifier readout data and multi-amplifier data stored as sections in a
+single two-dimensional image.  One example of this type of multi-amplifier
+data is the CTIO Arcon "dual" or "quad" readout format.  The format is
+described in the help topic \fBquadformat\fR.  This package is a
+combination of two earlier packages called \fBxccdred\fR and
+\fBared.quad\fR, each of which are variants of the original \fBccdred\fR
+package.
+
+The raw data contains overscan/prescan regions in the image.  For multi-amp
+data there are multiple overscan/prescan regions.  The first steps in
+processing the data is to use the overscan/prescan regions to determine
+the amplifier bias, subtract this bias, and trim the regions out of
+the data.  Once this is done the data are just simple images.  It is
+the special step of dealing with the overscan/prescan regions with
+the multi-amp format that is different from the standard \fBccdred\fR
+package.
+
+Two methods are provided for dealing with the special format.  One is a
+special version of \fBccdproc\fR which processes the sections directly.  If
+one uses this task then the reduction steps appear identical to using the
+\fBccdred\fR package.  The other method is to use the tasks \fBquadproc\fR,
+\fBqzerocombine\fR, \fBqdarkcombine\fR, and \fBqflatcombine\fR.  The latter
+calibration combining tasks are the same as the standard versions except
+they use \fBquadproc\fR instead of \fBccdproc\fR.  The task \fBquadproc\fR
+operates internally by splitting the multiple regions into temporary single
+amplifier images, processing them with \fBccdproc\fR, and then joining the
+pieces back together.
+
+The recommended method is to use \fBccdproc\fR.  However, the \fBquadproc\fR
+related tasks have a history of usage for CTIO data and so may also be
+used.
+
+The \fBquadred\fR package itself has several parameters which are
+common to many of the tasks in the package.  When images are processed or
+new image are created the output pixel datatype is that specified by the
+parameter \fBpixeltype\fR. Note that CCD processing replaces the original
+image by the processed image so the pixel type of the CCD images may change
+during processing.  It is unlikely that real images will be processed to
+short images but the reverse is quite likely.  Processing images from short
+to real pixel datatypes will generally increase the amount of disk space
+required (a factor of 2 on most computers).
+
+The tasks produce log output which may be printed on the standard
+output (the terminal unless redirected) and appended to a file.  The
+parameter \fIverbose\fR determines whether processing information
+is printed.  This may be desirable initially, but when using background
+jobs the verbose output should be turned off.  The user may look at
+the end of the log file (for example with \fBtail\fR) to determine
+the status of the processing.
+
+The package was designed to work with data from many different observatories
+and instruments.  In order to accomplish this an instrument translation
+file is used to define a mapping between the package parameters and
+the particular image header format.  The instrument translation file
+is specified to the package by the parameter \fIinstrument\fR.  This
+parameter is generally set by the task \fBsetinstrument\fR.  The other
+file used is a subset file.  This is generally created and maintained
+by the package and the user need not do anything.  For more sophisticated
+users see \fBinstruments\fR and \fBsubsets\fR.
+
+The package has very little graphics output.  The exception is the overscan
+bias subtraction.  The bias vector is logged in the metacode plot file if
+given.  The plot file may be examined with the tasks in the \fBplot\fR
+package such as \fBgkimosaic\fR.  When interactively fitting the overscan
+vector the graphics input and output devices must be specified.  The
+defaults should apply in most cases.
+
+Because processing replaces the input image by the processed image it may
+be desired to save the original image.  This may be done by specifying a
+backup prefix with the parameter \fIbackup\fR.  For example, if the prefix
+is "orig" and the image is "ccd001", the backup image will be
+"origccd001".  The prefix may be a directory but if so it must end with '/'
+or '$' (for logical directories) and the directory must already exist.
+.ih
+SEE ALSO
+quadformat, mscred
+.endhelp
diff --git a/noao/imred/quadred/doc/qhistogram.hlp b/noao/imred/quadred/doc/qhistogram.hlp
new file mode 100644
index 00000000..a34e412c
--- /dev/null
+++ b/noao/imred/quadred/doc/qhistogram.hlp
@@ -0,0 +1,37 @@
+.help qhistogram Aug01 noao.imred.quadred
+.ih
+NAME
+qhistogram -- Compute and print histogram for multi-amp data
+.ih
+USAGE
+qhistogram images
+.ih
+PARAMETERS
+.ls images
+List of image names in \fBquadformat\fR.
+.le
+.ls window = "datasec" (datasec|trimsec|biassec)
+Type of section to use for histogram.  The choices are "datasec" for the
+amplifier section which includes the bias if any is present, "trimsec" for
+the trim section, and "biassec" for the bias section.
+.le
+
+The remaining parameters come from the \fBimhistogram\fR task.
+.ih
+DESCRIPTION
+This script tasks uses the \fBquadsections\fR task to break the
+\fBquadformat\fR data into separate sections and runs the \fBimhistogram\fR
+task on the sections.  The graphics is collected onto a single page.
+.ih
+EXAMPLES
+
+1. To graph the histograms (default behavior).
+
+.nf
+    qu> qhist quad0072
+    [graph appears]
+.fi
+.ih
+SEE ALSO
+quadformat, quadsections, imhistogram
+.endhelp
diff --git a/noao/imred/quadred/doc/qstatistics.hlp b/noao/imred/quadred/doc/qstatistics.hlp
new file mode 100644
index 00000000..222ae778
--- /dev/null
+++ b/noao/imred/quadred/doc/qstatistics.hlp
@@ -0,0 +1,52 @@
+.help qstatistics Aug01 noao.imred.quadred
+.ih
+NAME
+qstatistics -- Compute and print statistics for multi-amp data
+.ih
+USAGE
+qstatistics images
+.ih
+PARAMETERS
+.ls images
+List of image names in \fBquadformat\fR.
+.le
+.ls window = "datasec" (datasec|trimsec|biassec)
+Type of section to output.  The choices are "datasec" for the amplifier
+section which includes the bias if any is present, "trimsec" for the trim
+section, and "biassec" for the bias section.
+.le
+
+The remaining parameters come from the \fBimstatistics\fR task.
+.ih
+DESCRIPTION
+This script tasks uses the \fBquadsections\fR task to break the
+\fBquadformat\fR data into separate sections and runs the \fBimstatistics\fR
+task on the sections.
+.ih
+EXAMPLES
+
+1. To compute the mean and stddev of the data section.
+
+.nf
+    qu> qstat quad0072 fields=image,mean,stddev
+    #               IMAGE      MEAN    STDDEV
+     quad0072[1:1034,1:1024]     5537.     2647.
+     quad0072[1163:2196,1:1024]     6210.     5439.
+     quad0072[1:1034,1025:2048]     5364.     2535.
+     quad0072[1163:2196,1025:2048]     5862.     1327.
+.fi
+
+2. To compute the mean and stdev of the bias section.
+
+.nf
+    qu> qstat quad0072 fields=image,mean,stddev window=biassec
+    #               IMAGE      MEAN    STDDEV
+     quad0072[1045:1098,1:1024]      713.     1.272
+     quad0072[1099:1152,1:1024]     516.2     1.425
+     quad0072[1045:1098,1025:2048]     554.3     1.347
+     quad0072[1099:1152,1025:2048]     530.3     1.377
+.fi
+.ih
+SEE ALSO
+quadformat, quadsections, imstatistics
+.endhelp
diff --git a/noao/imred/quadred/doc/quadformat.hlp b/noao/imred/quadred/doc/quadformat.hlp
new file mode 100644
index 00000000..eb5fbfbd
--- /dev/null
+++ b/noao/imred/quadred/doc/quadformat.hlp
@@ -0,0 +1,392 @@
+.help quadformat Aug01 imred.quadred
+.ih
+NAME
+quadformat - Description of the special multi-amplifier CCD format
+.ih
+DESCRIPTION
+CCDs may be readout from multiple amplifiers at the same time to increase
+the readout speed.  This produces multiple images of rectangular regions in
+the full CCD exposure.  The amplifier readout images may be recorded in
+various ways.  One way is as extensions in a multiextension FITS file.
+This type of format can be reduced using the MSCRED package.
+
+Another way is to paste the regions into a single two-dimensional image.
+This, along with specific keywords to describe the locations of the
+regions, constitutes the \fIquadformat\fR format described here and used by the
+QUADRED package.  The term "quad" originates from the possibility of using
+four amplifiers in quadrants but the format also includes any other
+number of amplifiers.
+
+It is important to realize that this is a special format only as long as
+the overscan or prescan data is included in the image data.  Once this
+information is used and removed as part of the processing the resulting
+image can be treated in the same way as a single amplifier CCD image.
+However, the image can still contain the format keywords allowing the
+regions from the different amplifiers to be identified and extracted as
+needed.
+
+The \fIquadformat\fR consists of a single 2D image for a single CCD
+exposure.  The image storage format may be any standard image type such
+as imh or fits.  Within the image are regions containing the CCD
+pixel data and regions containing overscan or prescan, which we will
+call bias regions, for each amplifier.  The \fIquadformat\fR requires
+the bias regions to be attached to the CCD regions such that a single
+rectangular region contains both.
+
+Generally the rectangular regions are of equal size in order to sequence
+the amplifiers simultaneously.  However, it is possible for the regions
+to be unequal in cases of subregion readouts with certain CCD controllers.
+The figure below illustrates a "dual" and "quad" readout with equal
+size regions.
+
+.nf
+        +-----+-+-+-----+   +-----+-+-+-----+   +----------+-+
+        |  D  !B|B!  D  |   |  D  !B|B!  D  |   |     D    !B|
+        |  3  !3|4!  4  |   |  1  !1|2!  2  |   |     2    !2|
+        |     ! | !     |   |     ! | !     |   |          ! |
+        +-----+-+-+-----+   |     ! | !     |   +----------+-+
+        |  D  !B|B!  D  |   |     ! | !     |   |     D    !B|
+        |  1  !1|2!  2  |   |     ! | !     |   |     1    !1|
+        |     ! | !     |   |     ! | !     |   |          ! |
+        +-----+-+-+-----+   +-----+-+-+-----+   +----------+-+
+.fi
+
+The areas labeled D are the data sections and those labeled B are the
+bias sections.  The data and biases are match by the amplifier labels
+which are 1-4 in these examples.  The combination of the data and
+bias sections are called the amplifier sections.
+
+The regions are identified in the header by various keywords.  There is
+a header translation facility which allows for alternative keyword names.
+Below we describe the default keyword names in the absence of a translation.
+The number of regions and the amplifier labels are described by the
+string keyword AMPLIST.  The value is a string of space separated
+amplifier labels.  For the above four amplifier example it would be
+
+.nf
+    AMPLIST = '1 2 3 4'
+.fi
+
+For CTIO data the labels are '11  12 21 22'.   Note that the labels
+are appended to rootnames so they should be relatively short. 
+
+The amplifier labels are appended to various root names.  The important
+ones define "section" keywords.  The values are image sections that
+describe regions in an raster such as the image or the CCD.  The format
+of a section follows the standard IRAF notation "[c1:c2,l1:l2]" where
+c1 and c2 are inclusive column endpoints and l1 and l2 are inclusive
+line endpoints.
+
+The various sections are defined below.  The labels again show the default
+untranslated keyword roots.
+
+.ls ASEC
+The section of the image containing the amplifier readout.  This is the
+combination of the data and bias regions as shown in the figures.
+.le
+.ls DSEC
+The section of the image containing the actual CCD data exclusive of
+bias data.  In the figures these are the D regions.
+.le
+.ls BSEC
+The section of the image containing the bias data.  In the figures these
+are the B regions.
+.le
+.ls TSEC
+The second of the image containing the useful CCD data.  This defines
+a "trimming" area and lies within the data section.  It may also be
+the same as the data region.  During trimming the final image will only
+include the regions in the trim sections.  Note that it generally does
+not make sense to trim between amplifier regions but does make sense to
+trim regions at the edges of the CCD.
+.le
+.ls CSEC
+The section of the CCD corresponding to the data section in the image.
+The CCD is considered an ideal raster (without bias regions) and a
+section corresponds to the pixels in the CCD.  The CCD section must be
+the same size as the data section.  It is the CCD sections that define
+how the amplifiers will be pieced together to form a single image
+after trimming the bias region.
+.le
+
+There may be other keyword root names for things such as gains which
+have the amplifier labels appended.  However, none of these are used
+by the current software.  Example image headers are given
+in the EXAMPLES section.
+
+There is a limitation in the current software that the regions be recorded
+without horizontal or vertical flips.  In other words, where amplifiers
+from opposite corners are used some of them must be flipped by the
+data acquisition system before recording then in this \fBquadformat\fR.
+
+.ih
+EXAMPLES
+
+1.  The following is an example of a full 2048x2048 CCD readout with
+four amplifiers at CTIO.
+
+.nf
+qu> imhad quad0020
+quad0020[2196,2048][ushort]: IC 1257 5290 180s
+No bad pixels, min=435., max=61973.
+Line storage mode, physdim [2304,2048], length of user area 3079 s.u.
+Created Thu 08:35:57 23-Aug-2001, Last modified Thu 08:35:57 23-Aug-2001
+Pixel file "HDR$pixels/quad0020.pix" [ok]
+'KPNO-IRAF'           /
+'06-07-99'            /
+IRAF-MAX=           6.197300E4  /  DATA MAX
+IRAF-MIN=           4.350000E2  /  DATA MIN
+IRAF-BPX=                   16  /  DATA BITS/PIXEL
+IRAFTYPE= 'USHORT  '            /  PIXEL TYPE
+OPICNUM =                  123 / Original picture number
+HDR_REV = '2.000  13Feb96     (add mode and group to hdrs)' /
+IMAGETYP= 'OBJECT  '           / Type of picture (object, dark, etc.)
+DETECTOR= 'Site2K_6'           / Detector (CCD type, photon counter, etc.)
+PREFLASH=             0.000000 / Preflash time in secs
+CCDSUM  = '1 1     '           / On chip summation (X,Y)
+DATE-OBS= '07/07/99'           / Date (dd/mm/yy) of observation
+UTSHUT  = '01:14:40.0'         / UT of shutter open
+UT      = ' 1:14:41.50'         /  UT of TCS coords
+OBSERVAT= 'CTIO    '           / Origin of data
+TELESCOP= 'CTIO 1.5 meter telescope' / Specific system
+NAMPSYX = '2 2     '           / Num amps in y & x (eg. '2 2'=quad)
+AMPLIST = '11 21 12 22'        / Readout order in y,x
+ASEC11  = '[1:1098,1:1024]'    / Section read with Amp11
+CSEC11  = '[1:1034,1:1024]'    / Section in full CCD for DSEC11
+DSEC11  = '[1:1034,1:1024]'    / Image area in raw frame for Amp11
+TSEC11  = '[11:1034,1:1024]'   / Trim section definition for Amp11
+BSEC11  = '[1045:1098,1:1024]' / Bias section definition for Amp11
+BSEC12  = '[1099:1152,1:1024]' / Bias section definition for Amp12
+ASEC12  = '[1099:2196,1:1024]' / Section read with Amp12
+CSEC12  = '[1035:2068,1:1024]' / Section in full CCD for DSEC12
+DSEC12  = '[1163:2196,1:1024]' / Image area in raw frame for Amp12
+TSEC12  = '[1163:2186,1:1024]' / Trim section definition for Amp12
+ASEC21  = '[1:1098,1025:2048]' / Section read with Amp21
+CSEC21  = '[1:1034,1025:2048]' / Section in full CCD for DSEC21
+DSEC21  = '[1:1034,1025:2048]' / Image area in raw frame for Amp21
+TSEC21  = '[11:1034,1025:2048]' / Trim section definition for Amp21
+BSEC21  = '[1045:1098,1025:2048]' / Bias section definition for Amp21
+BSEC22  = '[1099:1152,1025:2048]' / Bias section definition for Amp22
+ASEC22  = '[1099:2196,1025:2048]' / Section read with Amp22
+CSEC22  = '[1035:2068,1025:2048]' / Section in full CCD for DSEC22
+DSEC22  = '[1163:2196,1025:2048]' / Image area in raw frame for Amp22
+TSEC22  = '[1163:2186,1025:2048]' / Trim section definition for Amp22
+WAVEFILE= 'Obs Tue Jul  6 20:11:59 1999' /
+NOTE    = 'WARNING: Lower amps reaching full well before ADCs saturate' /
+WAVEMODE= 'MPP OverlapXmit EarlyReset' / Waveform mode switches on
+GTRON22 =                4.100 / (e-) predicted read noise, upper right
+GTRON21 =                3.900 / (e-) predicted read noise, upper left
+GTRON12 =                4.200 / (e-) predicted read noise, lower right
+GTRON11 =                4.200 / (e-) predicted read noise, lower left
+GTGAIN22=                2.800 / (e-/ADU), predicted gain, upper right
+GTGAIN21=                3.100 / (e-/ADU) predicted gain, upper left
+GTGAIN12=                2.900 / (e-/ADU) predicted gain, lower right
+GTGAIN11=                3.200 / (e-/ADU) predicted gain, lower left
+GTINDEX =                    2 / Gain selection (index into Gain Table)
+PIXELT  =                29520 / (ns) unbinned pixel read time
+DCS_TIME=                 7000 / (ns) Double Correlated Sample time
+RA      = '17:27:10.82'         /  right ascension (telescope)
+DEC     = '-7:06:35.40'         /  declination (telescope)
+EPOCH   =               2000.0 / epoch of RA & DEC
+ZD      =                 35.9 / zenith distance (degrees)
+HA      = '-01:57:23.7'        / hour angle (H:M:S)
+ST      = '15:29:46.00'         /  sidereal time
+AIRMASS =                1.234 / airmass
+EXPTIME =              180.000 / Exposure time in secs
+DARKTIME=              181.309 / Total elapsed time in secs
+OBSERVER= 'Jacoby'             / Observers
+PROPID  = '92'                 / Proposal Id
+COMMENT Globular PNe
+TELID   = 'ct60'               / CTIO 1.5-m Telescope
+ARCONVER= '17Oct97ver7_22'     / Arcon software version
+COMMENT INSTRUMENT PARAMETERS
+INSTRUME= 'cfccd'              / cassegrain direct imager
+FILTER1 = 'dia'                / Filter in wheel one
+FNAME1  = 'diaphragm'          / Full name of filter in  wheel1
+FILTER2 = 'ocon'               / Filter in wheel two
+FNAME2  = 'O cont'             / Full name of filter in  wheel2
+FILTERS = 'dia ocon'           / Filter positions
+TELFOCUS=                57550 / Telescope focus
+XPIXSIZE=                0.432 / Pixel size in X (arcsec/pix)
+YPIXSIZE=                0.432 / Pixel size in Y (arcsec/pix)
+RECID   = 'ct60.990707.011817' / NOAO Archive record ID
+.fi
+
+2.  The following is a more complex readout of a region where the
+full 2Kx2K CCD is not readout and where even the regions are not the
+same size.
+
+.nf
+qu> imhead quad0013
+quad0013[1686,1538][ushort]: R sky flat 7s
+No bad pixels, min=393., max=65535.
+Line storage mode, physdim [1792,1538], length of user area 3079 s.u.
+Created Thu 08:34:00 23-Aug-2001, Last modified Thu 08:34:00 23-Aug-2001
+Pixel file "HDR$pixels/quad0013.pix" [ok]
+'KPNO-IRAF'           /
+'06-07-99'            /
+IRAF-MAX=           6.553500E4  /  DATA MAX
+IRAF-MIN=           3.930000E2  /  DATA MIN
+IRAF-BPX=                   16  /  DATA BITS/PIXEL
+IRAFTYPE= 'USHORT  '            /  PIXEL TYPE
+OPICNUM =                   15 / Original picture number
+HDR_REV = '2.000  13Feb96     (add mode and group to hdrs)' /
+IMAGETYP= 'SKY FLAT'           / Type of picture (object, dark, etc.)
+DETECTOR= 'Site2K_6'           / Detector (CCD type, photon counter, etc.)
+PREFLASH=             0.000000 / Preflash time in secs
+CCDSUM  = '1 1     '           / On chip summation (X,Y)
+DATE-OBS= '06/07/99'           / Date (dd/mm/yy) of observation
+UTSHUT  = '22:25:22.0'         / UT of shutter open
+UT      = '22:25:34.00'         /  UT of TCS coords
+OBSERVAT= 'CTIO    '           / Origin of data
+TELESCOP= 'CTIO 1.5 meter telescope' / Specific system
+NAMPSYX = '2 2     '           / Num amps in y & x (eg. '2 2'=quad)
+AMPLIST = '11 21 12 22'        / Readout order in y,x
+ASEC11  = '[1:843,1:769]'      / Section read with Amp11
+CSEC11  = '[256:1034,256:1024]' / Section in full CCD for DSEC11
+DSEC11  = '[1:779,1:769]'      / Image area in raw frame for Amp11
+TSEC11  = '[11:779,1:769]'     / Trim section definition for Amp11
+BSEC11  = '[790:843,1:769]'    / Bias section definition for Amp11
+BSEC12  = '[844:897,1:769]'    / Bias section definition for Amp12
+ASEC12  = '[844:1686,1:769]'   / Section read with Amp12
+CSEC12  = '[1035:1813,256:1024]' / Section in full CCD for DSEC12
+DSEC12  = '[908:1686,1:769]'   / Image area in raw frame for Amp12
+TSEC12  = '[908:1418,1:769]'   / Trim section definition for Amp12
+ASEC21  = '[1:843,770:1538]'   / Section read with Amp21
+CSEC21  = '[256:1034,1025:1793]' / Section in full CCD for DSEC21
+DSEC21  = '[1:779,770:1538]'   / Image area in raw frame for Amp21
+TSEC21  = '[11:779,770:1280]'  / Trim section definition for Amp21
+BSEC21  = '[790:843,770:1538]' / Bias section definition for Amp21
+BSEC22  = '[844:897,770:1538]' / Bias section definition for Amp22
+ASEC22  = '[844:1686,770:1538]' / Section read with Amp22
+CSEC22  = '[1035:1813,1025:1793]' / Section in full CCD for DSEC22
+DSEC22  = '[908:1686,770:1538]' / Image area in raw frame for Amp22
+TSEC22  = '[908:1418,770:1280]' / Trim section definition for Amp22
+WAVEFILE= 'Obs Tue Jul  6 18:07:56 1999' /
+NOTE    = 'WARNING: Lower amps reaching full well before ADCs saturate' /
+WAVEMODE= 'MPP OverlapXmit EarlyReset' / Waveform mode switches on
+GTRON22 =                4.100 / (e-) predicted read noise, upper right
+GTRON21 =                3.900 / (e-) predicted read noise, upper left
+GTRON12 =                4.200 / (e-) predicted read noise, lower right
+GTRON11 =                4.200 / (e-) predicted read noise, lower left
+GTGAIN22=                2.800 / (e-/ADU), predicted gain, upper right
+GTGAIN21=                3.100 / (e-/ADU) predicted gain, upper left
+GTGAIN12=                2.900 / (e-/ADU) predicted gain, lower right
+GTGAIN11=                3.200 / (e-/ADU) predicted gain, lower left
+GTINDEX =                    2 / Gain selection (index into Gain Table)
+PIXELT  =                29520 / (ns) unbinned pixel read time
+DCS_TIME=                 7000 / (ns) Double Correlated Sample time
+RA      = '14:53:52.67'         /  right ascension (telescope)
+DEC     = '-19:20:10.70'        /  declination (telescope)
+EPOCH   =               2000.0 / epoch of RA & DEC
+ZD      =                 32.1 / zenith distance (degrees)
+HA      = '-02:13:40.3'        / hour angle (H:M:S)
+ST      = '12:40:10.80'         /  sidereal time
+AIRMASS =                1.180 / airmass
+EXPTIME =                7.000 / Exposure time in secs
+DARKTIME=                8.239 / Total elapsed time in secs
+OBSERVER= 'Jacoby'             / Observers
+PROPID  = '92'                 / Proposal Id
+COMMENT
+TELID   = 'ct60'               / CTIO 1.5-m Telescope
+ARCONVER= '17Oct97ver7_22'     / Arcon software version
+COMMENT INSTRUMENT PARAMETERS
+INSTRUME= 'cfccd'              / cassegrain direct imager
+FILTER1 = 'dia'                / Filter in wheel one
+FNAME1  = 'diaphragm'          / Full name of filter in  wheel1
+FILTER2 = 'r'                  / Filter in wheel two
+FNAME2  = 'R'                  / Full name of filter in  wheel2
+FILTERS = 'dia r'              / Filter positions
+TELFOCUS=                    0 / Telescope focus
+XPIXSIZE=                0.432 / Pixel size in X (arcsec/pix)
+YPIXSIZE=                0.432 / Pixel size in Y (arcsec/pix)
+RECID   = 'ct60.990706.222551' / NOAO Archive record ID
+.fi
+
+3.  The following is for the raw image of example 2 after it has been
+processed by CCDPROC.  Note that the various bias, trim, and CCD sections are
+removed.  The AMPLIST and ASEC keywords remain and may be used to split
+or evaluate the individual amplifier regions with tasks such as QUADSECTIONS,
+QUADSPLIT, and QSTATISTICS.
+
+.nf
+qu> imhead quad0013
+quad0013[1280,1280][real]: R sky flat 7s
+No bad pixels, min=unknown, max=unknown
+Line storage mode, physdim [1280,1280], length of user area 2795 s.u.
+Created Fri 13:29:40 24-Aug-2001, Last modified Fri 13:29:40 24-Aug-2001
+Pixel file "HDR$pixels/quad0013.pix" [ok]
+'KPNO-IRAF'           /
+'06-07-99'            /
+New copy of quad0013
+IRAF-MAX=           6.553500E4  /  DATA MAX
+IRAF-MIN=           3.930000E2  /  DATA MIN
+IRAF-BPX=                   16  /  DATA BITS/PIXEL
+IRAFTYPE= 'USHORT  '            /  PIXEL TYPE
+OPICNUM =                   15 / Original picture number
+HDR_REV = '2.000  13Feb96     (add mode and group to hdrs)' /
+IMAGETYP= 'SKY FLAT'           / Type of picture (object, dark, etc.)
+DETECTOR= 'Site2K_6'           / Detector (CCD type, photon counter, etc.)
+PREFLASH=             0.000000 / Preflash time in secs
+CCDSUM  = '1 1     '           / On chip summation (X,Y)
+DATE-OBS= '06/07/99'           / Date (dd/mm/yy) of observation
+UTSHUT  = '22:25:22.0'         / UT of shutter open
+UT      = '22:25:34.00'         /  UT of TCS coords
+OBSERVAT= 'CTIO    '           / Origin of data
+TELESCOP= 'CTIO 1.5 meter telescope' / Specific system
+NAMPSYX = '2 2     '           / Num amps in y & x (eg. '2 2'=quad)
+AMPLIST = '11 21 12 22'        / Readout order in y,x
+ASEC11  = '[1:769,1:769]'      / Section read with Amp11
+ASEC12  = '[770:1280,1:769]'   / Section read with Amp12
+ASEC21  = '[1:769,770:1280]'   / Section read with Amp21
+ASEC22  = '[770:1280,770:1280]' / Section read with Amp22
+WAVEFILE= 'Obs Tue Jul  6 18:07:56 1999' /
+NOTE    = 'WARNING: Lower amps reaching full well before ADCs saturate' /
+WAVEMODE= 'MPP OverlapXmit EarlyReset' / Waveform mode switches on
+GTRON22 =                4.100 / (e-) predicted read noise, upper right
+GTRON21 =                3.900 / (e-) predicted read noise, upper left
+GTRON12 =                4.200 / (e-) predicted read noise, lower right
+GTRON11 =                4.200 / (e-) predicted read noise, lower left
+GTGAIN22=                2.800 / (e-/ADU), predicted gain, upper right
+GTGAIN21=                3.100 / (e-/ADU) predicted gain, upper left
+GTGAIN12=                2.900 / (e-/ADU) predicted gain, lower right
+GTGAIN11=                3.200 / (e-/ADU) predicted gain, lower left
+GTINDEX =                    2 / Gain selection (index into Gain Table)
+PIXELT  =                29520 / (ns) unbinned pixel read time
+DCS_TIME=                 7000 / (ns) Double Correlated Sample time
+RA      = '14:53:52.67'         /  right ascension (telescope)
+DEC     = '-19:20:10.70'        /  declination (telescope)
+EPOCH   =               2000.0 / epoch of RA & DEC
+ZD      =                 32.1 / zenith distance (degrees)
+HA      = '-02:13:40.3'        / hour angle (H:M:S)
+ST      = '12:40:10.80'         /  sidereal time
+AIRMASS =                1.180 / airmass
+EXPTIME =                7.000 / Exposure time in secs
+DARKTIME=                8.239 / Total elapsed time in secs
+OBSERVER= 'Jacoby'             / Observers
+PROPID  = '92'                 / Proposal Id
+COMMENT
+TELID   = 'ct60'               / CTIO 1.5-m Telescope
+ARCONVER= '17Oct97ver7_22'     / Arcon software version
+COMMENT INSTRUMENT PARAMETERS
+INSTRUME= 'cfccd'              / cassegrain direct imager
+FILTER1 = 'dia'                / Filter in wheel one
+FNAME1  = 'diaphragm'          / Full name of filter in  wheel1
+FILTER2 = 'r'                  / Filter in wheel two
+FNAME2  = 'R'                  / Full name of filter in  wheel2
+FILTERS = 'dia r'              / Filter positions
+TELFOCUS=                    0 / Telescope focus
+XPIXSIZE=                0.432 / Pixel size in X (arcsec/pix)
+YPIXSIZE=                0.432 / Pixel size in Y (arcsec/pix)
+RECID   = 'ct60.990706.222551' / NOAO Archive record ID
+TRIM    = 'Aug 24 13:29 Trim multiple overscan sections'
+OVERSCAN= 'Aug 24 13:29 Overscan is [790:843,1:769] with mean=714.3438'
+OVRSCN2 = 'Aug 24 13:29 Overscan is [790:843,770:1538] with mean=554.01'
+OVRSCN3 = 'Aug 24 13:29 Overscan is [844:897,1:769] with mean=519.7755'
+OVRSCN4 = 'Aug 24 13:29 Overscan is [844:897,770:1538] with mean=531.69'
+CCDSEC  = '[266:1545,256:1535]'
+CCDMEAN =             9727.605
+CCDMEANT=            683126983
+CCDPROC = 'Aug 24 13:29 CCD processing done'
+.fi
+.endhelp
diff --git a/noao/imred/quadred/doc/quadjoin.hlp b/noao/imred/quadred/doc/quadjoin.hlp
new file mode 100644
index 00000000..2a3a075e
--- /dev/null
+++ b/noao/imred/quadred/doc/quadjoin.hlp
@@ -0,0 +1,43 @@
+.help quadjoin Aug01 noao.imred.quadred
+.ih
+NAME
+quadjoin -- Split quadformat data into single amplifier images
+.ih
+USAGE
+quadjoin input
+.ih
+PARAMETERS
+.ls input
+Root name of images to be joined.  Extensions based on the AMPLIST
+keyword are applied to the root name.  This task does not
+allow a list of input root names.
+.le
+.ls output = ""
+Output image name.  If one is not given then the input root name is used.
+.le
+.ls delete = no
+Delete subimages on completion?
+.le
+.ih
+DESCRIPTION
+Images in split "quadformat" (see help topic \fBquadformat\fR and
+\fBquadsplit\fR) are rejoined into "quadformat".  The input images
+have a common root name and then an extension given by the amplifier
+labels in the AMPLIST keyword are added.  The output name may be specified
+or the input root name may be used.
+.ih
+EXAMPLES
+1. To join a split set of images:
+
+.nf
+    qu> dir quad0072*
+    quad0072.11.imh     quad0072.21.imh
+    quad0072.12.imh     quad0072.22.imh     
+    qu> quadjoin quad0072 delete+
+    qu> dir quad0072*
+    quad0072.imh
+.fi
+.ih
+SEE ALSO
+quadformat, quadsplit
+.endhelp
diff --git a/noao/imred/quadred/doc/quadscale.hlp b/noao/imred/quadred/doc/quadscale.hlp
new file mode 100644
index 00000000..7493bcf6
--- /dev/null
+++ b/noao/imred/quadred/doc/quadscale.hlp
@@ -0,0 +1,37 @@
+.help quadscale Aug01 noao.imred.quadred
+.ih
+NAME
+quadscale -- Scale amplifier sections by separate gains
+.ih
+USAGE
+quadscale input output
+.ih
+PARAMETERS
+.ls input
+Input image in \fBquadformat\fR to be scaled.
+.le
+.ls output
+Output scaled image in \fBquadformat\fR.
+.le
+.ls gain11 = 1., gain12 = 1., gain21 = 1., gain22 = 1.
+Gain factors for each quadrant.
+.le
+.ls operation = "multiply" (multiply|divide)
+The operation to apply with the gains.
+.le
+.ih
+DESCRIPTION
+This task multiplies or divides by gain factors for each amplifier in
+\fBquadformat\fR.
+.ih
+EXAMPLES
+
+1. To multiply by different gain factors.
+
+.nf
+    qu> quadscale quad0072 test gain11=1.2 gain12=1.3 gain21=1.4
+.fi
+.ih
+SEE ALSO
+quadformat
+.endhelp
diff --git a/noao/imred/quadred/doc/quadsections.hlp b/noao/imred/quadred/doc/quadsections.hlp
new file mode 100644
index 00000000..2735e3d5
--- /dev/null
+++ b/noao/imred/quadred/doc/quadsections.hlp
@@ -0,0 +1,81 @@
+.help quadsections Aug01 noao.imred.quadred
+.ih
+NAME
+quadsections -- Create image sections
+.ih
+USAGE
+quadsplit images
+.ih
+PARAMETERS
+.ls images
+List of image names for images in \fBquadformat\fR.
+.le
+.ls window = "datasec" (datasec|trimsec|biassec)
+Type of section to output.  The choices are "datasec" for the amplifier
+section which includes the bias if any is present, "trimsec" for the trim
+section, and "biassec" for the bias section.
+.le
+.ls section = ""
+Section to be overlapped.  The output sections will be the parts of the
+amplifier windows which are included within this section.
+.le
+.ls template = ""
+Template for producing the output.  The template replaces occurs of
+$I with the image name, $S with the section, and $A with the amplifier
+label.  If none is specified then the default template "$I$S\\n" is
+used which produces the image name with section separated by new-lines.
+The special characters "\n" is the new-line and the extra "\" is
+required to pass the new-line through to the formatting routine.
+.le
+.ih
+DESCRIPTION
+Images in "quadformat" (see help topic \fBquadformat\fR) are broken down
+in sections and written to the standard output in a specified format.
+.ih
+EXAMPLES
+1. To print the default data sections.
+
+.nf
+    qu> quadsec quad0072
+    quad0072[1:1034,1:1024]
+    quad0072[1163:2196,1:1024]
+    quad0072[1:1034,1025:2048]
+    quad0072[1163:2196,1025:2048]
+.fi
+
+3. To apply an overlap section.
+
+.nf
+    qu> quadsec quad0072 section=[1000:2000,1000:2000]
+    quad0072[1000:1034,1000:1024]
+    quad0072[1163:2000,1000:1024]
+    quad0072[1000:1034,1025:2000]
+    quad0072[1163:2000,1025:2000]
+.fi
+
+2. To print the trim sections.
+
+.nf
+    qu> quadsec quad0072 window=trimsec
+    quad0072[11:1034,1:1024]
+    quad0072[1163:2186,1:1024]
+    quad0072[11:1034,1025:2048]
+    quad0072[1163:2186,1025:2048]
+.fi
+
+
+4.  To make a custom output.
+
+.nf
+    qu> quadsec quad0072 template="image=$I, section=$S, amplifier=$A\\n"
+    image=quad0072, section=[1:1034,1:1024], amplifier=11
+    image=quad0072, section=[1163:2196,1:1024], amplifier=12
+    image=quad0072, section=[1:1034,1025:2048], amplifier=21
+    image=quad0072, section=[1163:2196,1025:2048], amplifier=22
+    qu> quadsec quad0072 template="$I.$A,"
+    quad0072.11,quad0072.12,quad0072.21,quad0072.22,
+.fi
+.ih
+SEE ALSO
+quadformat
+.endhelp
diff --git a/noao/imred/quadred/doc/quadsplit.hlp b/noao/imred/quadred/doc/quadsplit.hlp
new file mode 100644
index 00000000..4a0adf66
--- /dev/null
+++ b/noao/imred/quadred/doc/quadsplit.hlp
@@ -0,0 +1,49 @@
+.help quadsplit Aug01 noao.imred.quadred
+.ih
+NAME
+quadsplit -- Split quadformat data into single amplifier images
+.ih
+USAGE
+quadsplit input
+.ih
+PARAMETERS
+.ls input
+Image name of \fIquadformat\fR image to be split.  This task does not
+allow a list of input names.
+.le
+.ls output = ""
+Output root name to which the AMPLIST amplifier identifiers will be
+appended to form the split images.  If no output name is given then
+the input name is used as the root name.
+.le
+.ls clobber = yes
+Clobber any existing images?
+.le
+.ih
+DESCRIPTION
+Images in "quadformat" (see help topic \fBquadformat\fR) are separated
+into images containing data from only one amplifier.  The output images
+have a common root name and then an extension given by the amplifier
+labels in the AMPLIST keyword.  The output root name may be specified
+or default to the input name.
+
+In addition to producing the individual images keywords, are added that
+are understood by the standard \fBccdproc\fR task for single amplifier
+CCD reductions.
+
+The task \fBquadjoin\fR may be used to rejoin images that were split
+by this task.
+.ih
+EXAMPLES
+1. To spit an image:
+
+.nf
+    qu> quadsplit quad0072
+    qu> dir quad0072*
+    quad0072.11.imh     quad0072.21.imh     quad0072.imh        
+    quad0072.12.imh     quad0072.22.imh     
+.fi
+.ih
+SEE ALSO
+quadformat, quadjoin
+.endhelp