Changeset 2725 in lmdz_wrf

Timestamp:

Oct 17, 2019, 1:50:41 PM (5 years ago)

Author:

lfita

Message:

Adding missing value in index_samevals1D_RK

File:

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

trunk/tools/module_generic.f90

@@ -1235 +1235 @@
   END SUBROUTINE rm_values_vecRK
 
-  SUBROUTINE index_samevals1D_RK(d1r, refv, d1v, vals, ii, indices, samevalues)
+  SUBROUTINE index_samevals1D_RK(d1r, refv, d1v, vals, ii, indices, missval, samevalues)
   ! Subroutine to search for the indices of the same values between 2 1D RK series of values allowing 
   !   repetitions
@@ -1242 +1242 @@
 
     INTEGER, INTENT(in)                                  :: d1r, d1v
+    REAL(r_k), INTENT(in)                                :: missval
     REAL(r_k), DIMENSION(d1r), INTENT(in)                :: refv
     REAL(r_k), DIMENSION(d1v), INTENT(in)                :: vals
@@ -1258 +1259 @@
 ! ii: quantity of same values found
 ! indices: output
+! missval: missing value
 ! samevalues: values where coincidence is found
 
@@ -1267 +1269 @@
     DO ir=1, d1r
       DO iv=1, d1v
-        IF (vals(iv) == refv(ir)) THEN
+        IF ( (refv(ir) /= missval) .AND. (vals(iv) == refv(ir)) )THEN
           ii = ii + 1
           indices(ii,:) = (/ ir, iv /)