import sys inp = open(sys.argv[1], 'r') out = open(sys.argv[2], 'w') def getentry(): line = inp.readline() while line.startswith('#'): line = inp.readline() for i in range(len(line)): if line[i] in [' ', '\t']: break else: raise Exception return float(line[:i]), float(line[i:]) def writeentry(x, y): out.write("%f\t%f\n" % (x, y)) entries = [] while True: try: entries.append(getentry()) except Exception: break if len(sys.argv) == 4: avg = float(sys.argv[3]) else: avg = sum([e[1] for e in entries]) / len(entries) for e in entries: writeentry(e[0], e[1] - avg)