Repatch (from linux) of OSX IRAF

1 files changed, 96 insertions, 0 deletions

diff --git a/pkg/utilities/nttools/tmatch/getweight.x b/pkg/utilities/nttools/tmatch/getweight.x
new file mode 100644
index 00000000..ef8b9204
--- /dev/null
+++ b/pkg/utilities/nttools/tmatch/getweight.x
@@ -0,0 +1,96 @@
+include	<math.h>
+include	<tbset.h>
+
+#* HISTORY *
+#* B.Simon	24-Aug-94	original
+
+# GETWEIGHT -- Get array of weights from list of factors or tables
+
+procedure getweight (ncol, col1, col2, factor, weight)
+
+int	ncol		# i: number of match columns
+pointer	col1[ARB]	# i: match columns from first table
+pointer	col2[ARB]	# i: match columns from second table
+char	factor[ARB]	# i: list of factors
+double	weight[ARB]	# o: array of weights
+#--
+double	unitval[6]
+int	invert[6]
+int	ic, jc, nc, icol, jcol, type1, type2, item
+pointer	sp, value, unit1, unit2, errmsg
+
+
+data	unitval	  / 1.0, 3600.0, 60.0, 1.0, 15.0, RADIAN /
+data	invert	  / NO,  YES,    YES,  NO,  NO,   NO     /
+
+string	unitlist  "|seconds|minutes|degrees|hours|radians|"
+string	badvalue  "Value in factor string is not a number (%s)"
+string	badunits  "Units mismatch in column %d of tables"
+
+int	ctod(), word_fetch(), strdic()
+
+begin
+	# Allocate memory for temporary strings
+
+	call smark (sp)
+	call salloc (value, SZ_FNAME, TY_CHAR)
+	call salloc (unit1, SZ_FNAME, TY_CHAR)
+	call salloc (unit2, SZ_FNAME, TY_CHAR)
+	call salloc (errmsg, SZ_LINE, TY_CHAR)
+
+	# Get each string from the list and convert to a number
+
+	ic = 1
+	icol = 0
+	while (word_fetch (factor, ic, Memc[value], SZ_FNAME) > 0) {
+	    icol = icol + 1
+
+	    jc = 1
+	    nc = ctod (Memc[value], jc, weight[icol])
+	    if (Memc[value+jc-1] != EOS) {
+		call sprintf (Memc[errmsg], SZ_LINE, badvalue)
+		call pargstr (Memc[value])
+
+		call error (1, Memc[errmsg])
+	    }
+	}
+
+	# Set remaining weights according to column units
+
+	do jcol = icol+1, ncol {
+	    # Read units from table
+
+	    call tbcigt (col1[jcol], TBL_COL_UNITS, Memc[unit1], SZ_FNAME)
+	    call tbcigt (col2[jcol], TBL_COL_UNITS, Memc[unit2], SZ_FNAME)
+
+	    # Search for units in dictionary
+
+	    call strlwr (Memc[unit1])
+	    call strlwr (Memc[unit2])
+
+	    type1 = strdic (Memc[unit1], Memc[unit1], SZ_FNAME, unitlist)
+	    type2 = strdic (Memc[unit2], Memc[unit2], SZ_FNAME, unitlist)
+
+	    # Take exit if units do not match
+
+	    if (type1 != type2) {
+		call sprintf (Memc[errmsg], SZ_LINE, badunits)
+		call pargi (jcol)
+
+		call error (1, Memc[errmsg])
+	    }
+
+	    # Read corresponding weight from unit value array
+	    # The first weight (1.0) is for missing or unknown units
+
+	    item = type1 + 1
+
+	    if (invert[item] == NO) {
+		weight[jcol] = unitval[item]
+	    } else {
+		weight[jcol] = 1.0 / unitval[item]
+	    }
+	}
+
+	call sfree (sp)
+end