diff --git a/flopy/mf6/utils/output_util.py b/flopy/mf6/utils/output_util.py
index e48526bf33..6424c17432 100644
--- a/flopy/mf6/utils/output_util.py
+++ b/flopy/mf6/utils/output_util.py
@@ -47,12 +47,14 @@ def __init__(self, obj):
         # capture the list file for Models and for OC packages
         if isinstance(obj, (ModelInterface, ModflowGwfoc, ModflowGwtoc)):
             if isinstance(obj, ModelInterface):
-                ml = obj
+                self._model = obj
             else:
-                ml = obj.parent
-            self._mtype = ml.model_type
-            nam_file = ml.model_nam_file[:-4]
-            self._lst = ml.name_file.blocks["options"].datasets["list"].array
+                self._model = obj.parent
+            self._mtype = self._model.model_type
+            nam_file = self._model.model_nam_file[:-4]
+            self._lst = (
+                self._model.name_file.blocks["options"].datasets["list"].array
+            )
             if self._lst is None:
                 self._lst = f"{nam_file}.lst"
             setattr(self, "list", self.__list)
@@ -60,6 +62,12 @@ def __init__(self, obj):
             if isinstance(obj, ModelInterface):
                 return
 
+        else:
+            if obj.model_or_sim.type == "Model":
+                self._model = obj.model_or_sim
+            else:
+                self._model = None
+
         obspkg = False
         if isinstance(obj, ModflowUtlobs):
             # this is a package
@@ -283,7 +291,11 @@ def __budget(self, precision="double"):
         if self._budget is not None:
             try:
                 budget_file = os.path.join(self._sim_ws, self._budget[0])
-                return CellBudgetFile(budget_file, precision=precision)
+                return CellBudgetFile(
+                    budget_file,
+                    precision=precision,
+                    modelgrid=self._model.modelgrid,
+                )
             except OSError:
                 return None
 
diff --git a/flopy/utils/binaryfile.py b/flopy/utils/binaryfile.py
index 471dd8b297..ae99bd2124 100755
--- a/flopy/utils/binaryfile.py
+++ b/flopy/utils/binaryfile.py
@@ -1666,7 +1666,11 @@ def get_record(self, idx, full3D=False):
                     s += f"a numpy recarray of size ({nlist}, 2)"
                 print(s)
             if full3D:
-                return self.__create3D(data)
+                data = self.__create3D(data)
+                if self.modelgrid is not None:
+                    return np.reshape(data, self.shape)
+                else:
+                    return data
             else:
                 return data.view(np.recarray)
         else: