Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/extrapol.f90

Timestamp:

Aug 2, 2024, 9:58:25 PM (3 months ago)

Author:

abarral

Message:

Put dimensions.h and paramet.h into modules

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3d_common/extrapol.f90 (modified) (14 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3d_common/extrapol.f90

@@ -2 +2 @@
 ! $Id$
 
-!
-!
+
+
 SUBROUTINE extrapol(pfild, kxlon, kylat, pmask, &
         norsud, ldper, knbor, pwork)
   IMPLICIT NONE
-  !
+
   ! OASIS routine (Adaptation: Laurent Li, le 14 mars 1997)
   ! Fill up missed values by using the neighbor points
-  !
+
   INTEGER :: kxlon, kylat ! longitude and latitude dimensions (Input)
   INTEGER :: knbor ! minimum neighbor number (Input)
@@ -18 +18 @@
   REAL :: pfild(kxlon,kylat) ! field to be extrapolated (Input/Output)
   REAL :: pwork(kxlon,kylat) ! working space
-  !
+
   REAL :: zwmsk
   INTEGER :: incre, idoit, i, j, k, inbor, ideb, ifin, ilon, jlat
   INTEGER :: ix(9), jy(9) ! index arrays for the neighbors coordinates
   REAL :: zmask(9)
-  !
+
   !  We search over the eight closest neighbors
-  !
+
   !        j+1  7  8  9
   !          j  4  5  6    Current point 5 --> (i,j)
   !        j-1  1  2  3
   !            i-1 i i+1
-  !
-  !
+
+
   IF (norsud) THEN
      DO j = 1, kylat
@@ -44 +44 @@
      ENDDO
   ENDIF
-  !
+
   incre = 0
-  !
+
   DO j = 1, kylat
   DO i = 1, kxlon
@@ -52 +52 @@
   ENDDO
   ENDDO
-  !
+
   !* To avoid problems in floating point tests
   zwmsk = pmask - 1.0
-  !
+
 200   CONTINUE
   incre = incre + 1
@@ -65 +65 @@
      ideb = 1
      ifin = 9
-  !
+
   !* Fill up ix array
      ix(1) = MAX (1,i-1)
@@ -76 +76 @@
      ix(8) = i
      ix(9) = MIN (kxlon,i+1)
-  !
+
   !* Fill up iy array
      jy(1) = MAX (1,j-1)
@@ -87 +87 @@
      jy(8) = MIN (kylat,j+1)
      jy(9) = MIN (kylat,j+1)
-  !
+
   !* Correct latitude bounds if southernmost or northernmost points
      IF (j == 1) ideb = 4
      IF (j == kylat) ifin = 6
-  !
+
   !* Account for periodicity in longitude
-  !
+
      IF (ldper) THEN
         IF (i == kxlon) THEN
@@ -121 +121 @@
            ix(6) = i
         ENDIF
-  !
+
         IF (i == 1 .OR. i == kxlon) THEN
            jy(1) = MAX (1,j-1)
@@ -129 +129 @@
            jy(5) = MIN (kylat,j+1)
            jy(6) = MIN (kylat,j+1)
-  !
+
            ideb = 1
            ifin = 6
@@ -136 +136 @@
         ENDIF
      ENDIF ! end for ldper test
-  !
+
   !* Find unmasked neighbors
-  !
+
      DO k = ideb, ifin
         zmask(k) = 0.
@@ -148 +148 @@
         ENDIF
   END DO
-  !
+
   !* Not enough points around point P are unmasked; interpolation on P
   !  will be done in a future CALL to extrap.
-  !
+
      IF (inbor >= knbor) THEN
         pwork(i,j) = 0.
@@ -161 +161 @@
         ENDDO
      ENDIF
-  !
+
   ENDIF
   END DO
   END DO
-  !
+
   !*    3. Writing back unmasked field in pfild
   !    ------------------------------------
-  !
+
   !* pfild then contains:
   ! - Values which were not masked
   ! - Interpolated values from the inbor neighbors
   ! - Values which are not yet interpolated
-  !
+
   idoit = 0
   DO j = 1, kylat
@@ -181 +181 @@
   ENDDO
   ENDDO
-  !
+
   IF (idoit /= 0) GOTO 200
   !cc      PRINT*, "Number of extrapolation steps incre =", incre
-  !
+
   IF (norsud) THEN
      DO j = 1, kylat
@@ -197 +197 @@
      ENDDO
   ENDIF
-  !
+
   RETURN
 END SUBROUTINE extrapol