--- ray/src/cal/cal/tregenza.cal 2009/06/14 20:07:59 1.2 +++ ray/src/cal/cal/tregenza.cal 2009/07/02 21:44:33 1.3 @@ -1,4 +1,4 @@ -{ RCSid $Id: tregenza.cal,v 1.2 2009/06/14 20:07:59 greg Exp $ } +{ RCSid $Id: tregenza.cal,v 1.3 2009/07/02 21:44:33 greg Exp $ } { Compute bin numbers corresponding to Tregenza sky discretization @@ -26,7 +26,7 @@ Atan2(y,x) : posangle(atan2(y,x)); alt = Asin(Dz)/DEGREE; azi = Atan2(Dx,Dy)/DEGREE; -tazi(inc) = if(azi-(360-.5*inc), 0, floor((azi+.5*inc)/inc)); +tazi(inc) = if(360-.5*inc - azi, floor((azi+.5*inc)/inc), 0); tbin = if(-alt, 0, select(floor(alt/12) + 1,