trunk/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	!
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