--- ray/src/cv/bsdf2rad.cal	2017/04/10 06:09:14	2.2
+++ ray/src/cv/bsdf2rad.cal	2017/04/11 02:34:10	2.5
@@ -1,4 +1,4 @@
-{ RCSid $Id: bsdf2rad.cal,v 2.2 2017/04/10 06:09:14 greg Exp $ }
+{ RCSid $Id: bsdf2rad.cal,v 2.5 2017/04/11 02:34:10 greg Exp $ }
 {
 	Calculation of false color scale for bsdf2rad.
 
@@ -12,9 +12,15 @@ LONGWIDTH : 1.*DEGREE;
 
 mylong = atan2(Ny,Nx);
 
-latlong = if(LATWIDTH-mod(Acos(Nz)+ LATWIDTH/2,LATSTEP), .05,
+latlong = if(LATWIDTH-mod(Acos(Nz)+LATWIDTH/2,LATSTEP), .05,
 		if(LONGWIDTH-mod(mylong+LONGWIDTH/2,90*DEGREE), .01,
-		if(LONGWIDTH-mod(mylong+LONGWIDTH/2,LONGSTEP), .2, 1)));
+		if(LONGWIDTH-mod(mylong+LONGWIDTH/2,LONGSTEP), .3, 1)));
+
+		{ Compute up vector corresponding to sphere position }
+ur2 = Nx*Nx + Ny*Ny;
+upx = Nx*Ny*(Nz - 1);
+upy = if(ur2 - FTINY, Nx*Nx + Ny*Ny*Nz, 1);
+upz = -Ny*ur2;
 
 interp_arr2`(i,x,f):(i+1-x)*f(i)+(x-i)*f(i+1);
 interp_arr`(x,f):if(x-1,if(f(0)-x,interp_arr2(floor(x),x,f),f(f(0))),f(1));