diff options
Diffstat (limited to 'src/cmd/astro/sun.c')
-rw-r--r-- | src/cmd/astro/sun.c | 94 |
1 files changed, 94 insertions, 0 deletions
diff --git a/src/cmd/astro/sun.c b/src/cmd/astro/sun.c new file mode 100644 index 00000000..8821121c --- /dev/null +++ b/src/cmd/astro/sun.c @@ -0,0 +1,94 @@ +#include "astro.h" + +void +sun(void) +{ + double mven, merth, mmars, mjup, msat; + double dmoon, mmoon, gmoon; + double pturbb, pturbl, pturbr, lograd; + + ecc = .01675104 - 4.180e-5 * capt - 1.26e-7*capt2; + incl = 0; + node = 0; + argp = 281.220833 + .0000470684*eday + .000453*capt2 + + .000003*capt3; + mrad = 1; + anom = 358.475845 + .9856002670*eday - .000150*capt2 + - .000003*capt3; + motion = .9856473354; + + dmoon = 350.737681+12.1907491914*eday-.001436*capt2; + gmoon = 11.250889 + 13.2293504490*eday - .003212*capt2; + mmoon = 296.104608 + 13.0649924465*eday + 9.192e-3*capt2; + mven = 212.448 + 1.602121635*eday; + merth = 358.476 + 0.985600267*eday; + mmars = 319.590 + .524024095*eday; + mjup = 225.269 + .083082362*eday; + msat = 175.593 + .033450794*eday; + + dmoon = fmod(dmoon, 360.)*radian; + gmoon = fmod(gmoon, 360.)*radian; + mmoon = fmod(mmoon, 360.)*radian; + mven *= radian; + merth *= radian; + mmars *= radian; + mjup *= radian; + msat *= radian; + + icosadd(sunfp, suncp); + anom += cosadd(4, mmars, merth, mven, mjup)/3600.; + anom += sinadd(5, mmars, merth, mven, mjup, .07884*capt)/3600.; + + incl *= radian; + node *= radian; + argp *= radian; + anom = fmod(anom, 360.)*radian; + +/* + * computation of elliptic orbit + */ + + lambda = anom + argp; + + pturbl = (6910.057 - 17.240*capt - 0.052*capt2)*sin(anom) + + (72.338 - 0.361*capt) * sin(2.*anom) + + (1.054 - 0.001*capt) * sin(3.*anom) + + 0.018 * sin(4.*anom); + + lambda += pturbl*radsec; + + beta = 0.; + + lograd = (30.57e-6 - 0.15e-6*capt) + - (7274.12e-6 - 18.14e-6*capt - 0.05e-6*capt2)*cos(anom) + - (91.38e-6 - 0.46e-6*capt) * cos(2.*anom) + - (1.45e-6 - 0.01e-6*capt) * cos(3.*anom) + - 0.02e-6 * cos(4.*anom); + + pturbl = cosadd(5, mmars, merth, mven, mjup, msat); + pturbl += sinadd(3, dmoon, mmoon, merth) + .9; + pturbl *= radsec; + + pturbb = cosadd(3, merth, mven, mjup); + pturbb += sinadd(3, gmoon, mmoon, dmoon); + pturbb *= radsec; + + pturbr = cosadd(5, mmars, merth, mven, mjup, msat); + pturbr += cosadd(3, dmoon, mmoon, merth); + + lambda += pturbl; + if(lambda > pipi) + lambda -= pipi; + + beta += pturbb; + + lograd = (lograd+pturbr) * 2.30258509; + rad = 1 + lograd * (1 + lograd * (.5 + lograd/6)); + + motion *= radian*mrad*mrad/(rad*rad); + + semi = 961.182; + if(flags['o']) + semi = 959.63; + mag = -26.5; +} |