diff options
author | wkj <devnull@localhost> | 2004-04-21 02:16:43 +0000 |
---|---|---|
committer | wkj <devnull@localhost> | 2004-04-21 02:16:43 +0000 |
commit | cd5bae7871bc0f0bc68b4d2a84703929a7a3c9d1 (patch) | |
tree | 336db54785d2b77113a6e570574be715c7eb7d1d /src/cmd/astro/mars.c | |
parent | 95f57b01e21feb457e79eaf52d593422c318024f (diff) | |
download | plan9port-cd5bae7871bc0f0bc68b4d2a84703929a7a3c9d1.tar.gz plan9port-cd5bae7871bc0f0bc68b4d2a84703929a7a3c9d1.tar.bz2 plan9port-cd5bae7871bc0f0bc68b4d2a84703929a7a3c9d1.zip |
Astro with some minor changes to placate Unix.
Diffstat (limited to 'src/cmd/astro/mars.c')
-rw-r--r-- | src/cmd/astro/mars.c | 67 |
1 files changed, 67 insertions, 0 deletions
diff --git a/src/cmd/astro/mars.c b/src/cmd/astro/mars.c new file mode 100644 index 00000000..327ab8d2 --- /dev/null +++ b/src/cmd/astro/mars.c @@ -0,0 +1,67 @@ +#include "astro.h" + +void +mars(void) +{ + double pturbl, pturbb, pturbr; + double lograd; + double dele, enom, vnom, nd, sl; + double lsun, elong, ci, dlong; + + + ecc = .09331290 + .000092064*capt; + incl = 1.850333 - 6.75e-4*capt; + node = 48.786442 + .770992*capt; + argp = 334.218203 + 1.840758*capt + 1.30e-4*capt2; + mrad = 1.5236915; + anom = 319.529425 + .5240207666*eday + 1.808e-4*capt2; + motion = 0.5240711638; + + + incl = incl*radian; + node = node*radian; + argp = argp*radian; + anom = fmod(anom,360.)*radian; + + enom = anom + ecc*sin(anom); + do { + dele = (anom - enom + ecc * sin(enom)) / + (1. - ecc*cos(enom)); + enom += dele; + } while(fabs(dele) > converge); + vnom = 2.*atan2(sqrt((1.+ecc)/(1.-ecc))*sin(enom/2.),cos(enom/2.)); + rad = mrad*(1. - ecc*cos(enom)); + + lambda = vnom + argp; + pturbl = 0.; + lambda = lambda + pturbl*radsec; + pturbb = 0.; + pturbr = 0.; + +/* + * reduce to the ecliptic + */ + + nd = lambda - node; + lambda = node + atan2(sin(nd)*cos(incl),cos(nd)); + + sl = sin(incl)*sin(nd) + pturbb*radsec; + beta = atan2(sl, pyth(sl)); + + lograd = pturbr*2.30258509; + rad *= 1. + lograd; + + + motion *= radian*mrad*mrad/(rad*rad); + semi = 4.68; + + lsun = 99.696678 + 0.9856473354*eday; + lsun *= radian; + elong = lambda - lsun; + ci = (rad - cos(elong))/sqrt(1. + rad*rad - 2.*rad*cos(elong)); + dlong = atan2(pyth(ci), ci)/radian; + mag = -1.30 + .01486*dlong; + + helio(); + geo(); +} |