13e12c5d1SDavid du Colombier #include "astro.h"
23e12c5d1SDavid du Colombier
33e12c5d1SDavid du Colombier void
jup(void)43e12c5d1SDavid du Colombier jup(void)
53e12c5d1SDavid du Colombier {
63e12c5d1SDavid du Colombier double pturbl, pturbb, pturbr;
73e12c5d1SDavid du Colombier double lograd;
83e12c5d1SDavid du Colombier double dele, enom, vnom, nd, sl;
93e12c5d1SDavid du Colombier
103e12c5d1SDavid du Colombier
113e12c5d1SDavid du Colombier ecc = .0483376 + 163.e-6*capt;
123e12c5d1SDavid du Colombier incl = 1.308660 - .0055*capt;
133e12c5d1SDavid du Colombier node = 99.43785 + 1.011*capt;
143e12c5d1SDavid du Colombier argp = 12.71165 + 1.611*capt;
153e12c5d1SDavid du Colombier mrad = 5.202803;
163e12c5d1SDavid du Colombier anom = 225.22165 + .0830912*eday - .0484*capt;
173e12c5d1SDavid du Colombier motion = 299.1284/3600.;
183e12c5d1SDavid du Colombier
193e12c5d1SDavid du Colombier
203e12c5d1SDavid du Colombier anom = anom;
213e12c5d1SDavid du Colombier incl *= radian;
223e12c5d1SDavid du Colombier node *= radian;
233e12c5d1SDavid du Colombier argp *= radian;
243e12c5d1SDavid du Colombier anom = fmod(anom,360.)*radian;
253e12c5d1SDavid du Colombier
263e12c5d1SDavid du Colombier enom = anom + ecc*sin(anom);
273e12c5d1SDavid du Colombier do {
283e12c5d1SDavid du Colombier dele = (anom - enom + ecc * sin(enom)) /
293e12c5d1SDavid du Colombier (1. - ecc*cos(enom));
303e12c5d1SDavid du Colombier enom += dele;
31*7dd7cddfSDavid du Colombier } while(fabs(dele) > converge);
323e12c5d1SDavid du Colombier vnom = 2.*atan2(sqrt((1.+ecc)/(1.-ecc))*sin(enom/2.),
333e12c5d1SDavid du Colombier cos(enom/2.));
343e12c5d1SDavid du Colombier rad = mrad*(1. - ecc*cos(enom));
353e12c5d1SDavid du Colombier
363e12c5d1SDavid du Colombier lambda = vnom + argp;
373e12c5d1SDavid du Colombier
383e12c5d1SDavid du Colombier pturbl = 0.;
393e12c5d1SDavid du Colombier
403e12c5d1SDavid du Colombier lambda += pturbl*radsec;
413e12c5d1SDavid du Colombier
423e12c5d1SDavid du Colombier pturbb = 0.;
433e12c5d1SDavid du Colombier
443e12c5d1SDavid du Colombier pturbr = 0.;
453e12c5d1SDavid du Colombier
463e12c5d1SDavid du Colombier /*
473e12c5d1SDavid du Colombier * reduce to the ecliptic
483e12c5d1SDavid du Colombier */
493e12c5d1SDavid du Colombier
503e12c5d1SDavid du Colombier nd = lambda - node;
513e12c5d1SDavid du Colombier lambda = node + atan2(sin(nd)*cos(incl),cos(nd));
523e12c5d1SDavid du Colombier
533e12c5d1SDavid du Colombier sl = sin(incl)*sin(nd) + pturbb*radsec;
543e12c5d1SDavid du Colombier beta = atan2(sl, pyth(sl));
553e12c5d1SDavid du Colombier
563e12c5d1SDavid du Colombier lograd = pturbr*2.30258509;
573e12c5d1SDavid du Colombier rad *= 1. + lograd;
583e12c5d1SDavid du Colombier
593e12c5d1SDavid du Colombier
603e12c5d1SDavid du Colombier lambda += 555.*radsec;
613e12c5d1SDavid du Colombier beta -= 51.*radsec;
623e12c5d1SDavid du Colombier motion *= radian*mrad*mrad/(rad*rad);
633e12c5d1SDavid du Colombier semi = 98.47;
643e12c5d1SDavid du Colombier
653e12c5d1SDavid du Colombier mag = -8.93;
663e12c5d1SDavid du Colombier helio();
673e12c5d1SDavid du Colombier geo();
683e12c5d1SDavid du Colombier }
69