!-----------------------------------------------------------------
      ! function  : sort
      !
      ! package   : li3 Potentials
      !
      ! Language  : Fortran 77
      !
      ! author    : V. Venturi (vanessa.venturi@nist.gov)
      !
      ! purpose   : sorting vector a and b in ascending order
      !-----------------------------------------------------------------
SUBROUTINE sort5 (a, b, n)
USE Numeric_Kinds_Module
IMPLICIT NONE
INTEGER :: n, i, j
REAL(KIND=DP_Kind) :: a (n), atem, b (n), btem
DO 10 j = 1, n
   DO 20 i = 1, n - j
      IF(a (i) .gt.a (i + 1) )THEN
         atem = a (i + 1)
         a (i + 1) = a (i)
         a (i) = atem
         btem = b (i + 1)
         b (i + 1) = b (i)
         b (i) = btem
      ENDIF
   20    ENDDO
   10 ENDDO
RETURN
END SUBROUTINE sort5