NOIRLAB: Re-add the vtel package

noao/imred/imred.cl

@@ -18,6 +18,7 @@ set kpnocoude = "imred$kpnocoude/"
 set kpnoslit = "imred$kpnoslit/"
 set quadred = "imred$quadred/"
 set specred = "imred$specred/"
+set vtel = "imred$vtel/"
 
 set apextract = "twodspec$apextract/"
 set doecslit = "imred$src/doecslit/"
@@ -49,5 +50,6 @@ task kpnocoude.pkg = "kpnocoude$kpnocoude.cl"
 task kpnoslit.pkg = "kpnoslit$kpnoslit.cl"
 task quadred.pkg = "quadred$quadred.cl"
 task specred.pkg = "specred$specred.cl"
+task vtel.pkg = "vtel$vtel.cl"
 
 clbye

noao/imred/imred.hd

@@ -18,6 +18,7 @@ $kpnocoude = "noao$imred/kpnocoude/"
 $kpnoslit = "noao$imred/kpnoslit/"
 $quadred = "noao$imred/quadred/"
 $specred = "noao$imred/specred/"
+$vtel = "noao$imred/vtel/"
 
 demos hlp=doc$demos.hlp
 #tutor hlp=doc$tutor.hlp, src=tutor.cl
@@ -110,3 +111,9 @@ specred men=specred$specred.men,
  sys=specred$specred.hlp,
  pkg=specred$specred.hd,
  src=specred$specred.cl
+
+vtel men=vtel$vtel.men,
+ hlp=..,
+ sys=vtel$vtel.hlp,
+ pkg=vtel$vtel.hd,
+ src=vtel$vtel.cl

noao/imred/imred.men

@@ -14,3 +14,4 @@
  kpnoslit - KPNO low/moderate dispersion slits (Goldcam, RCspec, Whitecam)
  quadred - CCD reductions for QUAD amplifier data
  specred - Generic slit and fiber spectral reduction package
+ vtel - NSO Solar vacuum telescope image reductions

noao/imred/mkpkg

@@ -15,4 +15,6 @@ update:
  $call update@irred
  $echo "------------------ IMRED.QUADRED -------------------"
  $call update@quadred
+ $echo "------------------ IMRED.VTEL ----------------------"
+ $call update@vtel
  ;

noao/imred/vtel/README

@@ -0,0 +1,81 @@
+This is the home directory for the Kitt Peak vacuum telescope
+reduction programs.
+
+README this file
+Revisions revisions file
+asciilook lookup table for ascii values into the pixelfont
+d1900.x calculate number of days since turn of century
+decodeheader.x decode/print vacuum telescope tape header
+destreak.par 
+destreak.x destreak 10830 full disk helium grams
+destreak5.cl script for processing 10830 tape containing 5 grams
+destreak5.par
+dicoplot.h header file containing defines for DICOPLOT
+dicoplot.par
+dicoplot.x program to make Carrington rotation mape on the Dicomed
+doc documentation directory
+ephem.x program to calculate solar ephemeris data
+fitslogr.cl script, make a log file of a fits tape (daily grams)
+fitslogr.par
+getsqib.par
+getsqib.x get the squibby brightness image from a full disk gram
+gryscl.dico greyscale lookup table for use with DICOPLOT
+imfglexr.x Get Line with EXtension Real for use with IMFilt
+imfilt.x convolve an image with gaussian kernel, used in destreak
+imratio.x find the ratio between two images, used in merge
+imtext.x subroutine to load text into an image by overwriting pixels
+lstsq.x least squares fitting subroutine
+makehelium.cl script to process a helium 10830 tape into daily grams (180)
+makehelium.par
+makeimages.cl script to process a magnetogram tape into daily grams (180)
+makeimages.par
+merge.par
+merge.x program to merge daily grams into Carrington rotation maps
+mkpkg make the package
+mrotlogr.cl script, make a log file of a fits tape (Carrington rotations)
+mrotlogr.par
+mscan.par
+mscan.x read vacuum telescope area scan tapes
+numeric.h header file for numeric subroutine
+numeric.x subroutine to calculate derivitives of latitude and longitude
+ with respect to x and y respectively (used in rmap)
+pimtext.par
+pimtext.x program to put text into images by overwriting pixels
+pixbit.x subroutine that looks up text pixel format in pixelfont
+pixelfont pixel font for use with pimtext (no lower case, no decenders)
+putsqib.par
+putsqib.x program to put the squibby brightness back in a full disk gram
+quickfit.par
+quickfit.x fit an ellipse to the limb of the sun
+readheader.x read a vacuum telescope header
+readss1.x subroutine to read a type 1 area scan
+readss2.x subroutine to read a type 2 area scan
+readss3.x subroutine to read a type 3 area scan
+readss4.x subroutine to read a type 4 area scan
+readsubswath.x subroutine to read a sub-swath
+readvt.par
+readvt.x read full disk grams from tape
+rmap.par
+rmap.x map full disk grams into daily grams (180x180)
+syndico.x Make dicomed print of daily grams 18 cm across.
+tcopy.par
+tcopy.x tape to tape copy program
+trim.par
+trim.x trim a full disk gram using squibby brightness info
+unwrap.par
+unwrap.x program to remove binary wrap-around from images
+vt.h
+vtblink.cl script to blink images on the IIS to check registration
+vtblink.par
+vtel.cl load the vacuum telescope package
+vtel.hd info about locations of various files
+vtel.men menu for package
+vtel.par
+vtexamine.par
+vtexamine.x program to examine a vacuum telescope tape (tell about record
+ lengths, header info, number of files, etc.)
+writetape.cl script to write five full disk grams to tape
+writetape.par
+writevt.par
+writevt.x program to write a full disk gram to tape in mountain format
+x_vtel.x package parent program

