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

LMDZ6/trunk/libf/dyn3d_common/coordij.f90

@@ -2 +2 @@
 ! $Header$
 !
-      SUBROUTINE coordij(lon,lat,ilon,jlat)
+SUBROUTINE coordij(lon,lat,ilon,jlat)
 
-c=======================================================================
-c
-c   calcul des coordonnees i et j de la maille scalaire dans
-c   laquelle se trouve le point (lon,lat) en radian
-c
-c=======================================================================
+  !=======================================================================
+  !
+  !   calcul des coordonnees i et j de la maille scalaire dans
+  !   laquelle se trouve le point (lon,lat) en radian
+  !
+  !=======================================================================
 
-      USE comconst_mod, ONLY: pi
-      
-      IMPLICIT NONE
-      REAL lon,lat
-      INTEGER ilon,jlat
-      INTEGER i,j
+  USE comconst_mod, ONLY: pi
 
-      include "dimensions.h"
-      include "paramet.h"
-      include "comgeom.h"
+  IMPLICIT NONE
+  REAL :: lon,lat
+  INTEGER :: ilon,jlat
+  INTEGER :: i,j
 
-      real zlon,zlat
+  include "dimensions.h"
+  include "paramet.h"
+  include "comgeom.h"
 
-      zlon=lon*pi/180.
-      zlat=lat*pi/180.
+  real :: zlon,zlat
 
-      DO i=1,iim+1
-         IF (rlonu(i).GT.zlon) THEN
-            ilon=i
-            GOTO 10
-         ENDIF
-      ENDDO
-10    CONTINUE
+  zlon=lon*pi/180.
+  zlat=lat*pi/180.
 
-      j=0
-      DO j=1,jjm
-         IF(rlatv(j).LT.zlat) THEN
-            jlat=j
-            GOTO 20
-         ENDIF
-      ENDDO
-20    CONTINUE
-      IF(j.EQ.0) j=jjm+1
+  DO i=1,iim+1
+     IF (rlonu(i).GT.zlon) THEN
+        ilon=i
+        GOTO 10
+     ENDIF
+  ENDDO
+10   CONTINUE
 
-      RETURN
-      END
+  j=0
+  DO j=1,jjm
+     IF(rlatv(j).LT.zlat) THEN
+        jlat=j
+        GOTO 20
+     ENDIF
+  ENDDO
+20   CONTINUE
+  IF(j.EQ.0) j=jjm+1
+
+  RETURN
+END SUBROUTINE coordij