--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/orig/optics2rad.csh	Sat Apr 09 17:55:16 2011 +0100
@@ -0,0 +1,97 @@
+#!/bin/csh -f
+# RCSid $Id: optics2rad.csh,v 2.3 2008/08/25 04:50:31 greg Exp $
+#
+# Convert Optics 5 output to correct Radiance input
+#
+if ($#argv < 1) then
+	echo "Usage: $0 optics.mat .."
+	exit 1
+endif
+set optf=`mktemp /tmp/optfmt.XXXXXXX`
+cat > $optf << '_EOF_'
+void glass $(name)_glass
+0
+0
+3 ${Rtn} ${Gtn} ${Btn}
+
+void BRTDfunc $(name)_front
+10
+	${fRrho} ${fGrho} ${fBrho}
+	${Rtau} ${Gtau} ${Btau}
+     0 0 0
+     .
+0
+9 0 0 0 0 0 0 0 0 0
+
+void BRTDfunc $(name)_back
+10
+	${bRrho} ${bGrho} ${bBrho}
+	${Rtau} ${Gtau} ${Btau}
+     0 0 0
+     .
+0
+9 0 0 0 0 0 0 0 0 0
+'_EOF_'
+
+echo "# Output generated by $0 from $*"
+sed -e '/^[^#]/d' -e '/^$/d' $*
+
+set outf=`mktemp /tmp/outfmt.XXXXXX`
+cat > $outf << '_EOF_'
+
+void glass $(name)
+0
+0
+3 ${Rtn} ${Gtn} ${Btn}
+'_EOF_'
+
+rcalc -l -e 'abs(x):if(x,x,-x);and(a,b):if(a,b,a);not(x):if(x,-1,1)' \
+	-e 'lum(r,g,b):.265*r+.670*g+.065*b' \
+	-e 'trans=lum(Rtau,Gtau,Btau)' \
+	-e 'rfront=lum(fRrho,fGrho,fBrho);rback=lum(bRrho,bGrho,bBrho)' \
+	-e 'cond=0.005-abs(rfront-rback)' \
+	-i $optf -o $outf $*
+
+cat > $outf << '_EOF_'
+
+void BRTDfunc $(name)
+10
+	rR_clear rG_clear rB_clear
+	${Rtau}*tR_clear ${Gtau}*tG_clear ${Btau}*tB_clear
+     0 0 0
+     window.cal
+0
+15 0 0 0 0 0 0 0 0 0
+	${fRrho} ${fGrho} ${fBrho}
+	${bRrho} ${bGrho} ${bBrho}
+'_EOF_'
+
+rcalc -l -e 'abs(x):if(x,x,-x);and(a,b):if(a,b,a);not(x):if(x,-1,1)' \
+	-e 'lum(r,g,b):.265*r+.670*g+.065*b' \
+	-e 'trans=lum(Rtau,Gtau,Btau)' \
+	-e 'rfront=lum(fRrho,fGrho,fBrho);rback=lum(bRrho,bGrho,bBrho)' \
+	-e 'cond=and(trans-.645, not(0.005-abs(rfront-rback)))' \
+	-i $optf -o $outf $*
+
+cat > $outf << '_EOF_'
+
+void BRTDfunc $(name)
+10
+	rR_bronze rG_bronze rB_bronze
+	${Rtau}*tR_bronze ${Gtau}*tG_bronze ${Btau}*tB_bronze
+     0 0 0
+     window.cal
+0
+15 0 0 0 0 0 0 0 0 0
+	${fRrho} ${fGrho} ${fBrho}
+	${bRrho} ${bGrho} ${bBrho}
+'_EOF_'
+
+rcalc -l -e 'abs(x):if(x,x,-x);and(a,b):if(a,b,a);not(x):if(x,-1,1)' \
+	-e 'lum(r,g,b):.265*r+.670*g+.065*b' \
+	-e 'trans=lum(Rtau,Gtau,Btau)' \
+	-e 'rfront=lum(fRrho,fGrho,fBrho);rback=lum(bRrho,bGrho,bBrho)' \
+	-e 'cond=and(not(trans-.645), not(0.005-abs(rfront-rback)))' \
+	-i $optf -o $outf $*
+
+rm -f $optf $outf