box geometry works. put slight tolerance on is_in functions, still not sure about this, should try and remove it and see what happens. coupling interface now expects a region, rather than a surface, and can be corrected for multiple crossings.