libf/filtrez/eigen_sort.f

!
! $Header: /home/cvsroot/LMDZ4/libf/filtrez/eigen_sort.F,v 1.1.1.1 2004/05/19 12:53:09 lmdzadmin Exp $
!
          SUBROUTINE eigen_sort(d,v,n,np)
          INTEGER n,np
          REAL d(np),v(np,np)
          INTEGER i,j,k
          REAL p

       DO i=1,n-1
          k=i
          p=d(i)
        DO j=i+1,n
           IF(d(j).ge.p) THEN
            k=j
            p=d(j)
           ENDIF
        ENDDO
          
        IF(k.ne.i) THEN
          d(k)=d(i)
          d(i)=p
         DO j=1,n
          p=v(j,i)
          v(j,i)=v(j,k)
          v(j,k)=p
         ENDDO
        ENDIF
       ENDDO

        RETURN
        END
1	guez	3	!
2			! $Header: /home/cvsroot/LMDZ4/libf/filtrez/eigen_sort.F,v 1.1.1.1 2004/05/19 12:53:09 lmdzadmin Exp $
3			!
4			SUBROUTINE eigen_sort(d,v,n,np)
5			INTEGER n,np
6			REAL d(np),v(np,np)
7			INTEGER i,j,k
8			REAL p
9
10			DO i=1,n-1
11			k=i
12			p=d(i)
13			DO j=i+1,n
14			IF(d(j).ge.p) THEN
15			k=j
16			p=d(j)
17			ENDIF
18			ENDDO
19
20			IF(k.ne.i) THEN
21			d(k)=d(i)
22			d(i)=p
23			DO j=1,n
24			p=v(j,i)
25			v(j,i)=v(j,k)
26			v(j,k)=p
27			ENDDO
28			ENDIF
29			ENDDO
30
31			RETURN
32			END