+#!/usr/bin/python
+
+""" meanstep.py
+ ...
+"""
+
+import sys
+import os
+
+results = {}
+
+if len(sys.argv) != 4:
+ print "usage: ./meanstep.py <datadir> <matching col 1>,...,<matching col n> <mean col 1>,...,<mean col n>"
+ exit(0)
+
+matchcols = sys.argv[2].split(",")
+
+for root, dirs, filenames in os.walk(sys.argv[1]):
+ for f in filenames:
+ if f[-4:] != ".dat":
+ continue
+
+ with open(root + "/" + f) as inf:
+ for line in inf:
+ if line[0:1] == "#":
+ continue
+
+ linecols = line.split()
+ matchkey = ""
+ for imatchcol in range(0,len(matchcols)):
+ matchkey = matchkey + linecols[int(matchcols[imatchcol])]
+
+
+
+ #results[f[:-4]] = "tefds"