DOUBLE PRECISION FUNCTION iau_DTDB ( DATE1, DATE2,
: UT, ELONG, U, V )
*+
* - - - - - - - - -
* i a u _ D T D B
* - - - - - - - - -
*
* An approximation to TDB-TT, the difference between barycentric
* dynamical time and terrestrial time, for an observer on the Earth.
*
* The different time scales - proper, coordinate and realized - are
* related to each other:
*
* TAI <- physically realized
* :
* offset <- observed (nominally +32.184s)
* :
* TT <- terrestrial time
* :
* rate adjustment (L_G) <- definition of TT
* :
* TCG <- time scale for GCRS
* :
* "periodic" terms <- iau_DTDB is an implementation
* :
* rate adjustment (L_C) <- function of solar-system ephemeris
* :
* TCB <- time scale for BCRS
* :
* rate adjustment (-L_B) <- definition of TDB
* :
* TDB <- TCB scaled to track TT
* :
* "periodic" terms <- -iau_DTDB is an approximation
* :
* TT <- terrestrial time
*
* Adopted values for the various constants can be found in the IERS
* Conventions (McCarthy & Petit 2003).
*
* This routine is part of the International Astronomical Union's
* SOFA (Standards of Fundamental Astronomy) software collection.
*
* Status: support routine.
*
* Given:
* DATE1,DATE2 d date, TDB (Notes 1-3)
* UT d universal time (UT1, fraction of one day)
* ELONG d longitude (east positive, 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 date DATE1+DATE2 is a Julian Date, apportioned in any
* convenient way between the arguments DATE1 and DATE2. For
* example, JD(TDB)=2450123.7 could be expressed in any of these
* ways, among others:
*
* DATE1 DATE2
*
* 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 date 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) TT can be regarded as a coordinate time that is realized as an
* offset of 32.184s from International Atomic Time, TAI. TT is a
* specific linear transformation of geocentric coordinate time TCG,
* which is the time scale for the Geocentric Celestial Reference
* System, GCRS.
*
* 3) TDB is a coordinate time, and is a specific linear transformation
* of barycentric coordinate time TCB, which is the time scale for
* the Barycentric Celestial Reference System, BCRS.
*
* 4) The difference TCG-TCB depends on the masses and positions of the
* bodies of the solar system and the velocity of the Earth. It is
* dominated by a rate difference, the residual being of a periodic
* character. The latter, which is modeled by the present routine,
* comprises 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.
* These effects come from the changing transverse Doppler effect
* and gravitational red-shift as the observer (on the Earth's
* surface) experiences variations in speed (with respect to the
* BCRS) and gravitational potential.
*
* 5) TDB can be regarded as the same as TCB but with a rate adjustment
* to keep it close to TT, which is convenient for many applications.
* The history of successive attempts to define TDB is set out in
* Resolution 3 adopted by the IAU General Assembly in 2006, which
* defines a fixed TDB(TCB) transformation that is consistent with
* contemporary solar-system ephemerides. Future ephemerides will
* imply slightly changed transformations between TCG and TCB, which
* could introduce a linear drift between TDB and TT; however, any
* such drift is unlikely to exceed 1 nanosecond per century.
*
* 6) The geocentric TDB-TT model used in the present routine is that of
* Fairhead & Bretagnon (1990), in its full form. It was originally
* supplied by Fairhead (private communications with P.T.Wallace,
* 1990) as 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 1D-20 s
* level.
*
* 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.
*
* 7) During the interval 1950-2050, the absolute accuracy is better
* than +/- 3 nanoseconds relative to time ephemerides obtained by
* direct numerical integrations based on the JPL DE405 solar system
* ephemeris.
*
* 8) 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 and hence between TT and TDB.
*
* References:
*
* Fairhead, L., & Bretagnon, P., Astron.Astrophys., 229, 240-247
* (1990).
*
* IAU 2006 Resolution 3.
*
* McCarthy, D. D., Petit, G. (eds.), IERS Conventions (2003),
* IERS Technical Note No. 32, BKG (2004)
*
* 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., Astron.Astrophys., 282, 663-683 (1994).
*
* This revision: 2010 July 29
*
* SOFA release 2012-03-01
*
* Copyright (C) 2012 IAU SOFA Board. See notes at end.
*
*-----------------------------------------------------------------------
IMPLICIT NONE
DOUBLE PRECISION DATE1, DATE2, 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.0), JD
DOUBLE PRECISION DJ00
PARAMETER ( DJ00 = 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.0),
* 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,
* T^2
: 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,
* T^3
: 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,
* T^4
: 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 = ( ( DATE1-DJ00 ) + DATE2 ) / 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) 2012
* Standards Of Fundamental Astronomy Board
* of the International Astronomical Union.
*
* =====================
* SOFA Software License
* =====================
*
* NOTICE TO USER:
*
* BY USING THIS SOFTWARE YOU ACCEPT THE FOLLOWING SIX TERMS AND
* CONDITIONS WHICH APPLY TO ITS USE.
*
* 1. The Software is owned by the IAU SOFA Board ("SOFA").
*
* 2. Permission is granted to anyone to use the SOFA software for any
* purpose, including commercial applications, free of charge and
* without payment of royalties, subject to the conditions and
* restrictions listed below.
*
* 3. You (the user) may copy and distribute SOFA source code to others,
* and use and adapt its code and algorithms in your own software,
* on a world-wide, royalty-free basis. That portion of your
* distribution that does not consist of intact and unchanged copies
* of SOFA source code files is a "derived work" that must comply
* with the following requirements:
*
* a) Your work shall be marked or carry a statement that it
* (i) uses routines and computations derived by you from
* software provided by SOFA under license to you; and
* (ii) does not itself constitute software provided by and/or
* endorsed by SOFA.
*
* b) The source code of your derived work must contain descriptions
* of how the derived work is based upon, contains and/or differs
* from the original SOFA software.
*
* c) The names of all routines in your derived work shall not
* include the prefix "iau" or "sofa" or trivial modifications
* thereof such as changes of case.
*
* d) The origin of the SOFA components of your derived work must
* not be misrepresented; you must not claim that you wrote the
* original software, nor file a patent application for SOFA
* software or algorithms embedded in the SOFA software.
*
* e) These requirements must be reproduced intact in any source
* distribution and shall apply to anyone to whom you have
* granted a further right to modify the source code of your
* derived work.
*
* Note that, as originally distributed, the SOFA software is
* intended to be a definitive implementation of the IAU standards,
* and consequently third-party modifications are discouraged. All
* variations, no matter how minor, must be explicitly marked as
* such, as explained above.
*
* 4. You shall not cause the SOFA software to be brought into
* disrepute, either by misuse, or use for inappropriate tasks, or
* by inappropriate modification.
*
* 5. The SOFA software is provided "as is" and SOFA makes no warranty
* as to its use or performance. SOFA does not and cannot warrant
* the performance or results which the user may obtain by using the
* SOFA software. SOFA 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 SOFA be
* liable to the user for any consequential, incidental, or special
* damages, including any lost profits or lost savings, even if a
* SOFA representative has been advised of such damages, or for any
* claim by any third party.
*
* 6. The provision of any version of the SOFA 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.
*
* In any published work or commercial product which uses the SOFA
* software directly, acknowledgement (see www.iausofa.org) is
* appreciated.
*
* Correspondence concerning SOFA software should be addressed as
* follows:
*
* By email: sofa@ukho.gov.uk
* By post: IAU SOFA Center
* HM Nautical Almanac Office
* UK Hydrographic Office
* Admiralty Way, Taunton
* Somerset, TA1 2DN
* United Kingdom
*
*-----------------------------------------------------------------------
END