old mode 100644
new mode 100755
copy from medcouple.py
copy to pymedcouple
--- a/medcouple.py
+++ b/pymedcouple
@@ -1,3 +1,5 @@
+#!/usr/bin/python
+
 import random
 import numpy as np
 
@@ -115,7 +117,7 @@
 
     # Scale inside [-0.5, 0.5], for greater numerical stability.
     Zden = -2*max(-Z[0], Z[-1])
-    Z = [z/zden for z in Zden]
+    Z = [z/Zden for z in Z]
     Zmed /= Zden
     
     Zeps = eps1*(eps1 + abs(Zmed))
@@ -203,3 +205,13 @@
         return -1
     else:
         return 0
+
+def main():
+    import sys
+    fname = sys.argv[1]
+    with open(fname) as f:
+        data = [float(x) for x in f.readlines() if x.strip() != ""]
+    print medcouple_1d(data)
+
+if __name__ == "__main__":
+    main()