program EXAMPLE c c -------------------------------------------------------------- c This is just an example for computing rectangular position and c velocities of Pluto for 36 dates by step of -30 days : c from January 1st 2000 12h, Julian Date 2451545.0 c to February 15th 1997 12h, Julian Date 2450495.0. c -------------------------------------------------------------- c implicit double precision (a-h,o-z) dimension r(6) c t0 = 2451545.d0 ndat = 36 step = -30 c do i=1,ndat tjd = t0 + (i-1)*step call PLUTO (tjd,r,ierr) write (*,1000) tjd, (r(k),k=1,6) enddo c 1000 format(/' Rectangular Coordinates of Pluto at Date: JD',F16.8 . /' X=',F18.14,' Y=',F18.14,' Z=',F18.14,' au' . /' X.=',F17.14,' Y.=',F17.14,' Z.=',F17.14,' au/day') c stop end * * * RESULTS * * Rectangular Coordinates of Pluto at Date: JD2451545.00000000 X= -9.87600250590294 Y=-27.97918968704569 Z= -5.75368219264631 au X.= .00302869050594 Y.= -.00112766304604 Z.= -.00126513753728 au/day Rectangular Coordinates of Pluto at Date: JD2451515.00000000 X= -9.96681919636484 Y=-27.94522823250150 Z= -5.71570178148893 au X.= .00302575617657 Y.= -.00113642808880 Z.= -.00126688576539 au/day Rectangular Coordinates of Pluto at Date: JD2451485.00000000 X=-10.05754818757939 Y=-27.91100417475840 Z= -5.67766918488472 au X.= .00302285005193 Y.= -.00114517784083 Z.= -.00126862085738 au/day Rectangular Coordinates of Pluto at Date: JD2451455.00000000 X=-10.14819041551763 Y=-27.87651719380927 Z= -5.63958442415157 au X.= .00301996437040 Y.= -.00115396191286 Z.= -.00127036524580 au/day Rectangular Coordinates of Pluto at Date: JD2451425.00000000 X=-10.23874586605841 Y=-27.84176606255038 Z= -5.60144718904665 au X.= .00301706233016 Y.= -.00116278606321 Z.= -.00127211814575 au/day Rectangular Coordinates of Pluto at Date: JD2451395.00000000 X=-10.32921393589269 Y=-27.80674959287793 Z= -5.56325722697165 au X.= .00301413620146 Y.= -.00117165441405 Z.= -.00127388191783 au/day Rectangular Coordinates of Pluto at Date: JD2451365.00000000 X=-10.41959337497062 Y=-27.77146629084686 Z= -5.52501416120429 au X.= .00301114465198 Y.= -.00118057061670 Z.= -.00127565568921 au/day Rectangular Coordinates of Pluto at Date: JD2451335.00000000 X=-10.50988175325351 Y=-27.73591538300988 Z= -5.48671803576693 au X.= .00300807004377 Y.= -.00118948556946 Z.= -.00127741544158 au/day Rectangular Coordinates of Pluto at Date: JD2451305.00000000 X=-10.60007721155225 Y=-27.70009741285859 Z= -5.44836943805137 au X.= .00300495667984 Y.= -.00119837652624 Z.= -.00127915537905 au/day Rectangular Coordinates of Pluto at Date: JD2451275.00000000 X=-10.69017888983058 Y=-27.66401289902811 Z= -5.40996883960385 au X.= .00300181715650 Y.= -.00120725508309 Z.= -.00128088187123 au/day Rectangular Coordinates of Pluto at Date: JD2451245.00000000 X=-10.78018611116573 Y=-27.62766233362997 Z= -5.37151670788036 au X.= .00299866595412 Y.= -.00121611300285 Z.= -.00128259121691 au/day Rectangular Coordinates of Pluto at Date: JD2451215.00000000 X=-10.87009905195423 Y=-27.59104604485563 Z= -5.33301338337801 au X.= .00299553184685 Y.= -.00122497877966 Z.= -.00128429877187 au/day Rectangular Coordinates of Pluto at Date: JD2451185.00000000 X=-10.95991782539640 Y=-27.55416315698773 Z= -5.29445866705867 au X.= .00299237872877 Y.= -.00123388667627 Z.= -.00128601649679 au/day Rectangular Coordinates of Pluto at Date: JD2451155.00000000 X=-11.04964128534106 Y=-27.51701263826687 Z= -5.25585242601208 au X.= .00298917863339 Y.= -.00124281643590 Z.= -.00128773181919 au/day Rectangular Coordinates of Pluto at Date: JD2451125.00000000 X=-11.13926819705775 Y=-27.47959400296240 Z= -5.21719478673037 au X.= .00298594181810 Y.= -.00125176313553 Z.= -.00128944412030 au/day Rectangular Coordinates of Pluto at Date: JD2451095.00000000 X=-11.22879723003038 Y=-27.44190670328977 Z= -5.17848583213434 au X.= .00298265034265 Y.= -.00126072269906 Z.= -.00129115038388 au/day Rectangular Coordinates of Pluto at Date: JD2451065.00000000 X=-11.31822683180701 Y=-27.40395085113462 Z= -5.13972596540844 au X.= .00297932164853 Y.= -.00126966386751 Z.= -.00129283785913 au/day Rectangular Coordinates of Pluto at Date: JD2451035.00000000 X=-11.40755659114986 Y=-27.36572686824480 Z= -5.10091561315320 au X.= .00297599584276 Y.= -.00127860571201 Z.= -.00129451990913 au/day Rectangular Coordinates of Pluto at Date: JD2451005.00000000 X=-11.49678638048157 Y=-27.32723418457539 Z= -5.06205469745719 au X.= .00297265072222 Y.= -.00128757785767 Z.= -.00129620840115 au/day Rectangular Coordinates of Pluto at Date: JD2450975.00000000 X=-11.58591528326430 Y=-27.28847197663094 Z= -5.02314309847962 au X.= .00296927112609 Y.= -.00129657363076 Z.= -.00129789881457 au/day Rectangular Coordinates of Pluto at Date: JD2450945.00000000 X=-11.67494224689124 Y=-27.24943937660896 Z= -4.98418067366769 au X.= .00296584898632 Y.= -.00130560734213 Z.= -.00129959796069 au/day Rectangular Coordinates of Pluto at Date: JD2450915.00000000 X=-11.76386531204222 Y=-27.21013524761163 Z= -4.94516721741911 au X.= .00296233854197 Y.= -.00131466778143 Z.= -.00130129669883 au/day Rectangular Coordinates of Pluto at Date: JD2450885.00000000 X=-11.85268175446358 Y=-27.17055957360837 Z= -4.90610311681482 au X.= .00295874908425 Y.= -.00132370441094 Z.= -.00130297228930 au/day Rectangular Coordinates of Pluto at Date: JD2450855.00000000 X=-11.94138992294454 Y=-27.13071326043074 Z= -4.86698908618445 au X.= .00295512391064 Y.= -.00133271333500 Z.= -.00130462701611 au/day Rectangular Coordinates of Pluto at Date: JD2450825.00000000 X=-12.02998889171437 Y=-27.09059699251057 Z= -4.82782568273014 au X.= .00295147074120 Y.= -.00134169960427 Z.= -.00130626286092 au/day Rectangular Coordinates of Pluto at Date: JD2450795.00000000 X=-12.11847821732684 Y=-27.05021153053962 Z= -4.78861349450148 au X.= .00294782187756 Y.= -.00135066391147 Z.= -.00130788194263 au/day Rectangular Coordinates of Pluto at Date: JD2450765.00000000 X=-12.20685838691573 Y=-27.00955687229443 Z= -4.74935268702747 au X.= .00294418798302 Y.= -.00135965537580 Z.= -.00130950796911 au/day Rectangular Coordinates of Pluto at Date: JD2450735.00000000 X=-12.29512908551758 Y=-26.96863169900186 Z= -4.71004288972583 au X.= .00294051519345 Y.= -.00136869570207 Z.= -.00131114594569 au/day Rectangular Coordinates of Pluto at Date: JD2450705.00000000 X=-12.38328871078026 Y=-26.92743487011751 Z= -4.67068392867023 au X.= .00293678397085 Y.= -.00137776370654 Z.= -.00131278493993 au/day Rectangular Coordinates of Pluto at Date: JD2450675.00000000 X=-12.47133546756759 Y=-26.88596559170406 Z= -4.63127576907705 au X.= .00293298575007 Y.= -.00138685866081 Z.= -.00131442546274 au/day Rectangular Coordinates of Pluto at Date: JD2450645.00000000 X=-12.55926692723893 Y=-26.84422336645171 Z= -4.59181853543207 au X.= .00292909790703 Y.= -.00139595146988 Z.= -.00131605257268 au/day Rectangular Coordinates of Pluto at Date: JD2450615.00000000 X=-12.64708084115579 Y=-26.80220889814741 Z= -4.55231288200450 au X.= .00292515947003 Y.= -.00140500785727 Z.= -.00131765433237 au/day Rectangular Coordinates of Pluto at Date: JD2450585.00000000 X=-12.73477632549175 Y=-26.75992300492028 Z= -4.51275937547623 au X.= .00292120151935 Y.= -.00141405214586 Z.= -.00131924528681 au/day Rectangular Coordinates of Pluto at Date: JD2450555.00000000 X=-12.82235254753751 Y=-26.71736577419868 Z= -4.47315824418078 au X.= .00291720708523 Y.= -.00142309533036 Z.= -.00132082821641 au/day Rectangular Coordinates of Pluto at Date: JD2450525.00000000 X=-12.90980849202301 Y=-26.67453730918752 Z= -4.43350976395465 au X.= .00291318549173 Y.= -.00143213790080 Z.= -.00132240374364 au/day Rectangular Coordinates of Pluto at Date: JD2450495.00000000 X=-12.99714325531504 Y=-26.63143722356541 Z= -4.39381397485578 au X.= .00290912027651 Y.= -.00144120570031 Z.= -.00132398243469 au/day