Changeset 1918 in lmdz_wrf

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

Location:

trunk/tools

Files:

• generic.py (modified) (3 diffs)
• generic_tools.py (modified) (2 diffs)

trunk/tools/generic.py

@@ -21 +21 @@
 import re
 
+# ASCIIfile_stats: Function to provide the statistics of a series of values from an ASCII file
 # coincident_CFtimes: Function to make coincident times for two different sets of CFtimes
 # count_cond: Function to count values of a variable which attain a condition
@@ -49 +50 @@
 
 # List of available operations
-operations=['coincident_CFtimes', 'count_cond', 'datetimeStr_conversion',            \
+operations=['ASCIIfile_stats', 'coincident_CFtimes', 'count_cond',                   \
+  'datetimeStr_conversion',                                                          \
   'grid_combinations',                                                               \
   'interpolate_locs', 'latex_fig_array', 'list_operations', 'PolyArea',              \
@@ -423 +425 @@
         print gen.WRFsetup(vals0, vals1, vals2)
 
+elif oper == 'ASCIIfile_stats':
+    Nvals = 3
+    vals = opts.values.split(cS)
+    if vals[0] == 'h':
+        print gen.ASCIIfile_stats.__doc__
+        quit(-1)
+    else:
+        if len(vals) != Nvals:
+            print errormsg
+            print '  ' + main + ": operation '" + oper + "' requires", Nvals, 'and', \
+              len(vals), ' has passed!!'
+            print gen.ASCIIfile_stats.__doc__
+            quit(-1)
+        vals0 = vals[0]
+        vals1 = vals[1]
+        vals2 = vals[2]
+        
+        print gen.ASCIIfile_stats(vals0, vals1, vals2)
+

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()