Changeset 1913 in lmdz_wrf for trunk

Timestamp:

May 14, 2018, 11:52:50 AM (7 years ago)

Author:

lfita

Message:

Adding:

• `NcountR': Subroutine to count real values

File:

• trunk/tools/module_scientific.f90 (modified) (2 diffs)

trunk/tools/module_scientific.f90

@@ -21 +21 @@
 ! look_clockwise_borders: Subroutine to look clock-wise for a next point within a collection of borders 
 !   (limits of a region)
+! NcountR: Subroutine to count real values
 ! paths_border: Subroutine to search the paths of a border field.
 ! path_properties: Subroutine to determine the properties of a path
@@ -3572 +3573 @@
 END SUBROUTINE StatsR_k
 
+  SUBROUTINE NcountR(values, d1, Ndiffvals, counts)
+! Subroutine to count real values
+
+    IMPLICIT NONE
+
+    INTEGER, INTENT(in)                                  :: d1
+    REAL(r_k), DIMENSION(d1), INTENT(in)                 :: values
+    INTEGER, INTENT(out)                                 :: Ndiffvals
+    REAL(r_k), DIMENSION(d1,2), INTENT(out)              :: counts
+! Local
+    INTEGER                                              :: i, ival
+    REAL(r_k), DIMENSION(d1)                             :: diffv
+
+!!!!!!! Variables
+! values: values to count
+! counts: counts of time for each value
+    
+    fname = 'NcountR'
+
+    counts = -1.
+
+    counts(1,1) = values(1)
+    counts(1,2) = 1
+    Ndiffvals = 1
+    DO i=2,d1
+      diffv(1:Ndiffvals) = counts(1:Ndiffvals,1) - values(i)
+      IF (ANY(diffv(1:Ndiffvals) == 0)) THEN
+        ival = Index1DArrayR(counts(1:Ndiffvals,1), Ndiffvals, values(i))
+        counts(ival,2) = counts(ival,2) + 1
+      ELSE
+        Ndiffvals = Ndiffvals + 1
+        counts(Ndiffvals,1) = values(i)
+        counts(Ndiffvals,2) = 1
+      END IF
+    END DO
+
+  END SUBROUTINE NcountR
+
 END MODULE module_scientific