Changeset 5246 for LMDZ6/trunk/libf/dynphy_lonlat/gr_dyn_fi.f90

Timestamp:

Oct 21, 2024, 2:58:45 PM (23 hours ago)

Author:

abarral

Message:

Convert fixed-form to free-form sources .F -> .{f,F}90
(WIP: some .F remain, will be handled in subsequent commits)

File:

• LMDZ6/trunk/libf/dynphy_lonlat/gr_dyn_fi.f90 (moved) (moved from LMDZ6/trunk/libf/dynphy_lonlat/gr_dyn_fi.F) (1 diff)

LMDZ6/trunk/libf/dynphy_lonlat/gr_dyn_fi.f90

@@ -2 +2 @@
 ! $Header$
 !
-      SUBROUTINE gr_dyn_fi(nfield,im,jm,ngrid,pdyn,pfi)
-      IMPLICIT NONE
-c=======================================================================
-c   passage d'un champ de la grille scalaire a la grille physique
-c=======================================================================
+SUBROUTINE gr_dyn_fi(nfield,im,jm,ngrid,pdyn,pfi)
+  IMPLICIT NONE
+  !=======================================================================
+  !   passage d'un champ de la grille scalaire a la grille physique
+  !=======================================================================
 
-c-----------------------------------------------------------------------
-c   declarations:
-c   -------------
+  !-----------------------------------------------------------------------
+  !   declarations:
+  !   -------------
 
-      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
 
-c-----------------------------------------------------------------------
-c   calcul:
-c   -------
+  !-----------------------------------------------------------------------
+  !   calcul:
+  !   -------
 
-      IF (ngrid.NE.2+(jm-2)*(im-1)) then
-         call abort_gcm("gr_dyn_fi", 'probleme de dim', 1)
-      end if
-c   traitement des poles
-      CALL SCOPY(nfield,pdyn,im*jm,pfi,ngrid)
-      CALL SCOPY(nfield,pdyn(1,jm,1),im*jm,pfi(ngrid,1),ngrid)
+  IF (ngrid.NE.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)
 
-c   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
-      ENDDO
+  !   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
+  ENDDO
 
-      RETURN
-      END
+  RETURN
+END SUBROUTINE gr_dyn_fi