subroutine rmysort(n, x, ipermut)
      implicit none
      integer n
      double precision x(n),xtemp
      integer ipermut(n)
      integer i,j,imax,ipm
      do i=1,n
        xtemp=x(i)
        imax=i
        do j=i+1,n
          if (x(ipermut(j)).lt.xtemp) then
            imax=j
            xtemp=x(ipermut(j))
          endif
        enddo
        ipm=ipermut(imax)
        ipermut(imax)=ipermut(i)
        ipermut(i)=ipm
      enddo
      return
      end