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

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

@@ -1 @@
-
 ! $Header$
 
-SUBROUTINE gr_dyn_fi(nfield,im,jm,ngrid,pdyn,pfi)
+SUBROUTINE gr_dyn_fi(nfield, im, jm, ngrid, pdyn, pfi)
+  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 :: j,ifield,ig
+  INTEGER :: j, ifield, ig
 
   !-----------------------------------------------------------------------
@@ -22 +23 @@
   !   -------
 
-  IF (ngrid/=2+(jm-2)*(im-1)) THEN
-     CALL abort_gcm("gr_dyn_fi", 'probleme de dim', 1)
+  IF (ngrid/=2 + (jm - 2) * (im - 1)) THEN
+    CALL abort_gcm("gr_dyn_fi", 'probleme de dim', 1)
   end if
   !   traitement des poles
-  CALL SCOPY(nfield,pdyn,im*jm,pfi,ngrid)
-  CALL SCOPY(nfield,pdyn(1,jm,1),im*jm,pfi(ngrid,1),ngrid)
+  CALL SCOPY(nfield, pdyn, im * jm, pfi, ngrid)
+  CALL SCOPY(nfield, pdyn(1, jm, 1), im * jm, pfi(ngrid, 1), ngrid)
 
   !   traitement des point normaux
-  DO ifield=1,nfield
-     DO j=2,jm-1
-        ig=2+(j-2)*(im-1)
-        CALL SCOPY(im-1,pdyn(1,j,ifield),1,pfi(ig,ifield),1)
-     ENDDO
+  DO ifield = 1, nfield
+    DO j = 2, jm - 1
+      ig = 2 + (j - 2) * (im - 1)
+      CALL SCOPY(im - 1, pdyn(1, j, ifield), 1, pfi(ig, ifield), 1)
+    ENDDO
   ENDDO