Schiller force model needs shifted equilibrium velocity

lbmpy/forcemodels.py

@@ -199,6 +199,9 @@ class Schiller:
     def macroscopic_momentum_density_shift(self, density):
         return default_momentum_density_shift(self._force)
 
+    def equilibrium_velocity_shift(self, density):
+        return default_velocity_shift(density, self._force)
+
 
 class Buick:
     r"""