--- a/jmedcouple.c++
+++ b/jmedcouple.c++
@@ -4,7 +4,6 @@
 #include <string>
 #include <fstream>
 #include <iomanip>
-#include <sstream>
 #include <algorithm>
 #include <utility>
 
@@ -68,8 +67,7 @@
 
 double medcouple(const std::vector<double>& X,
                  double eps1 = numeric_limits<double>::epsilon(), 
-                 double eps2 = numeric_limits<double>::min(),
-                 int trace_lev = 0)
+                 double eps2 = numeric_limits<double>::min())
 {
   size_t n = X.size(), n2 = (n - 1)/2;
   
@@ -226,13 +224,6 @@
     return 1;
   }
 
-  double trace_lev = 0;
-  if(argc > 2) {
-    stringstream ss;
-    ss << argv[2];
-    ss >> trace_lev;
-  }
-
   ifstream ifs(fname);
   if (not ifs) {
     cerr << "File " << fname << " not found." << endl;
@@ -246,7 +237,7 @@
   }
   
   cout << setprecision(16) 
-       << medcouple(data, eps1, eps2, trace_lev) << endl;
+       << medcouple(data, eps1, eps2) << endl;
 
   return 0;
 }