--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/orig/colorcal.csh	Sat Apr 09 17:55:16 2011 +0100
@@ -0,0 +1,38 @@
+#!/bin/csh -f
+# RCSid: $Id: colorcal.csh,v 1.1 2003/02/22 02:07:21 greg Exp $
+#
+# Compute CIE chromaticities from spectral reflectance data
+#
+if ( $#argv < 1 ) goto userr
+set cal = .
+if ( $argv[1] == "-i" ) then
+	if ( $#argv < 3 ) goto userr
+	shift argv
+	set illum=$argv[1]
+	shift argv
+	foreach r ( $argv[*] )
+		tabfunc -i rf < $r > /tmp/rf$$.cal
+		rcalc -f $cal/cieresp.cal -f /tmp/rf$$.cal \
+			-e 'r=rf($1);ty=$2*triy($1)' \
+			-e '$1=ty;$2=$2*r*trix($1);$3=r*ty' \
+			-e '$4=$2*r*triz($1)' \
+			-e 'cond=if($1-359,831-$1,-1)' \
+			$illum | total -m >> /tmp/rc$$.dat
+	end
+	rm -f /tmp/rf$$.cal
+else
+	foreach r ( $argv[*] )
+		rcalc -f $cal/cieresp.cal -e 'ty=triy($1);$1=ty' \
+		-e '$2=$2*trix($1);$3=$2*ty;$4=$2*triz($1)' \
+		-e 'cond=if($1-359,831-$1,-1)' $r \
+			| total -m >> /tmp/rc$$.dat
+	end
+endif
+rcalc -e 'X=$2/$1;Y=$3/$1;Z=$4/$1' \
+	-e 'x=X/(X+Y+Z);y=Y/(X+Y+Z);u=4*X/(X+15*Y+3*Z);v=9*Y/(X+15*Y+3*Z)' \
+	-o $cal/color.fmt /tmp/rc$$.dat
+rm -f /tmp/rc$$.dat
+exit 0
+userr:
+echo "Usage: $0 [-i illum.dat] refl.dat .."
+exit 1