Changeset 5105 for LMDZ6/branches/Amaury_dev/libf/filtrez/eigen_sort.f90

Timestamp:

Jul 23, 2024, 7:14:34 PM (8 weeks ago)

Author:

abarral

Message:

Replace 1DUTILS.h by module lmdz_1dutils.f90
Replace 1DConv.h by module lmdz_old_1dconv.f90 (it's only used by old_* files)
Convert *.F to *.f90
Fix gradsdef.h formatting
Remove unnecessary "RETURN" at the end of functions/subroutines

File:

• LMDZ6/branches/Amaury_dev/libf/filtrez/eigen_sort.f90 (moved) (moved from LMDZ6/branches/Amaury_dev/libf/filtrez/eigen_sort.F) (1 diff)

LMDZ6/branches/Amaury_dev/libf/filtrez/eigen_sort.f90

@@ -2 +2 @@
 ! $Header$
 
-          SUBROUTINE eigen_sort(d,v,n,np)
-          INTEGER n,np
-          REAL d(np),v(np,np)
-          INTEGER i,j,k
-          REAL p
+    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)>=p) THEN
-            k=j
-            p=d(j)
-           ENDIF
-        ENDDO
-          
-        IF(k/=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
+   DO i=1,n-1
+      k=i
+      p=d(i)
+    DO j=i+1,n
+       IF(d(j)>=p) THEN
+        k=j
+        p=d(j)
+       ENDIF
+    ENDDO
 
-        RETURN
-        END
+    IF(k/=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 SUBROUTINE eigen_sort