diff --git a/dirneucoupling.cc b/dirneucoupling.cc
index b9c16c015344202d474843b5e6709e08b439b486..cb89749a4a08fa3ff5d2fdb24e8344c8d9bc1e65 100644
--- a/dirneucoupling.cc
+++ b/dirneucoupling.cc
@@ -119,14 +119,10 @@ int main (int argc, char *argv[]) try
     for (int i=0; i<rodX.size(); i++) {
         rodX[i].r[0] = 0.5;
         rodX[i].r[1] = 0.5;
-        rodX[i].r[2] = 5 + (i* 5 /(rodX.size()-1));
+        rodX[i].r[2] = 5 + (i* 5.0 /(rodX.size()-1));
         rodX[i].q = Quaternion<double>::identity();
     }
 
-//     rodX[rodX.size()-1].r[0] = 0.5;
-//     rodX[rodX.size()-1].r[1] = 0.5;
-//     rodX[rodX.size()-1].r[2] = 11;
-
     // /////////////////////////////////////////
     //   Read Dirichlet values
     // /////////////////////////////////////////