Changeset 2210 in lmdz_wrf for trunk/tools

Timestamp:

Nov 2, 2018, 4:04:02 PM (6 years ago)

Author:

lfita

Message:

Adding:

• `runmean_F1D': Subroutine fo computing the running mean of a given set of float 1D values

File:

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

trunk/tools/module_scientific.f90

@@ -41 +41 @@
 ! reconstruct_matrix: Subroutine to reconstruct a 2D matrix from a pair of syncronized vectors with the positions on x 
 !   and y coordinates
+! runmean_F1D: Subroutine fo computing the running mean of a given set of float 1D values
 ! SortR_K*: Subroutine receives an array x() r_K and sorts it into ascending order.
 ! StatsR_K: Subroutine to provide the minmum, maximum, mean, the quadratic mean, and the standard deviation of a 
@@ -3611 +3612 @@
   END SUBROUTINE NcountR
 
+  SUBROUTINE runmean_F1D(d1, values, Nmean, headertail, runmean)
+! Subroutine fo computing the running mean of a given set of float 1D values
+
+  IMPLICIT NONE
+
+  INTEGER, INTENT(in)                                    :: d1, Nmean
+  REAL(r_k), DIMENSION(d1), INTENT(in)                   :: values
+  CHARACTER(len=*), INTENT(in)                           :: headertail
+  REAL(r_k), DIMENSION(d1), INTENT(out)                  :: runmean
+  
+! Local
+  INTEGER                                                :: i, j, Nmean2
+  CHARACTER(len=5)                                       :: NmeanS
+
+!!!!!!! Variables
+! values: values to compute the running mean
+! Nmean: number of odd points to use for the running mean
+! headertail: How to proceed for the grid points at the beginning of the values which are not 
+!   encompassed by the Nmean
+!   'miss': set as missing values (1.d20)
+!   'original': keep the original values
+!   'progressfill': mean the values as a progressive running filter (e.g. for Nmean=5):
+!     runmean[values(1)] = values(1)
+!     runmean[values(2)] = MEAN(values(1:3))
+!     runmean[values(3)] = MEAN(values(1:5))
+!     runmean[values(4)] = MEAN(values(2:6))
+!     (...)
+!     runmean[values(d1-2)] = MEAN(values(d1-5:d1))
+!     runmean[values(d1-1)] = MEAN(values(d1-2:d1))
+!     runmean[values(d1)] = MEAN(values(dd1))
+!   'zero': set as zero values
+! runmean: runnig mean values
+
+  fname = 'runmean_F1D'
+
+  IF (MOD(Nmean,2) == 0) THEN
+    WRITE(NmeanS,'(I5)')Nmean
+    msg="Nmean has to be odd!! value provided: "// NmeanS
+    CALL ErrMsg(msg, fname, -1)
+  END IF
+  Nmean2 = Nmean/2
+ 
+  SELECT CASE (TRIM(headertail))
+    CASE ('missing')
+      runmean = fillval64
+    CASE ('original')
+      runmean = values
+    CASE ('progressfill')
+      DO i=1, Nmean2
+        runmean(i) = SUM(values(1:2*(i-1)+1))/(2*(i-1)+1)
+      END DO
+      runmean(d1) = values(d1)
+      DO i=2, Nmean2
+        j = d1-(2*(i-1))
+        runmean(d1-(i-1)) = SUM(values(j:d1))/(2*(i-1)+1)
+      END DO
+    CASE ('zero')
+      runmean = zeroRK
+    CASE DEFAULT
+      msg = "'" // TRIM(headertail) // "' not available !!" //CHAR(44) // "  available ones: " //     &
+        "'missing', 'original', 'progressfill', 'zero'"
+      CALL ErrMsg(msg, fname, -1)
+  END SELECT
+
+  DO i= 1+Nmean2, d1 - Nmean2
+    runmean(i) = SUM(values(i-Nmean2:i+Nmean2))/Nmean
+  END DO
+
+  END SUBROUTINE runmean_F1D
+
 END MODULE module_scientific