# HG changeset patch
# User Jordi Gutiérrez Hermoso <jordigh@octave.org>
# Date 1421186449 18000
# Node ID cbe17f888c79e5122e1b6da5c6f53b7b5df87a11
# Parent  74d0d08dbc958baede29813c0a7d51e469f901f2
Make it executable

diff --git a/medcouple.py b/pymedcouple
old mode 100644
new mode 100755
rename from medcouple.py
rename 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()