| 213 | // the sidereal, and a couple of other things. |
| 214 | |
| 215 | real ProcessInput(void) |
| 216 | { |
| 217 | real Ln, Off; |
| 218 | |
| 219 | // Compute angle that the ecliptic is inclined to the Celestial Equator. |
| 220 | is.OB = 23.452294 - 0.0130125*is.T; |
| 221 | |
| 222 | Ln = Mod((933060-6962911*is.T+7.5*is.T*is.T)/3600.0); // Mean lunar node |
| 223 | Off = (259205536.0*is.T+2013816.0)/3600.0; // Mean Sun |
| 224 | Off = 17.23*RSin(RFromD(Ln)) + 1.27*RSin(RFromD(Off)) - |
| 225 | (5025.64+1.11*is.T)*is.T; |
| 226 | Off = (Off-84038.27)/3600.0; |
| 227 | is.rSid = (us.fSidereal ? Off + us.rZodiacOffset : 0.0) + |
| 228 | us.rZodiacOffsetAll; |
| 229 | return Off; |
| 230 | } |
| 231 | |
| 232 | |
| 233 | // Another modulus function, this time for the range of 0 to 2 Pi. |