| 13 |
|
mylong = atan2(Ny,Nx); |
| 14 |
|
|
| 15 |
|
latlong = if(LATWIDTH-mod(Acos(Nz)+LATWIDTH/2,LATSTEP), .05, |
| 16 |
< |
if(LONGWIDTH-mod(mylong+LONGWIDTH/2,90*DEGREE), .01, |
| 16 |
> |
if(LONGWIDTH-mod(mylong+LONGWIDTH/2,90*DEGREE), .02, |
| 17 |
|
if(LONGWIDTH-mod(mylong+LONGWIDTH/2,LONGSTEP), .3, 1))); |
| 18 |
|
|
| 19 |
|
{ Compute up vector corresponding to sphere position } |
| 20 |
|
ur2 = Nx*Nx + Ny*Ny; |
| 21 |
|
upx = Nx*Ny*(Nz - 1); |
| 22 |
< |
upy = if(ur2 - FTINY, Nx*Nx + Ny*Ny*Nz, 1); |
| 22 |
> |
upy = if(ur2 - FTINY, Nx*Nx + Ny*Ny*Nz, Nz); |
| 23 |
|
upz = -Ny*ur2; |
| 24 |
|
|
| 25 |
|
interp_arr2`(i,x,f):(i+1-x)*f(i)+(x-i)*f(i+1); |
