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

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

@@ -1 @@
-
 ! $Header$
 
-SUBROUTINE nxgrarot(klevel,xcov, ycov, lr, grx, gry )
+SUBROUTINE nxgrarot(klevel, xcov, ycov, lr, grx, gry)
   !   ***********************************************************
   !
@@ -17 +16 @@
   !
   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)
   !
-  REAL :: rot(ip1jm,llm)
+  REAL :: rot(ip1jm, llm)
 
-  INTEGER :: l,ij,iter,lr
+  INTEGER :: l, ij, iter, lr
   !
   !
   !
-  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)
   !
-  DO iter = 1,lr
-  CALL  rotat (klevel,grx, gry, rot )
-  CALL filtreg( rot, jjm, klevel, 2,1, .FALSE.,2)
-  CALL nxgrad (klevel,rot, grx, gry )
-  !
-  DO l = 1, klevel
-  DO ij = 1, ip1jm
-  gry( ij,l ) = - gry( ij,l ) * crot
-  END DO
-  DO ij = 1, ip1jmp1
-  grx( ij,l ) = - grx( ij,l ) * crot
-  END DO
-  END DO
-  !
+  DO iter = 1, lr
+    CALL  rotat (klevel, grx, gry, rot)
+    CALL filtreg(rot, jjm, klevel, 2, 1, .FALSE., 2)
+    CALL nxgrad (klevel, rot, grx, gry)
+    !
+    DO l = 1, klevel
+      DO ij = 1, ip1jm
+        gry(ij, l) = - gry(ij, l) * crot
+      END DO
+      DO ij = 1, ip1jmp1
+        grx(ij, l) = - grx(ij, l) * crot
+      END DO
+    END DO
+    !
   END DO
   RETURN