noao/imred/vtel/Revisions

@@ -0,0 +1,209 @@
+This is the vacuum telescope package revisions file.
+
+mkpkg
+ Added some missing file dependencies and removed unnecessary ones from
+ the mkpkg file. (9/30/99, Davis)
+
+doc/dicoplot.hlp
+doc/readvt.hlp
+doc/unwrap.hlp
+doc/pimtext.hlp
+ Fixed minor formating problems. (4/22/99, Valdes)
+
+=======
+V2.11.1
+=======
+
+May 16, 1989 by Dyer Lytle mods to 'readvt', 'syndico', and 'mscan'
+
+Fixed readvt to work with tape drives over the network [(if (mtfile(...].
+Modified syndico to take advantage of the disk-center info in the image
+header.
+
+Modified mscan to be much faster by taking out the geometrical correction.
+Also simplified it by removing the date/time pimtext call. Also made it
+create only the images needed. Also made it have a short file name option.
+Also made it work on tape drives over the net.
+
+
+June 5, 1988 by Dyer Lytle modification to PUTSQIB
+
+PUTSQIB had code in it for triming the limb as well as merging the two
+images. I simplified the program to just make the merge. The task TRIM
+can be used to trim the limb, and do a better job of it at that.
+
+September 29, 1987 by Dyer Lytle add SYNDICO to package
+
+Added this new program for makeing dicomed prints of daily
+grams 18 cm across.
+
+July 17, 1987 by Dyer Lytle fix bug in numeric.x
+
+There was a bug in the way an error flag was being set that made
+the program fail with a 'divide by zero' error on some data sets.
+
+June 8, 1987 by Dyer Lytle Overhaul of the package
+
+Major modifications were made to the code to make it conform to IRAF
+standards. Dynamic memory allocation replaced fixed memory allocation
+in many places. Readvt was modified to accept templates for input
+and output file names. New structures were provided for the vacuum
+telescope header, the tapeio buffer, and to reduce the argument count
+for the subroutine 'numeric'. Vtfix was dropped from the package
+since 'readvt' was modified to check for long records by doing its
+own buffering. Unwrap was updated to a new, more general and powerful
+version. A major bug was found and fixed in 'rmap' which was causing
+the total mapped pixel count to be off by about 20%.
+
+June 10, 1986 by Dyer Lytle Modification of PIMTEXT
+
+Pimtext was modified to allow the user to magnify the text in x and/or y
+and to get the date and/or time from a reference image if desired.
+
+May 21, 1986 by Dyer Lytle Addition of PIMTEXT to package
+
+Pimtext was added to the vacuum telescope package. This program allows
+the user to insert text directly into images. The default action of the
+program is to look up the date and time in the image headers and insert
+this information in the lower left corner of each image. The user can
+modify the parameters to write any text string.
+
+May 15, 1986 by Dyer Lytle Modification to Mscan
+
+Mscan was modified to write the date and time into the images using
+a pixel font. A hidden argument controls this option. The characters
+are written into the image itself to speed up the moviemaking process.
+Various hidden parameters were added to allow the user to specify
+things about the text such as postition, pixel value, background fill,
+and background value.
+
+May 7, 1986 by Dyer Lytle Modification to Makeimages and Destreak5
+
+Makeimages and Destreak5 were modified to accept as another argument
+the input scratch disk on which the input files are to be expected.
+
+February 19, 1986 by Dyer Lytle Modification to Fitslogr
+
+Rfits was changed to produce a short header by default instead of
+a long header. I changed fitslogr to force the long header it needs.
+
+February 6, 1986 by Dyer Lytle Modification to Dicoplot
+
+Dicoplot was plotting all of the dates in the input image header
+file. Sometimes, this list includes dates which should appear
+off the plot, before the zero or after the 360 degree marks.
+The modification involved teaching the program to dump these
+extra dates instead of putting them on the plots.
+
+January 30, 1986 by Dyer Lytle Modification to vtfix
+
+Vtfix was originally set up to correct extra long records on
+vacuum telescope tapes. It looked to record lengths of 10242
+bytes and truncated them to 10240 bytes. Today I found a tape
+with lots of different record lengths all larger than 10240 so
+I modified vtfix to look for records with lengths longer than
+10240 bytes and truncate them to 10240.
+
+January 29, 1986 by Dyer Lytle Modification to makehelium.
+
+Makehelium was modified to automatically delete the absolute
+value image output from RMAP since this image is junk anyway.
+
+January 29, 1986 by Dyer Lytle Bug fix and mods to dicoplot.
+
+Dicoplot had a bug which caused the Ratio (POLARITY) images to
+come out zero. This was corrected. Also some of the constants
+in GREYMAP were changed to increase the contrast in the weights
+image and in the abs. flux image. The greyscale as drawn on the
+images was modified to not have white boxes around each grey level
+and to have the number associated with each grey level printed on the
+plot.
+
+January 28, 1986 by Dyer Lytle Modifications to mscan.
+
+Mscan was using too much memory when processing large images.
+This was causing a lot of page fault errors on VMS. A modification
+was made to mscan to use fixed size subrasters, decreasing the
+memory needs drastically.
+
+January 20, 1986 by Dyer Lytle Modifications to readss4.x.
+
+Readss4, which is a subroutine called by mscan to read type 4
+sector scans was set up to add the average field to each pixel
+of the output image. This was found to be useful only in the
+special case of type 4 intensity scans and was removed.
+"It wasn't a BUG, it was a FEATURE!"
+
+January 20, 1986 by Dyer Lytle Modifications to destreak.x.
+
+Destreak was set up to use a temporary image for data storage
+between the two destreaking passes. The temporary image was
+hardwired into the name "tempim". This was found to unacceptable
+since two or more destreaking jobs run at the same time would have
+a collision at "tempim". The temporary image was made into an input
+parameter.
+
+January 20, 1986 by Dyer Lytle Modifications to CL scripts.
+
+The CL scripts makeimages.cl, makehelium.cl, destreak5.cl, and
+writetape.cl were modified to check for the existence of each file
+before it tries to use it. An error message is output if an image
+cannot be accessed.
+
+January 20, 1986 by Dyer Lytle Modification to vtblink.cl
+
+Vtblink was modified so that the command "stat" can be entered to the
+"next image" prompt and the script will list which images are loaded
+into which IIS memory plane.
+
+January 20, 1986 by Dyer Lytle Modification to merge.x
+
+Merge was not set up to handle the differences between the magnetogram
+reduction and the 10830 reduction. Magnetogram data has three(3) images
+per day and 10830 data has two(2) images per day. The extra image for
+magnetogram data is the absolute value immage. Merge was designed to
+expect all three images and to produce four(4) output images. When
+10830 data is input merge should expect two input images per day and
+only produce two output images. This modification was made.
+Also the output images were set up such that the data and absolute
+value images were output without being divided by the weight image.
+This was changed since no information is lost by doing this division
+since the weight image is also saved. Merge was also restructured
+quite a bit but is still a mess and needs rewriting, but it works.
+
+January 20, 1986 by Dyer Lytle Modification to rmap.x
+
+Rmap was changed to calculate the average field, the average absolute
+field, and the total number of pixels for each gram reduced.
+These parameters are stored in the reduced data image header as
+MEAN_FLD, MEANAFLD, and NUM_PIX.
+
+January 10, 1986 by Dyer Lytle Bug fix in tcopy.
+
+Tcopy was reporting errors incorrectly. The record number identified
+with the error was one less than the actual error record.
+
+January 10, 1986 by Dyer Lytle Modification to decodeheader.x.
+
+Changed the format used by decodeheader to print out the date and time,
+the format was of variable width depending on the size of the number printed.
+The new format has fixed length fields.
+
+January 9, 1986 by Dyer Lytle Modification to merge.
+
+Merge was modified to expect the images in the textfile 'mergelist' to be in the
+order (data, abs value, weights) instead of (data, weights, abs value).
+
+January 3, 1986 by Dyer Lytle Correction to dicoplot.
+
+Dicoplot had, for some integer expressions, TRUE/FALSE instead of YES/NO.
+This works fine on the UNIX system but was found to fail on VMS.
+
+January 3, 1986 by Dyer Lytle Correction to mscan.
+
+Mscan was not reading type one(1) area scans properly. The error occurred
+in readss1 where a temporary array was being salloced with the wrong length.
+The correction involved replacing "ny" by "2*ny".
+Also, readss1 and readss3 had a rather contrived error recovery mechanism built
+in, I removed this and will add a more general and reliable error procedure
+based on the fset(VALIDATE) call in a later revision.

noao/imred/vtel/asciilook.inc

@@ -0,0 +1,19 @@
+data (asciilook[i], i=1,7) / 449, 449, 449, 449, 449, 449, 449 /
+data (asciilook[i], i=8,14) / 449, 449, 449, 449, 449, 449, 449 /
+data (asciilook[i], i=15,21) / 449, 449, 449, 449, 449, 449, 449 /
+data (asciilook[i], i=22,28) / 449, 449, 449, 449, 449, 449, 449 /
+data (asciilook[i], i=29,35) / 449, 449, 449, 449, 001, 008, 015 /
+data (asciilook[i], i=36,42) / 022, 029, 036, 043, 050, 057, 064 /
+data (asciilook[i], i=43,49) / 071, 078, 085, 092, 099, 106, 113 /
+data (asciilook[i], i=50,56) / 120, 127, 134, 141, 148, 155, 162 /
+data (asciilook[i], i=57,63) / 169, 176, 183, 190, 197, 204, 211 /
+data (asciilook[i], i=64,70) / 218, 225, 232, 239, 246, 253, 260 /
+data (asciilook[i], i=71,77) / 267, 274, 281, 288, 295, 302, 309 /
+data (asciilook[i], i=78,84) / 316, 323, 330, 337, 344, 351, 358 /
+data (asciilook[i], i=85,91) / 365, 372, 379, 386, 393, 400, 407 /
+data (asciilook[i], i=92,98) / 414, 421, 428, 435, 442, 449, 232 /
+data (asciilook[i], i=99,105) / 239, 246, 253, 260, 267, 274, 281 /
+data (asciilook[i], i=106,112) / 288, 295, 302, 309, 316, 323, 330 /
+data (asciilook[i], i=113,119) / 337, 344, 351, 358, 365, 372, 379 /
+data (asciilook[i], i=120,126) / 386, 393, 400, 407, 449, 449, 449 /
+data (asciilook[i], i=127,128) / 449, 449/

noao/imred/vtel/d1900.x

@@ -0,0 +1,15 @@
+# D1900 -- Function to return the number of days since the turn of the
+# century.
+
+int procedure d1900 (month, day, year)
+
+int month, day, year # m,d,y of date
+
+int mac[12]
+data mac/0, 31, 59, 90, 120, 151, 181, 212, 243, 273, 304, 334/
+
+begin
+ d1900 = 365 * year + (year - 1) / 4 + mac[month] + day
+ if (month >= 3 && mod(year,4) == 0)
+ d1900 = d1900 + 1
+end