diff --git a/bilby/core/prior.py b/bilby/core/prior.py
index feea50a5c36be6f34affb0747bde4ea4be6da94c..20e3185cce0e82058a198e50a318de12791b1de0 100644
--- a/bilby/core/prior.py
+++ b/bilby/core/prior.py
@@ -2923,7 +2923,7 @@ class MultivariateGaussianDist(object):
                 # sample the multivariate Gaussian keys
                 vals = np.random.uniform(0, 1, len(self))
 
-                samp = self.rescale(vals, mode=mode)
+                samp = np.atleast_1d(self.rescale(vals, mode=mode))
                 samps[i, :] = samp
 
                 # check sample is in bounds (otherwise perform another draw)