--- ray/src/util/reinhartb.cal	2014/07/25 16:58:20	2.3
+++ ray/src/util/reinhartb.cal	2014/08/26 18:56:19	2.5
@@ -1,4 +1,4 @@
-{ RCSid $Id: reinhartb.cal,v 2.3 2014/07/25 16:58:20 greg Exp $ }
+{ RCSid $Id: reinhartb.cal,v 2.5 2014/08/26 18:56:19 greg Exp $ }
 {
 	Compute Reinhart high-density sky patch bin (rbin)
 
@@ -20,14 +20,14 @@ Atan2(y,x) : posangle(atan2(y,x));
 				{ Default to Tregenza sky }
 MF = 1;
 				{ Default axis orientation (+Y == north) }
-rNx = 0; rNy = 0; rNz = 1;
+rNx = 0; rNy = 0; rNz = -1;
 Ux = 0; Uy = 1; Uz = 0;
 				{ Compute oriented axis angles }
 inc_dz = -Dx*rNx-Dy*rNy-Dz*rNz;
 inc_rx = -Dx*(Uy*rNz-Uz*rNy) - Dy*(Uz*rNx-Ux*rNz) - Dz*(Ux*rNy-Uy*rNx);
-inc_ry = -Dx*Ux-Dy*Uy-Dz*Uz - inc_dz*(rNx*Ux+rNy*Uy+rNz*Uz);
+inc_ry = Dx*Ux+Dy*Uy+Dz*Uz + inc_dz*(rNx*Ux+rNy*Uy+rNz*Uz);
 
-r_alt = Asin(inc_dz)/DEGREE;
+r_alt = Asin(inc_dz);
 r_azi = Atan2(inc_rx,inc_ry);
 
 alpha = 90/(MF*7 + .5);		{ Separation between rows in degrees }