1 files changed, 127 insertions, 0 deletions
diff --git a/noao/mtlocal/pds/t_pdsread.x b/noao/mtlocal/pds/t_pdsread.x
new file mode 100644
index 00000000..94b13320
--- /dev/null
+++ b/noao/mtlocal/pds/t_pdsread.x
@@ -0,0 +1,127 @@
+include <error.h>
+include <fset.h>
+
+# T_PDSREAD -- Read PDS format data.  Further documentation given
+# in pds.hlp
+
+define	MAX_RANGES	100
+
+procedure t_pdsread()
+
+char	infile[SZ_FNAME]		# the input file name list
+char	outfile[SZ_FNAME]		# the output file name list
+char	file_list[SZ_FNAME]		# the input file number list
+int	offset				# the output file name offset
+
+char	in_fname[SZ_FNAME], out_fname[SZ_FNAME]
+int	range[MAX_RANGES*2+1], nfiles, file_number, stat, junk
+int	lenlist
+pointer	list
+
+bool	clgetb()
+char	clgetc()
+int	fstati(), btoi(), clgeti(), fntlenb(), fntgfnb(), mtfile()
+int	mtneedfileno()
+int	pds_read(), decode_ranges(), get_next_number(), pds_get_image_type()
+pointer	fntopnb()
+
+include	"rpds.com"
+
+begin
+	# Set up the standard output to flush on a newline.
+	if (fstati (STDOUT, F_REDIR) == NO)
+	    call fseti (STDOUT, F_FLUSHNL, YES)
+
+	# Get the input file name(s).
+	call clgstr ("pds_file", infile, SZ_FNAME)
+	if (mtfile (infile) == YES) {
+	    list = NULL
+	    if (mtneedfileno (infile) == YES)
+	        call clgstr ("file_list", file_list, SZ_LINE)
+	    else
+	        call strcpy ("1", file_list, SZ_LINE)
+	} else {
+	    list = fntopnb (infile, YES)
+	    lenlist = fntlenb (list)
+	    call sprintf (file_list, SZ_LINE, "1-%d")
+		call pargi (lenlist)
+	}
+
+	# Decode the ranges.
+	if (decode_ranges (file_list, range, MAX_RANGES, nfiles) == ERR)
+	    call error (1, "Illegal file number list")
+
+	# Setup the output options.
+	long_header = btoi (clgetb ("long_header"))
+	short_header = btoi (clgetb ("short_header"))
+	make_image = btoi (clgetb ("make_image"))
+	tenbit = btoi (clgetb ("tenbit"))
+	ninetrack = btoi (clgetb ("ninetrack"))
+	offset = clgeti ("offset")
+
+	# Set the output image data type.
+	if (make_image == YES) {
+	    data_type = pds_get_image_type (clgetc ("datatype"))
+	    call clgstr ("iraf_file", outfile, SZ_FNAME)
+	} else
+	    outfile[1] = EOS
+
+	# Read successive PDS files, convert and write into a numbered
+	# succession of output IRAF files.
+
+	file_number = 0
+	while (get_next_number (range, file_number) != EOF) {
+
+	    # Get the input file name.
+	    if (list != NULL)
+		junk = fntgfnb (list, in_fname, SZ_FNAME)
+	    else {
+	        if (mtneedfileno (infile) == YES)
+		    call mtfname (infile, file_number, in_fname, SZ_FNAME)
+		else
+	            call strcpy (infile, in_fname, SZ_FNAME)
+	    }
+
+	    # Get the output file name.
+	    if (nfiles > 1) {
+		call sprintf (out_fname[1], SZ_FNAME, "%s%03d")
+		    call pargstr (outfile)
+		    call pargi (file_number + offset)
+	    } else
+	        call strcpy (outfile, out_fname, SZ_FNAME)
+
+	    # Convert PDS file to the output IRAF file. If EOT is reached
+	    # then exit. If an error is detected then print a warning and
+	    # continue with the next file.
+
+	    iferr (stat = pds_read (in_fname, out_fname))
+		call erract (EA_FATAL)
+	    if (stat == EOF)	# EOT found
+		break
+	}
+
+	if (list != NULL)
+	    call fntclsb (list)
+end
+
+# GET_IMAGE_TYPE -- Convert a character to an IRAF image type.
+
+define	NTYPES	7
+
+int procedure pds_get_image_type(c)
+
+char	c
+int	i, type_codes[NTYPES]
+int	stridx()
+
+string	types "usilrds"			# supported image data types
+data	type_codes /TY_USHORT, TY_SHORT, TY_INT, TY_LONG, TY_REAL,
+		    TY_DOUBLE, TY_COMPLEX/
+
+begin
+	i = stridx (c, types)
+	if (i == 0)
+	    return (ERR)
+	else
+	    return (type_codes[stridx (c, types)])
+end