diff --git a/fact/VERSION b/fact/VERSION
index 04a373e..2a0970c 100644
--- a/fact/VERSION
+++ b/fact/VERSION
@@ -1 +1 @@
-0.16.0
+0.16.1
diff --git a/fact/coordinates/representation.py b/fact/coordinates/representation.py
index 986d04b..4be9c0c 100644
--- a/fact/coordinates/representation.py
+++ b/fact/coordinates/representation.py
@@ -1,7 +1,7 @@
 from astropy.coordinates import BaseRepresentation, CartesianRepresentation
 import astropy.units as u
 from collections import OrderedDict
-from astropy.utils.compat.numpy import broadcast_arrays
+import numpy as np
 
 
 class PlanarRepresentation(BaseRepresentation):
@@ -42,7 +42,7 @@ def __init__(self, x, y, copy=True, **kwargs):
             raise u.UnitsError("x and y should have matching physical types")
 
         try:
-            x, y = broadcast_arrays(x, y, subok=True)
+            x, y = np.broadcast_arrays(x, y, subok=True)
         except ValueError:
             raise ValueError("Input parameters x and y cannot be broadcast")