MODULE REAL_PRECISION
C THIS IS FOR 64-BIT IEEE ARITHMETIC.
      INTEGER, PARAMETER:: R8=SELECTED_REAL_KIND(15,307)
      END MODULE REAL_PRECISION