| 282 | // astrological vertex is also calculated here. |
| 283 | |
| 284 | void ComputeVariables(real *vtx) |
| 285 | { |
| 286 | real B, L, G, tim = TT; |
| 287 | |
| 288 | if (us.fProgress && us.nProgress != ptSolarArc) { |
| 289 | tim = JulianDayFromTime(is.T) + 0.5; |
| 290 | tim = RFract(tim)*24.0; |
| 291 | } |
| 292 | is.RA = Mod((6.6460656 + 2400.0513*is.T + 2.58E-5*is.T*is.T + |
| 293 | tim)*15.0 - OO); |
| 294 | L = RFromD(is.RA) + rPi; B = rPiHalf - RAbs(RFromD(AA)); |
| 295 | if (AA < 0.0) |
| 296 | B = -B; |
| 297 | G = RecToSph(B, L, -is.OB); |
| 298 | *vtx = Mod(is.rSid + DFromR(G+rPiHalf)); // Vertex |
| 299 | } |
| 300 | |
| 301 | |
| 302 | /* |