DOUBLE PRECISION FUNCTION iau_DTDB ( EPOCH1, EPOCH2, UT, : ELONG, U, V ) *+ * - - - - - - - - - * i a u _ D T D B * - - - - - - - - - * * The periodic part of the offset between coordinate time in the * solar system barycenter frame of reference and proper time for a * clock on the Earth. * * The routine is an implementation of TDB-TT, where TDB is barycentric * dynamical time and TT is terrestrial time. It also models the * periodic component of the conversion from geocentric coordinate time * (TCG) to barycentric coordinate time (TCB). * * The different proper/coordinate timescales are related to each other * and to TAI as shown: * * TAI <- physically realized * : * offset <- observed (currently +32.184s) * : * TT <- proper time on the geoid * : * rate adjustment (L_G) <- definition of TT * : * TCG <- coordinate time at the geocenter * : * periodic terms <- iau_DTDB is an implementation * : * rate adjustment (L_C) <- modeled * : * TCB <- coordinate time at the barycenter * : * rate adjustment (-L_B) <- modeled * : * TDB <- superseded * : * periodic terms <- -iau_DTDB is an implementation * : * TT * * Adopted values for the various constants can be found in the IERS * Conventions (McCarthy 1996). * * This routine is part of the International Astronomical Union's * SOFA (Standards of Fundamental Astronomy) software collection. * * Status: canonical model. * * Given: * EPOCH1,EPOCH2 d epoch, TDB (Notes 1-3) * UT d universal time (UT1, fraction of one day) * ELONG d longitude (east +ve, radians) * U d distance from Earth spin axis (km) * V d distance north of equatorial plane (km) * * Returned: * iau_DTDB d TDB-TT (seconds) * * Notes: * * 1) The epoch EPOCH1+EPOCH2 is a Julian Date, apportioned in any * convenient way between the arguments EPOCH1 and EPOCH2. For * example, JD(TDB)=2450123.7 could be expressed in any of these * ways, among others: * * EPOCH1 EPOCH2 * * 2450123.7D0 0D0 (JD method) * 2451545D0 -1421.3D0 (J2000 method) * 2400000.5D0 50123.2D0 (MJD method) * 2450123.5D0 0.2D0 (date & time method) * * The JD method is the most natural and convenient to use in * cases where the loss of several decimal digits of resolution * is acceptable. The J2000 method is best matched to the way * the argument is handled internally and will deliver the * optimum resolution. The MJD method and the date & time methods * are both good compromises between resolution and convenience. * * Although the epoch is, formally, barycentric dynamical time (TDB), * the terrestrial dynamical time (TT) can be used with no practical * effect on the accuracy of the prediction. * * 2) TDB is a canonical implementation of the coordinate time in the * Solar System barycenter frame of reference, and TT is the proper * time given by clocks at mean sea level on the Earth. The present * routine provides a TDB-TT implementation. The result, which is * in seconds, has a main (annual) sinusoidal term of amplitude * approximately 0.00166 seconds, plus planetary terms up to about * 20 microseconds, and lunar and diurnal terms up to 2 microseconds. * The variation arises from the transverse Doppler effect and the * gravitational red-shift as the observer varies in speed and * experiences different gravitational potentials. * * 3) The IAU 1976 definition of TDB was that it must differ from TT * (then called TDT) only by periodic terms. Though practical, this * is an imprecise definition which ignores the existence of very * long-period and secular effects in the dynamics of the solar * system. Consequently, different implementations of TDB will, in * general, have different zero-points and will drift linearly * relative to one other. Subsequent IAU resolutions (1991) * introduced new timescales to overcome these objections, namely * TCG and TCB (see Seidelmann et al. 1992 for details). * * 4) The geocentric model is that of Fairhead & Bretagnon (1990), in * its full form. It was originally supplied by Fairhead (private * communications with P.T.Wallace, 1990) in the form of a Fortran * subroutine. The present routine contains an adaptation of the * Fairhead code. The numerical results are essentially unaffected * by the changes, the differences with respect to the Fairhead & * Bretagnon original being at the 1E-20 sec level. This model has * been endorsed by the IERS Conventions (1996). * * The topocentric part of the model is from Moyer (1981) and * Murray (1983), with fundamental arguments adapted from * Simon et al. 1994. It is an approximation to the expression * ( v / c ) . ( r / c ), where v is the barycentric velocity of * the Earth, r is the geocentric position of the observer and * c is the speed of light. * * By supplying zeroes for U and V, the topocentric part of the * model can be nullified, and the routine will return the Fairhead * & Bretagnon result alone. * * 5) During the interval 1950-2050, the absolute accuracy is better * than +/- 3 nanoseconds relative to direct numerical integrations * using the JPL DE200/LE200 solar system ephemeris. * * 6) It must be stressed that the present routine is merely a model, * and that numerical integration of solar-system ephemerides is the * definitive method for predicting the relationship between TCG and * TCB. In particular, published values of L_C, the constant * component of the relationship, may not be optimally matched with * the periodic component modeled here. * * References: * * Fairhead,L., & Bretagnon,P., Astron.Astrophys., 229, 240-247 * (1990). * * McCarthy,D.D., IERS Conventions (1996), IERS Technical Note 21, * Observatoire de Paris (1996). * * Moyer,T.D., Cel.Mech., 23, 33 (1981). * * Murray,C.A., Vectorial Astrometry, Adam Hilger (1983). * * Seidelmann,P.K. et al., Explanatory Supplement to the * Astronomical Almanac, Chapter 2, University Science Books * (1992). * * Simon J.L., Bretagnon P., Chapront J., Chapront-Touze M., * Francou G. & Laskar J., 1994, Astron.Astrophys., 282, 663-683. * * This revision: 2003 March 11 * * Copyright (C) 2003 IAU SOFA Review Board. See notes at end. * *----------------------------------------------------------------------- IMPLICIT NONE DOUBLE PRECISION EPOCH1, EPOCH2, UT, ELONG, U, V * 2Pi DOUBLE PRECISION D2PI PARAMETER ( D2PI = 6.283185307179586476925287D0 ) * Degrees to radians DOUBLE PRECISION DD2R PARAMETER ( DD2R = 1.745329251994329576923691D-2 ) * Reference epoch (J2000), JD DOUBLE PRECISION DJ0 PARAMETER ( DJ0 = 2451545D0 ) * Days per Julian millennium DOUBLE PRECISION DJM PARAMETER ( DJM = 365250D0 ) DOUBLE PRECISION T, TSOL, W, ELSUN, EMSUN, D, ELJ, ELS, : WT, W0, W1, W2, W3, W4, WF, WJ * * ===================== * Fairhead et al. model * ===================== * * 787 sets of three coefficients. * * Each set is amplitude (microseconds) * frequency (radians per Julian millennium since J2000), * phase (radians). * * Sets 1-474 are the T**0 terms, * " 475-679 " " T**1 " * " 680-764 " " T**2 " * " 765-784 " " T**3 " * " 785-787 " " T**4 " . * DOUBLE PRECISION FAIRHD(3,787) INTEGER I,J DATA ((FAIRHD(I,J),I=1,3),J= 1, 10) / : 1656.674564D-6, 6283.075849991D0, 6.240054195D0, : 22.417471D-6, 5753.384884897D0, 4.296977442D0, : 13.839792D-6, 12566.151699983D0, 6.196904410D0, : 4.770086D-6, 529.690965095D0, 0.444401603D0, : 4.676740D-6, 6069.776754553D0, 4.021195093D0, : 2.256707D-6, 213.299095438D0, 5.543113262D0, : 1.694205D-6, -3.523118349D0, 5.025132748D0, : 1.554905D-6, 77713.771467920D0, 5.198467090D0, : 1.276839D-6, 7860.419392439D0, 5.988822341D0, : 1.193379D-6, 5223.693919802D0, 3.649823730D0 / DATA ((FAIRHD(I,J),I=1,3),J= 11, 20) / : 1.115322D-6, 3930.209696220D0, 1.422745069D0, : 0.794185D-6, 11506.769769794D0, 2.322313077D0, : 0.447061D-6, 26.298319800D0, 3.615796498D0, : 0.435206D-6, -398.149003408D0, 4.349338347D0, : 0.600309D-6, 1577.343542448D0, 2.678271909D0, : 0.496817D-6, 6208.294251424D0, 5.696701824D0, : 0.486306D-6, 5884.926846583D0, 0.520007179D0, : 0.432392D-6, 74.781598567D0, 2.435898309D0, : 0.468597D-6, 6244.942814354D0, 5.866398759D0, : 0.375510D-6, 5507.553238667D0, 4.103476804D0 / DATA ((FAIRHD(I,J),I=1,3),J= 21, 30) / : 0.243085D-6, -775.522611324D0, 3.651837925D0, : 0.173435D-6, 18849.227549974D0, 6.153743485D0, : 0.230685D-6, 5856.477659115D0, 4.773852582D0, : 0.203747D-6, 12036.460734888D0, 4.333987818D0, : 0.143935D-6, -796.298006816D0, 5.957517795D0, : 0.159080D-6, 10977.078804699D0, 1.890075226D0, : 0.119979D-6, 38.133035638D0, 4.551585768D0, : 0.118971D-6, 5486.777843175D0, 1.914547226D0, : 0.116120D-6, 1059.381930189D0, 0.873504123D0, : 0.137927D-6, 11790.629088659D0, 1.135934669D0 / DATA ((FAIRHD(I,J),I=1,3),J= 31, 40) / : 0.098358D-6, 2544.314419883D0, 0.092793886D0, : 0.101868D-6, -5573.142801634D0, 5.984503847D0, : 0.080164D-6, 206.185548437D0, 2.095377709D0, : 0.079645D-6, 4694.002954708D0, 2.949233637D0, : 0.062617D-6, 20.775395492D0, 2.654394814D0, : 0.075019D-6, 2942.463423292D0, 4.980931759D0, : 0.064397D-6, 5746.271337896D0, 1.280308748D0, : 0.063814D-6, 5760.498431898D0, 4.167901731D0, : 0.048042D-6, 2146.165416475D0, 1.495846011D0, : 0.048373D-6, 155.420399434D0, 2.251573730D0 / DATA ((FAIRHD(I,J),I=1,3),J= 41, 50) / : 0.058844D-6, 426.598190876D0, 4.839650148D0, : 0.046551D-6, -0.980321068D0, 0.921573539D0, : 0.054139D-6, 17260.154654690D0, 3.411091093D0, : 0.042411D-6, 6275.962302991D0, 2.869567043D0, : 0.040184D-6, -7.113547001D0, 3.565975565D0, : 0.036564D-6, 5088.628839767D0, 3.324679049D0, : 0.040759D-6, 12352.852604545D0, 3.981496998D0, : 0.036507D-6, 801.820931124D0, 6.248866009D0, : 0.036955D-6, 3154.687084896D0, 5.071801441D0, : 0.042732D-6, 632.783739313D0, 5.720622217D0 / DATA ((FAIRHD(I,J),I=1,3),J= 51, 60) / : 0.042560D-6, 161000.685737473D0, 1.270837679D0, : 0.040480D-6, 15720.838784878D0, 2.546610123D0, : 0.028244D-6, -6286.598968340D0, 5.069663519D0, : 0.033477D-6, 6062.663207553D0, 4.144987272D0, : 0.034867D-6, 522.577418094D0, 5.210064075D0, : 0.032438D-6, 6076.890301554D0, 0.749317412D0, : 0.030215D-6, 7084.896781115D0, 3.389610345D0, : 0.029247D-6, -71430.695617928D0, 4.183178762D0, : 0.033529D-6, 9437.762934887D0, 2.404714239D0, : 0.032423D-6, 8827.390269875D0, 5.541473556D0 / DATA ((FAIRHD(I,J),I=1,3),J= 61, 70) / : 0.027567D-6, 6279.552731642D0, 5.040846034D0, : 0.029862D-6, 12139.553509107D0, 1.770181024D0, : 0.022509D-6, 10447.387839604D0, 1.460726241D0, : 0.020937D-6, 8429.241266467D0, 0.652303414D0, : 0.020322D-6, 419.484643875D0, 3.735430632D0, : 0.024816D-6, -1194.447010225D0, 1.087136918D0, : 0.025196D-6, 1748.016413067D0, 2.901883301D0, : 0.021691D-6, 14143.495242431D0, 5.952658009D0, : 0.017673D-6, 6812.766815086D0, 3.186129845D0, : 0.022567D-6, 6133.512652857D0, 3.307984806D0 / DATA ((FAIRHD(I,J),I=1,3),J= 71, 80) / : 0.016155D-6, 10213.285546211D0, 1.331103168D0, : 0.014751D-6, 1349.867409659D0, 4.308933301D0, : 0.015949D-6, -220.412642439D0, 4.005298270D0, : 0.015974D-6, -2352.866153772D0, 6.145309371D0, : 0.014223D-6, 17789.845619785D0, 2.104551349D0, : 0.017806D-6, 73.297125859D0, 3.475975097D0, : 0.013671D-6, -536.804512095D0, 5.971672571D0, : 0.011942D-6, 8031.092263058D0, 2.053414715D0, : 0.014318D-6, 16730.463689596D0, 3.016058075D0, : 0.012462D-6, 103.092774219D0, 1.737438797D0 / DATA ((FAIRHD(I,J),I=1,3),J= 81, 90) / : 0.010962D-6, 3.590428652D0, 2.196567739D0, : 0.015078D-6, 19651.048481098D0, 3.969480770D0, : 0.010396D-6, 951.718406251D0, 5.717799605D0, : 0.011707D-6, -4705.732307544D0, 2.654125618D0, : 0.010453D-6, 5863.591206116D0, 1.913704550D0, : 0.012420D-6, 4690.479836359D0, 4.734090399D0, : 0.011847D-6, 5643.178563677D0, 5.489005403D0, : 0.008610D-6, 3340.612426700D0, 3.661698944D0, : 0.011622D-6, 5120.601145584D0, 4.863931876D0, : 0.010825D-6, 553.569402842D0, 0.842715011D0 / DATA ((FAIRHD(I,J),I=1,3),J= 91,100) / : 0.008666D-6, -135.065080035D0, 3.293406547D0, : 0.009963D-6, 149.563197135D0, 4.870690598D0, : 0.009858D-6, 6309.374169791D0, 1.061816410D0, : 0.007959D-6, 316.391869657D0, 2.465042647D0, : 0.010099D-6, 283.859318865D0, 1.942176992D0, : 0.007147D-6, -242.728603974D0, 3.661486981D0, : 0.007505D-6, 5230.807466803D0, 4.920937029D0, : 0.008323D-6, 11769.853693166D0, 1.229392026D0, : 0.007490D-6, -6256.777530192D0, 3.658444681D0, : 0.009370D-6, 149854.400134205D0, 0.673880395D0 / DATA ((FAIRHD(I,J),I=1,3),J=101,110) / : 0.007117D-6, 38.027672636D0, 5.294249518D0, : 0.007857D-6, 12168.002696575D0, 0.525733528D0, : 0.007019D-6, 6206.809778716D0, 0.837688810D0, : 0.006056D-6, 955.599741609D0, 4.194535082D0, : 0.008107D-6, 13367.972631107D0, 3.793235253D0, : 0.006731D-6, 5650.292110678D0, 5.639906583D0, : 0.007332D-6, 36.648562930D0, 0.114858677D0, : 0.006366D-6, 4164.311989613D0, 2.262081818D0, : 0.006858D-6, 5216.580372801D0, 0.642063318D0, : 0.006919D-6, 6681.224853400D0, 6.018501522D0 / DATA ((FAIRHD(I,J),I=1,3),J=111,120) / : 0.006826D-6, 7632.943259650D0, 3.458654112D0, : 0.005308D-6, -1592.596013633D0, 2.500382359D0, : 0.005096D-6, 11371.704689758D0, 2.547107806D0, : 0.004841D-6, 5333.900241022D0, 0.437078094D0, : 0.005582D-6, 5966.683980335D0, 2.246174308D0, : 0.006304D-6, 11926.254413669D0, 2.512929171D0, : 0.006603D-6, 23581.258177318D0, 5.393136889D0, : 0.005123D-6, -1.484472708D0, 2.999641028D0, : 0.004648D-6, 1589.072895284D0, 1.275847090D0, : 0.005119D-6, 6438.496249426D0, 1.486539246D0 / DATA ((FAIRHD(I,J),I=1,3),J=121,130) / : 0.004521D-6, 4292.330832950D0, 6.140635794D0, : 0.005680D-6, 23013.539539587D0, 4.557814849D0, : 0.005488D-6, -3.455808046D0, 0.090675389D0, : 0.004193D-6, 7234.794256242D0, 4.869091389D0, : 0.003742D-6, 7238.675591600D0, 4.691976180D0, : 0.004148D-6, -110.206321219D0, 3.016173439D0, : 0.004553D-6, 11499.656222793D0, 5.554998314D0, : 0.004892D-6, 5436.993015240D0, 1.475415597D0, : 0.004044D-6, 4732.030627343D0, 1.398784824D0, : 0.004164D-6, 12491.370101415D0, 5.650931916D0 / DATA ((FAIRHD(I,J),I=1,3),J=131,140) / : 0.004349D-6, 11513.883316794D0, 2.181745369D0, : 0.003919D-6, 12528.018664345D0, 5.823319737D0, : 0.003129D-6, 6836.645252834D0, 0.003844094D0, : 0.004080D-6, -7058.598461315D0, 3.690360123D0, : 0.003270D-6, 76.266071276D0, 1.517189902D0, : 0.002954D-6, 6283.143160294D0, 4.447203799D0, : 0.002872D-6, 28.449187468D0, 1.158692983D0, : 0.002881D-6, 735.876513532D0, 0.349250250D0, : 0.003279D-6, 5849.364112115D0, 4.893384368D0, : 0.003625D-6, 6209.778724132D0, 1.473760578D0 / DATA ((FAIRHD(I,J),I=1,3),J=141,150) / : 0.003074D-6, 949.175608970D0, 5.185878737D0, : 0.002775D-6, 9917.696874510D0, 1.030026325D0, : 0.002646D-6, 10973.555686350D0, 3.918259169D0, : 0.002575D-6, 25132.303399966D0, 6.109659023D0, : 0.003500D-6, 263.083923373D0, 1.892100742D0, : 0.002740D-6, 18319.536584880D0, 4.320519510D0, : 0.002464D-6, 202.253395174D0, 4.698203059D0, : 0.002409D-6, 2.542797281D0, 5.325009315D0, : 0.003354D-6, -90955.551694697D0, 1.942656623D0, : 0.002296D-6, 6496.374945429D0, 5.061810696D0 / DATA ((FAIRHD(I,J),I=1,3),J=151,160) / : 0.003002D-6, 6172.869528772D0, 2.797822767D0, : 0.003202D-6, 27511.467873537D0, 0.531673101D0, : 0.002954D-6, -6283.008539689D0, 4.533471191D0, : 0.002353D-6, 639.897286314D0, 3.734548088D0, : 0.002401D-6, 16200.772724501D0, 2.605547070D0, : 0.003053D-6, 233141.314403759D0, 3.029030662D0, : 0.003024D-6, 83286.914269554D0, 2.355556099D0, : 0.002863D-6, 17298.182327326D0, 5.240963796D0, : 0.002103D-6, -7079.373856808D0, 5.756641637D0, : 0.002303D-6, 83996.847317911D0, 2.013686814D0 / DATA ((FAIRHD(I,J),I=1,3),J=161,170) / : 0.002303D-6, 18073.704938650D0, 1.089100410D0, : 0.002381D-6, 63.735898303D0, 0.759188178D0, : 0.002493D-6, 6386.168624210D0, 0.645026535D0, : 0.002366D-6, 3.932153263D0, 6.215885448D0, : 0.002169D-6, 11015.106477335D0, 4.845297676D0, : 0.002397D-6, 6243.458341645D0, 3.809290043D0, : 0.002183D-6, 1162.474704408D0, 6.179611691D0, : 0.002353D-6, 6246.427287062D0, 4.781719760D0, : 0.002199D-6, -245.831646229D0, 5.956152284D0, : 0.001729D-6, 3894.181829542D0, 1.264976635D0 / DATA ((FAIRHD(I,J),I=1,3),J=171,180) / : 0.001896D-6, -3128.388765096D0, 4.914231596D0, : 0.002085D-6, 35.164090221D0, 1.405158503D0, : 0.002024D-6, 14712.317116458D0, 2.752035928D0, : 0.001737D-6, 6290.189396992D0, 5.280820144D0, : 0.002229D-6, 491.557929457D0, 1.571007057D0, : 0.001602D-6, 14314.168113050D0, 4.203664806D0, : 0.002186D-6, 454.909366527D0, 1.402101526D0, : 0.001897D-6, 22483.848574493D0, 4.167932508D0, : 0.001825D-6, -3738.761430108D0, 0.545828785D0, : 0.001894D-6, 1052.268383188D0, 5.817167450D0 / DATA ((FAIRHD(I,J),I=1,3),J=181,190) / : 0.001421D-6, 20.355319399D0, 2.419886601D0, : 0.001408D-6, 10984.192351700D0, 2.732084787D0, : 0.001847D-6, 10873.986030480D0, 2.903477885D0, : 0.001391D-6, -8635.942003763D0, 0.593891500D0, : 0.001388D-6, -7.046236698D0, 1.166145902D0, : 0.001810D-6, -88860.057071188D0, 0.487355242D0, : 0.001288D-6, -1990.745017041D0, 3.913022880D0, : 0.001297D-6, 23543.230504682D0, 3.063805171D0, : 0.001335D-6, -266.607041722D0, 3.995764039D0, : 0.001376D-6, 10969.965257698D0, 5.152914309D0 / DATA ((FAIRHD(I,J),I=1,3),J=191,200) / : 0.001745D-6, 244287.600007027D0, 3.626395673D0, : 0.001649D-6, 31441.677569757D0, 1.952049260D0, : 0.001416D-6, 9225.539273283D0, 4.996408389D0, : 0.001238D-6, 4804.209275927D0, 5.503379738D0, : 0.001472D-6, 4590.910180489D0, 4.164913291D0, : 0.001169D-6, 6040.347246017D0, 5.841719038D0, : 0.001039D-6, 5540.085789459D0, 2.769753519D0, : 0.001004D-6, -170.672870619D0, 0.755008103D0, : 0.001284D-6, 10575.406682942D0, 5.306538209D0, : 0.001278D-6, 71.812653151D0, 4.713486491D0 / DATA ((FAIRHD(I,J),I=1,3),J=201,210) / : 0.001321D-6, 18209.330263660D0, 2.624866359D0, : 0.001297D-6, 21228.392023546D0, 0.382603541D0, : 0.000954D-6, 6282.095528923D0, 0.882213514D0, : 0.001145D-6, 6058.731054289D0, 1.169483931D0, : 0.000979D-6, 5547.199336460D0, 5.448375984D0, : 0.000987D-6, -6262.300454499D0, 2.656486959D0, : 0.001070D-6, -154717.609887482D0, 1.827624012D0, : 0.000991D-6, 4701.116501708D0, 4.387001801D0, : 0.001155D-6, -14.227094002D0, 3.042700750D0, : 0.001176D-6, 277.034993741D0, 3.335519004D0 / DATA ((FAIRHD(I,J),I=1,3),J=211,220) / : 0.000890D-6, 13916.019109642D0, 5.601498297D0, : 0.000884D-6, -1551.045222648D0, 1.088831705D0, : 0.000876D-6, 5017.508371365D0, 3.969902609D0, : 0.000806D-6, 15110.466119866D0, 5.142876744D0, : 0.000773D-6, -4136.910433516D0, 0.022067765D0, : 0.001077D-6, 175.166059800D0, 1.844913056D0, : 0.000954D-6, -6284.056171060D0, 0.968480906D0, : 0.000737D-6, 5326.786694021D0, 4.923831588D0, : 0.000845D-6, -433.711737877D0, 4.749245231D0, : 0.000819D-6, 8662.240323563D0, 5.991247817D0 / DATA ((FAIRHD(I,J),I=1,3),J=221,230) / : 0.000852D-6, 199.072001436D0, 2.189604979D0, : 0.000723D-6, 17256.631536341D0, 6.068719637D0, : 0.000940D-6, 6037.244203762D0, 6.197428148D0, : 0.000885D-6, 11712.955318231D0, 3.280414875D0, : 0.000706D-6, 12559.038152982D0, 2.824848947D0, : 0.000732D-6, 2379.164473572D0, 2.501813417D0, : 0.000764D-6, -6127.655450557D0, 2.236346329D0, : 0.000908D-6, 131.541961686D0, 2.521257490D0, : 0.000907D-6, 35371.887265976D0, 3.370195967D0, : 0.000673D-6, 1066.495477190D0, 3.876512374D0 / DATA ((FAIRHD(I,J),I=1,3),J=231,240) / : 0.000814D-6, 17654.780539750D0, 4.627122566D0, : 0.000630D-6, 36.027866677D0, 0.156368499D0, : 0.000798D-6, 515.463871093D0, 5.151962502D0, : 0.000798D-6, 148.078724426D0, 5.909225055D0, : 0.000806D-6, 309.278322656D0, 6.054064447D0, : 0.000607D-6, -39.617508346D0, 2.839021623D0, : 0.000601D-6, 412.371096874D0, 3.984225404D0, : 0.000646D-6, 11403.676995575D0, 3.852959484D0, : 0.000704D-6, 13521.751441591D0, 2.300991267D0, : 0.000603D-6, -65147.619767937D0, 4.140083146D0 / DATA ((FAIRHD(I,J),I=1,3),J=241,250) / : 0.000609D-6, 10177.257679534D0, 0.437122327D0, : 0.000631D-6, 5767.611978898D0, 4.026532329D0, : 0.000576D-6, 11087.285125918D0, 4.760293101D0, : 0.000674D-6, 14945.316173554D0, 6.270510511D0, : 0.000726D-6, 5429.879468239D0, 6.039606892D0, : 0.000710D-6, 28766.924424484D0, 5.672617711D0, : 0.000647D-6, 11856.218651625D0, 3.397132627D0, : 0.000678D-6, -5481.254918868D0, 6.249666675D0, : 0.000618D-6, 22003.914634870D0, 2.466427018D0, : 0.000738D-6, 6134.997125565D0, 2.242668890D0 / DATA ((FAIRHD(I,J),I=1,3),J=251,260) / : 0.000660D-6, 625.670192312D0, 5.864091907D0, : 0.000694D-6, 3496.032826134D0, 2.668309141D0, : 0.000531D-6, 6489.261398429D0, 1.681888780D0, : 0.000611D-6, -143571.324284214D0, 2.424978312D0, : 0.000575D-6, 12043.574281889D0, 4.216492400D0, : 0.000553D-6, 12416.588502848D0, 4.772158039D0, : 0.000689D-6, 4686.889407707D0, 6.224271088D0, : 0.000495D-6, 7342.457780181D0, 3.817285811D0, : 0.000567D-6, 3634.621024518D0, 1.649264690D0, : 0.000515D-6, 18635.928454536D0, 3.945345892D0 / DATA ((FAIRHD(I,J),I=1,3),J=261,270) / : 0.000486D-6, -323.505416657D0, 4.061673868D0, : 0.000662D-6, 25158.601719765D0, 1.794058369D0, : 0.000509D-6, 846.082834751D0, 3.053874588D0, : 0.000472D-6, -12569.674818332D0, 5.112133338D0, : 0.000461D-6, 6179.983075773D0, 0.513669325D0, : 0.000641D-6, 83467.156352816D0, 3.210727723D0, : 0.000520D-6, 10344.295065386D0, 2.445597761D0, : 0.000493D-6, 18422.629359098D0, 1.676939306D0, : 0.000478D-6, 1265.567478626D0, 5.487314569D0, : 0.000472D-6, -18.159247265D0, 1.999707589D0 / DATA ((FAIRHD(I,J),I=1,3),J=271,280) / : 0.000559D-6, 11190.377900137D0, 5.783236356D0, : 0.000494D-6, 9623.688276691D0, 3.022645053D0, : 0.000463D-6, 5739.157790895D0, 1.411223013D0, : 0.000432D-6, 16858.482532933D0, 1.179256434D0, : 0.000574D-6, 72140.628666286D0, 1.758191830D0, : 0.000484D-6, 17267.268201691D0, 3.290589143D0, : 0.000550D-6, 4907.302050146D0, 0.864024298D0, : 0.000399D-6, 14.977853527D0, 2.094441910D0, : 0.000491D-6, 224.344795702D0, 0.878372791D0, : 0.000432D-6, 20426.571092422D0, 6.003829241D0 / DATA ((FAIRHD(I,J),I=1,3),J=281,290) / : 0.000481D-6, 5749.452731634D0, 4.309591964D0, : 0.000480D-6, 5757.317038160D0, 1.142348571D0, : 0.000485D-6, 6702.560493867D0, 0.210580917D0, : 0.000426D-6, 6055.549660552D0, 4.274476529D0, : 0.000480D-6, 5959.570433334D0, 5.031351030D0, : 0.000466D-6, 12562.628581634D0, 4.959581597D0, : 0.000520D-6, 39302.096962196D0, 4.788002889D0, : 0.000458D-6, 12132.439962106D0, 1.880103788D0, : 0.000470D-6, 12029.347187887D0, 1.405611197D0, : 0.000416D-6, -7477.522860216D0, 1.082356330D0 / DATA ((FAIRHD(I,J),I=1,3),J=291,300) / : 0.000449D-6, 11609.862544012D0, 4.179989585D0, : 0.000465D-6, 17253.041107690D0, 0.353496295D0, : 0.000362D-6, -4535.059436924D0, 1.583849576D0, : 0.000383D-6, 21954.157609398D0, 3.747376371D0, : 0.000389D-6, 17.252277143D0, 1.395753179D0, : 0.000331D-6, 18052.929543158D0, 0.566790582D0, : 0.000430D-6, 13517.870106233D0, 0.685827538D0, : 0.000368D-6, -5756.908003246D0, 0.731374317D0, : 0.000330D-6, 10557.594160824D0, 3.710043680D0, : 0.000332D-6, 20199.094959633D0, 1.652901407D0 / DATA ((FAIRHD(I,J),I=1,3),J=301,310) / : 0.000384D-6, 11933.367960670D0, 5.827781531D0, : 0.000387D-6, 10454.501386605D0, 2.541182564D0, : 0.000325D-6, 15671.081759407D0, 2.178850542D0, : 0.000318D-6, 138.517496871D0, 2.253253037D0, : 0.000305D-6, 9388.005909415D0, 0.578340206D0, : 0.000352D-6, 5749.861766548D0, 3.000297967D0, : 0.000311D-6, 6915.859589305D0, 1.693574249D0, : 0.000297D-6, 24072.921469776D0, 1.997249392D0, : 0.000363D-6, -640.877607382D0, 5.071820966D0, : 0.000323D-6, 12592.450019783D0, 1.072262823D0 / DATA ((FAIRHD(I,J),I=1,3),J=311,320) / : 0.000341D-6, 12146.667056108D0, 4.700657997D0, : 0.000290D-6, 9779.108676125D0, 1.812320441D0, : 0.000342D-6, 6132.028180148D0, 4.322238614D0, : 0.000329D-6, 6268.848755990D0, 3.033827743D0, : 0.000374D-6, 17996.031168222D0, 3.388716544D0, : 0.000285D-6, -533.214083444D0, 4.687313233D0, : 0.000338D-6, 6065.844601290D0, 0.877776108D0, : 0.000276D-6, 24.298513841D0, 0.770299429D0, : 0.000336D-6, -2388.894020449D0, 5.353796034D0, : 0.000290D-6, 3097.883822726D0, 4.075291557D0 / DATA ((FAIRHD(I,J),I=1,3),J=321,330) / : 0.000318D-6, 709.933048357D0, 5.941207518D0, : 0.000271D-6, 13095.842665077D0, 3.208912203D0, : 0.000331D-6, 6073.708907816D0, 4.007881169D0, : 0.000292D-6, 742.990060533D0, 2.714333592D0, : 0.000362D-6, 29088.811415985D0, 3.215977013D0, : 0.000280D-6, 12359.966151546D0, 0.710872502D0, : 0.000267D-6, 10440.274292604D0, 4.730108488D0, : 0.000262D-6, 838.969287750D0, 1.327720272D0, : 0.000250D-6, 16496.361396202D0, 0.898769761D0, : 0.000325D-6, 20597.243963041D0, 0.180044365D0 / DATA ((FAIRHD(I,J),I=1,3),J=331,340) / : 0.000268D-6, 6148.010769956D0, 5.152666276D0, : 0.000284D-6, 5636.065016677D0, 5.655385808D0, : 0.000301D-6, 6080.822454817D0, 2.135396205D0, : 0.000294D-6, -377.373607916D0, 3.708784168D0, : 0.000236D-6, 2118.763860378D0, 1.733578756D0, : 0.000234D-6, 5867.523359379D0, 5.575209112D0, : 0.000268D-6, -226858.238553767D0, 0.069432392D0, : 0.000265D-6, 167283.761587465D0, 4.369302826D0, : 0.000280D-6, 28237.233459389D0, 5.304829118D0, : 0.000292D-6, 12345.739057544D0, 4.096094132D0 / DATA ((FAIRHD(I,J),I=1,3),J=341,350) / : 0.000223D-6, 19800.945956225D0, 3.069327406D0, : 0.000301D-6, 43232.306658416D0, 6.205311188D0, : 0.000264D-6, 18875.525869774D0, 1.417263408D0, : 0.000304D-6, -1823.175188677D0, 3.409035232D0, : 0.000301D-6, 109.945688789D0, 0.510922054D0, : 0.000260D-6, 813.550283960D0, 2.389438934D0, : 0.000299D-6, 316428.228673312D0, 5.384595078D0, : 0.000211D-6, 5756.566278634D0, 3.789392838D0, : 0.000209D-6, 5750.203491159D0, 1.661943545D0, : 0.000240D-6, 12489.885628707D0, 5.684549045D0 / DATA ((FAIRHD(I,J),I=1,3),J=351,360) / : 0.000216D-6, 6303.851245484D0, 3.862942261D0, : 0.000203D-6, 1581.959348283D0, 5.549853589D0, : 0.000200D-6, 5642.198242609D0, 1.016115785D0, : 0.000197D-6, -70.849445304D0, 4.690702525D0, : 0.000227D-6, 6287.008003254D0, 2.911891613D0, : 0.000197D-6, 533.623118358D0, 1.048982898D0, : 0.000205D-6, -6279.485421340D0, 1.829362730D0, : 0.000209D-6, -10988.808157535D0, 2.636140084D0, : 0.000208D-6, -227.526189440D0, 4.127883842D0, : 0.000191D-6, 415.552490612D0, 4.401165650D0 / DATA ((FAIRHD(I,J),I=1,3),J=361,370) / : 0.000190D-6, 29296.615389579D0, 4.175658539D0, : 0.000264D-6, 66567.485864652D0, 4.601102551D0, : 0.000256D-6, -3646.350377354D0, 0.506364778D0, : 0.000188D-6, 13119.721102825D0, 2.032195842D0, : 0.000185D-6, -209.366942175D0, 4.694756586D0, : 0.000198D-6, 25934.124331089D0, 3.832703118D0, : 0.000195D-6, 4061.219215394D0, 3.308463427D0, : 0.000234D-6, 5113.487598583D0, 1.716090661D0, : 0.000188D-6, 1478.866574064D0, 5.686865780D0, : 0.000222D-6, 11823.161639450D0, 1.942386641D0 / DATA ((FAIRHD(I,J),I=1,3),J=371,380) / : 0.000181D-6, 10770.893256262D0, 1.999482059D0, : 0.000171D-6, 6546.159773364D0, 1.182807992D0, : 0.000206D-6, 70.328180442D0, 5.934076062D0, : 0.000169D-6, 20995.392966449D0, 2.169080622D0, : 0.000191D-6, 10660.686935042D0, 5.405515999D0, : 0.000228D-6, 33019.021112205D0, 4.656985514D0, : 0.000184D-6, -4933.208440333D0, 3.327476868D0, : 0.000220D-6, -135.625325010D0, 1.765430262D0, : 0.000166D-6, 23141.558382925D0, 3.454132746D0, : 0.000191D-6, 6144.558353121D0, 5.020393445D0 / DATA ((FAIRHD(I,J),I=1,3),J=381,390) / : 0.000180D-6, 6084.003848555D0, 0.602182191D0, : 0.000163D-6, 17782.732072784D0, 4.960593133D0, : 0.000225D-6, 16460.333529525D0, 2.596451817D0, : 0.000222D-6, 5905.702242076D0, 3.731990323D0, : 0.000204D-6, 227.476132789D0, 5.636192701D0, : 0.000159D-6, 16737.577236597D0, 3.600691544D0, : 0.000200D-6, 6805.653268085D0, 0.868220961D0, : 0.000187D-6, 11919.140866668D0, 2.629456641D0, : 0.000161D-6, 127.471796607D0, 2.862574720D0, : 0.000205D-6, 6286.666278643D0, 1.742882331D0 / DATA ((FAIRHD(I,J),I=1,3),J=391,400) / : 0.000189D-6, 153.778810485D0, 4.812372643D0, : 0.000168D-6, 16723.350142595D0, 0.027860588D0, : 0.000149D-6, 11720.068865232D0, 0.659721876D0, : 0.000189D-6, 5237.921013804D0, 5.245313000D0, : 0.000143D-6, 6709.674040867D0, 4.317625647D0, : 0.000146D-6, 4487.817406270D0, 4.815297007D0, : 0.000144D-6, -664.756045130D0, 5.381366880D0, : 0.000175D-6, 5127.714692584D0, 4.728443327D0, : 0.000162D-6, 6254.626662524D0, 1.435132069D0, : 0.000187D-6, 47162.516354635D0, 1.354371923D0 / DATA ((FAIRHD(I,J),I=1,3),J=401,410) / : 0.000146D-6, 11080.171578918D0, 3.369695406D0, : 0.000180D-6, -348.924420448D0, 2.490902145D0, : 0.000148D-6, 151.047669843D0, 3.799109588D0, : 0.000157D-6, 6197.248551160D0, 1.284375887D0, : 0.000167D-6, 146.594251718D0, 0.759969109D0, : 0.000133D-6, -5331.357443741D0, 5.409701889D0, : 0.000154D-6, 95.979227218D0, 3.366890614D0, : 0.000148D-6, -6418.140930027D0, 3.384104996D0, : 0.000128D-6, -6525.804453965D0, 3.803419985D0, : 0.000130D-6, 11293.470674356D0, 0.939039445D0 / DATA ((FAIRHD(I,J),I=1,3),J=411,420) / : 0.000152D-6, -5729.506447149D0, 0.734117523D0, : 0.000138D-6, 210.117701700D0, 2.564216078D0, : 0.000123D-6, 6066.595360816D0, 4.517099537D0, : 0.000140D-6, 18451.078546566D0, 0.642049130D0, : 0.000126D-6, 11300.584221356D0, 3.485280663D0, : 0.000119D-6, 10027.903195729D0, 3.217431161D0, : 0.000151D-6, 4274.518310832D0, 4.404359108D0, : 0.000117D-6, 6072.958148291D0, 0.366324650D0, : 0.000165D-6, -7668.637425143D0, 4.298212528D0, : 0.000117D-6, -6245.048177356D0, 5.379518958D0 / DATA ((FAIRHD(I,J),I=1,3),J=421,430) / : 0.000130D-6, -5888.449964932D0, 4.527681115D0, : 0.000121D-6, -543.918059096D0, 6.109429504D0, : 0.000162D-6, 9683.594581116D0, 5.720092446D0, : 0.000141D-6, 6219.339951688D0, 0.679068671D0, : 0.000118D-6, 22743.409379516D0, 4.881123092D0, : 0.000129D-6, 1692.165669502D0, 0.351407289D0, : 0.000126D-6, 5657.405657679D0, 5.146592349D0, : 0.000114D-6, 728.762966531D0, 0.520791814D0, : 0.000120D-6, 52.596639600D0, 0.948516300D0, : 0.000115D-6, 65.220371012D0, 3.504914846D0 / DATA ((FAIRHD(I,J),I=1,3),J=431,440) / : 0.000126D-6, 5881.403728234D0, 5.577502482D0, : 0.000158D-6, 163096.180360983D0, 2.957128968D0, : 0.000134D-6, 12341.806904281D0, 2.598576764D0, : 0.000151D-6, 16627.370915377D0, 3.985702050D0, : 0.000109D-6, 1368.660252845D0, 0.014730471D0, : 0.000131D-6, 6211.263196841D0, 0.085077024D0, : 0.000146D-6, 5792.741760812D0, 0.708426604D0, : 0.000146D-6, -77.750543984D0, 3.121576600D0, : 0.000107D-6, 5341.013788022D0, 0.288231904D0, : 0.000138D-6, 6281.591377283D0, 2.797450317D0 / DATA ((FAIRHD(I,J),I=1,3),J=441,450) / : 0.000113D-6, -6277.552925684D0, 2.788904128D0, : 0.000115D-6, -525.758811831D0, 5.895222200D0, : 0.000138D-6, 6016.468808270D0, 6.096188999D0, : 0.000139D-6, 23539.707386333D0, 2.028195445D0, : 0.000146D-6, -4176.041342449D0, 4.660008502D0, : 0.000107D-6, 16062.184526117D0, 4.066520001D0, : 0.000142D-6, 83783.548222473D0, 2.936315115D0, : 0.000128D-6, 9380.959672717D0, 3.223844306D0, : 0.000135D-6, 6205.325306007D0, 1.638054048D0, : 0.000101D-6, 2699.734819318D0, 5.481603249D0 / DATA ((FAIRHD(I,J),I=1,3),J=451,460) / : 0.000104D-6, -568.821874027D0, 2.205734493D0, : 0.000103D-6, 6321.103522627D0, 2.440421099D0, : 0.000119D-6, 6321.208885629D0, 2.547496264D0, : 0.000138D-6, 1975.492545856D0, 2.314608466D0, : 0.000121D-6, 137.033024162D0, 4.539108237D0, : 0.000123D-6, 19402.796952817D0, 4.538074405D0, : 0.000119D-6, 22805.735565994D0, 2.869040566D0, : 0.000133D-6, 64471.991241142D0, 6.056405489D0, : 0.000129D-6, -85.827298831D0, 2.540635083D0, : 0.000131D-6, 13613.804277336D0, 4.005732868D0 / DATA ((FAIRHD(I,J),I=1,3),J=461,470) / : 0.000104D-6, 9814.604100291D0, 1.959967212D0, : 0.000112D-6, 16097.679950283D0, 3.589026260D0, : 0.000123D-6, 2107.034507542D0, 1.728627253D0, : 0.000121D-6, 36949.230808424D0, 6.072332087D0, : 0.000108D-6, -12539.853380183D0, 3.716133846D0, : 0.000113D-6, -7875.671863624D0, 2.725771122D0, : 0.000109D-6, 4171.425536614D0, 4.033338079D0, : 0.000101D-6, 6247.911759770D0, 3.441347021D0, : 0.000113D-6, 7330.728427345D0, 0.656372122D0, : 0.000113D-6, 51092.726050855D0, 2.791483066D0 / DATA ((FAIRHD(I,J),I=1,3),J=471,480) / : 0.000106D-6, 5621.842923210D0, 1.815323326D0, : 0.000101D-6, 111.430161497D0, 5.711033677D0, : 0.000103D-6, 909.818733055D0, 2.812745443D0, : 0.000101D-6, 1790.642637886D0, 1.965746028D0, * T : 102.156724D-6, 6283.075849991D0, 4.249032005D0, : 1.706807D-6, 12566.151699983D0, 4.205904248D0, : 0.269668D-6, 213.299095438D0, 3.400290479D0, : 0.265919D-6, 529.690965095D0, 5.836047367D0, : 0.210568D-6, -3.523118349D0, 6.262738348D0, : 0.077996D-6, 5223.693919802D0, 4.670344204D0 / DATA ((FAIRHD(I,J),I=1,3),J=481,490) / : 0.054764D-6, 1577.343542448D0, 4.534800170D0, : 0.059146D-6, 26.298319800D0, 1.083044735D0, : 0.034420D-6, -398.149003408D0, 5.980077351D0, : 0.032088D-6, 18849.227549974D0, 4.162913471D0, : 0.033595D-6, 5507.553238667D0, 5.980162321D0, : 0.029198D-6, 5856.477659115D0, 0.623811863D0, : 0.027764D-6, 155.420399434D0, 3.745318113D0, : 0.025190D-6, 5746.271337896D0, 2.980330535D0, : 0.022997D-6, -796.298006816D0, 1.174411803D0, : 0.024976D-6, 5760.498431898D0, 2.467913690D0 / DATA ((FAIRHD(I,J),I=1,3),J=491,500) / : 0.021774D-6, 206.185548437D0, 3.854787540D0, : 0.017925D-6, -775.522611324D0, 1.092065955D0, : 0.013794D-6, 426.598190876D0, 2.699831988D0, : 0.013276D-6, 6062.663207553D0, 5.845801920D0, : 0.011774D-6, 12036.460734888D0, 2.292832062D0, : 0.012869D-6, 6076.890301554D0, 5.333425680D0, : 0.012152D-6, 1059.381930189D0, 6.222874454D0, : 0.011081D-6, -7.113547001D0, 5.154724984D0, : 0.010143D-6, 4694.002954708D0, 4.044013795D0, : 0.009357D-6, 5486.777843175D0, 3.416081409D0 / DATA ((FAIRHD(I,J),I=1,3),J=501,510) / : 0.010084D-6, 522.577418094D0, 0.749320262D0, : 0.008587D-6, 10977.078804699D0, 2.777152598D0, : 0.008628D-6, 6275.962302991D0, 4.562060226D0, : 0.008158D-6, -220.412642439D0, 5.806891533D0, : 0.007746D-6, 2544.314419883D0, 1.603197066D0, : 0.007670D-6, 2146.165416475D0, 3.000200440D0, : 0.007098D-6, 74.781598567D0, 0.443725817D0, : 0.006180D-6, -536.804512095D0, 1.302642751D0, : 0.005818D-6, 5088.628839767D0, 4.827723531D0, : 0.004945D-6, -6286.598968340D0, 0.268305170D0 / DATA ((FAIRHD(I,J),I=1,3),J=511,520) / : 0.004774D-6, 1349.867409659D0, 5.808636673D0, : 0.004687D-6, -242.728603974D0, 5.154890570D0, : 0.006089D-6, 1748.016413067D0, 4.403765209D0, : 0.005975D-6, -1194.447010225D0, 2.583472591D0, : 0.004229D-6, 951.718406251D0, 0.931172179D0, : 0.005264D-6, 553.569402842D0, 2.336107252D0, : 0.003049D-6, 5643.178563677D0, 1.362634430D0, : 0.002974D-6, 6812.766815086D0, 1.583012668D0, : 0.003403D-6, -2352.866153772D0, 2.552189886D0, : 0.003030D-6, 419.484643875D0, 5.286473844D0 / DATA ((FAIRHD(I,J),I=1,3),J=521,530) / : 0.003210D-6, -7.046236698D0, 1.863796539D0, : 0.003058D-6, 9437.762934887D0, 4.226420633D0, : 0.002589D-6, 12352.852604545D0, 1.991935820D0, : 0.002927D-6, 5216.580372801D0, 2.319951253D0, : 0.002425D-6, 5230.807466803D0, 3.084752833D0, : 0.002656D-6, 3154.687084896D0, 2.487447866D0, : 0.002445D-6, 10447.387839604D0, 2.347139160D0, : 0.002990D-6, 4690.479836359D0, 6.235872050D0, : 0.002890D-6, 5863.591206116D0, 0.095197563D0, : 0.002498D-6, 6438.496249426D0, 2.994779800D0 / DATA ((FAIRHD(I,J),I=1,3),J=531,540) / : 0.001889D-6, 8031.092263058D0, 3.569003717D0, : 0.002567D-6, 801.820931124D0, 3.425611498D0, : 0.001803D-6, -71430.695617928D0, 2.192295512D0, : 0.001782D-6, 3.932153263D0, 5.180433689D0, : 0.001694D-6, -4705.732307544D0, 4.641779174D0, : 0.001704D-6, -1592.596013633D0, 3.997097652D0, : 0.001735D-6, 5849.364112115D0, 0.417558428D0, : 0.001643D-6, 8429.241266467D0, 2.180619584D0, : 0.001680D-6, 38.133035638D0, 4.164529426D0, : 0.002045D-6, 7084.896781115D0, 0.526323854D0 / DATA ((FAIRHD(I,J),I=1,3),J=541,550) / : 0.001458D-6, 4292.330832950D0, 1.356098141D0, : 0.001437D-6, 20.355319399D0, 3.895439360D0, : 0.001738D-6, 6279.552731642D0, 0.087484036D0, : 0.001367D-6, 14143.495242431D0, 3.987576591D0, : 0.001344D-6, 7234.794256242D0, 0.090454338D0, : 0.001438D-6, 11499.656222793D0, 0.974387904D0, : 0.001257D-6, 6836.645252834D0, 1.509069366D0, : 0.001358D-6, 11513.883316794D0, 0.495572260D0, : 0.001628D-6, 7632.943259650D0, 4.968445721D0, : 0.001169D-6, 103.092774219D0, 2.838496795D0 / DATA ((FAIRHD(I,J),I=1,3),J=551,560) / : 0.001162D-6, 4164.311989613D0, 3.408387778D0, : 0.001092D-6, 6069.776754553D0, 3.617942651D0, : 0.001008D-6, 17789.845619785D0, 0.286350174D0, : 0.001008D-6, 639.897286314D0, 1.610762073D0, : 0.000918D-6, 10213.285546211D0, 5.532798067D0, : 0.001011D-6, -6256.777530192D0, 0.661826484D0, : 0.000753D-6, 16730.463689596D0, 3.905030235D0, : 0.000737D-6, 11926.254413669D0, 4.641956361D0, : 0.000694D-6, 3340.612426700D0, 2.111120332D0, : 0.000701D-6, 3894.181829542D0, 2.760823491D0 / DATA ((FAIRHD(I,J),I=1,3),J=561,570) / : 0.000689D-6, -135.065080035D0, 4.768800780D0, : 0.000700D-6, 13367.972631107D0, 5.760439898D0, : 0.000664D-6, 6040.347246017D0, 1.051215840D0, : 0.000654D-6, 5650.292110678D0, 4.911332503D0, : 0.000788D-6, 6681.224853400D0, 4.699648011D0, : 0.000628D-6, 5333.900241022D0, 5.024608847D0, : 0.000755D-6, -110.206321219D0, 4.370971253D0, : 0.000628D-6, 6290.189396992D0, 3.660478857D0, : 0.000635D-6, 25132.303399966D0, 4.121051532D0, : 0.000534D-6, 5966.683980335D0, 1.173284524D0 / DATA ((FAIRHD(I,J),I=1,3),J=571,580) / : 0.000543D-6, -433.711737877D0, 0.345585464D0, : 0.000517D-6, -1990.745017041D0, 5.414571768D0, : 0.000504D-6, 5767.611978898D0, 2.328281115D0, : 0.000485D-6, 5753.384884897D0, 1.685874771D0, : 0.000463D-6, 7860.419392439D0, 5.297703006D0, : 0.000604D-6, 515.463871093D0, 0.591998446D0, : 0.000443D-6, 12168.002696575D0, 4.830881244D0, : 0.000570D-6, 199.072001436D0, 3.899190272D0, : 0.000465D-6, 10969.965257698D0, 0.476681802D0, : 0.000424D-6, -7079.373856808D0, 1.112242763D0 / DATA ((FAIRHD(I,J),I=1,3),J=581,590) / : 0.000427D-6, 735.876513532D0, 1.994214480D0, : 0.000478D-6, -6127.655450557D0, 3.778025483D0, : 0.000414D-6, 10973.555686350D0, 5.441088327D0, : 0.000512D-6, 1589.072895284D0, 0.107123853D0, : 0.000378D-6, 10984.192351700D0, 0.915087231D0, : 0.000402D-6, 11371.704689758D0, 4.107281715D0, : 0.000453D-6, 9917.696874510D0, 1.917490952D0, : 0.000395D-6, 149.563197135D0, 2.763124165D0, : 0.000371D-6, 5739.157790895D0, 3.112111866D0, : 0.000350D-6, 11790.629088659D0, 0.440639857D0 / DATA ((FAIRHD(I,J),I=1,3),J=591,600) / : 0.000356D-6, 6133.512652857D0, 5.444568842D0, : 0.000344D-6, 412.371096874D0, 5.676832684D0, : 0.000383D-6, 955.599741609D0, 5.559734846D0, : 0.000333D-6, 6496.374945429D0, 0.261537984D0, : 0.000340D-6, 6055.549660552D0, 5.975534987D0, : 0.000334D-6, 1066.495477190D0, 2.335063907D0, : 0.000399D-6, 11506.769769794D0, 5.321230910D0, : 0.000314D-6, 18319.536584880D0, 2.313312404D0, : 0.000424D-6, 1052.268383188D0, 1.211961766D0, : 0.000307D-6, 63.735898303D0, 3.169551388D0 / DATA ((FAIRHD(I,J),I=1,3),J=601,610) / : 0.000329D-6, 29.821438149D0, 6.106912080D0, : 0.000357D-6, 6309.374169791D0, 4.223760346D0, : 0.000312D-6, -3738.761430108D0, 2.180556645D0, : 0.000301D-6, 309.278322656D0, 1.499984572D0, : 0.000268D-6, 12043.574281889D0, 2.447520648D0, : 0.000257D-6, 12491.370101415D0, 3.662331761D0, : 0.000290D-6, 625.670192312D0, 1.272834584D0, : 0.000256D-6, 5429.879468239D0, 1.913426912D0, : 0.000339D-6, 3496.032826134D0, 4.165930011D0, : 0.000283D-6, 3930.209696220D0, 4.325565754D0 / DATA ((FAIRHD(I,J),I=1,3),J=611,620) / : 0.000241D-6, 12528.018664345D0, 3.832324536D0, : 0.000304D-6, 4686.889407707D0, 1.612348468D0, : 0.000259D-6, 16200.772724501D0, 3.470173146D0, : 0.000238D-6, 12139.553509107D0, 1.147977842D0, : 0.000236D-6, 6172.869528772D0, 3.776271728D0, : 0.000296D-6, -7058.598461315D0, 0.460368852D0, : 0.000306D-6, 10575.406682942D0, 0.554749016D0, : 0.000251D-6, 17298.182327326D0, 0.834332510D0, : 0.000290D-6, 4732.030627343D0, 4.759564091D0, : 0.000261D-6, 5884.926846583D0, 0.298259862D0 / DATA ((FAIRHD(I,J),I=1,3),J=621,630) / : 0.000249D-6, 5547.199336460D0, 3.749366406D0, : 0.000213D-6, 11712.955318231D0, 5.415666119D0, : 0.000223D-6, 4701.116501708D0, 2.703203558D0, : 0.000268D-6, -640.877607382D0, 0.283670793D0, : 0.000209D-6, 5636.065016677D0, 1.238477199D0, : 0.000193D-6, 10177.257679534D0, 1.943251340D0, : 0.000182D-6, 6283.143160294D0, 2.456157599D0, : 0.000184D-6, -227.526189440D0, 5.888038582D0, : 0.000182D-6, -6283.008539689D0, 0.241332086D0, : 0.000228D-6, -6284.056171060D0, 2.657323816D0 / DATA ((FAIRHD(I,J),I=1,3),J=631,640) / : 0.000166D-6, 7238.675591600D0, 5.930629110D0, : 0.000167D-6, 3097.883822726D0, 5.570955333D0, : 0.000159D-6, -323.505416657D0, 5.786670700D0, : 0.000154D-6, -4136.910433516D0, 1.517805532D0, : 0.000176D-6, 12029.347187887D0, 3.139266834D0, : 0.000167D-6, 12132.439962106D0, 3.556352289D0, : 0.000153D-6, 202.253395174D0, 1.463313961D0, : 0.000157D-6, 17267.268201691D0, 1.586837396D0, : 0.000142D-6, 83996.847317911D0, 0.022670115D0, : 0.000152D-6, 17260.154654690D0, 0.708528947D0 / DATA ((FAIRHD(I,J),I=1,3),J=641,650) / : 0.000144D-6, 6084.003848555D0, 5.187075177D0, : 0.000135D-6, 5756.566278634D0, 1.993229262D0, : 0.000134D-6, 5750.203491159D0, 3.457197134D0, : 0.000144D-6, 5326.786694021D0, 6.066193291D0, : 0.000160D-6, 11015.106477335D0, 1.710431974D0, : 0.000133D-6, 3634.621024518D0, 2.836451652D0, : 0.000134D-6, 18073.704938650D0, 5.453106665D0, : 0.000134D-6, 1162.474704408D0, 5.326898811D0, : 0.000128D-6, 5642.198242609D0, 2.511652591D0, : 0.000160D-6, 632.783739313D0, 5.628785365D0 / DATA ((FAIRHD(I,J),I=1,3),J=651,660) / : 0.000132D-6, 13916.019109642D0, 0.819294053D0, : 0.000122D-6, 14314.168113050D0, 5.677408071D0, : 0.000125D-6, 12359.966151546D0, 5.251984735D0, : 0.000121D-6, 5749.452731634D0, 2.210924603D0, : 0.000136D-6, -245.831646229D0, 1.646502367D0, : 0.000120D-6, 5757.317038160D0, 3.240883049D0, : 0.000134D-6, 12146.667056108D0, 3.059480037D0, : 0.000137D-6, 6206.809778716D0, 1.867105418D0, : 0.000141D-6, 17253.041107690D0, 2.069217456D0, : 0.000129D-6, -7477.522860216D0, 2.781469314D0 / DATA ((FAIRHD(I,J),I=1,3),J=661,670) / : 0.000116D-6, 5540.085789459D0, 4.281176991D0, : 0.000116D-6, 9779.108676125D0, 3.320925381D0, : 0.000129D-6, 5237.921013804D0, 3.497704076D0, : 0.000113D-6, 5959.570433334D0, 0.983210840D0, : 0.000122D-6, 6282.095528923D0, 2.674938860D0, : 0.000140D-6, -11.045700264D0, 4.957936982D0, : 0.000108D-6, 23543.230504682D0, 1.390113589D0, : 0.000106D-6, -12569.674818332D0, 0.429631317D0, : 0.000110D-6, -266.607041722D0, 5.501340197D0, : 0.000115D-6, 12559.038152982D0, 4.691456618D0 / DATA ((FAIRHD(I,J),I=1,3),J=671,680) / : 0.000134D-6, -2388.894020449D0, 0.577313584D0, : 0.000109D-6, 10440.274292604D0, 6.218148717D0, : 0.000102D-6, -543.918059096D0, 1.477842615D0, : 0.000108D-6, 21228.392023546D0, 2.237753948D0, : 0.000101D-6, -4535.059436924D0, 3.100492232D0, : 0.000103D-6, 76.266071276D0, 5.594294322D0, : 0.000104D-6, 949.175608970D0, 5.674287810D0, : 0.000101D-6, 13517.870106233D0, 2.196632348D0, : 0.000100D-6, 11933.367960670D0, 4.056084160D0, : 4.322990D-6, 6283.075849991D0, 2.642893748D0 / DATA ((FAIRHD(I,J),I=1,3),J=681,690) / : 0.406495D-6, 0.000000000D0, 4.712388980D0, : 0.122605D-6, 12566.151699983D0, 2.438140634D0, : 0.019476D-6, 213.299095438D0, 1.642186981D0, : 0.016916D-6, 529.690965095D0, 4.510959344D0, : 0.013374D-6, -3.523118349D0, 1.502210314D0, : 0.008042D-6, 26.298319800D0, 0.478549024D0, : 0.007824D-6, 155.420399434D0, 5.254710405D0, : 0.004894D-6, 5746.271337896D0, 4.683210850D0, : 0.004875D-6, 5760.498431898D0, 0.759507698D0, : 0.004416D-6, 5223.693919802D0, 6.028853166D0 / DATA ((FAIRHD(I,J),I=1,3),J=691,700) / : 0.004088D-6, -7.113547001D0, 0.060926389D0, : 0.004433D-6, 77713.771467920D0, 3.627734103D0, : 0.003277D-6, 18849.227549974D0, 2.327912542D0, : 0.002703D-6, 6062.663207553D0, 1.271941729D0, : 0.003435D-6, -775.522611324D0, 0.747446224D0, : 0.002618D-6, 6076.890301554D0, 3.633715689D0, : 0.003146D-6, 206.185548437D0, 5.647874613D0, : 0.002544D-6, 1577.343542448D0, 6.232904270D0, : 0.002218D-6, -220.412642439D0, 1.309509946D0, : 0.002197D-6, 5856.477659115D0, 2.407212349D0 / DATA ((FAIRHD(I,J),I=1,3),J=701,710) / : 0.002897D-6, 5753.384884897D0, 5.863842246D0, : 0.001766D-6, 426.598190876D0, 0.754113147D0, : 0.001738D-6, -796.298006816D0, 2.714942671D0, : 0.001695D-6, 522.577418094D0, 2.629369842D0, : 0.001584D-6, 5507.553238667D0, 1.341138229D0, : 0.001503D-6, -242.728603974D0, 0.377699736D0, : 0.001552D-6, -536.804512095D0, 2.904684667D0, : 0.001370D-6, -398.149003408D0, 1.265599125D0, : 0.001889D-6, -5573.142801634D0, 4.413514859D0, : 0.001722D-6, 6069.776754553D0, 2.445966339D0 / DATA ((FAIRHD(I,J),I=1,3),J=711,720) / : 0.001124D-6, 1059.381930189D0, 5.041799657D0, : 0.001258D-6, 553.569402842D0, 3.849557278D0, : 0.000831D-6, 951.718406251D0, 2.471094709D0, : 0.000767D-6, 4694.002954708D0, 5.363125422D0, : 0.000756D-6, 1349.867409659D0, 1.046195744D0, : 0.000775D-6, -11.045700264D0, 0.245548001D0, : 0.000597D-6, 2146.165416475D0, 4.543268798D0, : 0.000568D-6, 5216.580372801D0, 4.178853144D0, : 0.000711D-6, 1748.016413067D0, 5.934271972D0, : 0.000499D-6, 12036.460734888D0, 0.624434410D0 / DATA ((FAIRHD(I,J),I=1,3),J=721,730) / : 0.000671D-6, -1194.447010225D0, 4.136047594D0, : 0.000488D-6, 5849.364112115D0, 2.209679987D0, : 0.000621D-6, 6438.496249426D0, 4.518860804D0, : 0.000495D-6, -6286.598968340D0, 1.868201275D0, : 0.000456D-6, 5230.807466803D0, 1.271231591D0, : 0.000451D-6, 5088.628839767D0, 0.084060889D0, : 0.000435D-6, 5643.178563677D0, 3.324456609D0, : 0.000387D-6, 10977.078804699D0, 4.052488477D0, : 0.000547D-6, 161000.685737473D0, 2.841633844D0, : 0.000522D-6, 3154.687084896D0, 2.171979966D0 / DATA ((FAIRHD(I,J),I=1,3),J=731,740) / : 0.000375D-6, 5486.777843175D0, 4.983027306D0, : 0.000421D-6, 5863.591206116D0, 4.546432249D0, : 0.000439D-6, 7084.896781115D0, 0.522967921D0, : 0.000309D-6, 2544.314419883D0, 3.172606705D0, : 0.000347D-6, 4690.479836359D0, 1.479586566D0, : 0.000317D-6, 801.820931124D0, 3.553088096D0, : 0.000262D-6, 419.484643875D0, 0.606635550D0, : 0.000248D-6, 6836.645252834D0, 3.014082064D0, : 0.000245D-6, -1592.596013633D0, 5.519526220D0, : 0.000225D-6, 4292.330832950D0, 2.877956536D0 / DATA ((FAIRHD(I,J),I=1,3),J=741,750) / : 0.000214D-6, 7234.794256242D0, 1.605227587D0, : 0.000205D-6, 5767.611978898D0, 0.625804796D0, : 0.000180D-6, 10447.387839604D0, 3.499954526D0, : 0.000229D-6, 199.072001436D0, 5.632304604D0, : 0.000214D-6, 639.897286314D0, 5.960227667D0, : 0.000175D-6, -433.711737877D0, 2.162417992D0, : 0.000209D-6, 515.463871093D0, 2.322150893D0, : 0.000173D-6, 6040.347246017D0, 2.556183691D0, : 0.000184D-6, 6309.374169791D0, 4.732296790D0, : 0.000227D-6, 149854.400134205D0, 5.385812217D0 / DATA ((FAIRHD(I,J),I=1,3),J=751,760) / : 0.000154D-6, 8031.092263058D0, 5.120720920D0, : 0.000151D-6, 5739.157790895D0, 4.815000443D0, : 0.000197D-6, 7632.943259650D0, 0.222827271D0, : 0.000197D-6, 74.781598567D0, 3.910456770D0, : 0.000138D-6, 6055.549660552D0, 1.397484253D0, : 0.000149D-6, -6127.655450557D0, 5.333727496D0, : 0.000137D-6, 3894.181829542D0, 4.281749907D0, : 0.000135D-6, 9437.762934887D0, 5.979971885D0, : 0.000139D-6, -2352.866153772D0, 4.715630782D0, : 0.000142D-6, 6812.766815086D0, 0.513330157D0 / DATA ((FAIRHD(I,J),I=1,3),J=761,770) / : 0.000120D-6, -4705.732307544D0, 0.194160689D0, : 0.000131D-6, -71430.695617928D0, 0.000379226D0, : 0.000124D-6, 6279.552731642D0, 2.122264908D0, : 0.000108D-6, -6256.777530192D0, 0.883445696D0, : 0.143388D-6, 6283.075849991D0, 1.131453581D0, : 0.006671D-6, 12566.151699983D0, 0.775148887D0, : 0.001480D-6, 155.420399434D0, 0.480016880D0, : 0.000934D-6, 213.299095438D0, 6.144453084D0, : 0.000795D-6, 529.690965095D0, 2.941595619D0, : 0.000673D-6, 5746.271337896D0, 0.120415406D0 / DATA ((FAIRHD(I,J),I=1,3),J=771,780) / : 0.000672D-6, 5760.498431898D0, 5.317009738D0, : 0.000389D-6, -220.412642439D0, 3.090323467D0, : 0.000373D-6, 6062.663207553D0, 3.003551964D0, : 0.000360D-6, 6076.890301554D0, 1.918913041D0, : 0.000316D-6, -21.340641002D0, 5.545798121D0, : 0.000315D-6, -242.728603974D0, 1.884932563D0, : 0.000278D-6, 206.185548437D0, 1.266254859D0, : 0.000238D-6, -536.804512095D0, 4.532664830D0, : 0.000185D-6, 522.577418094D0, 4.578313856D0, : 0.000245D-6, 18849.227549974D0, 0.587467082D0 / DATA ((FAIRHD(I,J),I=1,3),J=781,787) / : 0.000180D-6, 426.598190876D0, 5.151178553D0, : 0.000200D-6, 553.569402842D0, 5.355983739D0, : 0.000141D-6, 5223.693919802D0, 1.336556009D0, : 0.000104D-6, 5856.477659115D0, 4.239842759D0, : 0.003826D-6, 6283.075849991D0, 5.705257275D0, : 0.000303D-6, 12566.151699983D0, 5.407132842D0, : 0.000209D-6, 155.420399434D0, 1.989815753D0 / * - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - * Time since J2000.0 in Julian millennia. T = ( ( EPOCH1-DJ0 ) + EPOCH2 ) / DJM * ================= * Topocentric terms * ================= * Convert UT to local solar time in radians. TSOL = MOD(UT,1D0) * D2PI + ELONG * FUNDAMENTAL ARGUMENTS: Simon et al. 1994. * Combine time argument (millennia) with deg/arcsec factor. W = T / 3600D0 * Sun Mean Longitude. ELSUN = MOD(280.46645683D0 + 1296027711.03429D0 * W, 360D0) * DD2R * Sun Mean Anomaly. EMSUN = MOD(357.52910918D0 + 1295965810.481D0 * W, 360D0) * DD2R * Mean Elongation of Moon from Sun. D = MOD(297.85019547D0 + 16029616012.090D0 * W, 360D0) * DD2R * Mean Longitude of Jupiter. ELJ = MOD(34.35151874D0 + 109306899.89453D0 * W, 360D0) * DD2R * Mean Longitude of Saturn. ELS = MOD(50.07744430D0 + 44046398.47038D0 * W, 360D0) * DD2R * TOPOCENTRIC TERMS: Moyer 1981 and Murray 1983. WT = + 0.00029D-10 * U * SIN(TSOL + ELSUN - ELS) : + 0.00100D-10 * U * SIN(TSOL - 2D0*EMSUN) : + 0.00133D-10 * U * SIN(TSOL - D) : + 0.00133D-10 * U * SIN(TSOL + ELSUN - ELJ) : - 0.00229D-10 * U * SIN(TSOL + 2D0*ELSUN + EMSUN) : - 0.0220 D-10 * V * COS(ELSUN + EMSUN) : + 0.05312D-10 * U * SIN(TSOL - EMSUN) : - 0.13677D-10 * U * SIN(TSOL + 2D0*ELSUN) : - 1.3184 D-10 * V * COS(ELSUN) : + 3.17679D-10 * U * SIN(TSOL) * ===================== * Fairhead et al. model * ===================== * T**0 W0 = 0D0 DO 10 J=474,1,-1 W0 = W0 + FAIRHD(1,J) * SIN(FAIRHD(2,J)*T + FAIRHD(3,J)) 10 CONTINUE * T**1 W1 = 0D0 DO 11 J=679,475,-1 W1 = W1 + FAIRHD(1,J) * SIN(FAIRHD(2,J)*T + FAIRHD(3,J)) 11 CONTINUE * T**2 W2 = 0D0 DO 12 J=764,680,-1 W2 = W2 + FAIRHD(1,J) * SIN(FAIRHD(2,J)*T + FAIRHD(3,J)) 12 CONTINUE * T**3 W3 = 0D0 DO 13 J=784,765,-1 W3 = W3 + FAIRHD(1,J) * SIN(FAIRHD(2,J)*T + FAIRHD(3,J)) 13 CONTINUE * T**4 W4 = 0D0 DO 14 J=787,785,-1 W4 = W4 + FAIRHD(1,J) * SIN(FAIRHD(2,J)*T + FAIRHD(3,J)) 14 CONTINUE * Multiply by powers of T and combine. WF = T * ( T * ( T * ( T * W4 + W3 ) + W2 ) + W1 ) + W0 * Adjustments to use JPL planetary masses instead of IAU. WJ = 0.00065D-6 * SIN(6069.776754D0*T + 4.021194D0) + : 0.00033D-6 * SIN( 213.299095D0*T + 5.543132D0) + : ( -0.00196D-6 * SIN(6208.294251D0*T + 5.696701D0) ) + : ( -0.00173D-6 * SIN( 74.781599D0*T + 2.435900D0) ) + : 0.03638D-6 * T * T * ============ * Final result * ============ * TDB-TT in seconds. iau_DTDB = WT + WF + WJ * Finished. *+---------------------------------------------------------------------- * * Copyright (C) 2003 * Standards Of Fundamental Astronomy Review Board * of the International Astronomical Union. * * ===================== * SOFA Software License * ===================== * * NOTICE TO USER: * * BY USING THIS SOFTWARE YOU ACCEPT THE FOLLOWING TERMS AND CONDITIONS * WHICH APPLY TO ITS USE. * * 1. The Software is owned by the IAU SOFA Review Board ("the Board"). * * 2. The Software is made available free of charge for use by: * * a) private individuals for non-profit research; and * * b) non-profit educational, academic and research institutions. * * 3. Commercial use of the Software is specifically excluded from the * terms and conditions of this license. Commercial use of the * Software is subject to the prior written agreement of the Board on * terms to be agreed. * * 4. The provision of any version of the Software under the terms and * conditions specified herein does not imply that future versions * will also be made available under the same terms and conditions. * * 5. The user may modify the Software for his/her own purposes. The * user may distribute the modified software provided that the Board * is informed and that a copy of the modified software is made * available to the Board on request. All modifications made by the * user shall be clearly identified to show how the modified software * differs from the original Software, and the name(s) of the * affected routine(s) shall be changed. The original SOFA Software * License text must be present. * * 6. In any published work produced by the user and which includes * results achieved by using the Software, the user shall acknowledge * that the Software was used in producing the information contained * in such publication. * * 7. The user may incorporate or embed the Software into other software * products which he/she may then give away free of charge but not * sell provided the user makes due acknowledgement of the use which * he/she has made of the Software in creating such software * products. Any redistribution of the Software in this way shall be * made under the same terms and conditions under which the user * received it from the SOFA Center. * * 8. The user shall not cause the Software to be brought into * disrepute, either by misuse, or use for inappropriate tasks, or by * inappropriate modification. * * 9. The Software is provided to the user "as is" and the Board makes * no warranty as to its use or performance. The Board does not and * cannot warrant the performance or results which the user may * obtain by using the Software. The Board makes no warranties, * express or implied, as to non-infringement of third party rights, * merchantability, or fitness for any particular purpose. In no * event will the Board be liable to the user for any consequential, * incidental, or special damages, including any lost profits or lost * savings, even if a Board representative has been advised of such * damages, or for any claim by any third party. * * Correspondence concerning SOFA software should be addressed as * follows: * * Internet email: sofa@rl.ac.uk * Postal address: IAU SOFA Center * Rutherford Appleton Laboratory * Chilton, Didcot, Oxon OX11 0QX * United Kingdom * * *----------------------------------------------------------------------- END