Changeset 5119 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divgrad.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/divgrad.f90 (modified) (2 diffs)

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

@@ -1 @@
-
 ! $Header$
 
-SUBROUTINE divgrad(klevel,h, lh, divgra )
+SUBROUTINE divgrad(klevel, h, lh, divgra)
   USE lmdz_filtreg, ONLY: filtreg
+  USE lmdz_ssum_scopy, ONLY: scopy
   IMPLICIT NONE
   !
@@ -27 +27 @@
   !
   INTEGER :: klevel
-  REAL :: h( ip1jmp1,klevel ), divgra( ip1jmp1,klevel )
+  REAL :: h(ip1jmp1, klevel), divgra(ip1jmp1, klevel)
   !
-  REAL :: ghy(ip1jm,llm), ghx(ip1jmp1,llm)
+  REAL :: ghy(ip1jm, llm), ghx(ip1jmp1, llm)
 
-  INTEGER :: l,ij,iter,lh
+  INTEGER :: l, ij, iter, lh
   !
   !
   !
-  CALL SCOPY ( ip1jmp1*klevel,h,1,divgra,1 )
+  CALL SCOPY (ip1jmp1 * klevel, h, 1, divgra, 1)
   !
-  DO iter = 1,lh
+  DO iter = 1, lh
 
-  CALL filtreg ( divgra,jjp1,klevel,2,1,.TRUE.,1  )
+    CALL filtreg (divgra, jjp1, klevel, 2, 1, .TRUE., 1)
 
-  CALL    grad (klevel,divgra, ghx  , ghy          )
-  CALL  diverg (klevel,  ghx , ghy  , divgra       )
+    CALL    grad (klevel, divgra, ghx, ghy)
+    CALL  diverg (klevel, ghx, ghy, divgra)
 
-  CALL filtreg ( divgra,jjp1,klevel,2,1,.TRUE.,1)
+    CALL filtreg (divgra, jjp1, klevel, 2, 1, .TRUE., 1)
 
-  DO l = 1,klevel
-  DO ij = 1, ip1jmp1
-  divgra( ij,l ) = - cdivh * divgra( ij,l )
-  END DO
-  END DO
-  !
+    DO l = 1, klevel
+      DO ij = 1, ip1jmp1
+        divgra(ij, l) = - cdivh * divgra(ij, l)
+      END DO
+    END DO
+    !
   END DO
   RETURN