PDS_VERSION_ID = PDS3 RECORD_TYPE = STREAM OBJECT = TEXT NOTE = "EPGIOTTO.TXT file for GEM Data Set Collection Volume" PUBLICATION_DATE = 1995-12-31 END_OBJECT END ***** File EPGIOTTO.TXT GIOTTO EXPERIMENTER SUPPORT DATA. I. Abstract: A new experimenter support tape (EXSUP5) was made available to investigating teams and was provided for this disk by Trevor Morley. The Giotto Ephemeris on that tape represented the most recent predicted orbit. The documentation enclosed is adapted from the documentation for tape EXSUP5, which updated EXSUP4 with the following amendments:- a detailed tape layout, a sample listing of part of the position data file, and a sample of position/velocities. In addition, for EXSUP5, the subroutine PETRA, used to compute Kepler orbits at the home sites of experimenters, was contained on the tape and is provided here as a source listing found in appendix A4. II. Table of Contents: 1. Introduction. 2. Description of Tape EXSUP5. 3. The Position Data File. 3.1 The Control Record. 3.2 The Position Data Record. 4. The Interface/Auxiliary Software. 4.1 Subroutine GIODAT ( interface routine to the sequential formatted Position Data File ). 4.2 Subroutine TRAVEL ( generation of experimenter support data ). 4.3 Subroutine PETRA ( compute Kepler Orbits ). 5. Station Coordinates. 5.1 General. 5.2 Table of Station Coordinates. Appendix: A1 The Position Data File A2 Source listing of subroutine GIODAT. A3 Source listing of subroutine TRAVEL. A4 Source listing of subroutine PETRA. 1. Introduction. In this document, the Position Data, the interface routine to this data and the software for generation of experimenter support data, will be described. To ensure a safe 'transport' of the Position Data and to limit potential conversion problems between different computers it has been decided to deliver this as a sequential access file containing formatted data. In addition to this data file the following interface/auxiliary software can be found on the tape:- Software for:- - Acquisition of Position Data from the sequential file and determination of the approximate velocities derived from the position data. - Generation of Giotto experimenter support data ( signal travel time, etc. ). In order to reduce the processing time, the use of a binary direct access file for the Position Data is recommended. Software for generation of a direct access file from the formatted sequential file, and for acquisition of Position Data from this, is not included in this export package. Please note that the subroutine PETRA is now enclosed in this export package. See section 4.3 and Appendix A4. 2. Description of tape EXSUP5 Tape EXSUP5 is not described, since the Position data file and subroutine listings are included in the Appendices. 3. The Position Data File. In this section the data layout for the position data and control data for this will be described. 3.1 The Control Record. The logical Control Record consists of two physical records each totaling 80 characters per line as shown below by :- Variable Format Description -------------------------------------------------------- IY I5 Year in which the tape was generated. MO I5 Ditto for month. DA D16.10 Ditto for day. NS I5 Size of Position Record in characters. NR I5 Number of logical Position Records. --------------------------------------------------------- ET1 D16.10 Earliest Epoch - MJD ( E.T ) ET2 D16.10 Latest Epoch - MJD ( E.T ) SUTC D16.10 MJD ( E.T ) at which UTC jumps. ETUTC1 D16.10 E.T - UTC before SUTC. ETUTC2 D16.10 E.T - UTC after SUTC. ----------------------------------------------------------- Note: SUTC and ETUTC2 are T.B.D at present. SUTC is preset to mjd 13230 ( after latest epoch on file ). ETUTC2 is preset to ETUTC1. 3.2 The Position Data Record. The logical Position Data Record consists of three physical records each totaling 80 characters. Each logical Position Record is shown below by -: Variable Format Description ------------------------------------------------------------------------ ET D16.10 Epoch = MJD ( E.T ) PG(1) D16.10 Component of position vector of the Giotto probe at ET. PG(2) D16.10 Component of position vector of the Giotto probe at ET. PG(3) D16.10 Component of position vector of the Giotto probe at ET. PH(1) D16.10 Component of position vector of P/Grigg- Skjellerup at ET. ------------------------------------------------------------------------ PH(2) D16.10 Component of position vector of P/Grigg- Skjellerup at ET. PH(3) D16.10 Component of position vector of P/Grigg- Skjellerup at ET. PE(1) D16.10 Component of position vector of the Earth at ET. PE(2) D16.10 Component of position vector of the Earth at ET. PE(3) D16.10 Component of position vector of the Earth at ET. ------------------------------------------------------------------------ PS(1) D16.10 Component of position vector of the Sun at ET. PS(2) D16.10 Component of position vector of the Sun at ET. PS(3) D16.10 Component of position vector of the Sun at ET. ------------------------------------------------------------------------ Note: Reference frame: M.E.S of 1950. Origin = solar system barycenter. Position units: km Time units: days The Position Data File is stored in Appendix A1. 4. Interface/Auxiliary Software. In this section the auxiliary software found on the tape will be described. 4.1 Subroutine GIODAT ( interface routine to the sequential formatted Position Data File). Input variable: TUTC : Universal coordinated time since 1950.0 in days. Output variables: PVG(1-3) : Position of the Giotto probe at time TUTC ( units: km ). PVG(4-6) : Approximate velocity of the Giotto probe at time TUTC (units: km/s) PVH(1-3) : PVH(4-6) : Ditto for Component of position vector of P/Grigg-Skjellerup. PVE(1-3) : PVE(4-6) : Ditto for Earth. PVS(1-3) : PVS(4-6) : Ditto for Sun. The positions will be acquired from the sequential file using quadratic interpolation. For the Position Records on the file a time step size was selected which will give an accuracy better than 1 sigma-error limit. The approximate velocities are derived from the position data. A Source Listing Of Subroutine Giodat Can Be Seen In Appendix A2 Subroutine TRAVEL ( generation of Experimenter Support Data ). Subroutine TRAVEL provides the signal travel time between a station of Earth and the Giotto S/C as well as position of station and S/C at the times of signal reception and transmission. The states of the Giotto S/C, the Earth and the sun are acquired from the formatted Position Data File on unit 7 ( described in section 3 ) via subroutine GIODAT (described in section 4.1 ). The station coordinates are computed without precession. Thus an accuracy of approx. 1 msec for the travel times computed can be expected. Input variables: TUTC : Signal reception time ( Universal Coordinated Time since 1950.0 in days ). S(3) : Cartesian Coordinates of station in the body fixed system at Epoch TUTC ( units: km ). Output variables: DTG : Signal travel time ( sec ) Giotto S/C - Station P(1-3): Position of station at TUTC in Heliocentric EME50 ( units: km ). G(1-3): Ditto for the Giotto S/C at TUTC-DTG. See Appendix A3 for a source listing of subroutine TRAVEL. 4.3 Subroutine PETRA ( computing of Kepler Orbits ). About one day prior to the encounter, position and velocity data for the relevant bodies as well as misc. data, was sent to the experimenters. Subroutine PETRA can be used to compute a Kepler orbit, thus obtaining position/velocities at the Encounter Epoch. See Appendix A4 for a source listing. 5. The Station Coordinates. In this section some general comments on station coordinates will be given as well as a table of coordinates for stations relevant to the Giotto mission. 5.1 General. Section 5.2 contains a table with the coordinates of some stations that are involved in the Giotto mission. These coordinates are needed in the computation of the signal travel times between the Giotto S/C and the station ( see also subroutine TRAVEL ). The coordinates are Cartesian position vectors ( X,Y,Z ) in the Body Fixed System. The X-Y plane of this system is the equatorial plane of the Earth with the X-axis going through the Greenwich Meridian. 5.2 Table of Station Coordinates. Station name X ( km ) Y ( km ) Z ( km ) Carnarvon -2328 5302 -2666 Tidbinbilla -4461 2682 -3675 Madrid 4847 -370 4117 Goldstone -2350 -4652 3666 ********************************************************************** The format is the same as was used for the P/Halley encounter. The data set consists of 1 header record and 280 logical data records. The total number of physical records is 842, each 80 bytes long. The data covers the period 1992 July 1.0 to July 14.0. The routines GIODAT, TRAVEL and PETRA which were supplied to the experimenters for P/Halley are still applicable. Use of interpolation routine GIODAT will result in small differences compared with the data on the original ESOC files. For P/G-S, the average (maximum) position difference over the interval from July 1.1 until July 14.0 is 8.2 km (16.4 km). The position difference at the encounter time is 12.2 km. For Giotto, the average (maximum) position difference over the interval from July 1.1 until July 14.0 is 6.6 km (13.1 km) except for around July 2.0. At this time, the orbit file was updated from a new orbit determination and there is a position discontinuity of 406.6 km. The interpolation smoothes this out so the maximum position difference between July 1.9 and July 2.1 is 239 km. The position difference at the encounter time is 9.8 km. Apart from around July 2.0 for Giotto, the differences detailed above are of course, far smaller than the uncertainties in the orbits of the two bodies. From the data, the experimenters will find that the time of closest approach is 10 July at 15:18:43.6 UT or l5:l9:42.8 ET, both on-board times, with a miss-distance of l4.7 km. Again, both the time and distance are subject to uncertainties as reflected by the error ellipses. As we all know, OPE estimates a miss-distance of circa 200 km. To test that they are using the orbital data correctly, the experimenters can compare their results with those given below. They result from interpolating the data using GIODAT at the encounter time and again precisely 10 minutes later, and are barycentric Bl950.0. EPOCH = 15531.638004815 (UTC) ---------------------------------------------------- GIOTTO POSITION/VELOCITY = PVG(l-6):- ---------------------------------------------------- -.1421241922D+09 -.3700354006D+08 -.3170794771D+08 0.1078429300D+02 -.2653479808D+02 -.1086507984D+02 COMET P/G-S POSITION/VELOCITY = PVH(1-6):- ---------------------------------------------------- -.1421241817D+09 -.3700353133D+08 -.3170794224D+08 0.1364733772D+02 -.3569294824D+02 -.6828218254D+00 EARTH POSITION/VELOCITY = PVE(1-6):- ---------------------------------------------------- 0.4738010968D+08 -.1322267079D+09 -.5734347600D+08 0.2783258405D+02 0.8367409578D+01 0.3626565039D+01 SUN POSITION/VELOCITY = PVS(1-6):- ---------------------------------------------------- 0.3662957289D+06 0.4704742983D+06 0.1908138597D+06 0.6631734562D-04 0.8841270648D-02 0.3848191013D-02 EPOCH = 15531.644949259 (UTC) ---------------------------------------------------- GIOTTO POSITION/VELOCITY = PVG(1-6):- ---------------------------------------------------- -.1421177160D+09 -.3701947210D+08 -.3171447121D+08 0.1078759287D+02 -.2653393328D+02 -.1086434143D+02 COMET P/G-S POSITION/VELOCITY = PVH(1-6):- ---------------------------------------------------- -.1421159865D+09 -.3702496220D+08 -.3170835200D+08 0.1365062171D+02 -.3569208831D+02 -.6820797674D+00 EARTH POSITION/VELOCITY = PVE(1-6):- ---------------------------------------------------- 0.4739682089D+08 -.1322216829D+09 -.5734129811D+08 0.2783151782D+02 0.8370415169D+01 0.3627863388D+01 SUN POSITION/VELOCITY = PVS(1-6):- ---------------------------------------------------- 0.3662957687D+06 0.4704796069D+06 0.1908161703D+06 0.6619033189D-04 0.8841297634D-02 0.3848198422D-02 The subroutine TRAVEL may be used to compute the one-way light time from the spacecraft to the ground station. The relevant ground station co-ordinates are given below. They are cartesian in the body-fixed system. The x-y plane is the equator and the x-axis passes through the Greenwich meridian. Ground Station x (km) y (km) z (km) Weilheim +4206 + 824 +4709 Madrid DSS 63 +4849 - 360 +4115 Canberra DSS 43 -4461 +2682 -3675 Goldstone DSS 14 -2354 -4641 +3677 ************************************************************************* Appendix A1 The Position Data File 1992 80.1800000000D+02 208 280 0.1552200000D+050.1553500000D+050.1552200000D+050.6734259260D-030.6850000000D-03 0.1552200000D+05-.1491694787D+09-.1450565992D+08-.2227386924D+08-.1515783709D+09 -.6917125671D+07-.3072484224D+080.2367921871D+08-.1374351530D+09-.5960159019D+08 0.3661752418D+060.4631192870D+060.1876144192D+06 0.1552204659D+05-.1491448067D+09-.1461585598D+08-.2232110997D+08-.1515417892D+09 -.7063649257D+07-.3073151132D+080.2379580926D+08-.1374185850D+09-.5959440280D+08 0.3661761370D+060.4631548013D+060.1876298592D+06 0.1552209319D+05-.1491200430D+09-.1472604279D+08-.2236833690D+08-.1515051202D+09 -.7210168511D+07-.3073816261D+080.2391238533D+08-.1374019320D+09-.5958717865D+08 0.3661770291D+060.4631903161D+060.1876452994D+06 0.1552213978D+05-.1490951877D+09-.1483622028D+08-.2241554999D+08-.1514683637D+09 -.7356683347D+07-.3074479610D+080.2402894682D+08-.1373851943D+09-.5957991774D+08 0.3661779182D+060.4632258313D+060.1876607398D+06 0.1552218638D+05-.1490702407D+09-.1494638839D+08-.2246274922D+08-.1514315198D+09 -.7503193679D+07-.3075141177D+080.2414549366D+08-.1373683718D+09-.5957262008D+08 0.3661788043D+060.4632613469D+060.1876761806D+06 0.1552223297D+05-.1490452020D+09-.1505654705D+08-.2250993456D+08-.1513945884D+09 -.7649699420D+07-.3075800962D+080.2426202577D+08-.1373514645D+09-.5956528567D+08 0.3661796873D+060.4632968630D+060.1876916216D+06 0.1552227957D+05-.1490200718D+09-.1516669619D+08-.2255710597D+08-.1513575695D+09 -.7796200483D+07-.3076458962D+080.2437854307D+08-.1373344724D+09-.5955791453D+08 0.3661805674D+060.4633323795D+060.1877070629D+06 0.1552232616D+05-.1489948498D+09-.1527683573D+08-.2260426343D+08-.1513204632D+09 -.7942696780D+07-.3077115177D+080.2449504549D+08-.1373173955D+09-.5955050666D+08 0.3661814444D+060.4633678964D+060.1877225045D+06 0.1552237276D+05-.1489695362D+09-.1538696562D+08-.2265140691D+08-.1512832692D+09 -.8089188225D+07-.3077769605D+080.2461153294D+08-.1373002339D+09-.5954306206D+08 0.3661823183D+060.4634034137D+060.1877379463D+06 0.1552241935D+05-.1489441310D+09-.1549708578D+08-.2269853636D+08-.1512459878D+09 -.8235674731D+07-.3078422245D+080.2472800535D+08-.1372829875D+09-.5953558075D+08 0.3661831893D+060.4634389315D+060.1877533884D+06 0.1552246595D+05-.1489186342D+09-.1560719615D+08-.2274565177D+08-.1512086187D+09 -.8382156210D+07-.3079073096D+080.2484446264D+08-.1372656564D+09-.5952806272D+08 0.3661840572D+060.4634744497D+060.1877688308D+06 0.1552251254D+05-.1488930457D+09-.1571729664D+08-.2279275311D+08-.1511711620D+09 -.8528632575D+07-.3079722156D+080.2496090474D+08-.1372482406D+09-.5952050798D+08 0.3661849221D+060.4635099683D+060.1877842735D+06 0.1552255914D+05-.1488673656D+09-.1582738721D+08-.2283984033D+08-.1511336177D+09 -.8675103738D+07-.3080369424D+080.2507733155D+08-.1372307401D+09-.5951291654D+08 0.3661857840D+060.4635454873D+060.1877997164D+06 0.1552260573D+05-.1488415938D+09-.1593746777D+08-.2288691341D+08-.1510959857D+09 -.8821569613D+07-.3081014898D+080.2519374301D+08-.1372131549D+09-.5950528841D+08 0.3661866428D+060.4635810068D+060.1878151596D+06 0.1552265233D+05-.1488157305D+09-.1604753826D+08-.2293397233D+08-.1510582660D+09 -.8968030112D+07-.3081658578D+080.2531013904D+08-.1371954851D+09-.5949762359D+08 0.3661874986D+060.4636165267D+060.1878306030D+06 0.1552269892D+05-.1487897755D+09-.1615759861D+08-.2298101705D+08-.1510204587D+09 -.9114485146D+07-.3082300462D+080.2542651956D+08-.1371777306D+09-.5948992209D+08 0.3661883514D+060.4636520470D+060.1878460468D+06 0.1552274552D+05-.1487637290D+09-.1626764874D+08-.2302804753D+08-.1509825636D+09 -.9260934629D+07-.3082940548D+080.2554288449D+08-.1371598914D+09-.5948218392D+08 0.3661892012D+060.4636875677D+060.1878614908D+06 0.1552279211D+05-.1487375908D+09-.1637768860D+08-.2307506375D+08-.1509445807D+09 -.9407378472D+07-.3083578836D+080.2565923376D+08-.1371419676D+09-.5947440908D+08 0.3661900479D+060.4637230888D+060.1878769350D+06 0.1552283871D+05-.1487113610D+09-.1648771812D+08-.2312206568D+08-.1509065101D+09 -.9553816588D+07-.3084215325D+080.2577556728D+08-.1371239591D+09-.5946659757D+08 0.3661908916D+060.4637586104D+060.1878923796D+06 0.1552288530D+05-.1486850397D+09-.1659773721D+08-.2316905329D+08-.1508683517D+09 -.9700248888D+07-.3084850012D+080.2589188498D+08-.1371058661D+09-.5945874941D+08 0.3661917323D+060.4637941324D+060.1879078244D+06 0.1552293190D+05-.1486586267D+09-.1670774582D+08-.2321602654D+08-.1508301054D+09 -.9846675285D+07-.3085482896D+080.2600818679D+08-.1370876885D+09-.5945086461D+08 0.3661925699D+060.4638296548D+060.1879232695D+06 0.1552297849D+05-.1486321222D+09-.1681774388D+08-.2326298541D+08-.1507917713D+09 -.9993095690D+07-.3086113977D+080.2612447261D+08-.1370694263D+09-.5944294316D+08 0.3661934045D+060.4638651777D+060.1879387148D+06 0.1552302509D+05-.1486054854D+09-.1692758951D+08-.2331030867D+08-.1507533494D+09 -.1013951002D+08-.3086743253D+080.2624074239D+08-.1370510795D+09-.5943498508D+08 0.3661942361D+060.4639007009D+060.1879541604D+06 0.1552307168D+05-.1485787976D+09-.1703756628D+08-.2335723851D+08-.1507148395D+09 -.1028591817D+08-.3087370722D+080.2635699603D+08-.1370326482D+09-.5942699037D+08 0.3661950647D+060.4639362246D+060.1879696063D+06 0.1552311828D+05-.1485520182D+09-.1714753230D+08-.2340415386D+08-.1506762418D+09 -.1043232007D+08-.3087996384D+080.2647323346D+08-.1370141323D+09-.5941895905D+08 0.3661958902D+060.4639717487D+060.1879850525D+06 0.1552316487D+05-.1485251473D+09-.1725748749D+08-.2345105472D+08-.1506375561D+09 -.1057871563D+08-.3088620237D+080.2658945461D+08-.1369955319D+09-.5941089111D+08 0.3661967127D+060.4640072732D+060.1880004989D+06 0.1552321147D+05-.1484981848D+09-.1736743177D+08-.2349794103D+08-.1505987824D+09 -.1072510475D+08-.3089242280D+080.2670565939D+08-.1369768471D+09-.5940278656D+08 0.3661975322D+060.4640427982D+060.1880159456D+06 0.1552325806D+05-.1484711307D+09-.1747736509D+08-.2354481278D+08-.1505599208D+09 -.1087148735D+08-.3089862511D+080.2682184774D+08-.1369580777D+09-.5939464542D+08 0.3661983486D+060.4640783235D+060.1880313925D+06 0.1552330466D+05-.1484439852D+09-.1758728737D+08-.2359166992D+08-.1505209711D+09 -.1101786334D+08-.3090480930D+080.2693801956D+08-.1369392239D+09-.5938646768D+08 0.3661991620D+060.4641138493D+060.1880468397D+06 0.1552335125D+05-.1484167480D+09-.1769719855D+08-.2363851244D+08-.1504819335D+09 -.1116423264D+08-.3091097535D+080.2705417479D+08-.1369202856D+09-.5937825336D+08 0.3661999724D+060.4641493755D+060.1880622872D+06 0.1552339785D+05-.1483894193D+09-.1780709855D+08-.2368534029D+08-.1504428078D+09 -.1131059514D+08-.3091712324D+080.2717031335D+08-.1369012628D+09-.5937000247D+08 0.3662007797D+060.4641849021D+060.1880777350D+06 0.1552344444D+05-.1483619991D+09-.1791698731D+08-.2373215346D+08-.1504035940D+09 -.1145695076D+08-.3092325297D+080.2728643517D+08-.1368821556D+09-.5936171500D+08 0.3662015840D+060.4642204291D+060.1880931830D+06 0.1552349104D+05-.1483344874D+09-.1802686475D+08-.2377895190D+08-.1503642922D+09 -.1160329942D+08-.3092936452D+080.2740254016D+08-.1368629641D+09-.5935339097D+08 0.3662023853D+060.4642559565D+060.1881086313D+06 0.1552353763D+05-.1483068841D+09-.1813673082D+08-.2382573559D+08-.1503249022D+09 -.1174964103D+08-.3093545788D+080.2751862824D+08-.1368436881D+09-.5934503039D+08 0.3662031836D+060.4642914844D+060.1881240799D+06 0.1552358423D+05-.1482791894D+09-.1824658543D+08-.2387250450D+08-.1502854241D+09 -.1189597549D+08-.3094153303D+080.2763469935D+08-.1368243277D+09-.5933663326D+08 0.3662039788D+060.4643270126D+060.1881395287D+06 0.1552363082D+05-.1482514031D+09-.1835642852D+08-.2391925860D+08-.1502458579D+09 -.1204230271D+08-.3094758997D+080.2775075340D+08-.1368048830D+09-.5932819959D+08 0.3662047710D+060.4643625413D+060.1881549778D+06 0.1552367742D+05-.1482235253D+09-.1846626002D+08-.2396599785D+08-.1502062035D+09 -.1218862262D+08-.3095362868D+080.2786679032D+08-.1367853539D+09-.5931972938D+08 0.3662055601D+060.4643980704D+060.1881704271D+06 0.1552372401D+05-.1481955560D+09-.1857607994D+08-.2401272221D+08-.1501664609D+09 -.1233493511D+08-.3095964915D+080.2798281003D+08-.1367657405D+09-.5931122266D+08 0.3662063463D+060.4644335999D+060.1881858767D+06 0.1552377061D+05-.1481674952D+09-.1868588820D+08-.2405943167D+08-.1501266301D+09 -.1248124010D+08-.3096565136D+080.2809881246D+08-.1367460428D+09-.5930267941D+08 0.3662071293D+060.4644691298D+060.1882013266D+06 0.1552381720D+05-.1481393429D+09-.1879568466D+08-.2410612619D+08-.1500867110D+09 -.1262753750D+08-.3097163531D+080.2821479753D+08-.1367262607D+09-.5929409966D+08 0.3662079094D+060.4645046602D+060.1882167768D+06 0.1552386380D+05-.1481110992D+09-.1890546926D+08-.2415280574D+08-.1500467037D+09 -.1277382721D+08-.3097760098D+080.2833076516D+08-.1367063944D+09-.5928548341D+08 0.3662086864D+060.4645401909D+060.1882322272D+06 0.1552391039D+05-.1480827639D+09-.1901524193D+08-.2419947029D+08-.1500066081D+09 -.1292010916D+08-.3098354836D+080.2844671528D+08-.1366864439D+09-.5927683066D+08 0.3662094604D+060.4645757221D+060.1882476779D+06 0.1552395699D+05-.1480543372D+09-.1912500259D+08-.2424611982D+08-.1499664243D+09 -.1306638325D+08-.3098947743D+080.2856264781D+08-.1366664090D+09-.5926814143D+08 0.3662102314D+060.4646112536D+060.1882631288D+06 0.1552400358D+05-.1480258191D+09-.1923475117D+08-.2429275429D+08-.1499261521D+09 -.1321264939D+08-.3099538819D+080.2867856267D+08-.1366462900D+09-.5925941572D+08 0.3662109993D+060.4646467856D+060.1882785800D+06 0.1552405018D+05-.1479972095D+09-.1934448762D+08-.2433937367D+08-.1498857916D+09 -.1335890748D+08-.3100128061D+080.2879445979D+08-.1366260867D+09-.5925065353D+08 0.3662117642D+060.4646823180D+060.1882940315D+06 0.1552409677D+05-.1479685084D+09-.1945421185D+08-.2438597793D+08-.1498453428D+09 -.1350515745D+08-.3100715470D+080.2891033910D+08-.1366057992D+09-.5924185489D+08 0.3662125261D+060.4647178508D+060.1883094832D+06 0.1552414337D+05-.1479397159D+09-.1956392381D+08-.2443256704D+08-.1498048055D+09 -.1365139920D+08-.3101301043D+080.2902620052D+08-.1365854276D+09-.5923301979D+08 0.3662132849D+060.4647533840D+060.1883249352D+06 0.1552418996D+05-.1479108320D+09-.1967362341D+08-.2447914097D+08-.1497641799D+09 -.1379763264D+08-.3101884780D+080.2914204396D+08-.1365649718D+09-.5922414825D+08 0.3662140407D+060.4647889176D+060.1883403875D+06 0.1552423656D+05-.1478818566D+09-.1978331060D+08-.2452569969D+08-.1497234659D+09 -.1394385768D+08-.3102466678D+080.2925786937D+08-.1365444318D+09-.5921524027D+08 0.3662147934D+060.4648244516D+060.1883558400D+06 0.1552428315D+05-.1478527898D+09-.1989298529D+08-.2457224316D+08-.1496826635D+09 -.1409007422D+08-.3103046738D+080.2937367665D+08-.1365238077D+09-.5920629586D+08 0.3662155432D+060.4648599861D+060.1883712928D+06 0.1552432975D+05-.1478236316D+09-.2000264743D+08-.2461877137D+08-.1496417726D+09 -.1423628219D+08-.3103624957D+080.2948946574D+08-.1365030996D+09-.5919731503D+08 0.3662162899D+060.4648955209D+060.1883867459D+06 0.1552437634D+05-.1477943820D+09-.2011229695D+08-.2466528427D+08-.1496007932D+09 -.1438248149D+08-.3104201335D+080.2960523657D+08-.1364823073D+09-.5918829778D+08 0.3662170335D+060.4649310562D+060.1884021992D+06 0.1552442294D+05-.1477650410D+09-.2022193376D+08-.2471178184D+08-.1495597254D+09 -.1452867202D+08-.3104775870D+080.2972098904D+08-.1364614309D+09-.5917924413D+08 0.3662177741D+060.4649665918D+060.1884176528D+06 0.1552446953D+05-.1477356086D+09-.2033155782D+08-.2475826404D+08-.1495185690D+09 -.1467485370D+08-.3105348560D+080.2983672310D+08-.1364404705D+09-.5917015408D+08 0.3662185117D+060.4650021279D+060.1884331066D+06 0.1552451613D+05-.1477060849D+09-.2044116904D+08-.2480473085D+08-.1494773242D+09 -.1482102644D+08-.3105919406D+080.2995243867D+08-.1364194260D+09-.5916102764D+08 0.3662192463D+060.4650376643D+060.1884485607D+06 0.1552456272D+05-.1476764697D+09-.2055076735D+08-.2485118223D+08-.1494359908D+09 -.1496719015D+08-.3106488405D+080.3006813566D+08-.1363982976D+09-.5915186483D+08 0.3662199778D+060.4650732012D+060.1884640151D+06 0.1552460932D+05-.1476467632D+09-.2066035270D+08-.2489761816D+08-.1493945688D+09 -.1511334473D+08-.3107055556D+080.3018381401D+08-.1363770851D+09-.5914266564D+08 0.3662207063D+060.4651087385D+060.1884794697D+06 0.1552465591D+05-.1476169653D+09-.2076992500D+08-.2494403860D+08-.1493530583D+09 -.1525949010D+08-.3107620859D+080.3029947365D+08-.1363557886D+09-.5913343008D+08 0.3662214317D+060.4651442762D+060.1884949246D+06 0.1552470251D+05-.1475870760D+09-.2087948419D+08-.2499044353D+08-.1493114592D+09 -.1540562616D+08-.3108184311D+080.3041511449D+08-.1363344082D+09-.5912415817D+08 0.3662221541D+060.4651798143D+060.1885103798D+06 0.1552474910D+05-.1475570954D+09-.2098903020D+08-.2503683291D+08-.1492697714D+09 -.1555175283D+08-.3108745912D+080.3053073646D+08-.1363129438D+09-.5911484992D+08 0.3662228735D+060.4652153527D+060.1885258352D+06 0.1552479570D+05-.1475270235D+09-.2109856296D+08-.2508320671D+08-.1492279951D+09 -.1569787000D+08-.3109305661D+080.3064633949D+08-.1362913955D+09-.5910550533D+08 0.3662235898D+060.4652508916D+060.1885412908D+06 0.1552484229D+05-.1474968602D+09-.2120808241D+08-.2512956490D+08-.1491861301D+09 -.1584397760D+08-.3109863555D+080.3076192351D+08-.1362697633D+09-.5909612440D+08 0.3662243031D+060.4652864309D+060.1885567468D+06 0.1552488889D+05-.1474666056D+09-.2131758846D+08-.2517590745D+08-.1491441764D+09 -.1599007553D+08-.3110419595D+080.3087748843D+08-.1362480472D+09-.5908670716D+08 0.3662250134D+060.4653219706D+060.1885722030D+06 0.1552493548D+05-.1474362597D+09-.2142708106D+08-.2522223433D+08-.1491021341D+09 -.1613616370D+08-.3110973779D+080.3099303419D+08-.1362262472D+09-.5907725360D+08 0.3662257206D+060.4653575107D+060.1885876594D+06 0.1552498208D+05-.1474058225D+09-.2153656013D+08-.2526854552D+08-.1490600031D+09 -.1628224201D+08-.3111526105D+080.3110856071D+08-.1362043634D+09-.5906776374D+08 0.3662264248D+060.4653930512D+060.1886031162D+06 0.1552502867D+05-.1473752939D+09-.2164602561D+08-.2531484097D+08-.1490177833D+09 -.1642831038D+08-.3112076573D+080.3122406793D+08-.1361823958D+09-.5905823758D+08 0.3662271259D+060.4654285921D+060.1886185731D+06 0.1552507527D+05-.1473446741D+09-.2175547742D+08-.2536112067D+08-.1489754749D+09 -.1657436872D+08-.3112625181D+080.3133955575D+08-.1361603443D+09-.5904867513D+08 0.3662278241D+060.4654641334D+060.1886340304D+06 0.1552512186D+05-.1473139630D+09-.2186491550D+08-.2540738457D+08-.1489330777D+09 -.1672041693D+08-.3113171927D+080.3145502411D+08-.1361382090D+09-.5903907641D+08 0.3662285191D+060.4654996751D+060.1886494879D+06 0.1552516846D+05-.1472831606D+09-.2197433977D+08-.2545363265D+08-.1488905917D+09 -.1686645492D+08-.3113716812D+080.3157047295D+08-.1361159900D+09-.5902944142D+08 0.3662292112D+060.4655352172D+060.1886649456D+06 0.1552521505D+05-.1472522670D+09-.2208375017D+08-.2549986487D+08-.1488480170D+09 -.1701248260D+08-.3114259834D+080.3168590217D+08-.1360936872D+09-.5901977016D+08 0.3662299002D+060.4655707597D+060.1886804037D+06 0.1552526165D+05-.1472212820D+09-.2219314663D+08-.2554608122D+08-.1488053535D+09 -.1715849988D+08-.3114800991D+080.3180131172D+08-.1360713007D+09-.5901006266D+08 0.3662305861D+060.4656063026D+060.1886958619D+06 0.1552530824D+05-.1471902059D+09-.2230252908D+08-.2559228165D+08-.1487626011D+09 -.1730450666D+08-.3115340282D+080.3191670151D+08-.1360488304D+09-.5900031891D+08 0.3662312690D+060.4656418459D+060.1887113205D+06 0.1552535484D+05-.1471590384D+09-.2241189744D+08-.2563846614D+08-.1487197600D+09 -.1745050286D+08-.3115877706D+080.3203207148D+08-.1360262765D+09-.5899053892D+08 0.3662319489D+060.4656773896D+060.1887267793D+06 0.1552540143D+05-.1471277798D+09-.2252125166D+08-.2568463465D+08-.1486768300D+09 -.1759648838D+08-.3116413262D+080.3214742154D+08-.1360036389D+09-.5898072271D+08 0.3662326258D+060.4657129337D+060.1887422383D+06 0.1552544803D+05-.1470964299D+09-.2263059166D+08-.2573078716D+08-.1486338111D+09 -.1774246313D+08-.3116946949D+080.3226275164D+08-.1359809176D+09-.5897087029D+08 0.3662332996D+060.4657484782D+060.1887576976D+06 0.1552549462D+05-.1470649888D+09-.2273991736D+08-.2577692363D+08-.1485907034D+09 -.1788842703D+08-.3117478765D+080.3237806169D+08-.1359581127D+09-.5896098165D+08 0.3662339703D+060.4657840230D+060.1887731572D+06 0.1552554122D+05-.1470334565D+09-.2284922871D+08-.2582304404D+08-.1485475068D+09 -.1803437996D+08-.3118008710D+080.3249335162D+08-.1359352241D+09-.5895105682D+08 0.3662346381D+060.4658195683D+060.1887886170D+06 0.1552558781D+05-.1470018329D+09-.2295852563D+08-.2586914835D+08-.1485042213D+09 -.1818032186D+08-.3118536782D+080.3260862136D+08-.1359122520D+09-.5894109579D+08 0.3662353028D+060.4658551140D+060.1888040771D+06 0.1552563441D+05-.1469701182D+09-.2306780806D+08-.2591523654D+08-.1484608469D+09 -.1832625261D+08-.3119062980D+080.3272387084D+08-.1358891963D+09-.5893109859D+08 0.3662359644D+060.4658906601D+060.1888195375D+06 0.1552568100D+05-.1469383123D+09-.2317707592D+08-.2596130856D+08-.1484173836D+09 -.1847217214D+08-.3119587303D+080.3283909998D+08-.1358660570D+09-.5892106521D+08 0.3662366230D+060.4659262065D+060.1888349981D+06 0.1552572760D+05-.1469064152D+09-.2328632914D+08-.2600736441D+08-.1483738313D+09 -.1861808034D+08-.3120109750D+080.3295430872D+08-.1358428342D+09-.5891099567D+08 0.3662372786D+060.4659617534D+060.1888504589D+06 0.1552577419D+05-.1468744270D+09-.2339556765D+08-.2605340403D+08-.1483301901D+09 -.1876397713D+08-.3120630319D+080.3306949698D+08-.1358195279D+09-.5890088997D+08 0.3662379311D+060.4659973006D+060.1888659200D+06 0.1552582079D+05-.1468423476D+09-.2350479140D+08-.2609942740D+08-.1482864599D+09 -.1890986241D+08-.3121149010D+080.3318466468D+08-.1357961381D+09-.5889074813D+08 0.3662385806D+060.4660328483D+060.1888813814D+06 0.1552586738D+05-.1468101770D+09-.2361400030D+08-.2614543450D+08-.1482426407D+09 -.1905573609D+08-.3121665821D+080.3329981176D+08-.1357726648D+09-.5888057015D+08 0.3662392271D+060.4660683963D+060.1888968431D+06 0.1552591398D+05-.1467779153D+09-.2372319428D+08-.2619142529D+08-.1481987326D+09 -.1920159808D+08-.3122180751D+080.3341493814D+08-.1357491080D+09-.5887035604D+08 0.3662398705D+060.4661039447D+060.1889123049D+06 0.1552596057D+05-.1467455625D+09-.2383237341D+08-.2623739972D+08-.1481547354D+09 -.1934744828D+08-.3122693799D+080.3353004376D+08-.1357254678D+09-.5886010581D+08 0.3662405108D+060.4661394935D+060.1889277671D+06 0.1552600717D+05-.1467131185D+09-.2394153753D+08-.2628335778D+08-.1481106492D+09 -.1949328661D+08-.3123204963D+080.3364512853D+08-.1357017442D+09-.5884981948D+08 0.3662411482D+060.4661750428D+060.1889432295D+06 0.1552605376D+05-.1466805834D+09-.2405068654D+08-.2632929944D+08-.1480664739D+09 -.1963911296D+08-.3123714244D+080.3376019239D+08-.1356779372D+09-.5883949704D+08 0.3662417824D+060.4662105924D+060.1889586921D+06 0.1552610036D+05-.1466479572D+09-.2415982036D+08-.2637522467D+08-.1480222096D+09 -.1978492725D+08-.3124221638D+080.3387523527D+08-.1356540469D+09-.5882913851D+08 0.3662424137D+060.4662461423D+060.1889741551D+06 0.1552614695D+05-.1466152400D+09-.2426893891D+08-.2642113344D+08-.1479778563D+09 -.1993072939D+08-.3124727147D+080.3399025710D+08-.1356300731D+09-.5881874390D+08 0.3662430419D+060.4662816927D+060.1889896182D+06 0.1552619355D+05-.1465824316D+09-.2437804214D+08-.2646702571D+08-.1479334139D+09 -.2007651928D+08-.3125230767D+080.3410525780D+08-.1356060161D+09-.5880831322D+08 0.3662436670D+060.4663172435D+060.1890050817D+06 0.1552624014D+05-.1465495322D+09-.2448712997D+08-.2651290147D+08-.1478888823D+09 -.2022229682D+08-.3125732498D+080.3422023730D+08-.1355818757D+09-.5879784647D+08 0.3662442892D+060.4663527947D+060.1890205453D+06 0.1552628674D+05-.1465165417D+09-.2459620233D+08-.2655876067D+08-.1478442617D+09 -.2036806193D+08-.3126232340D+080.3433519554D+08-.1355576521D+09-.5878734367D+08 0.3662449082D+060.4663883462D+060.1890360093D+06 0.1552633333D+05-.1464834601D+09-.2470525915D+08-.2660460329D+08-.1477995520D+09 -.2051381451D+08-.3126730290D+080.3445013244D+08-.1355333451D+09-.5877680482D+08 0.3662455243D+060.4664238981D+060.1890514735D+06 0.1552637993D+05-.1464502875D+09-.2481430036D+08-.2665042929D+08-.1477547531D+09 -.2065955447D+08-.3127226347D+080.3456504792D+08-.1355089550D+09-.5876622994D+08 0.3662461372D+060.4664594505D+060.1890669379D+06 0.1552642652D+05-.1464170239D+09-.2492332589D+08-.2669623866D+08-.1477098651D+09 -.2080528171D+08-.3127720512D+080.3467994193D+08-.1354844815D+09-.5875561902D+08 0.3662467472D+060.4664950032D+060.1890824026D+06 0.1552647312D+05-.1463836692D+09-.2503233568D+08-.2674203134D+08-.1476648880D+09 -.2095099615D+08-.3128212781D+080.3479481438D+08-.1354599249D+09-.5874497209D+08 0.3662473541D+060.4665305563D+060.1890978676D+06 0.1552651971D+05-.1463502236D+09-.2514132966D+08-.2678780733D+08-.1476198217D+09 -.2109669768D+08-.3128703155D+080.3490966521D+08-.1354352851D+09-.5873428915D+08 0.3662479580D+060.4665661098D+060.1891133328D+06 0.1552656631D+05-.1463166869D+09-.2525030774D+08-.2683356658D+08-.1475746662D+09 -.2124238622D+08-.3129191632D+080.3502449435D+08-.1354105622D+09-.5872357021D+08 0.3662485588D+060.4666016636D+060.1891287982D+06 0.1552661290D+05-.1462830592D+09-.2535926988D+08-.2687930907D+08-.1475294215D+09 -.2138806167D+08-.3129678212D+080.3513930173D+08-.1353857561D+09-.5871281528D+08 0.3662491565D+060.4666372179D+060.1891442639D+06 0.1552665950D+05-.1462493406D+09-.2546821599D+08-.2692503476D+08-.1474840877D+09 -.2153372394D+08-.3130162892D+080.3525408727D+08-.1353608668D+09-.5870202436D+08 0.3662497513D+060.4666727725D+060.1891597299D+06 0.1552670609D+05-.1462155309D+09-.2557714600D+08-.2697074362D+08-.1474386646D+09 -.2167937293D+08-.3130645672D+080.3536885091D+08-.1353358945D+09-.5869119748D+08 0.3662503430D+060.4667083275D+060.1891751961D+06 0.1552675269D+05-.1461816304D+09-.2568605985D+08-.2701643563D+08-.1473931523D+09 -.2182500856D+08-.3131126550D+080.3548359257D+08-.1353108391D+09-.5868033463D+08 0.3662509316D+060.4667438829D+060.1891906626D+06 0.1552679928D+05-.1461476388D+09-.2579495747D+08-.2706211076D+08-.1473475508D+09 -.2197063072D+08-.3131605527D+080.3559831219D+08-.1352857006D+09-.5866943582D+08 0.3662515172D+060.4667794387D+060.1892061293D+06 0.1552684588D+05-.1461135563D+09-.2590383878D+08-.2710776896D+08-.1473018600D+09 -.2211623932D+08-.3132082599D+080.3571300969D+08-.1352604791D+09-.5865850108D+08 0.3662520998D+060.4668149949D+060.1892215963D+06 0.1552689247D+05-.1460793829D+09-.2601270373D+08-.2715341022D+08-.1472560800D+09 -.2226183427D+08-.3132557768D+080.3582768501D+08-.1352351746D+09-.5864753039D+08 0.3662526793D+060.4668505515D+060.1892370635D+06 0.1552693907D+05-.1460451186D+09-.2612155223D+08-.2719903451D+08-.1472102108D+09 -.2240741548D+08-.3133031030D+080.3594233808D+08-.1352097871D+09-.5863652378D+08 0.3662532557D+060.4668861084D+060.1892525310D+06 0.1552698566D+05-.1460107634D+09-.2623038422D+08-.2724464179D+08-.1471642522D+09 -.2255298285D+08-.3133502386D+080.3605696882D+08-.1351843166D+09-.5862548126D+08 0.3662538292D+060.4669216657D+060.1892679988D+06 0.1552703226D+05-.1459763172D+09-.2633919963D+08-.2729023202D+08-.1471182044D+09 -.2269853629D+08-.3133971833D+080.3617157717D+08-.1351587632D+09-.5861440282D+08 0.3662543995D+060.4669572234D+060.1892834668D+06 0.1552707885D+05-.1459417802D+09-.2644799840D+08-.2733580520D+08-.1470720673D+09 -.2284407570D+08-.3134439372D+080.3628616305D+08-.1351331268D+09-.5860328849D+08 0.3662549669D+060.4669927815D+060.1892989350D+06 0.1552712545D+05-.1459071523D+09-.2655678044D+08-.2738136127D+08-.1470258409D+09 -.2298960099D+08-.3134905001D+080.3640072641D+08-.1351074075D+09-.5859213826D+08 0.3662555311D+060.4670283400D+060.1893144035D+06 0.1552717204D+05-.1458724335D+09-.2666554569D+08-.2742690021D+08-.1469795252D+09 -.2313511207D+08-.3135368718D+080.3651526716D+08-.1350816054D+09-.5858095216D+08 0.3662560924D+060.4670638988D+060.1893298722D+06 0.1552721864D+05-.1458376239D+09-.2677429409D+08-.2747242200D+08-.1469331202D+09 -.2328060884D+08-.3135830524D+080.3662978524D+08-.1350557203D+09-.5856973019D+08 0.3662566506D+060.4670994580D+060.1893453412D+06 0.1552726523D+05-.1458027235D+09-.2688302555D+08-.2751792660D+08-.1468866258D+09 -.2342609120D+08-.3136290416D+080.3674428058D+08-.1350297525D+09-.5855847235D+08 0.3662572057D+060.4671350176D+060.1893608105D+06 0.1552731183D+05-.1457677322D+09-.2699174002D+08-.2756341398D+08-.1468400421D+09 -.2357155907D+08-.3136748393D+080.3685875311D+08-.1350037018D+09-.5854717867D+08 0.3662577578D+060.4671705776D+060.1893762800D+06 0.1552735842D+05-.1457326500D+09-.2710043742D+08-.2760888411D+08-.1467933691D+09 -.2371701234D+08-.3137204455D+080.3697320276D+08-.1349775683D+09-.5853584914D+08 0.3662583069D+060.4672061380D+060.1893917497D+06 0.1552740502D+05-.1456974771D+09-.2720911769D+08-.2765433696D+08-.1467466067D+09 -.2386245093D+08-.3137658601D+080.3708762947D+08-.1349513520D+09-.5852448378D+08 0.3662588529D+060.4672416987D+060.1894072197D+06 0.1552745161D+05-.1456622134D+09-.2731778075D+08-.2769977249D+08-.1466997549D+09 -.2400787474D+08-.3138110829D+080.3720203316D+08-.1349250530D+09-.5851308259D+08 0.3662593959D+060.4672772598D+060.1894226900D+06 0.1552749821D+05-.1456268589D+09-.2742642653D+08-.2774519069D+08-.1466528138D+09 -.2415328367D+08-.3138561138D+080.3731641376D+08-.1348986712D+09-.5850164559D+08 0.3662599358D+060.4673128213D+060.1894381605D+06 0.1552754480D+05-.1455914136D+09-.2753505497D+08-.2779059152D+08-.1466057832D+09 -.2429867764D+08-.3139009527D+080.3743077121D+08-.1348722067D+09-.5849017278D+08 0.3662604726D+060.4673483831D+060.1894536312D+06 0.1552759140D+05-.1455558776D+09-.2764366600D+08-.2783597495D+08-.1465586633D+09 -.2444405654D+08-.3139455996D+080.3754510544D+08-.1348456596D+09-.5847866418D+08 0.3662610065D+060.4673839454D+060.1894691022D+06 0.1552763799D+05-.1455202508D+09-.2775225954D+08-.2788134095D+08-.1465114540D+09 -.2458942028D+08-.3139900543D+080.3765941638D+08-.1348190297D+09-.5846711979D+08 0.3662615372D+060.4674195080D+060.1894845735D+06 0.1552768459D+05-.1454845333D+09-.2786083552D+08-.2792668949D+08-.1464641553D+09 -.2473476876D+08-.3140343167D+080.3777370395D+08-.1347923172D+09-.5845553962D+08 0.3662620650D+060.4674550710D+060.1895000450D+06 0.1552773118D+05-.1454487250D+09-.2796939398D+08-.2797202052D+08-.1464167672D+09 -.2488010190D+08-.3140783866D+080.3788796810D+08-.1347655221D+09-.5844392368D+08 0.3662625896D+060.4674906343D+060.1895155167D+06 0.1552777778D+05-.1454128260D+09-.2807793479D+08-.2801733403D+08-.1463692897D+09 -.2502541959D+08-.3141222641D+080.3800220876D+08-.1347386444D+09-.5843227199D+08 0.3662631113D+060.4675261981D+060.1895309887D+06 0.1552782437D+05-.1453768363D+09-.2818645784D+08-.2806262999D+08-.1463217227D+09 -.2517072175D+08-.3141659490D+080.3811642585D+08-.1347116841D+09-.5842058454D+08 0.3662636299D+060.4675617622D+060.1895464610D+06 0.1552787097D+05-.1453407559D+09-.2829496305D+08-.2810790836D+08-.1462740663D+09 -.2531600828D+08-.3142094411D+080.3823061931D+08-.1346846413D+09-.5840886136D+08 0.3662641454D+060.4675973267D+060.1895619335D+06 0.1552791756D+05-.1453045848D+09-.2840345037D+08-.2815316913D+08-.1462263205D+09 -.2546127907D+08-.3142527404D+080.3834478907D+08-.1346575159D+09-.5839710244D+08 0.3662646579D+060.4676328915D+060.1895774062D+06 0.1552796416D+05-.1452683231D+09-.2851191971D+08-.2819841225D+08-.1461784852D+09 -.2560653404D+08-.3142958469D+080.3845893506D+08-.1346303080D+09-.5838530781D+08 0.3662651673D+060.4676684567D+060.1895928792D+06 0.1552801075D+05-.1452319707D+09-.2862037101D+08-.2824363770D+08-.1461305605D+09 -.2575177310D+08-.3143387602D+080.3857305721D+08-.1346030176D+09-.5837347746D+08 0.3662656737D+060.4677040223D+060.1896083524D+06 0.1552805735D+05-.1451955276D+09-.2872880420D+08-.2828884544D+08-.1460825463D+09 -.2589699614D+08-.3143814805D+080.3868715547D+08-.1345756447D+09-.5836161141D+08 0.3662661770D+060.4677395883D+060.1896238259D+06 0.1552810394D+05-.1451589940D+09-.2883721921D+08-.2833403546D+08-.1460344426D+09 -.2604220307D+08-.3144240075D+080.3880122975D+08-.1345481893D+09-.5834970966D+08 0.3662666773D+060.4677751546D+060.1896392996D+06 0.1552815054D+05-.1451223697D+09-.2894561597D+08-.2837920771D+08-.1459862494D+09 -.2618739380D+08-.3144663412D+080.3891527999D+08-.1345206516D+09-.5833777223D+08 0.3662671746D+060.4678107214D+060.1896547736D+06 0.1552819713D+05-.1450856548D+09-.2905399442D+08-.2842436217D+08-.1459379668D+09 -.2633256824D+08-.3145084814D+080.3902930612D+08-.1344930314D+09-.5832579913D+08 0.3662676688D+060.4678462884D+060.1896702478D+06 0.1552824373D+05-.1450488493D+09-.2916235447D+08-.2846949880D+08-.1458895947D+09 -.2647772627D+08-.3145504281D+080.3914330808D+08-.1344653289D+09-.5831379036D+08 0.3662681599D+060.4678818559D+060.1896857223D+06 0.1552829032D+05-.1450119533D+09-.2927069606D+08-.2851461758D+08-.1458411331D+09 -.2662286782D+08-.3145921812D+080.3925728580D+08-.1344375440D+09-.5830174593D+08 0.3662686480D+060.4679174237D+060.1897011970D+06 0.1552833692D+05-.1449749667D+09-.2937901913D+08-.2855971848D+08-.1457925820D+09 -.2676799279D+08-.3146337405D+080.3937123921D+08-.1344096767D+09-.5828966586D+08 0.3662691331D+060.4679529919D+060.1897166720D+06 0.1552838351D+05-.1449378895D+09-.2948732360D+08-.2860480147D+08-.1457439414D+09 -.2691310107D+08-.3146751060D+080.3948516824D+08-.1343817272D+09-.5827755015D+08 0.3662696151D+060.4679885604D+060.1897321472D+06 0.1552843011D+05-.1449007218D+09-.2959560941D+08-.2864986651D+08-.1456952113D+09 -.2705819258D+08-.3147162775D+080.3959907283D+08-.1343536953D+09-.5826539881D+08 0.3662700940D+060.4680241293D+060.1897476227D+06 0.1552847670D+05-.1448634635D+09-.2970387647D+08-.2869491358D+08-.1456463917D+09 -.2720326722D+08-.3147572550D+080.3971295291D+08-.1343255812D+09-.5825321185D+08 0.3662705699D+060.4680596986D+060.1897630984D+06 0.1552852330D+05-.1448261148D+09-.2981212473D+08-.2873994265D+08-.1455974826D+09 -.2734832489D+08-.3147980383D+080.3982680841D+08-.1342973848D+09-.5824098929D+08 0.3662710427D+060.4680952683D+060.1897785743D+06 0.1552856989D+05-.1447886755D+09-.2992035412D+08-.2878495368D+08-.1455484839D+09 -.2749336550D+08-.3148386274D+080.3994063927D+08-.1342691062D+09-.5822873112D+08 0.3662715125D+060.4681308383D+060.1897940505D+06 0.1552861649D+05-.1447511458D+09-.3002856456D+08-.2882994665D+08-.1454993957D+09 -.2763838895D+08-.3148790221D+080.4005444541D+08-.1342407454D+09-.5821643736D+08 0.3662719793D+060.4681664087D+060.1898095270D+06 0.1552866308D+05-.1447135256D+09-.3013675598D+08-.2887492152D+08-.1454502180D+09 -.2778339515D+08-.3149192223D+080.4016822678D+08-.1342123024D+09-.5820410802D+08 0.3662724430D+060.4682019794D+060.1898250037D+06 0.1552870968D+05-.1446758149D+09-.3024492832D+08-.2891987827D+08-.1454009508D+09 -.2792838400D+08-.3149592281D+080.4028198330D+08-.1341837772D+09-.5819174311D+08 0.3662729036D+060.4682375505D+060.1898404806D+06 0.1552875627D+05-.1446380138D+09-.3035308151D+08-.2896481687D+08-.1453515940D+09 -.2807335540D+08-.3149990391D+080.4039571491D+08-.1341551699D+09-.5817934264D+08 0.3662733612D+060.4682731220D+060.1898559578D+06 0.1552880287D+05-.1446001223D+09-.3046121547D+08-.2900973728D+08-.1453021477D+09 -.2821830927D+08-.3150386554D+080.4050942154D+08-.1341264805D+09-.5816690661D+08 0.3662738157D+060.4683086939D+060.1898714352D+06 0.1552884946D+05-.1445621403D+09-.3056933014D+08-.2905463948D+08-.1452526118D+09 -.2836324550D+08-.3150780769D+080.4062310313D+08-.1340977089D+09-.5815443504D+08 0.3662742672D+060.4683442661D+060.1898869129D+06 0.1552889606D+05-.1445240679D+09-.3067742544D+08-.2909952343D+08-.1452029863D+09 -.2850816399D+08-.3151173034D+080.4073675960D+08-.1340688553D+09-.5814192793D+08 0.3662747157D+060.4683798386D+060.1899023908D+06 0.1552894265D+05-.1444859051D+09-.3078550131D+08-.2914438911D+08-.1451532714D+09 -.2865306466D+08-.3151563349D+080.4085039090D+08-.1340399197D+09-.5812938530D+08 0.3662751610D+060.4684154116D+060.1899178689D+06 0.1552898925D+05-.1444476520D+09-.3089355768D+08-.2918923649D+08-.1451034668D+09 -.2879794741D+08-.3151951713D+080.4096399695D+08-.1340109020D+09-.5811680715D+08 0.3662756034D+060.4684509848D+060.1899333473D+06 0.1552903584D+05-.1444093085D+09-.3100159448D+08-.2923406553D+08-.1450535727D+09 -.2894281213D+08-.3152338124D+080.4107757769D+08-.1339818023D+09-.5810419350D+08 0.3662760427D+060.4684865585D+060.1899488260D+06 0.1552908244D+05-.1443708746D+09-.3110961164D+08-.2927887621D+08-.1450035891D+09 -.2908765874D+08-.3152722582D+080.4119113305D+08-.1339526207D+09-.5809154434D+08 0.3662764789D+060.4685221325D+060.1899643049D+06 0.1552912903D+05-.1443323504D+09-.3121760909D+08-.2932366849D+08-.1449535158D+09 -.2923248713D+08-.3153105085D+080.4130466297D+08-.1339233570D+09-.5807885970D+08 0.3662769121D+060.4685577069D+060.1899797840D+06 0.1552917563D+05-.1442937359D+09-.3132558675D+08-.2936844235D+08-.1449033530D+09 -.2937729722D+08-.3153485633D+080.4141816738D+08-.1338940115D+09-.5806613958D+08 0.3662773422D+060.4685932816D+060.1899952634D+06 0.1552922222D+05-.1442550311D+09-.3143354457D+08-.2941319776D+08-.1448531007D+09 -.2952208890D+08-.3153864225D+080.4153164622D+08-.1338645840D+09-.5805338399D+08 0.3662777692D+060.4686288567D+060.1900107430D+06 0.1552926882D+05-.1442162360D+09-.3154148246D+08-.2945793468D+08-.1448027587D+09 -.2966686209D+08-.3154240859D+080.4164509941D+08-.1338350746D+09-.5804059294D+08 0.3662781933D+060.4686644321D+060.1900262228D+06 0.1552931541D+05-.1441773506D+09-.3164940037D+08-.2950265308D+08-.1447523272D+09 -.2981161667D+08-.3154615536D+080.4175852689D+08-.1338054833D+09-.5802776644D+08 0.3662786142D+060.4687000079D+060.1900417029D+06 0.1552936201D+05-.1441383749D+09-.3175729822D+08-.2954735295D+08-.1447018061D+09 -.2995635257D+08-.3154988253D+080.4187192861D+08-.1337758102D+09-.5801490450D+08 0.3662790321D+060.4687355841D+060.1900571833D+06 0.1552940860D+05-.1440993090D+09-.3186517593D+08-.2959203424D+08-.1446511954D+09 -.3010106967D+08-.3155359010D+080.4198530448D+08-.1337460553D+09-.5800200713D+08 0.3662794470D+060.4687711606D+060.1900726639D+06 0.1552945520D+05-.1440601528D+09-.3197303345D+08-.2963669692D+08-.1446004952D+09 -.3024576789D+08-.3155727806D+080.4209865444D+08-.1337162186D+09-.5798907433D+08 0.3662798588D+060.4688067375D+060.1900881447D+06 0.1552950179D+05-.1440209065D+09-.3208087071D+08-.2968134097D+08-.1445497053D+09 -.3039044713D+08-.3156094640D+080.4221197844D+08-.1336863000D+09-.5797610612D+08 0.3662802675D+060.4688423147D+060.1901036258D+06 0.1552954839D+05-.1439815699D+09-.3218868762D+08-.2972596636D+08-.1444988259D+09 -.3053510729D+08-.3156459512D+080.4232527639D+08-.1336562998D+09-.5796310251D+08 0.3662806732D+060.4688778923D+060.1901191071D+06 0.1552959498D+05-.1439421431D+09-.3229648413D+08-.2977057305D+08-.1444478569D+09 -.3067974827D+08-.3156822419D+080.4243854825D+08-.1336262177D+09-.5795006350D+08 0.3662810759D+060.4689134703D+060.1901345886D+06 0.1552964158D+05-.1439026261D+09-.3240426016D+08-.2981516102D+08-.1443967983D+09 -.3082436998D+08-.3157183362D+080.4255179393D+08-.1335960540D+09-.5793698910D+08 0.3662814755D+060.4689490486D+060.1901500704D+06 0.1552968817D+05-.1438630189D+09-.3251201572D+08-.2985973004D+08-.1443456501D+09 -.3096897233D+08-.3157542339D+080.4266501339D+08-.1335658085D+09-.5792387933D+08 0.3662818720D+060.4689846272D+060.1901655524D+06 0.1552973477D+05-.1438233186D+09-.3261975268D+08-.2990427510D+08-.1442944124D+09 -.3111355521D+08-.3157899349D+080.4277820654D+08-.1335354814D+09-.5791073419D+08 0.3662822655D+060.4690202062D+060.1901810347D+06 0.1552978136D+05-.1437835274D+09-.3272746934D+08-.2994879994D+08-.1442430850D+09 -.3125811853D+08-.3158254391D+080.4289137333D+08-.1335050727D+09-.5789755369D+08 0.3662826559D+060.4690557856D+060.1901965172D+06 0.1552982796D+05-.1437436460D+09-.3283516577D+08-.2999330606D+08-.1441916681D+09 -.3140266220D+08-.3158607465D+080.4300451369D+08-.1334745823D+09-.5788433785D+08 0.3662830432D+060.4690913653D+060.1902120000D+06 0.1552987455D+05-.1437036744D+09-.3294284136D+08-.3003779327D+08-.1441401616D+09 -.3154718611D+08-.3158958570D+080.4311762755D+08-.1334440103D+09-.5787108666D+08 0.3662834276D+060.4691269454D+060.1902274830D+06 0.1552992115D+05-.1436636128D+09-.3305049579D+08-.3008226146D+08-.1440885655D+09 -.3169169017D+08-.3159307704D+080.4323071486D+08-.1334133567D+09-.5785780015D+08 0.3662838088D+060.4691625258D+060.1902429662D+06 0.1552996774D+05-.1436234612D+09-.3315812904D+08-.3012671062D+08-.1440368798D+09 -.3183617428D+08-.3159654867D+080.4334377553D+08-.1333826216D+09-.5784447831D+08 0.3662841870D+060.4691981066D+060.1902584497D+06 0.1553001434D+05-.1435832202D+09-.3326574006D+08-.3017114048D+08-.1439851045D+09 -.3198063835D+08-.3160000058D+080.4345680952D+08-.1333518049D+09-.5783112116D+08 0.3662845622D+060.4692336877D+060.1902739334D+06 0.1553006093D+05-.1435428887D+09-.3337333144D+08-.3021555187D+08-.1439332396D+09 -.3212508227D+08-.3160343276D+080.4356981675D+08-.1333209067D+09-.5781772871D+08 0.3662849342D+060.4692692692D+060.1902894173D+06 0.1553010753D+05-.1435024672D+09-.3348090164D+08-.3025994424D+08-.1438812851D+09 -.3226950597D+08-.3160684520D+080.4368279716D+08-.1332899270D+09-.5780430097D+08 0.3662853033D+060.4693048510D+060.1903049015D+06 0.1553015412D+05-.1434619556D+09-.3358845060D+08-.3030431756D+08-.1438292411D+09 -.3241390932D+08-.3161023789D+080.4379575068D+08-.1332588659D+09-.5779083795D+08 0.3662856693D+060.4693404332D+060.1903203860D+06 0.1553020072D+05-.1434213541D+09-.3369597824D+08-.3034867180D+08-.1437771075D+09 -.3255829225D+08-.3161361082D+080.4390867725D+08-.1332277233D+09-.5777733965D+08 0.3662860322D+060.4693760157D+060.1903358706D+06 0.1553024731D+05-.1433806626D+09-.3380348451D+08-.3039300692D+08-.1437248843D+09 -.3270265464D+08-.3161696398D+080.4402157681D+08-.1331964992D+09-.5776380608D+08 0.3662863920D+060.4694115986D+060.1903513555D+06 0.1553029391D+05-.1433398811D+09-.3391096932D+08-.3043732290D+08-.1436725715D+09 -.3284699642D+08-.3162029737D+080.4413444928D+08-.1331651938D+09-.5775023726D+08 0.3662867488D+060.4694471818D+060.1903668407D+06 0.1553034050D+05-.1432990096D+09-.3401843261D+08-.3048161971D+08-.1436201691D+09 -.3299131747D+08-.3162361097D+080.4424729461D+08-.1331338070D+09-.5773663319D+08 0.3662871026D+060.4694827654D+060.1903823261D+06 0.1553038710D+05-.1432580483D+09-.3412587432D+08-.3052589732D+08-.1435676772D+09 -.3313561770D+08-.3162690478D+080.4436011273D+08-.1331023388D+09-.5772299389D+08 0.3662874533D+060.4695183493D+060.1903978117D+06 0.1553043369D+05-.1432169970D+09-.3423329436D+08-.3057015570D+08-.1435150957D+09 -.3327989702D+08-.3163017879D+080.4447290357D+08-.1330707893D+09-.5770931935D+08 0.3662878009D+060.4695539336D+060.1904132975D+06 0.1553048029D+05-.1431758558D+09-.3434069267D+08-.3061439482D+08-.1434624246D+09 -.3342415532D+08-.3163343298D+080.4458566707D+08-.1330391585D+09-.5769560960D+08 0.3662881455D+060.4695895182D+060.1904287836D+06 0.1553052688D+05-.1431346248D+09-.3444806918D+08-.3065861464D+08-.1434096639D+09 -.3356839252D+08-.3163666736D+080.4469840316D+08-.1330074463D+09-.5768186464D+08 0.3662884870D+060.4696251031D+060.1904442700D+06 0.1553057348D+05-.1430933038D+09-.3455542383D+08-.3070281515D+08-.1433568137D+09 -.3371260851D+08-.3163988190D+080.4481111179D+08-.1329756529D+09-.5766808448D+08 0.3662888255D+060.4696606884D+060.1904597565D+06 0.1553062007D+05-.1430518931D+09-.3466275654D+08-.3074699631D+08-.1433038739D+09 -.3385680319D+08-.3164307661D+080.4492379288D+08-.1329437783D+09-.5765426913D+08 0.3662891609D+060.4696962741D+060.1904752433D+06 0.1553066667D+05-.1430103925D+09-.3477006723D+08-.3079115808D+08-.1432508445D+09 -.3400097648D+08-.3164625147D+080.4503644637D+08-.1329118224D+09-.5764041859D+08 0.3662894933D+060.4697318601D+060.1904907304D+06 0.1553071326D+05-.1429688020D+09-.3487735586D+08-.3083530045D+08-.1431977256D+09 -.3414512826D+08-.3164940648D+080.4514907220D+08-.1328797854D+09-.5762653289D+08 0.3662898225D+060.4697674464D+060.1905062177D+06 0.1553075986D+05-.1429271218D+09-.3498462233D+08-.3087942338D+08-.1431445171D+09 -.3428925846D+08-.3165254163D+080.4526167031D+08-.1328476671D+09-.5761261202D+08 0.3662901488D+060.4698030331D+060.1905217052D+06 0.1553080645D+05-.1428853518D+09-.3509186659D+08-.3092352683D+08-.1430912191D+09 -.3443336696D+08-.3165565690D+080.4537424061D+08-.1328154677D+09-.5759865599D+08 0.3662904720D+060.4698386201D+060.1905371929D+06 0.1553085305D+05-.1428434920D+09-.3519908856D+08-.3096761079D+08-.1430378315D+09 -.3457745367D+08-.3165875229D+080.4548678307D+08-.1327831872D+09-.5758466482D+08 0.3662907921D+060.4698742075D+060.1905526809D+06 0.1553089964D+05-.1428015425D+09-.3530628817D+08-.3101167522D+08-.1429843544D+09 -.3472151849D+08-.3166182779D+080.4559929760D+08-.1327508255D+09-.5757063852D+08 0.3662911091D+060.4699097952D+060.1905681691D+06 0.1553094624D+05-.1427595033D+09-.3541346536D+08-.3105572010D+08-.1429307877D+09 -.3486556133D+08-.3166488339D+080.4571178414D+08-.1327183828D+09-.5755657709D+08 0.3662914231D+060.4699453832D+060.1905836576D+06 0.1553099283D+05-.1427173743D+09-.3552062006D+08-.3109974538D+08-.1428771315D+09 -.3500958209D+08-.3166791909D+080.4582424264D+08-.1326858589D+09-.5754248054D+08 0.3662917341D+060.4699809716D+060.1905991463D+06 0.1553103943D+05-.1426751556D+09-.3562775219D+08-.3114375105D+08-.1428233858D+09 -.3515358067D+08-.3167093488D+080.4593667302D+08-.1326532541D+09-.5752834889D+08 0.3662920420D+060.4700165603D+060.1906146352D+06 0.1553108602D+05-.1426328472D+09-.3573486168D+08-.3118773707D+08-.1427695505D+09 -.3529755697D+08-.3167393074D+080.4604907523D+08-.1326205682D+09-.5751418213D+08 0.3662923468D+060.4700521494D+060.1906301244D+06 0.1553113262D+05-.1425904492D+09-.3584194847D+08-.3123170341D+08-.1427156257D+09 -.3544151090D+08-.3167690666D+080.4616144919D+08-.1325878013D+09-.5749998029D+08 0.3662926485D+060.4700877388D+060.1906456138D+06 0.1553117921D+05-.1425479616D+09-.3594901249D+08-.3127565005D+08-.1426616114D+09 -.3558544237D+08-.3167986265D+080.4627379485D+08-.1325549534D+09-.5748574337D+08 0.3662929472D+060.4701233286D+060.1906611034D+06 0.1553122581D+05-.1425053842D+09-.3605605367D+08-.3131957695D+08-.1426075076D+09 -.3572935126D+08-.3168279869D+080.4638611213D+08-.1325220246D+09-.5747147137D+08 0.3662932429D+060.4701589187D+060.1906765933D+06 0.1553127240D+05-.1424627173D+09-.3616307193D+08-.3136348408D+08-.1425533143D+09 -.3587323749D+08-.3168571478D+080.4649840098D+08-.1324890149D+09-.5745716432D+08 0.3662935355D+060.4701945091D+060.1906920834D+06 0.1553131900D+05-.1424199608D+09-.3627006721D+08-.3140737142D+08-.1424990314D+09 -.3601710096D+08-.3168861090D+080.4661066133D+08-.1324559242D+09-.5744282221D+08 0.3662938250D+060.4702300999D+060.1907075737D+06 0.1553136559D+05-.1423771147D+09-.3637703944D+08-.3145123893D+08-.1424446591D+09 -.3616094156D+08-.3169148705D+080.4672289312D+08-.1324227526D+09-.5742844506D+08 0.3662941115D+060.4702656910D+060.1907230643D+06 0.1553141219D+05-.1423341791D+09-.3648398855D+08-.3149508659D+08-.1423901972D+09 -.3630475921D+08-.3169434321D+080.4683509628D+08-.1323895002D+09-.5741403287D+08 0.3662943949D+060.4703012824D+060.1907385551D+06 0.1553145878D+05-.1422911539D+09-.3659091447D+08-.3153891437D+08-.1423356459D+09 -.3644855380D+08-.3169717939D+080.4694727076D+08-.1323561669D+09-.5739958566D+08 0.3662946752D+060.4703368742D+060.1907540461D+06 0.1553150538D+05-.1422480392D+09-.3669781713D+08-.3158272223D+08-.1422810051D+09 -.3659232525D+08-.3169999557D+080.4705941647D+08-.1323227528D+09-.5738510344D+08 0.3662949525D+060.4703724663D+060.1907695374D+06 0.1553155197D+05-.1422048349D+09-.3680469645D+08-.3162651015D+08-.1422262748D+09 -.3673607344D+08-.3170279175D+080.4717153337D+08-.1322892580D+09-.5737058621D+08 0.3662952267D+060.4704080587D+060.1907850289D+06 0.1553159857D+05-.1421615412D+09-.3691155238D+08-.3167027809D+08-.1421714551D+09 -.3687979828D+08-.3170556791D+080.4728362139D+08-.1322556823D+09-.5735603398D+08 0.3662954979D+060.4704436515D+060.1908005206D+06 0.1553164516D+05-.1421181580D+09-.3701838484D+08-.3171402603D+08-.1421165458D+09 -.3702349968D+08-.3170832405D+080.4739568046D+08-.1322220259D+09-.5734144676D+08 0.3662957660D+060.4704792446D+060.1908160126D+06 0.1553169176D+05-.1420746854D+09-.3712519376D+08-.3175775394D+08-.1420615472D+09 -.3716717754D+08-.3171106016D+080.4750771052D+08-.1321882887D+09-.5732682456D+08 0.3662960310D+060.4705148381D+060.1908315048D+06 0.1553173835D+05-.1420311233D+09-.3723197908D+08-.3180146179D+08-.1420064591D+09 -.3731083176D+08-.3171377624D+080.4761971150D+08-.1321544708D+09-.5731216740D+08 0.3662962930D+060.4705504319D+060.1908469972D+06 0.1553178495D+05-.1419874718D+09-.3733874071D+08-.3184514955D+08-.1419512815D+09 -.3745446224D+08-.3171647227D+080.4773168335D+08-.1321205723D+09-.5729747527D+08 0.3662965519D+060.4705860260D+060.1908624899D+06 0.1553183154D+05-.1419437309D+09-.3744547860D+08-.3188881719D+08-.1418960145D+09 -.3759806889D+08-.3171914824D+080.4784362600D+08-.1320865930D+09-.5728274820D+08 0.3662968078D+060.4706216204D+060.1908779828D+06 0.1553187814D+05-.1418999006D+09-.3755219268D+08-.3193246467D+08-.1418406580D+09 -.3774165160D+08-.3172180416D+080.4795553939D+08-.1320525332D+09-.5726798618D+08 0.3662970606D+060.4706572152D+060.1908934759D+06 0.1553192473D+05-.1418559809D+09-.3765888286D+08-.3197609198D+08-.1417852122D+09 -.3788521028D+08-.3172444001D+080.4806742344D+08-.1320183927D+09-.5725318922D+08 0.3662973103D+060.4706928103D+060.1909089693D+06 0.1553197133D+05-.1418119719D+09-.3776554910D+08-.3201969908D+08-.1417296769D+09 -.3802874483D+08-.3172705578D+080.4817927810D+08-.1319841716D+09-.5723835735D+08 0.3662975570D+060.4707284058D+060.1909244628D+06 0.1553201792D+05-.1417678736D+09-.3787219130D+08-.3206328593D+08-.1416740523D+09 -.3817225516D+08-.3172965147D+080.4829110331D+08-.1319498699D+09-.5722349056D+08 0.3662978006D+060.4707640015D+060.1909399567D+06 0.1553206452D+05-.1417236860D+09-.3797880941D+08-.3210685252D+08-.1416183382D+09 -.3831574117D+08-.3173222706D+080.4840289900D+08-.1319154877D+09-.5720858887D+08 0.3662980412D+060.4707995977D+060.1909554507D+06 0.1553211111D+05-.1416794090D+09-.3808540336D+08-.3215039882D+08-.1415625347D+09 -.3845920275D+08-.3173478256D+080.4851466510D+08-.1318810250D+09-.5719365228D+08 0.3662982787D+060.4708351941D+060.1909709450D+06 0.1553215771D+05-.1416350428D+09-.3819197308D+08-.3219392478D+08-.1415066419D+09 -.3860263982D+08-.3173731794D+080.4862640156D+08-.1318464817D+09-.5717868080D+08 0.3662985131D+060.4708707909D+060.1909864395D+06 0.1553220430D+05-.1415905874D+09-.3829851849D+08-.3223743039D+08-.1414506597D+09 -.3874605227D+08-.3173983322D+080.4873810831D+08-.1318118580D+09-.5716367445D+08 0.3662987444D+060.4709063879D+060.1910019342D+06 0.1553225090D+05-.1415460427D+09-.3840503953D+08-.3228091561D+08-.1413945881D+09 -.3888944000D+08-.3174232837D+080.4884978529D+08-.1317771538D+09-.5714863322D+08 0.3662989727D+060.4709419854D+060.1910174292D+06 0.1553229749D+05-.1415014088D+09-.3851153613D+08-.3232438042D+08-.1413384272D+09 -.3903280293D+08-.3174480339D+080.4896143243D+08-.1317423691D+09-.5713355714D+08 0.3662991980D+060.4709775831D+060.1910329244D+06 0.1553234409D+05-.1414566857D+09-.3861800821D+08-.3236782478D+08-.1412821769D+09 -.3917614094D+08-.3174725828D+080.4907304967D+08-.1317075040D+09-.5711844621D+08 0.3662994202D+060.4710131812D+060.1910484198D+06 0.1553239068D+05-.1414118734D+09-.3872445572D+08-.3241124867D+08-.1412258373D+09 -.3931945395D+08-.3174969302D+080.4918463694D+08-.1316725586D+09-.5710330044D+08 0.3662996393D+060.4710487796D+060.1910639155D+06 0.1553243728D+05-.1413669720D+09-.3883087857D+08-.3245465206D+08-.1411694083D+09 -.3946274185D+08-.3175210760D+080.4929619419D+08-.1316375327D+09-.5708811983D+08 0.3662998553D+060.4710843783D+060.1910794114D+06 0.1553248387D+05-.1413219814D+09-.3893727671D+08-.3249803491D+08-.1411128901D+09 -.3960600455D+08-.3175450203D+080.4940772135D+08-.1316024265D+09-.5707290441D+08 0.3663000683D+060.4711199774D+060.1910949075D+06 0.1553253047D+05-.1412769017D+09-.3904365005D+08-.3254139720D+08-.1410562825D+09 -.3974924195D+08-.3175687629D+080.4951921836D+08-.1315672400D+09-.5705765417D+08 0.3663002782D+060.4711555768D+060.1911104038D+06 0.1553257706D+05-.1412317330D+09-.3914999854D+08-.3258473890D+08-.1409995856D+09 -.3989245395D+08-.3175923038D+080.4963068514D+08-.1315319732D+09-.5704236913D+08 0.3663004851D+060.4711911765D+060.1911259004D+06 0.1553262366D+05-.1411864751D+09-.3925632210D+08-.3262805997D+08-.1409427995D+09 -.4003564046D+08-.3176156428D+080.4974212165D+08-.1314966260D+09-.5702704929D+08 0.3663006889D+060.4712267765D+060.1911413972D+06 0.1553267025D+05-.1411411282D+09-.3936262066D+08-.3267136040D+08-.1408859240D+09 -.4017880137D+08-.3176387800D+080.4985352782D+08-.1314611987D+09-.5701169467D+08 0.3663008896D+060.4712623768D+060.1911568942D+06 0.1553271685D+05-.1410956922D+09-.3946889415D+08-.3271464014D+08-.1408289594D+09 -.4032193659D+08-.3176617152D+080.4996490358D+08-.1314256910D+09-.5699630527D+08 0.3663010873D+060.4712979775D+060.1911723915D+06 0.1553276344D+05-.1410501672D+09-.3957514256D+08-.3275789917D+08-.1407719054D+09 -.4046504602D+08-.3176844484D+080.5007624886D+08-.1313901032D+09-.5698088110D+08 0.3663012819D+060.4713335785D+060.1911878890D+06 0.1553281004D+05-.1410045531D+09-.3968136587D+08-.3280113745D+08-.1407147622D+09 -.4060812956D+08-.3177069794D+080.5018756362D+08-.1313544352D+09-.5696542218D+08 0.3663014734D+060.4713691798D+060.1912033867D+06 0.1553285663D+05-.1409588501D+09-.3978756392D+08-.3284435496D+08-.1406575298D+09 -.4075118712D+08-.3177293083D+080.5029884778D+08-.1313186870D+09-.5694992850D+08 0.3663016619D+060.4714047815D+060.1912188846D+06 0.1553290323D+05-.1409130581D+09-.3989373661D+08-.3288755167D+08-.1406002081D+09 -.4089421859D+08-.3177514350D+080.5041010128D+08-.1312828586D+09-.5693440009D+08 0.3663018473D+060.4714403834D+060.1912343828D+06 0.1553294982D+05-.1408671772D+09-.3999988389D+08-.3293072756D+08-.1405427972D+09 -.4103722389D+08-.3177733593D+080.5052132405D+08-.1312469501D+09-.5691883695D+08 0.3663020296D+060.4714759857D+060.1912498812D+06 0.1553299642D+05-.1408212073D+09-.4010600569D+08-.3297388259D+08-.1404852972D+09 -.4118020290D+08-.3177950812D+080.5063251604D+08-.1312109615D+09-.5690323909D+08 0.3663022089D+060.4715115883D+060.1912653798D+06 0.1553304301D+05-.1407751486D+09-.4021210194D+08-.3301701674D+08-.1404277079D+09 -.4132315554D+08-.3178166007D+080.5074367719D+08-.1311748928D+09-.5688760651D+08 0.3663023851D+060.4715471913D+060.1912808786D+06 0.1553308961D+05-.1407290009D+09-.4031817256D+08-.3306012997D+08-.1403700295D+09 -.4146608170D+08-.3178379177D+080.5085480742D+08-.1311387441D+09-.5687193923D+08 0.3663025583D+060.4715827945D+060.1912963777D+06 0.1553313620D+05-.1406827644D+09-.4042421749D+08-.3310322226D+08-.1403122619D+09 -.4160898130D+08-.3178590320D+080.5096590667D+08-.1311025153D+09-.5685623726D+08 0.3663027284D+060.4716183981D+060.1913118770D+06 0.1553318280D+05-.1406364391D+09-.4053023666D+08-.3314629357D+08-.1402544051D+09 -.4175185422D+08-.3178799437D+080.5107697489D+08-.1310662065D+09-.5684050060D+08 0.3663028954D+060.4716540020D+060.1913273765D+06 0.1553322939D+05-.1405900249D+09-.4063623000D+08-.3318934388D+08-.1401964592D+09 -.4189470038D+08-.3179006527D+080.5118801200D+08-.1310298177D+09-.5682472927D+08 0.3663030593D+060.4716896062D+060.1913428762D+06 0.1553327599D+05-.1405435220D+09-.4074219744D+08-.3323237316D+08-.1401384242D+09 -.4203751967D+08-.3179211588D+080.5129901795D+08-.1309933490D+09-.5680892328D+08 0.3663032202D+060.4717252107D+060.1913583762D+06 0.1553332258D+05-.1404969302D+09-.4084813891D+08-.3327538139D+08-.1400803001D+09 -.4218031199D+08-.3179414621D+080.5140999268D+08-.1309568002D+09-.5679308262D+08 0.3663033780D+060.4717608155D+060.1913738764D+06 0.1553336918D+05-.1404502497D+09-.4095405434D+08-.3331836852D+08-.1400220869D+09 -.4232307726D+08-.3179615625D+080.5152093611D+08-.1309201716D+09-.5677720732D+08 0.3663035328D+060.4717964207D+060.1913893768D+06 0.1553341577D+05-.1404034804D+09-.4105994367D+08-.3336133453D+08-.1399637845D+09 -.4246581536D+08-.3179814598D+080.5163184819D+08-.1308834631D+09-.5676129738D+08 0.3663036845D+060.4718320261D+060.1914048775D+06 0.1553346237D+05-.1403566225D+09-.4116580681D+08-.3340427939D+08-.1399053931D+09 -.4260852621D+08-.3180011541D+080.5174272886D+08-.1308466746D+09-.5674535281D+08 0.3663038331D+060.4718676319D+060.1914203783D+06 0.1553350896D+05-.1403096758D+09-.4127164371D+08-.3344720308D+08-.1398469127D+09 -.4275120970D+08-.3180206453D+080.5185357804D+08-.1308098064D+09-.5672937363D+08 0.3663039787D+060.4719032380D+060.1914358794D+06 0.1553355556D+05-.1402626404D+09-.4137745430D+08-.3349010556D+08-.1397883432D+09 -.4289386575D+08-.3180399332D+080.5196439568D+08-.1307728582D+09-.5671335983D+08 0.3663041211D+060.4719388444D+060.1914513807D+06 0.1553360215D+05-.1402155164D+09-.4148323849D+08-.3353298680D+08-.1397296846D+09 -.4303649424D+08-.3180590179D+080.5207518172D+08-.1307358303D+09-.5669731143D+08 0.3663042606D+060.4719744511D+060.1914668823D+06 0.1553364875D+05-.1401683038D+09-.4158899624D+08-.3357584678D+08-.1396709371D+09 -.4317909508D+08-.3180778992D+080.5218593609D+08-.1306987226D+09-.5668122844D+08 0.3663043969D+060.4720100582D+060.1914823840D+06 0.1553369534D+05-.1401210025D+09-.4169472745D+08-.3361868547D+08-.1396121005D+09 -.4332166818D+08-.3180965772D+080.5229665873D+08-.1306615351D+09-.5666511086D+08 0.3663045302D+060.4720456655D+060.1914978860D+06 0.1553374194D+05-.1400736127D+09-.4180043208D+08-.3366150283D+08-.1395531750D+09 -.4346421343D+08-.3181150517D+080.5240734957D+08-.1306242678D+09-.5664895871D+08 0.3663046604D+060.4720812732D+060.1915133882D+06 0.1553378853D+05-.1400261342D+09-.4190611004D+08-.3370429884D+08-.1394941604D+09 -.4360673074D+08-.3181333226D+080.5251800856D+08-.1305869208D+09-.5663277200D+08 0.3663047876D+060.4721168812D+060.1915288907D+06 0.1553383513D+05-.1399785673D+09-.4201176127D+08-.3374707347D+08-.1394350569D+09 -.4374922001D+08-.3181513900D+080.5262863563D+08-.1305494942D+09-.5661655073D+08 0.3663049117D+060.4721524895D+060.1915443933D+06 0.1553388172D+05-.1399309118D+09-.4211738570D+08-.3378982669D+08-.1393758645D+09 -.4389168114D+08-.3181692537D+080.5273923071D+08-.1305119878D+09-.5660029492D+08 0.3663050327D+060.4721880981D+060.1915598962D+06 0.1553392832D+05-.1398831677D+09-.4222298326D+08-.3383255847D+08-.1393165831D+09 -.4403411404D+08-.3181869137D+080.5284979375D+08-.1304744018D+09-.5658400456D+08 0.3663051507D+060.4722237070D+060.1915753993D+06 0.1553397491D+05-.1398353352D+09-.4232855388D+08-.3387526878D+08-.1392572128D+09 -.4417651860D+08-.3182043699D+080.5296032469D+08-.1304367362D+09-.5656767969D+08 0.3663052656D+060.4722593162D+060.1915909026D+06 0.1553402151D+05-.1397874143D+09-.4243409749D+08-.3391795760D+08-.1391977536D+09 -.4431889474D+08-.3182216223D+080.5307082345D+08-.1303989909D+09-.5655132029D+08 0.3663053774D+060.4722949257D+060.1916064062D+06 0.1553406810D+05-.1397394048D+09-.4253961402D+08-.3396062488D+08-.1391382056D+09 -.4446124234D+08-.3182386707D+080.5318128997D+08-.1303611660D+09-.5653492638D+08 0.3663054861D+060.4723305356D+060.1916219099D+06 0.1553411470D+05-.1396913070D+09-.4264510340D+08-.3400327062D+08-.1390785686D+09 -.4460356132D+08-.3182555152D+080.5329172420D+08-.1303232616D+09-.5651849798D+08 0.3663055918D+060.4723661457D+060.1916374139D+06 0.1553416129D+05-.1396431208D+09-.4275056557D+08-.3404589476D+08-.1390188428D+09 -.4474585157D+08-.3182721557D+080.5340212607D+08-.1302852777D+09-.5650203508D+08 0.3663056944D+060.4724017562D+060.1916529181D+06 0.1553420789D+05-.1395948461D+09-.4285600045D+08-.3408849730D+08-.1389590282D+09 -.4488811300D+08-.3182885921D+080.5351249552D+08-.1302472142D+09-.5648553770D+08 0.3663057940D+060.4724373669D+060.1916684226D+06 0.1553425448D+05-.1395464832D+09-.4296140797D+08-.3413107819D+08-.1388991247D+09 -.4503034552D+08-.3183048243D+080.5362283248D+08-.1302090712D+09-.5646900585D+08 0.3663058905D+060.4724729780D+060.1916839272D+06 0.1553430108D+05-.1394980318D+09-.4306678807D+08-.3417363740D+08-.1388391324D+09 -.4517254901D+08-.3183208523D+080.5373313690D+08-.1301708487D+09-.5645243954D+08 0.3663059839D+060.4725085894D+060.1916994321D+06 0.1553434767D+05-.1394494922D+09-.4317214068D+08-.3421617492D+08-.1387790514D+09 -.4531472339D+08-.3183366761D+080.5384340870D+08-.1301325468D+09-.5643583877D+08 0.3663060742D+060.4725442011D+060.1917149372D+06 0.1553439427D+05-.1394008643D+09-.4327746572D+08-.3425869071D+08-.1387188816D+09 -.4545686855D+08-.3183522955D+080.5395364783D+08-.1300941654D+09-.5641920356D+08 0.3663061615D+060.4725798130D+060.1917304425D+06 0.1553444086D+05-.1393521481D+09-.4338276313D+08-.3430118474D+08-.1386586230D+09 -.4559898441D+08-.3183677106D+080.5406385422D+08-.1300557047D+09-.5640253391D+08 0.3663062457D+060.4726154253D+060.1917459480D+06 0.1553448746D+05-.1393033437D+09-.4348803284D+08-.3434365698D+08-.1385982757D+09 -.4574107085D+08-.3183829211D+080.5417402782D+08-.1300171645D+09-.5638582984D+08 0.3663063268D+060.4726510379D+060.1917614538D+06 0.1553453405D+05-.1392544510D+09-.4359327478D+08-.3438610741D+08-.1385378396D+09 -.4588312779D+08-.3183979272D+080.5428416855D+08-.1299785450D+09-.5636909135D+08 0.3663064049D+060.4726866508D+060.1917769597D+06 0.1553458065D+05-.1392054701D+09-.4369848888D+08-.3442853599D+08-.1384773149D+09 -.4602515512D+08-.3184127288D+080.5439427635D+08-.1299398461D+09-.5635231846D+08 0.3663064799D+060.4727222640D+060.1917924659D+06 0.1553462724D+05-.1391564011D+09-.4380367522D+08-.3447094268D+08-.1384167015D+09 -.4616715275D+08-.3184273257D+080.5450435117D+08-.1299010680D+09-.5633551116D+08 0.3663065518D+060.4727578776D+060.1918079723D+06 0.1553467384D+05-.1391072438D+09-.4390883359D+08-.3451332746D+08-.1383559994D+09 -.4630912059D+08-.3184417179D+080.5461439294D+08-.1298622105D+09-.5631866948D+08 0.3663066207D+060.4727934914D+060.1918234790D+06 0.1553472043D+05-.1390579984D+09-.4401396391D+08-.3455569031D+08-.1382952087D+09 -.4645105852D+08-.3184559054D+080.5472440160D+08-.1298232737D+09-.5630179342D+08 0.3663066865D+060.4728291055D+060.1918389858D+06 0.1553476703D+05-.1390086649D+09-.4411906611D+08-.3459803120D+08-.1382343293D+09 -.4659296646D+08-.3184698881D+080.5483437708D+08-.1297842577D+09-.5628488299D+08 0.3663067492D+060.4728647199D+060.1918544929D+06 0.1553481362D+05-.1389592433D+09-.4422414013D+08-.3464035009D+08-.1381733613D+09 -.4673484431D+08-.3184836659D+080.5494431932D+08-.1297451625D+09-.5626793820D+08 0.3663068089D+060.4729003346D+060.1918700002D+06 0.1553486022D+05-.1389097337D+09-.4432918590D+08-.3468264697D+08-.1381123048D+09 -.4687669197D+08-.3184972389D+080.5505422826D+08-.1297059880D+09-.5625095905D+08 0.3663068655D+060.4729359496D+060.1918855077D+06 0.1553490681D+05-.1388601360D+09-.4443420334D+08-.3472492181D+08-.1380511596D+09 -.4701850934D+08-.3185106069D+080.5516410384D+08-.1296667344D+09-.5623394557D+08 0.3663069190D+060.4729715649D+060.1919010154D+06 0.1553495341D+05-.1388104502D+09-.4453919239D+08-.3476717456D+08-.1379899259D+09 -.4716029632D+08-.3185237698D+080.5527394599D+08-.1296274016D+09-.5621689775D+08 0.3663069694D+060.4730071806D+060.1919165233D+06 0.1553500000D+05-.1387606765D+09-.4464415297D+08-.3480940521D+08-.1379286037D+09 -.4730205282D+08-.3185367277D+080.5538375465D+08-.1295879897D+09-.5619981560D+08 0.3663070168D+060.4730427965D+060.1919320315D+06 Appendix A2 Source Module for Subroutine GIODAT C----------------------------------------------------------------------- SUBROUTINE GIODAT(TUTC,PVG,PVH,PVE,PVS,IERR) C------------------------------------------------- C S.J ISSUE 04 - 14.01.85 C------------------------ C ************ SUBROUTINE IDENTIFICATION ************ * * * GIODAT IS AN INTERFACE ROUTINE TO THE EXPERIMENTER SUPPORT DATA * * FILE ON LOGICAL UNIT 7 WHICH CONTAINS POSITION DATA FOR THE * * GIOTTO PROBE, TARGET COMET, THE EARTH AND THE SUN GIVEN IN THE * * BARYCENTRIC REFERENCE FRAME OF 1950 ( UNITS = KM ) * * THE APPROXIMATE VELOCITIES FOR THE FOUR BODIES WILL BE DERIVED * * FROM THE POSITION DATA AND GIVEN IN KM/S. * * * ************************************************************************ * * ************ VARIABLE IDENTIFICATIONS ************* * * ***** INPUT VARIABLES ( ARGUMENT LIST ) *** * * * TUTC - UNIVERSAL COORDINATED TIME * * * ***** OUTPUT VARIABLES ( ARGUMENT LIST ) *** * * * PVG(1-3) - POSITION OF THE GIOTTO PROBE AT TIME TUTC * * PVG(4-6) - VELOCITY OF THE GIOTTO PROBE AT TIME TUTC * * PVH(1-3) * * PVH(4-6) - DITTO FOR COMET * * PVE(1-3) * * PVE(4-6) - DITTO FOR THE EARTH * * PVS(1-3) * * PVS(4-6) - DITTO FOR THE SUN * * IERR - ERROR CODE * * =0 : NO ERROR * * =1 : TUTC TOO EARLY * * =2 : TUTC TOO LATE * * * ***** LOCAL VARIABLES *** * * * ET - MODIFIED JULIAN DAY ( 1950.0 ) * * PGI(3,3) - TO STORE POSITION DATA FOR THREE SAMPLES FOR THE * * GIOTTO PROBE. USED IN THE QUADRATIC INTERPOLATION * * PHI(3,3) - DITTO FOR COMET * * PEI(3,3) - DITTO FOR THE EARTH * * PSI(3,3) - DITTO FOR THE SUN * * T(3) - TO STORE REFERENCE TIME T0 FOR THE INTERPOLATION. * * SKIP(3) - USED FOR SKIPPING RECORDS ON THE DATA FILE * * * ***** CONSTANT IDENTIFICATION *** * * * HSEC - (STEP SIZE)*86400 * * LUN - LOGICAL UNIT NUMBER OF FILE/TAPE CONTAINING * * THE POSITION DATA. * * STPUTC - STEP IN UTC AT TIME SUTC * * * ***** TEMPORARY VARIABLES *** * * * (1) * * TMP(1) - 0.5*(P-1)*(P-2) * * TMP(2) - P*(P-2) * * TMP(3) - 0.5*P*(P-1) * * (2) * * TMP(1) - (P-1.5)/HSEC * * TMP(2) - (2*P-2)/HSEC * * TMP(3) - (P-0.5)/HSEC * * WHERE ET(P) IS THE EPOCH FOR WHICH THE POSITIONS ARE REQUIRED. * * ( RANGE OF P : 0:2 ) * * * ************************************************************************ * * ************ TYPE DECLARATION AND STORAGE ALLOCATION ************ * * IMPLICIT REAL*8(A-H,O-Z) LOGICAL FIRST CHARACTER*1 SKIP(3) DIMENSION PVG(6),PVH(6),PVE(6),PVS(6),PGI(3,3), & PHI(3,3),PEI(3,3),PSI(3,3),T(3),TMP(3) DATA FIRST/.TRUE./ DATA LUN/7/ * * ************************************************************************ * * ************ INITIALISE AND VERIFY INPUT DATA ***BLOCK0000 * * IERR=0 C IF FIRST CALL READ CONTROL RECORD AND SET CONSTANTS IF(FIRST)THEN READ(LUN,10) IYGEN,MOGEN,DGEN,NS,NR,ET1,ET2,SUTC,ETUTC1,ETUTC2 10 FORMAT(2I5,D16.10,2I5,/,5D16.10) C DETERMINE STEP IN UTC AT TIME SUTC STPUTC=ETUTC2-ETUTC1 C DETERMINE STEP SIZE H=(ET2-ET1)/(DFLOAT(NR)-1.D0) HSEC=H*86400.D0 FIRST=.FALSE. ENDIF C DETERMINE ET ( MODIFIED JULIAN DAY ) ET=TUTC+ETUTC1 IF(ET.GE.SUTC)ET=ET+STPUTC C VERIFY EPOCH ON FILE IF(ET.LT.ET1-1.D-4)GOTO 910 IF(ET.GT.ET2+1.D-4)GOTO 920 DO 20 I=1,6 PVG(I)=0.D0 PVH(I)=0.D0 PVE(I)=0.D0 PVS(I)=0.D0 20 CONTINUE * * ************************************************************************ * * ******* DETERMINE THE POSITION VECTORS FOR EPOCH ET **BLOCK0100 * * C FIND RECORD NUMBER IREC=DINT((ET-ET1)/H)+1 110 CONTINUE C POSITION THE T0 SAMPLE TO ENSURE T1 AND T2 AVAILABLE IF(NR-IREC.LT.2)THEN IREC=IREC-1 GOTO 110 ENDIF IF(FIRST)THEN IF(IREC.EQ.1)GOTO 135 DO 130 I=1,IREC-1 READ(LUN,120)SKIP 120 FORMAT(A1) 130 CONTINUE GOTO 135 ENDIF IDIF=IREC-IRECOL-1 IF(IDIF.EQ.0)GOTO 135 ISKP=ABS(IDIF) IF(IDIF.LT.0)THEN DO 131 I=1,ISKP*3 BACKSPACE LUN 131 CONTINUE ELSE DO 132 I=1,ISKP READ(LUN,120)SKIP 132 CONTINUE ENDIF 135 CONTINUE IRECOL=IREC+2 DO 150 I=1,3 READ(LUN,140)T(I),(PGI(I,J),J=1,3),(PHI(I,J),J=1,3), & (PEI(I,J),J=1,3),(PSI(I,J),J=1,3) 140 FORMAT(5D16.10) 150 CONTINUE * * * * ************************************************************************ * * ******* DETERMINE POSITIONS USING QUADRATIC INTERPOLATION **BLOCK0200 * * P=(ET-T(1))/H TMP(1)=0.5*(P-1.D0)*(P-2.D0) TMP(2)=P*(2.D0-P) TMP(3)=0.5D0*P*(P-1.D0) DO 220 I=1,3 DO 210 J=1,3 PVG(I)=PVG(I)+PGI(J,I)*TMP(J) PVH(I)=PVH(I)+PHI(J,I)*TMP(J) PVE(I)=PVE(I)+PEI(J,I)*TMP(J) PVS(I)=PVS(I)+PSI(J,I)*TMP(J) 210 CONTINUE 220 CONTINUE * * ************************************************************************ * * ******* DETERMINE VELOCITIES FROM DERIVATIVES OF POSITIONS **BLOCK0300 * * TMP(1)=(P-1.5D0)/HSEC TMP(2)=(2.D0-2.D0*P)/HSEC TMP(3)=(P-0.5D0)/HSEC DO 320 I=1,3 I1=I+3 DO 310 J=1,3 PVG(I1)=PVG(I1)+PGI(J,I)*TMP(J) PVH(I1)=PVH(I1)+PHI(J,I)*TMP(J) PVE(I1)=PVE(I1)+PEI(J,I)*TMP(J) PVS(I1)=PVS(I1)+PSI(J,I)*TMP(J) 310 CONTINUE 320 CONTINUE * * ************************************************************************ * * ************ ERROR MESSAGE BLOCK ***BLOCK0900 * * GOTO 1000 910 CONTINUE IERR=1 WRITE(6,911) 911 FORMAT(/' GIOPOS : ET TOO EARLY ') GOTO 1000 920 CONTINUE IERR=2 WRITE(6,921) 921 FORMAT(/' GIOPOS : ET TOO LATE ') 1000 CONTINUE RETURN END Appendix A3 Source Module for Subroutine TRAVEL C----------------------------------------------------------------------- SUBROUTINE TRAVEL(TUTC,S,DTG,P,G,IER) C------------------------------------------------- C ISSUE 01 - 18.01.85 C---------------------- C ************ SUBROUTINE IDENTIFICATION ************ * * * TRAVEL PROVIDES THE FOLLOWING GIOTTO EXPERIMENTER SUPPORT DATA: * * - SIGNAL TRAVEL TIME BETWEEN THE GIOTTO S/C AND A STATION ON EARTH * * - POSITIONS OF STATION AND S/C IN THE HELIOCENTRIC EME50 * * AT SIGNAL RECEPTION AND TRANSMISSION TIMES INVOLVED. * * STATION COORDINATES ARE COMPUTED WITHOUT PRECESSION, I.E. * * THE TRAVEL TIMES GIVEN ARE ACCURATE TO ABOUT 1 MSEC. * * TRAVEL RETRIEVES THE STATES OF THE PARTICIPATING BODIES FROM * * THE POSITION DATA FILE BY MEANS OF THE INTERFACE ROUTINE GIODAT. * * * ************************************************************************ * * ************ VARIABLE IDENTIFICATIONS ************* * * ***** INPUT VARIABLES ( ARGUMENT LIST ) *** * * * TUTC - SIGNAL RECEPTION TIME (UTC DAYS SINCE 1950.0) * * S(1-3) - CARTESIAN COORDINATES OF STATION IN BODY-FIXED SYSTEM * * * ***** OUTPUT VARIABLES ( ARGUMENT LIST ) *** * * * DTG - SIGNAL TRAVEL TIME(SEC) GIOTTO S/C - STATION * * P(1-3) - POSITION OF STATION AT TUTC IN HELIOCENTRIC EME50 * * G(1-3) - POSITION OF GIOTTO AT TUTC-DTG IN HELIOCENTRIC EME50 * * IER - ERROR CODE * * =0 : NO ERROR * * =1 : TUTC TOO EARLY, NO DATA ON POSITION DATA FILE * * =2 : TUTC TOO LATE , NO DATA ON POSITION DATA FILE * * =3 : NO CONVERGENCE IN SOLUTION OF LIGHT TIME EQUATION* * * ***** LOCAL VARIABLES OF SPECIAL INTEREST *** * * * SG(6) GIOTTO STATE IN BARYCENTRIC EME50 * * SH(6) - DITTO FOR COMET * * SU(6) - DITTO FOR SUN * * * ***** CONSTANT IDENTIFICATION *** * * * RJ,RK,RL - COEFFICIENTS IN FORMULA FOR SIDERIAL ANGLE * * C - VELOCITY OF LIGHT * * TWOPI - TWO PI * * * ************************************************************************ * * ************ TYPE DECLARATION AND STORAGE ALLOCATION ************ * * IMPLICIT REAL*8(A-H,O-Z) DIMENSION S(3),P(3),G(3),SS(6),SG(6),SH(6),SU(6) * * ************************************************************************ * * ************ SET CONSTANTS ***BLOCK0000 DATA RJ,RK,RL /0.2769194D0,100.0021359028D0,0.107523D-5/ DATA C,TWOPI /299792.458D0,6.283185307179587D0/ * * ************************************************************************ * * ******* STATION STATE AT TUTC **BLOCK0100 * * CALL GIODAT(TUTC,SG,SH,SS,SU,IER) IF(IER.NE.0) RETURN Y=0.5D0+TUTC/36525.D0 IT=TUTC X=TUTC-IT X=(X+RJ+(RK+RL*Y)*Y)*TWOPI Y=DCOS(X) Z=DSIN(X) P(1)=S(1)*Y-S(2)*Z+SS(1) P(2)=S(1)*Z+S(2)*Y+SS(2) P(3)=S(3)+SS(3) G(1)=SU(1) G(2)=SU(2) G(3)=SU(3) * * ************************************************************************ * * ******* SOLUTION OF LIGHT TIME EQUATION **BLOCK0200 * * DTG=0.D0 DO 210 I=1,5 X=SG(1)-P(1) Y=SG(2)-P(2) Z=SG(3)-P(3) DGB=DSQRT(X*X+Y*Y+Z*Z) Z=(X*SG(4)+Y*SG(5)+Z*SG(6))/DGB Y=(DTG-DGB/C)/(1.D0-Z/C) DTG=DTG-Y IF(DABS(Y).LT.1.D-4) GOTO 1000 X=TUTC-DTG/864.D2 CALL GIODAT(X,SG,SH,SS,SU,IER) IF(IER.NE.0) RETURN 210 CONTINUE * * ************************************************************************ * * ******* ERROR MESSAGES **BLOCK0900 * * IER=3 WRITE(6,910) IER 910 FORMAT(' TRAVEL ERROR MESSAGE',I3/ &' = NO CONVERGENCE AFTER 5 ITERATIONS'/ &' IN SOLUTION OF LIGHT TIME EQUATION') * * ************************************************************************ * * ******* SET OUTPUT VARIABLES AND RETURN **BLOCK1000 * * 1000 P(1)=P(1)-G(1) P(2)=P(2)-G(2) P(3)=P(3)-G(3) G(1)=SG(1)-SU(1) G(2)=SG(2)-SU(2) G(3)=SG(3)-SU(3) RETURN END Appendix A4 Source Module for Suroutine PETRA SUBROUTINE PETRA(XMUP,X,INIT,KODE,T,Y,RY,AUT) CP GIVES STATE VECTOR OF AN ELLIPTICAL, PARABOLICAL OR CP HYPERBOLICAL TRAJECTORY CENTRED AT THE ORIGIN OF THE COORD. SYST. C SAME AS 'TRAJEC', EXCEPT FOR ADDITIONAL INPUT: 'XMUP'. C CI XMUP = GRAVITY POTENTIAL OF CENTRAL BODY CI X( ) = POSITION (X(1)..X(3)) AND VELOCITY (X(4)..X(6)), (KM,KM/SEC) CI INIT = 0 FOR EACH NEW X. OLD VALUE SHOULD HAVE INIT.GT.0 CI KODE = 1 ONLY POSITION IN Y(1)..Y(3) CI = 2 POSITION + VELOCITY IN Y(1)..Y(6) CI E = 14: AS 2 + OUTPUT VECTOR AUT(14) FOR PARTIAL DERIVATIVES. CI T = TIME DIFFERENCE (SEC) BETWEEN X AND Y C CO Y( ) = STATE VECTOR AFTER TIME T CO RY = SQRT(Y(1)**2+Y(2)**2+Y(3)**2) ONLY IF KODE.GE.2 CO AUT(1) = INVERSE SEMI-MAJ.AXIS; NEGATIVE FOR HYPERBOLAE. CO AUT(2)...AUT(14) = INPUT TO SUBROUTINE 'PEDER' FOR PARTIAL CO DERIVATIVES. ONLY IF KODE.GE.3 C IMPLICIT REAL*8(A-H,O-Z) DIMENSION X(6),Y(6),AUT(KODE),C2(7),C3(7) DATA KNIT/0/ COEFFICIENTS FOR STUMPFF-FUNCTIONS DATA C3/0.7650283228592385 D-12,-0.1605953765319026 D-9, & 0.2505210785999854 D-7,-0.2755731919939401 D-5, & 0.1984126984129264 D-3, -0.8333333333333568D-2, & 0.166666666666666665 D0/,C2/0.1147134108311665 D-10, & -0.2087759200397967 D-8,0.2755731917059028 D-6, & -0.2480158725995993 D-4,0.1388888888888999 D-2, & -0.4166666666667176 D-1,0.4999999999999998 D0/ IF(INIT*KNIT .GT. 0) GOTO 115 KNIT = 1 WMUP = DSQRT(XMUP) RX = DSQRT(X(1)**2+X(2)**2+X(3)**2) A = 2.D0/RX -(X(4)**2+X(5)**2+X(6)**2)/XMUP QA = DSQRT(DABS(A)) RVWM= (X(1)*X(4)+X(2)*X(5)+X(3)*X(6))/WMUP 115 TM = T*WMUP C FIND SUITABLE INITIAL VALUE FOR B AND FOR MAXIMUM STEP LENGTH: DBMAX = DMAX1(DABS(A*TM),DABS(TM)/RX) C IF WE ARE FAR OUT ON A PARABOLA: ATM6 = 6.D0*DABS(TM) IF(DBMAX**3.LT.ATM6) DBMAX = ATM6**(1.D0/3.D0) IF(A.GE.-1.D-30) GOTO 135 C IF WE ARE FAR OUT ON A HYPERBOLA: UP = DABS(TM*QA*A*2.D0) DO = 1.D0 - A*RX IF(UP.LE.DO) GOTO 135 DBMAX = DMIN1(DBMAX,DLOG(UP/DO)/QA) 135 B = DSIGN(DBMAX,TM) C MAXIMUM STEP = 1 RAD FOR A PARABOLA: IF(A*DBMAX.GT.QA) DBMAX = 1.D0/QA DB = 0.D0 C END OF INITIALISATION, START NEWTON'S ITERATION FOR ECC. ANOM. 155 W = B*QA C USE STUMPFF-SERIES WHEN CLOSE TO 0: IF(DABS(W).GE.0.3D0) GOTO 242 BB = B*B WW = BB*A G2 = C2(1) G3 = C3(1) DO 221 I=2,7 G2 = G2*WW + C2(I) 221 G3 = G3*WW + C3(I) G2 = G2*BB G3 = G3*BB*B GO = 1.D0-A*G2 G1 = B-A*G3 GOTO 595 242 IF(A.LT.0.D0) GOTO 511 C ONLY ELLIPSES GO = DCOS(W) G1 = DSIN(W)/QA GOTO 553 C ONLY HYPERBOLAE: 511 GO = DCOSH(W) G1 = DSINH(W)/QA 553 G2 = (1.D0-GO)/A G3 = (B-G1)/A 595 DOLD = DB DB = (TM-RX*G1-RVWM*G2-G3 & )/(RX*GO+RVWM*G1+G2) C DECREASE MAX STEP IF WE HAVE ALTERNATING SIGNS: IF(DB*DOLD.LT.-1.D-20) DBMAX = DBMAX*0.5D0 C PUT A LIMIT TO THE STEP LENGTH: DB = DMAX1(-DBMAX,DMIN1(DBMAX,DB)) B = B + DB CHECK IF WE HAVE CONVERGED: IF(DABS(DB).GT.DABS(B)*1.D-8) GOTO 155 C ITERATION END; UPDATE THE G'S LINEARLY FROM LAST STEP. GO = GO-DB*G1*A G1 = G1+DB*GO G2 = G2+DB*G1 G3 = G3+DB*G2 FX = 1.D0-G2/RX GX = T-G3/WMUP DO 635 J = 1,3 635 Y(J) = FX*X(J)+GX*X(J+3) AUT(1) = A IF(KODE .LE. 1) RETURN RY = DSQRT(Y(1)**2+Y(2)**2+Y(3)**2) FT = -WMUP*G1/(RX*RY) GT = 1.D0-G2/RY DO 660 J=4,6 660 Y(J) = FT*X(J-3)+GT*X(J) IF(KODE .LE. 2) RETURN C PREPARATION OF PARTIAL DERIVATIVES AUT(2)=B AUT(3)=RVWM AUT(4)=GO AUT(5)=G1 AUT(6)=G2 AUT(7)=G3 AUT(8)=RX AUT(9)=FX AUT(10)=GX AUT(11)=FT AUT(12)=GT AUT(13)=XMUP AUT(14)=WMUP RETURN END