See #108. Currently only topboundary='confined' is implemented.
See #108. Currently only topboundary='confined' is implemented.