fix incompressible advection diffusion kernel

src/gpu/core/Kernel/Incompressible/AdvectionDiffusion/F16/F16IncompressibleAdvectionDiffusion.cu

@@ -21,7 +21,8 @@ void F16IncompressibleAdvectionDiffusion::run()
         para->getParD(level)->neighborZ,
         para->getParD(level)->distributions.f[0], 
         para->getParD(level)->distributionsAD.f[0], 
-        para->getParD(level)->numberOfNodes,
+        para->getParD(level)->numberOfNodes, 
+        para->getParD(level)->forcing,
         para->getParD(level)->isEvenTimestep);
     getLastCudaError("F16IncompressibleAdvectionDiffusion_Device execution failed");
 }

src/gpu/core/Kernel/Incompressible/AdvectionDiffusion/F16/F16IncompressibleAdvectionDiffusion_Device.cuh

@@ -5,13 +5,14 @@
 #include <curand.h>
 
 __global__ void F16IncompressibleAdvectionDiffusion_Device(real diffusivity,
-	unsigned int* bcMatD,
-	unsigned int* neighborX,
-	unsigned int* neighborY,
-	unsigned int* neighborZ,
-	real* DDStart,
-	real* DD27,
-	unsigned long long numberOfLBnodes,
-	bool EvenOrOdd);
+    unsigned int* bcMatD,
+    unsigned int* neighborX,
+    unsigned int* neighborY,
+    unsigned int* neighborZ,
+    real* DDStart,
+    real* DD27, 
+    unsigned long long numberOfLBnodes, 
+    real* forces,
+    bool EvenOrOdd);
 
 #endif
\ No newline at end of file