Execute LBM scenario with OpenCL

It compiles, it runs, so it's correct I guess...

lbmpy/lbstep.py

@@ -29,11 +29,16 @@ class LatticeBoltzmannStep:
             if domain_size is not None:
                 raise ValueError("When passing a data_handling, the domain_size parameter can not be specified")
 
+        target = optimization['target']
         if data_handling is None:
             if domain_size is None:
                 raise ValueError("Specify either domain_size or data_handling")
-            data_handling = create_data_handling(domain_size, default_ghost_layers=1,
-                                                 periodicity=periodicity, parallel=False)
+            __import__('pdb').set_trace()
+            data_handling = create_data_handling(domain_size,
+                                                 default_ghost_layers=1,
+                                                 periodicity=periodicity,
+                                                 default_target=target,
+                                                 parallel=False)
 
         if 'stencil' not in method_parameters:
             method_parameters['stencil'] = 'D2Q9' if data_handling.dim == 2 else 'D3Q27'
@@ -47,7 +52,6 @@ class LatticeBoltzmannStep:
             q = len(lbm_kernel.method.stencil)
         else:
             q = len(get_stencil(method_parameters['stencil']))
-        target = optimization['target']
 
         self.name = name
         self._data_handling = data_handling
@@ -58,7 +62,7 @@ class LatticeBoltzmannStep:
         self.density_data_index = density_data_index
         self._optimization = optimization
 
-        self._gpu = target == 'gpu'
+        self._gpu = target == 'gpu' or target == 'opencl'
         layout = optimization['field_layout']
 
         alignment = False