Merge pull request #529 from paulromano/bounding-box-bugfix

Fix bug in Complement.bounding_box

Author: Will Boyd

openmc/region.py

@@ -355,6 +355,8 @@ def bounding_box(self):
             temp_region = Intersection(*[~n for n in self.node.nodes])
         elif isinstance(self.node, Intersection):
             temp_region = Union(*[~n for n in self.node.nodes])
+        elif isinstance(self.node, Complement):
+            temp_region = self.node.node
         else:
-            temp_region = ~n
+            temp_region = ~self.node
         return temp_region.bounding_box

src/state_point.F90

@@ -870,7 +870,6 @@ subroutine write_source_bank(group_id)
     integer(HSIZE_T) :: dims(1)
     type(c_ptr) :: f_ptr
 #ifdef PHDF5
-    integer :: data_xfer_mode
     integer(HID_T) :: plist    ! property list
 #else
     integer :: i
@@ -989,7 +988,6 @@ subroutine read_source_bank(group_id)
     integer(HSIZE_T) :: offset(1) ! offset of data
     type(c_ptr) :: f_ptr
 #ifdef PHDF5
-    integer :: data_xfer_mode
     integer(HID_T) :: plist    ! property list
 #endif