Bug raus von interpolation

program/source/interpol/interpol.cc

@@ -740,7 +740,7 @@ if(boxWSD.z < 0 && boxENT.z > 0.0 && boxWSD.y < 0.5 && boxENT.y > 0.5 && boxWSD.
             }
 
             #pragma omp critical
-            if(stop==false && typ_tet[ind_global] != 6) {
+            if(stop==false) { // && typ_tet[ind_global] != 6) {  // Phillip:   das darf man so nicht machen!
               ids_hex[ind_global] = id_hex;
               ids_i[ind_global] = i;
               ids_j[ind_global] = j;
@@ -1108,9 +1108,9 @@ Interpolate_on_block_grid::Interpolate_on_block_grid(int nx_, int ny_, int nz_,
     hz = (pENT.z - pWSD.z) / (nz-1);
     
     
-    /*
     
-    // test GGGG
+    /*
+    // test
     cout << "\n WSD: " ; pWSD.Print();
     cout << "\n ENT: " ; pENT.Print();
     cout << "nx: " << nx << " ny: " << ny << " nz: " << nz << endl;

program/source/interpol/interpol.h

@@ -194,6 +194,7 @@ void Interpolate_on_structured_grid::interpolate(Variable<DTyp>& u, DTyp* data,
 	         << endl;
   	  */
 
+
       if(typ==0) data[ind_global] = interpolate_in_tet(lambda[ind_global],
 							   du.WND(),du.WNT(),du.WST(),du.EST());
       if(typ==1) data[ind_global] = interpolate_in_tet(lambda[ind_global],