rsort.f

Go to the documentation of this file.

      SUBROUTINE rsort(R,SWITCH,IND)
*
*
*       Sorting of mutual distances.
*       ----------------------------
*
      REAL*8  r(6)
      INTEGER  ind(6)
      LOGICAL  switch
*
*
*       Sort particle separations in increasing order.
    1 nsw = 0
      DO 2 k = 1,5
          IF (r(ind(k+1)).GE.r(ind(k))) go to 2
          ind1 = ind(k)
          ind(k) = ind(k+1)
          ind(k+1) = ind1
          nsw = 1
    2 CONTINUE
*
      IF (nsw.GT.0) go to 1
*
*       Check the switching conditions.
      switch = .false.
      IF (ind(1).GT.3.OR.ind(2).EQ.6) switch = .true.
*
      RETURN
*
      END