diff --git a/src/gpu/VirtualFluids_GPU/GPU/PressBCs27.cu b/src/gpu/VirtualFluids_GPU/GPU/PressBCs27.cu
index c32dd9117b80ab2876625343c7cd16a4803b3c6f..02cfd2bce3723162b645cef568c87ca3b1dd2720 100644
--- a/src/gpu/VirtualFluids_GPU/GPU/PressBCs27.cu
+++ b/src/gpu/VirtualFluids_GPU/GPU/PressBCs27.cu
@@ -852,28 +852,16 @@ __global__ void QPressDeviceNEQ27(
     unsigned long long numberOfLBnodes,
     bool isEvenTimestep)
 {
+   ////////////////////////////////////////////////////////////////////////////////
+   //! The pressure boundary condition is executed in the following steps
+   //!
+
    ////////////////////////////////////////////////////////////////////////////////
    //! - Get node index coordinates from threadIdx, blockIdx, blockDim and gridDim.
    //!
    const unsigned nodeIndex = getNodeIndex();
 
-    
-   // //////////////////////////////////////////////////////////////////////////
-   ////! The pressure boundary condition is executed in the following steps
-   ////!
-   //////////////////////////////////////////////////////////////////////////////////
-   ////! - Get node index coordinates from threadIdx, blockIdx, blockDim and gridDim.
-   ////!
-   //const unsigned x = threadIdx.x;    // global x-index
-   //const unsigned y = blockIdx.x;     // global y-index
-   //const unsigned z = blockIdx.y;     // global z-index
-
-   //const unsigned nx = blockDim.x;
-   //const unsigned ny = gridDim.x;
-
-   //const unsigned kThread = nx*(ny*z + y) + x;
-
-   //////////////////////////////////////////////////////////////////////////
+   ////////////////////////////////////////////////////////////////////////////////
    //! - Run for all indices in size of boundary condition (numberOfBCnodes)
    //!
    if(nodeIndex < numberOfBCnodes)
diff --git a/src/gpu/VirtualFluids_GPU/GPU/SlipBCs27.cu b/src/gpu/VirtualFluids_GPU/GPU/SlipBCs27.cu
index e189463e0fb764a68c52aa6c0cebcab71b817c22..632712746e1101771dc651861c36ab94eee04b63 100644
--- a/src/gpu/VirtualFluids_GPU/GPU/SlipBCs27.cu
+++ b/src/gpu/VirtualFluids_GPU/GPU/SlipBCs27.cu
@@ -703,6 +703,9 @@ __global__ void QSlipDeviceComp27(
     unsigned long long numberOfLBnodes, 
     bool isEvenTimestep)
 {
+   //! The slip boundary condition is executed in the following steps
+   //!
+
    ////////////////////////////////////////////////////////////////////////////////
    //! - Get node index coordinates from threadIdx, blockIdx, blockDim and gridDim.
    //!
@@ -1232,6 +1235,9 @@ __global__ void BBSlipDeviceComp27(
     unsigned long long numberOfLBnodes, 
     bool isEvenTimestep)
 {
+   //! The slip boundary condition is executed in the following steps
+   //!
+
    ////////////////////////////////////////////////////////////////////////////////
    //! - Get node index coordinates from threadIdx, blockIdx, blockDim and gridDim.
    //!
@@ -1719,6 +1725,9 @@ __global__ void QSlipDeviceComp27TurbViscosity(
     unsigned long long numberOfLBnodes, 
     bool isEvenTimestep)
 {
+   //! The slip boundary condition is executed in the following steps
+   //!
+
    ////////////////////////////////////////////////////////////////////////////////
    //! - Get node index coordinates from threadIdx, blockIdx, blockDim and gridDim.
    //!
@@ -2257,6 +2266,8 @@ __global__ void QSlipPressureDeviceComp27TurbViscosity(
     unsigned long long numberOfLBnodes, 
     bool isEvenTimestep)
 {
+   //! The slip boundary condition is executed in the following steps
+   //!
    ////////////////////////////////////////////////////////////////////////////////
    //! - Get node index coordinates from threadIdx, blockIdx, blockDim and gridDim.
    //!