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

Timestamp:

Jul 24, 2024, 6:46:45 PM (4 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/nxgraro2.f90 (modified) (3 diffs)

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

@@ -1 @@
-
 ! $Header$
 
- SUBROUTINE nxgraro2(klevel,xcov, ycov, lr, grx, gry )
+SUBROUTINE nxgraro2(klevel, xcov, ycov, lr, grx, gry)
   !
   !  P.Le Van .
@@ -16 +15 @@
   !
   USE lmdz_filtreg, ONLY: filtreg
+  USE lmdz_ssum_scopy, ONLY: scopy
+
   IMPLICIT NONE
   !
@@ -25 +26 @@
   !
   INTEGER :: klevel
-  REAL :: xcov( ip1jmp1,klevel ), ycov( ip1jm,klevel )
-  REAL :: grx( ip1jmp1,klevel ),  gry( ip1jm,klevel )
+  REAL :: xcov(ip1jmp1, klevel), ycov(ip1jm, klevel)
+  REAL :: grx(ip1jmp1, klevel), gry(ip1jm, klevel)
   !
   !    ......   variables locales     ........
   !
-  REAL :: rot(ip1jm,llm) , signe, nugradrs
-  INTEGER :: l,ij,iter,lr
+  REAL :: rot(ip1jm, llm), signe, nugradrs
+  INTEGER :: l, ij, iter, lr
   !    ........................................................
   !
   !
   !
-  signe    = (-1.)**lr
+  signe = (-1.)**lr
   nugradrs = signe * crot
   !
-  CALL SCOPY ( ip1jmp1* klevel, xcov, 1, grx, 1 )
-  CALL SCOPY (  ip1jm * klevel, ycov, 1, gry, 1 )
+  CALL SCOPY (ip1jmp1 * klevel, xcov, 1, grx, 1)
+  CALL SCOPY (ip1jm * klevel, ycov, 1, gry, 1)
   !
-  CALL     rotatf     ( klevel, grx, gry, rot )
+  CALL     rotatf     (klevel, grx, gry, rot)
   !
-  CALL laplacien_rot ( klevel, rot, rot,grx,gry      )
+  CALL laplacien_rot (klevel, rot, rot, grx, gry)
 
   !
   !    .....   Iteration de l'operateur laplacien_rotgam  .....
   !
-  DO  iter = 1, lr -2
-    CALL laplacien_rotgam ( klevel, rot, rot )
+  DO  iter = 1, lr - 2
+    CALL laplacien_rotgam (klevel, rot, rot)
   ENDDO
   !
   !
-  CALL filtreg( rot, jjm, klevel, 2,1, .FALSE.,1)
-  CALL nxgrad ( klevel, rot, grx, gry )
+  CALL filtreg(rot, jjm, klevel, 2, 1, .FALSE., 1)
+  CALL nxgrad (klevel, rot, grx, gry)
   !
   DO    l = 1, klevel
-     DO  ij = 1, ip1jm
-      gry( ij,l ) = gry( ij,l ) * nugradrs
-     ENDDO
-     DO  ij = 1, ip1jmp1
-      grx( ij,l ) = grx( ij,l ) * nugradrs
-     ENDDO
+    DO  ij = 1, ip1jm
+      gry(ij, l) = gry(ij, l) * nugradrs
+    ENDDO
+    DO  ij = 1, ip1jmp1
+      grx(ij, l) = grx(ij, l) * nugradrs
+    ENDDO
   ENDDO
   !