PageRenderTime 29ms CodeModel.GetById 17ms RepoModel.GetById 0ms app.codeStats 0ms

/measures/fortran/sofa/gst00a.f

http://casacore.googlecode.com/
FORTRAN Legacy | 173 lines | 7 code | 8 blank | 158 comment | 0 complexity | 4b67c08c004b6e3de65e3a8648edfbf8 MD5 | raw file
Possible License(s): GPL-2.0
  1. DOUBLE PRECISION FUNCTION iau_GST00A ( UTA, UTB, TTA, TTB )
  2. *+
  3. * - - - - - - - - - - -
  4. * i a u _ G S T 0 0 A
  5. * - - - - - - - - - - -
  6. *
  7. * Greenwich Apparent Sidereal Time (consistent with IAU 2000
  8. * resolutions).
  9. *
  10. * This routine is part of the International Astronomical Union's
  11. * SOFA (Standards of Fundamental Astronomy) software collection.
  12. *
  13. * Status: canonical model.
  14. *
  15. * Given:
  16. * UTA, UTB d UT1 as a 2-part Julian Date (Notes 1,2)
  17. * TTA, TTB d TT as a 2-part Julian Date (Notes 1,2)
  18. *
  19. * The result is the Greenwich Apparent Sidereal Time (radians), in the
  20. * range 0 to 2pi.
  21. *
  22. * Notes:
  23. *
  24. * 1) The UT1 and TT dates UTA+UTB and TTA+TTB respectively, are both
  25. * Julian Dates, apportioned in any convenient way between the
  26. * argument pairs. For example, JD=2450123.7 could be expressed in
  27. * any of these ways, among others:
  28. *
  29. * Part A Part B
  30. *
  31. * 2450123.7D0 0D0 (JD method)
  32. * 2451545D0 -1421.3D0 (J2000 method)
  33. * 2400000.5D0 50123.2D0 (MJD method)
  34. * 2450123.5D0 0.2D0 (date & time method)
  35. *
  36. * The JD method is the most natural and convenient to use in
  37. * cases where the loss of several decimal digits of resolution
  38. * is acceptable (in the case of UT; the TT is not at all critical
  39. * in this respect). The J2000 and MJD methods are good compromises
  40. * between resolution and convenience. For UT, the date & time
  41. * method is best matched to the algorithm that is used by the Earth
  42. * Rotation Angle routine, called internally: maximum accuracy (or,
  43. * at least, minimum noise) is delivered when the UTA argument is for
  44. * 0hrs UT1 on the day in question and the UTB argument lies in the
  45. * range 0 to 1, or vice versa.
  46. *
  47. * 2) Both UT1 and TT are required, UT1 to predict the Earth rotation
  48. * and TT to predict the effects of precession-nutation. If UT1 is
  49. * used for both purposes, errors of order 100 microarcseconds
  50. * result.
  51. *
  52. * 3) This GAST is compatible with the IAU 2000 resolutions and must be
  53. * used only in conjunction with other IAU 2000 compatible components
  54. * such as precession-nutation.
  55. *
  56. * 4) The algorithm is from Capitaine et al. (2003) and IERS Conventions
  57. * 2003.
  58. *
  59. * Called:
  60. * iau_GMST00 Greenwich Mean Sidereal Time, IAU 2000
  61. * iau_EE00A equation of the equinoxes, IAU 2000A
  62. * iau_ANP normalize angle into range 0 to 2pi
  63. *
  64. * References:
  65. *
  66. * Capitaine, N., Wallace, P.T. and McCarthy, D.D., "Expressions to
  67. * implement the IAU 2000 definition of UT1", Astronomy &
  68. * Astrophysics, 406, 1135-1149 (2003)
  69. *
  70. * McCarthy, D. D., Petit, G. (eds.), IERS Conventions (2003),
  71. * IERS Technical Note No. 32, BKG (2004)
  72. *
  73. * This revision: 2005 August 24
  74. *
  75. * Copyright (C) 2005 IAU SOFA Review Board. See notes at end.
  76. *
  77. *-----------------------------------------------------------------------
  78. IMPLICIT NONE
  79. DOUBLE PRECISION UTA, UTB, TTA, TTB
  80. DOUBLE PRECISION iau_ANP, iau_GMST00, iau_EE00A
  81. * - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
  82. iau_GST00A = iau_ANP ( iau_GMST00 ( UTA,UTB, TTA,TTB ) +
  83. : iau_EE00A ( TTA,TTB ) )
  84. * Finished.
  85. *+----------------------------------------------------------------------
  86. *
  87. * Copyright (C) 2005
  88. * Standards Of Fundamental Astronomy Review Board
  89. * of the International Astronomical Union.
  90. *
  91. * =====================
  92. * SOFA Software License
  93. * =====================
  94. *
  95. * NOTICE TO USER:
  96. *
  97. * BY USING THIS SOFTWARE YOU ACCEPT THE FOLLOWING TERMS AND CONDITIONS
  98. * WHICH APPLY TO ITS USE.
  99. *
  100. * 1. The Software is owned by the IAU SOFA Review Board ("the Board").
  101. *
  102. * 2. The Software is made available free of charge for use by:
  103. *
  104. * a) private individuals for non-profit research; and
  105. *
  106. * b) non-profit educational, academic and research institutions.
  107. *
  108. * 3. Commercial use of the Software is specifically excluded from the
  109. * terms and conditions of this license. Commercial use of the
  110. * Software is subject to the prior written agreement of the Board on
  111. * terms to be agreed.
  112. *
  113. * 4. The provision of any version of the Software under the terms and
  114. * conditions specified herein does not imply that future versions
  115. * will also be made available under the same terms and conditions.
  116. *
  117. * 5. The user may modify the Software for his/her own purposes. The
  118. * user may distribute the modified software provided that the Board
  119. * is informed and that a copy of the modified software is made
  120. * available to the Board on request. All modifications made by the
  121. * user shall be clearly identified to show how the modified software
  122. * differs from the original Software, and the name(s) of the
  123. * affected routine(s) shall be changed. The original SOFA Software
  124. * License text must be present.
  125. *
  126. * 6. In any published work produced by the user and which includes
  127. * results achieved by using the Software, the user shall acknowledge
  128. * that the Software was used in producing the information contained
  129. * in such publication.
  130. *
  131. * 7. The user may incorporate or embed the Software into other software
  132. * products which he/she may then give away free of charge but not
  133. * sell provided the user makes due acknowledgement of the use which
  134. * he/she has made of the Software in creating such software
  135. * products. Any redistribution of the Software in this way shall be
  136. * made under the same terms and conditions under which the user
  137. * received it from the SOFA Center.
  138. *
  139. * 8. The user shall not cause the Software to be brought into
  140. * disrepute, either by misuse, or use for inappropriate tasks, or by
  141. * inappropriate modification.
  142. *
  143. * 9. The Software is provided to the user "as is" and the Board makes
  144. * no warranty as to its use or performance. The Board does not and
  145. * cannot warrant the performance or results which the user may
  146. * obtain by using the Software. The Board makes no warranties,
  147. * express or implied, as to non-infringement of third party rights,
  148. * merchantability, or fitness for any particular purpose. In no
  149. * event will the Board be liable to the user for any consequential,
  150. * incidental, or special damages, including any lost profits or lost
  151. * savings, even if a Board representative has been advised of such
  152. * damages, or for any claim by any third party.
  153. *
  154. * Correspondence concerning SOFA software should be addressed as
  155. * follows:
  156. *
  157. * Internet email: sofa@rl.ac.uk
  158. * Postal address: IAU SOFA Center
  159. * Rutherford Appleton Laboratory
  160. * Chilton, Didcot, Oxon OX11 0QX
  161. * United Kingdom
  162. *
  163. *
  164. *-----------------------------------------------------------------------
  165. END