5 |
|
|
6 |
|
Modified from Tregenza.cal |
7 |
|
|
8 |
+ |
Need to specify MF constant before loading (usu. power of two) |
9 |
+ |
|
10 |
|
June 2009 G. Ward |
11 |
|
} |
12 |
|
DEGREE : PI/180; |
17 |
|
alt = Asin(Dz)/DEGREE; |
18 |
|
azi = Atan2(Dx,Dy)/DEGREE; |
19 |
|
|
18 |
– |
MF : 2^2; { Multiplication factor } |
20 |
|
alpha : 90/(MF*7 + .5); { Separation between rows in degrees } |
21 |
|
{ Number of patches per row } |
22 |
|
tnaz(r) : select(r, 30, 30, 24, 24, 18, 12, 6); |
28 |
|
|
29 |
|
r_inc = 360/rnaz(r_row); |
30 |
|
|
31 |
< |
r_azn = if(azi-(360-.5*r_inc), 0, floor((azi+.5*r_inc)/r_inc)); |
31 |
> |
r_azn = if(359.9999-.5*r_inc - azi, floor((azi+.5*r_inc)/r_inc), 0); |
32 |
|
|
33 |
|
rbin = if (-alt, 0, raccum(r_row) + r_azn + 1); |
34 |
|
|