lbmpy boundary handling works on pure flag fields, geometry setup
 classes could only set up boundary handling, not pure flag fields.