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

LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/gr_fi_dyn.f90

@@ -1 @@
-
 ! $Header$
 
-SUBROUTINE gr_fi_dyn(nfield,ngrid,im,jm,pfi,pdyn)
+SUBROUTINE gr_fi_dyn(nfield, ngrid, im, jm, pfi, pdyn)
+  USE lmdz_ssum_scopy, ONLY: scopy
+
   IMPLICIT NONE
   !=======================================================================
@@ -12 +13 @@
   !   -------------
 
-  INTEGER :: im,jm,ngrid,nfield
-  REAL :: pdyn(im,jm,nfield)
-  REAL :: pfi(ngrid,nfield)
+  INTEGER :: im, jm, ngrid, nfield
+  REAL :: pdyn(im, jm, nfield)
+  REAL :: pfi(ngrid, nfield)
 
-  INTEGER :: i,j,ifield,ig
+  INTEGER :: i, j, ifield, ig
 
   !-----------------------------------------------------------------------
@@ -22 +23 @@
   !   -------
 
-  DO ifield=1,nfield
-  !   traitement des poles
-     DO i=1,im
-        pdyn(i,1,ifield)=pfi(1,ifield)
-        pdyn(i,jm,ifield)=pfi(ngrid,ifield)
-     ENDDO
+  DO ifield = 1, nfield
+    !   traitement des poles
+    DO i = 1, im
+      pdyn(i, 1, ifield) = pfi(1, ifield)
+      pdyn(i, jm, ifield) = pfi(ngrid, ifield)
+    ENDDO
 
-  !   traitement des point normaux
-     DO j=2,jm-1
-        ig=2+(j-2)*(im-1)
-        CALL SCOPY(im-1,pfi(ig,ifield),1,pdyn(1,j,ifield),1)
-        pdyn(im,j,ifield)=pdyn(1,j,ifield)
-     ENDDO
+    !   traitement des point normaux
+    DO j = 2, jm - 1
+      ig = 2 + (j - 2) * (im - 1)
+      CALL SCOPY(im - 1, pfi(ig, ifield), 1, pdyn(1, j, ifield), 1)
+      pdyn(im, j, ifield) = pdyn(1, j, ifield)
+    ENDDO
   ENDDO