Corrected gallilean correction

lbmpy/turbulence_models.py

@@ -33,6 +33,9 @@ def add_smagorinsky_model(collision_rule, smagorinsky_constant, omega_output_fie
     if sp.Symbol('rho') not in symbols and sp.Symbol('m_00') in symbols:
         collision_rule.subexpressions += [Assignment(sp.Symbol('rho'), sp.Symbol('m_00'))]
 
+    if sp.Symbol('rho') not in symbols and sp.Symbol('m_000') in symbols:
+        collision_rule.subexpressions += [Assignment(sp.Symbol('rho'), sp.Symbol('m_000'))]
+
     adapted_omega = sp.Symbol("smagorinsky_omega")
 
     collision_rule = collision_rule.new_with_substitutions({omega_s: adapted_omega})