Changeset 1653 in lmdz_wrf

Timestamp:

Sep 21, 2017, 8:20:40 PM (7 years ago)

Author:

lfita

Message:

Adding:

• `Index2DArrayR_K': Function to provide the first index of a given value inside a 2D real(r_k) array

File:

• trunk/tools/module_generic.f90 (modified) (4 diffs)

trunk/tools/module_generic.f90

@@ -10 +10 @@
 ! Index1DArrayR_K: Function to provide the first index of a given value inside a 1D real(r_k) array
 ! Index2DArrayR: Function to provide the first index of a given value inside a 2D real array
+! Index2DArrayR_K: Function to provide the first index of a given value inside a 2D real(r_k) array
 ! mat2DPosition: Function to provide the i, j indices of a given value inside a 2D matrix
 ! Nstrings: Function to repeat a number of times a given string
@@ -139 +140 @@
     IMPLICIT NONE
 
-    INTEGER, PARAMETER                                   :: r_k = KIND(1.d0)
-    INTEGER, INTENT(in)                                  :: d1
-    REAL(r_k), INTENT(in)                                :: val
-    REAL(r_k), DIMENSION(d1), INTENT(in)                 :: array1D
+    INTEGER, INTENT(in)                                  :: d1
+    REAL, INTENT(in)                                     :: val
+    REAL, DIMENSION(d1), INTENT(in)                      :: array1D
 
 ! Local
@@ -190 +190 @@
     IMPLICIT NONE
 
-    INTEGER, PARAMETER                                   :: r_k = KIND(1.d0)
     INTEGER, INTENT(in)                                  :: d1, d2
-    REAL(r_k), INTENT(in)                                :: val
-    REAL(r_k), DIMENSION(d1,d2), INTENT(in)              :: array2D
+    REAL, INTENT(in)                                     :: val
+    REAL, DIMENSION(d1,d2), INTENT(in)                   :: array2D
     INTEGER, DIMENSION(2)                                :: Index2DArrayR
 
@@ -214 +213 @@
 
   END FUNCTION Index2DArrayR
+
+  FUNCTION Index2DArrayR_K(array2D, d1, d2, val)
+! Function to provide the first index of a given value inside a 2D real array
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: d1, d2
+    REAL(r_k), INTENT(in)                                :: val
+    REAL(r_k), DIMENSION(d1,d2), INTENT(in)              :: array2D
+    INTEGER, DIMENSION(2)                                :: Index2DArrayR_K
+
+! Local
+    INTEGER                                              :: i, j
+
+    fname = 'Index2DArrayR_K'
+
+    Index2DArrayR_K = -1
+
+    DO i=1,d1
+      DO j=1,d2
+        IF (array2d(i,j) == val) THEN
+          Index2DArrayR_K(1) = i
+          Index2DArrayR_K(2) = j
+          EXIT
+        END IF
+      END DO
+    END DO
+
+  END FUNCTION Index2DArrayR_K
 
   FUNCTION RangeI(d1, iniv, endv)