12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273747576777879808182838485868788899091929394 |
- #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;
- }
|