Changeset 1265 in lmdz_wrf

Timestamp:

Nov 3, 2016, 3:09:35 PM (8 years ago)

Author:

lfita

Message:

Making use of the numpy variance' calculation in field_stats' and `field_stats_dim'
Fixing documentation of both functions

File:

• trunk/tools/nc_var_tools.py (modified) (4 diffs)

trunk/tools/nc_var_tools.py

@@ -10936 +10936 @@
     """ Function to retrieve statistics from a field
     field_stats(values, ncfile, varn)
-      [values]= [stats],[fillVals]
+      [values]= [stats],[fillVals],[countVals]
         [stats]: kind of statistics 
           'full': all statistics given variable
@@ -11022 +11022 @@
             meanv = np.mean(field)
             mean2v = np.mean(field*field)
-            varv = np.sqrt(mean2v - meanv*meanv)
+            #varv = np.sqrt(mean2v - meanv*meanv)
+            varv = np.var(field)
 
             print '    Fstats shape:',field.shape
@@ -11063 +11064 @@
     """ Function to retrieve statistics from a field along dimensions
     field_stats_dim(values, ncfile, varn)
-      [values]= [stats],[fillVals],[dimns],[vardims],[stdout]
+      [values]= [stats],[fillVals],[countVals],[dimns],[vardims],[stdout]
         [stats]: kind of statistics 
           'full': all statistics given variable (min, max, mean, mean2, variability, [countVals])
@@ -11255 +11256 @@
                 meanv = np.mean(field)
                 mean2v = np.mean(field*field)
-                varv = np.sqrt(mean2v - meanv*meanv)
+                #varv = np.sqrt(mean2v - meanv*meanv)
+                varv = np.var(field)
 
                 if stdout: