diff --git a/manimlib/mobject/types/dot_cloud.py b/manimlib/mobject/types/dot_cloud.py
index 28b1713b..92c20ac1 100644
--- a/manimlib/mobject/types/dot_cloud.py
+++ b/manimlib/mobject/types/dot_cloud.py
@@ -8,6 +8,8 @@ from manimlib.utils.iterables import resize_preserving_order
 
 
 DEFAULT_DOT_CLOUD_RADIUS = 0.05
+DEFAULT_GRID_HEIGHT = 6
+DEFAULT_BUFF_RATIO = 0.5
 
 
 class DotCloud(PMobject):
@@ -34,36 +36,57 @@ class DotCloud(PMobject):
         self.data["radii"] = np.zeros((1, 1))
         self.set_radii(self.radii)
 
-    def set_points_by_grid(self, n_rows, n_cols, height=None, width=None):
-        # TODO, add buff/hbuff/vbuff args...
-        new_points = np.zeros((n_rows * n_cols, 3))
-        new_points[:, 0] = np.tile(range(n_cols), n_rows)
-        new_points[:, 1] = np.repeat(range(n_rows), n_cols)
-        new_points[:, 2] = 0
-        self.set_points(new_points)
+    def to_grid(self, n_rows, n_cols, n_layers=1,
+                buff_ratio=None,
+                h_buff_ratio=1.0,
+                v_buff_ratio=1.0,
+                d_buff_ratio=1.0,
+                height=DEFAULT_GRID_HEIGHT,
+                ):
+        n_points = n_rows * n_cols * n_layers
+        points = np.repeat(range(n_points), 3, axis=0).reshape((n_points, 3))
+        points[:, 0] = points[:, 0] % n_cols
+        points[:, 1] = (points[:, 1] // n_cols) % n_rows
+        points[:, 2] = points[:, 2] // (n_rows * n_cols)
+        self.set_points(points.astype(float))
 
-        radius = self.data["radii"].max()
-        if height is None:
-            height = n_rows * 3 * radius
-        if width is None:
-            width = n_cols * 3 * radius
+        if buff_ratio is not None:
+            v_buff_ratio = buff_ratio
+            h_buff_ratio = buff_ratio
+            d_buff_ratio = buff_ratio
 
-        self.set_height(height, stretch=True)
-        self.set_width(width, stretch=True)
+        radius = self.get_radius()
+        ns = [n_cols, n_rows, n_layers]
+        brs = [h_buff_ratio, v_buff_ratio, d_buff_ratio]
+        for n, br, dim in zip(ns, brs, range(3)):
+            self.rescale_to_fit(2 * radius * (1 + br) * (n - 1), dim, stretch=True)
+        if height is not None:
+            self.set_height(height)
         self.center()
-
         return self
 
     def set_radii(self, radii):
         if not isinstance(radii, numbers.Number):
-            radii = resize_preserving_order(radii, self.get_num_points())
-        self.data["radii"][:, 0] = radii
+            radii = resize_preserving_order(radii, len(self.data["radii"]))
+        self.data["radii"][:] = radii
         return self
 
+    def get_radii(self):
+        return self.data["radii"]
+
+    def get_radius(self):
+        return self.get_radii().max()
+
     def scale(self, scale_factor, scale_radii=True, **kwargs):
         super().scale(scale_factor, **kwargs)
         if scale_radii:
-            self.data["radii"] *= scale_factor
+            self.set_radii(scale_factor * self.get_radii())
+        return self
+
+    def style_for_3d(self, gloss=0.5, shadow=0.2):
+        self.set_gloss(gloss)
+        self.set_shadow(shadow)
+        self.apply_depth_test()
         return self
 
     def get_shader_data(self):