Changeset 5246 for LMDZ6/trunk/libf/dyn3d_common/gr_int_dyn.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/dyn3d_common/gr_int_dyn.f90 (moved) (moved from LMDZ6/trunk/libf/dyn3d_common/gr_int_dyn.F) (1 diff)

LMDZ6/trunk/libf/dyn3d_common/gr_int_dyn.f90

@@ -2 +2 @@
 ! $Header$
 !
-      subroutine gr_int_dyn(champin,champdyn,iim,jp1)
-      implicit none
-c=======================================================================
-c   passage d'un champ interpole a un champ sur grille scalaire
-c=======================================================================
-c-----------------------------------------------------------------------
-c   declarations:
-c   -------------
+subroutine gr_int_dyn(champin,champdyn,iim,jp1)
+  implicit none
+  !=======================================================================
+  !   passage d'un champ interpole a un champ sur grille scalaire
+  !=======================================================================
+  !-----------------------------------------------------------------------
+  !   declarations:
+  !   -------------
 
-      INTEGER iim
-      integer ip1, jp1
-      REAL champin(iim, jp1)
-      REAL champdyn(iim+1, jp1)
+  INTEGER :: iim
+  integer :: ip1, jp1
+  REAL :: champin(iim, jp1)
+  REAL :: champdyn(iim+1, jp1)
 
-      INTEGER i, j
-      real polenord, polesud
+  INTEGER :: i, j
+  real :: polenord, polesud
 
-c-----------------------------------------------------------------------
-c   calcul:
-c   -------
+  !-----------------------------------------------------------------------
+  !   calcul:
+  !   -------
 
-      ip1 = iim + 1
-      polenord = 0.
-      polesud = 0.
-      do i = 1, iim
-        polenord = polenord + champin (i, 1)
-        polesud = polesud + champin (i, jp1)
-      enddo
-      polenord = polenord / iim
-      polesud = polesud / iim
-      do j = 1, jp1
-        do i = 1, iim
-          if (j .eq. 1) then
-            champdyn(i, j) = polenord
-          else if (j .eq. jp1) then
-            champdyn(i, j) = polesud
-          else
-            champdyn(i, j) = champin (i, j)
-          endif
-        enddo
-        champdyn(ip1, j) = champdyn(1, j)
-      enddo
+  ip1 = iim + 1
+  polenord = 0.
+  polesud = 0.
+  do i = 1, iim
+    polenord = polenord + champin (i, 1)
+    polesud = polesud + champin (i, jp1)
+  enddo
+  polenord = polenord / iim
+  polesud = polesud / iim
+  do j = 1, jp1
+    do i = 1, iim
+      if (j .eq. 1) then
+        champdyn(i, j) = polenord
+      else if (j .eq. jp1) then
+        champdyn(i, j) = polesud
+      else
+        champdyn(i, j) = champin (i, j)
+      endif
+    enddo
+    champdyn(ip1, j) = champdyn(1, j)
+  enddo
 
-      RETURN
-      END
+  RETURN
+END SUBROUTINE gr_int_dyn