Changeset 2170 in lmdz_wrf for trunk/tools

Timestamp:

Oct 8, 2018, 3:19:17 PM (7 years ago)

Author:

lfita

Message:

Adding percentiles to `stats_time2D'

File:

• trunk/tools/generic_tools.py (modified) (5 diffs)

trunk/tools/generic_tools.py

@@ -5616 +5616 @@
         self.stdv[t]: temporal standard deviation at each grid point
         self.anomv[t]: temporal anomaly from the whole mean at each grid point
+        self.precen[t]: temporal percentiles at each grid point
     """
   
@@ -5633 +5634 @@
             self.mean2v = None
             self.stdv = None
+            self.percen = None
         else:
             dimx = vals.shape[2]
@@ -5638 +5640 @@
             stats=np.zeros((dimy,dimx,5), dtype=np.float)
             absmean = np.mean(vals,axis=0)
+            pecents = np.zeros((20,dimy,dimx), dtype=np.float)
 
             stats[:,:,0]=np.min(vals, axis=0)
@@ -5649 +5652 @@
             stats = np.where(stats is np.inf, fillValue, stats)
             stats = np.where(stats is None, fillValue, stats)
+
+            for j in range(dimy):
+                for i in range(dimx):
+                    percents[:,j,i] = Quantiles(vals[:,j,i],20).quantilesv
 
             self.minv=stats[:,:,0]
@@ -5656 +5663 @@
             self.stdv=np.sqrt(stats[:,:,3]-stats[:,:,2]*stats[:,:,2])
             self.anomv=stats[:,:,4]
+            self.percen=percents
 
         return