Changeset 1918 in lmdz_wrf for trunk/tools/generic_tools.py

Timestamp:

Jun 18, 2018, 4:09:24 PM (7 years ago)

Author:

lfita

Message:

Adding:

• 'ASCIIfile_stats': Function to provide the statistics of a series of values from an ASCII file

File:

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

trunk/tools/generic_tools.py

@@ -52 +52 @@
 # advance_matDate: Function to advance matrix-date with a matrix-increment
 # angle_DegMinSec: Function to transform an angle to Degrees Minutes Seconds
+# ASCIIfile_stats: Function to provide the statistics of a series of values from an ASCII file
 # ASCII_HTML: Function to transform from an ASCII line to HTML codification
 # ASCII_LaTeX: Function to transform from an ASCII character to LaTeX codification
@@ -13264 +13265 @@
     return minNCvalues
 
+def ASCIIfile_stats(filen, comment='#', kind='R'):
+    """ Function to provide the statistics of a series of values from an ASCII file
+      filen= name of the file
+      comment= ':' separated list of comment characters
+      kind= kind of values
+        'I': integer
+        'R': float
+        'D': double
+    """
+    fname = 'ACIIfile_stats'
+
+    if not os.path.isfile(filen):
+        print errormsg
+        print '  ' + fname + ": ASCII file '" + filen + "' does not exist !!"
+        quit(-1)
+
+    comments = str_list(comment, ':')
+
+    objf = open(filen, 'r')
+    values = []
+    for line in objf:
+        linev = line.replace('\n','').replace('\t','').replace('\r','')
+        if len(linev) >= 1:
+            ch1 = linev[0:1]
+            if not searchInlist(comments, ch1):
+                print "'" + linev + "'"
+                values.append(typemod(linev,kind))
+
+    objf.close()
+
+    matv = np.array(values)
+    # Statistics
+    minv = np.min(matv)
+    maxv = np.max(matv)
+    meanv = np.mean(matv)
+    stdv = np.std(matv)
+    quantv = Quantiles(matv, 20)
+    
+    print "# Statistics of '" + filen + "' _______"
+    print '#@# min:', minv
+    print '#@# max:', maxv
+    print '#@# mean:', meanv
+    print '#@# stddev:', stdv
+    print '#@# quantiles:', quantv.quantilesv
+
+    quantS = []
+    for iq in quantv.quantilesv: quantS.append(str(iq))
+
+    print '#@# allvalues; min:' + str(minv) + ';max:' + str(maxv) + ';mean:' +       \
+      str(meanv) + ';stddev:' + str(stdv) + ';quant:' + ','.join(quantS)
+
+    return
+
+#ASCIIfile_stats('times_NOCDXWRF2.dat')
+
 #quit()