Initial commit

1 files changed, 208 insertions, 0 deletions

diff --git a/noao/twodspec/multispec/msget.x b/noao/twodspec/multispec/msget.x
new file mode 100644
index 00000000..e187015a
--- /dev/null
+++ b/noao/twodspec/multispec/msget.x
@@ -0,0 +1,208 @@
+include <imhdr.h>
+include "ms.h"
+
+# MSGET -- Allocate memory and get data from the MULTISPEC database
+# and associated image.
+#
+# MSGHDR	-- Allocate memory and get MULTISPEC header information.
+# MSGCOMMENTS	-- Allocate memory and get MULTISPEC comments.
+# MSGPARAM	-- Allocate memory and get a line of MULTISPEC parameter data.
+# MSGSAMPLE	-- Allocate memory and get SAMPLE line array.
+# MSGFIT	-- Get parameter fit for a spectrum.
+# MSGFITS	-- Get parameter fit for all spectra.
+# MSGGAUSS5	-- Get a line of GAUSS5 parameter data.
+# MSGIMAGE	-- Get a line of the image with possible averaging.
+
+
+# MSGHDR -- Allocate memory and get MULTISPEC header information.
+
+procedure msghdr (ms)
+
+pointer	ms				# MULTISPEC data structure
+
+int	i
+
+int	dbread()
+
+begin
+	if (MS_DATA(ms, HDR) == NULL)
+	    call calloc (MS_DATA(ms, HDR), LEN_MS_HDR, TY_STRUCT)
+	i = dbread (MS_DB(ms), NAME(ms, HDR), HEADER(ms), 1)
+end
+
+# MSGCOMMENTS -- Allocate memory and get MULTISPEC comments.
+
+procedure msgcomments (ms)
+
+pointer	ms				# MULTISPEC data structure
+
+int	i
+
+int	dbread()
+
+begin
+	if (MS_DATA(ms, COMMENTS) == NULL)
+	    call calloc (MS_DATA(ms, COMMENTS), SZ_MS_COMMENTS, TY_CHAR)
+	i = dbread (MS_DB(ms), NAME(ms, COMMENTS), COMMENT(ms, 1), 1)
+end
+
+# MSGPARAM -- Allocate memory and get a line of MULTISPEC parameter data.
+
+procedure msgparam (ms, parameter, line)
+
+pointer	ms				# MULTISPEC data structure
+int	parameter			# Parameter ID
+int	line				# Sample line to be obtained
+
+int	i
+char	reference[SZ_MS_KEY]
+
+bool	is_param_id()
+int	dbread()
+
+begin
+	# Check if the the requested parameter is valid.
+	if (!is_param_id (parameter))
+	    call error (MS_ERROR, "Bad parameter identifier")
+
+	if (MS_DATA(ms, parameter) == NULL)
+	    call calloc (MS_DATA(ms, parameter), MS_NSPECTRA(ms), TY_REAL)
+
+	# Make reference to the desired database record.
+	call sprintf (reference, SZ_MS_KEY, "%s[%d]")
+	    call pargstr (NAME(ms, parameter))
+	    call pargi (line)
+
+	i = dbread (MS_DB(ms), reference, PARAMETER(ms, parameter, 1), 1)
+end
+
+# MSGSAMPLE -- Allocate memory and get SAMPLE line array.
+
+procedure msgsample (ms)
+
+pointer	ms				# MULTISPEC data structure
+
+int	i
+
+int	dbread()
+
+begin
+	if (MS_DATA(ms, SAMPLE) == NULL)
+	    call malloc (MS_DATA(ms, SAMPLE), MS_NSAMPLES(ms), TY_INT)
+	i = dbread (MS_DB(ms), NAME(ms, SAMPLE), LINE(ms,1), 1)
+end
+
+
+# MSGFIT -- Get parameter fit for a spectrum.
+
+procedure msgfit (ms, parameter, spectrum)
+
+pointer	ms				# MULTISPEC data structure
+int	parameter			# Parameter ID for desired fit
+int	spectrum			# Spectrum
+
+int	i
+char	reference[SZ_MS_KEY]
+pointer	sp, fit
+
+bool	is_fit_id()
+int	dbread()
+
+errchk	cvrestore
+
+begin
+	# Check if for valid parameter id.
+	if (!is_fit_id (parameter))
+	    call error (MS_ERROR, "Bad fit identifier")
+
+	# Allocate memory for the curfit pointers.
+	if (MS_DATA(ms, parameter) == NULL)
+	    call malloc (MS_DATA(ms, parameter), MS_NSPECTRA(ms), TY_INT)
+
+	# Allocate memory for the curfit coefficients.
+	call smark (sp)
+	call salloc (fit, 7 + MS_NSAMPLES(ms), TY_REAL)
+
+	# Reference appropriate data.
+	call sprintf (reference, SZ_MS_KEY, "%s[%d]")
+	    call pargstr (NAME(ms, parameter))
+	    call pargi (spectrum)
+
+	i = dbread (MS_DB(ms), reference, Memr[fit], 1)
+	iferr (call cvrestore (CV(ms, parameter, spectrum), Memr[fit]))
+	    ;
+
+	call sfree (sp)
+end
+
+
+# MSGFITS -- Get parameter fits.
+
+procedure msgfits (ms, parameter)
+
+pointer	ms				# MULTISPEC data structure
+int	parameter			# Parameter ID for desired fit
+
+int	i
+
+begin
+	do i = 1, MS_NSPECTRA(ms)
+	    call msgfit (ms, parameter, i)
+end
+
+
+# MSGGAUSS5 -- Get a line of GAUSS5 parameter data.
+
+procedure msggauss5 (ms, line)
+
+pointer	ms				# MULTISPEC data structure
+int	line				# Sample line to be obtained
+
+begin
+	call msgparam (ms, I0, line)
+	call msgparam (ms, X0, line)
+	call msgparam (ms, S0, line)
+	call msgparam (ms, S1, line)
+	call msgparam (ms, S2, line)
+end
+
+
+# MSGIMAGE -- Get a line of the image with possible averaging.
+
+procedure msgimage (im, line, naverage, data)
+
+pointer	im				# Image descriptor
+int	line				# Line to be gotten from the image
+int	naverage			# Number of line to use in average
+real	data[ARB]			# The output data array
+
+int	i, line_start, line_end
+real	nlines
+pointer	buf
+
+pointer	imgl2r()
+
+begin
+	# If naverage is <= 1 copy the image line to the data array
+	# Else average the several lines.
+
+	if (naverage <= 1) {
+	    call amovr (Memr[imgl2r (im, line)], data, IM_LEN(im,1))
+	} else {
+	    # Determine starting and ending lines for the average.
+	    line_start = max (1, line - naverage / 2)
+	    line_end = min (IM_LEN(im, 2), line_start + naverage - 1)
+
+	    # Clear data array for accumulating sum and then vector
+	    # add the image lines.
+	    call aclrr (data, IM_LEN(im, 1))
+	    do i = line_start, line_end {
+	        buf = imgl2r (im, i)
+	        call aaddr (Memr[buf], data, data, IM_LEN(im, 1))
+	    }
+
+	    # Vector divide by the number of lines to form average.
+	    nlines = line_end - line_start + 1
+	    call adivkr (data, nlines, data, IM_LEN(im, 1))
+	}
+end