source: PSPA/parmelaPSPA/trunk/sort.f @ 436

      subroutine sort(n,x,y)
c-----------------------------------------------------------------------
      save
c
      include 'param_sz.h'
      include 'ucom.h'
c
      dimension x(imaa),y(imaa)
c-----------------------------------------------------------------------
c*
      if(x(2).ge.x(1))go to 2
      r=x(2)
      s=y(2)
      x(2)=x(1)
      y(2)=y(1)
      x(1)=r
      y(1)=s
    2 continue
      do 10 i=3,n
      l=i-1
      if(x(i).ge.x(l))go to 10
      r=x(i)
      s=y(i)
      x(i)=x(l)
      y(i)=y(l)
      do 4 j=2,l
      k=i-j
      if(r.ge.x(k))go to 6
      x(k+1)=x(k)
      y(k+1)=y(k)
    4 continue
      k=0
    6 continue
      j=k+1
      x(j)=r
      y(j)=s
   10 continue
      return
      end
c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*