--- ray/src/cv/pabopto2bsdf.c	2016/02/02 00:42:11	2.29
+++ ray/src/cv/pabopto2bsdf.c	2017/05/16 20:41:03	2.30
@@ -1,5 +1,5 @@
 #ifndef lint
-static const char RCSid[] = "$Id: pabopto2bsdf.c,v 2.29 2016/02/02 00:42:11 greg Exp $";
+static const char RCSid[] = "$Id: pabopto2bsdf.c,v 2.30 2017/05/16 20:41:03 greg Exp $";
 #endif
 /*
  * Load measured BSDF data in PAB-Opto format.
@@ -128,7 +128,10 @@ init_pabopto_inp(const int i, const char *fname)
 	dv[0] = cos(M_PI/180.*inpfile[i].phi)*dv[2];
 	dv[1] = sin(M_PI/180.*inpfile[i].phi)*dv[2];
 	dv[2] = sqrt(1. - dv[2]*dv[2]);
-	pos_from_vec(inpfile[i].igp, dv);
+	if (inpfile[i].theta <= FTINY)
+		inpfile[i].igp[0] = inpfile[i].igp[1] = grid_res/2 - 1;
+	else
+		pos_from_vec(inpfile[i].igp, dv);
 	return(1);
 }