SUBROUTINE HEAPSORT(N, RA)
C======================================================================*
C heapsort -  sort data using the HEAPSORT algorithm
C
C This is the routine SORT from "Numerical Recipies".  It takes
C as input the array RA of N real numbers and returns the array
C sorted in ascending order.
C
C From "Numerical Recipes" by W.H. Press, B.P. Flannery, S.A. Teukolsky 
C and W.T. Vetterling, (Cambridge University Press, 1986), pg. 231
C======================================================================*

    Due to copyright restrictions we cannot distribute a copy
    of this file on the web.  However, you are perfectly free to 
    get the book and type in your own copy for private use.


C==EOF==