Changeset 5119 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/gradiv.f90

Timestamp:

Jul 24, 2024, 6:46:45 PM (2 months ago)

Author:

abarral

Message:

enforce PRIVATE by default in several modules, expose PUBLIC as needed
move eigen.f90 to obsolete/
(lint) aslong the way

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3d_common/gradiv.f90 (modified) (3 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3d_common/gradiv.f90

@@ -1 @@
-
 ! $Header$
 
-SUBROUTINE gradiv(klevel, xcov, ycov, ld, gdx, gdy )
+SUBROUTINE gradiv(klevel, xcov, ycov, ld, gdx, gdy)
   !
   !    Auteur :   P. Le Van
@@ -18 +17 @@
   !
   USE lmdz_filtreg, ONLY: filtreg
+  USE lmdz_ssum_scopy, ONLY: scopy
+
   IMPLICIT NONE
   !
@@ -26 +27 @@
   INTEGER :: klevel
   !
-  REAL :: xcov( ip1jmp1,klevel ), ycov( ip1jm,klevel )
-  REAL :: gdx( ip1jmp1,klevel ),   gdy( ip1jm,klevel )
+  REAL :: xcov(ip1jmp1, klevel), ycov(ip1jm, klevel)
+  REAL :: gdx(ip1jmp1, klevel), gdy(ip1jm, klevel)
 
-  REAL :: div(ip1jmp1,llm)
+  REAL :: div(ip1jmp1, llm)
 
-  INTEGER :: l,ij,iter,ld
+  INTEGER :: l, ij, iter, ld
   !
   !
   !
-  CALL SCOPY( ip1jmp1*klevel,xcov,1,gdx,1 )
-  CALL SCOPY( ip1jm*klevel,  ycov,1,gdy,1 )
+  CALL SCOPY(ip1jmp1 * klevel, xcov, 1, gdx, 1)
+  CALL SCOPY(ip1jm * klevel, ycov, 1, gdy, 1)
   !
-  DO iter = 1,ld
-  !
-  CALL  diverg( klevel,  gdx , gdy, div          )
-  CALL filtreg( div, jjp1, klevel, 2,1, .TRUE.,2 )
-  CALL    grad( klevel,  div, gdx, gdy           )
-  !
-  DO l = 1, klevel
-  DO ij = 1, ip1jmp1
-  gdx( ij,l ) = - gdx( ij,l ) * cdivu
-  END DO
-  DO ij = 1, ip1jm
-  gdy( ij,l ) = - gdy( ij,l ) * cdivu
-  END DO
-  END DO
-  !
+  DO iter = 1, ld
+    !
+    CALL  diverg(klevel, gdx, gdy, div)
+    CALL filtreg(div, jjp1, klevel, 2, 1, .TRUE., 2)
+    CALL    grad(klevel, div, gdx, gdy)
+    !
+    DO l = 1, klevel
+      DO ij = 1, ip1jmp1
+        gdx(ij, l) = - gdx(ij, l) * cdivu
+      END DO
+      DO ij = 1, ip1jm
+        gdy(ij, l) = - gdy(ij, l) * cdivu
+      END DO
+    END DO
+    !
   END DO
   RETURN