diff --git a/.github/workflows/tests.yml b/.github/workflows/tests.yml
index 2fd881c..991a311 100644
--- a/.github/workflows/tests.yml
+++ b/.github/workflows/tests.yml
@@ -15,7 +15,7 @@ jobs:
       fail-fast: false
       matrix:
         os: ["ubuntu-latest"]
-        python-version: ["3.8", "3.10"]
+        python-version: ["3.10", "3.11"]
         test-file: ["tests/test_dcip.py", "tests/test_em.py", "tests/test_gpr.py tests/test_seismic.py", "tests/test_inversion.py", "tests/test_gravity.py tests/test_mag.py"]
 
     steps:
diff --git a/geoscilabs/dcip/DC_Pseudosections.py b/geoscilabs/dcip/DC_Pseudosections.py
index aee3982..a588901 100644
--- a/geoscilabs/dcip/DC_Pseudosections.py
+++ b/geoscilabs/dcip/DC_Pseudosections.py
@@ -53,7 +53,7 @@ def _transform(self, m):
             sig1, sig2, sig3 = np.exp(sig1), np.exp(sig2), np.exp(sig3)
         sigma = np.ones(mesh.nC) * sig1
         sigma[mesh.gridCC[:, 1] < zh] = sig2
-        blkind = utils.model_builder.getIndicesSphere(np.r_[x, zc], r, mesh.gridCC)
+        blkind = utils.model_builder.get_indices_sphere(np.r_[x, zc], r, mesh.gridCC)
         sigma[blkind] = sig3
         return sigma
 
diff --git a/geoscilabs/em/TDEMGroundedSource.py b/geoscilabs/em/TDEMGroundedSource.py
index ca87c24..948a45c 100644
--- a/geoscilabs/em/TDEMGroundedSource.py
+++ b/geoscilabs/em/TDEMGroundedSource.py
@@ -100,7 +100,7 @@ def run_simulation(fname="tdem_gs_half.h5", sigma_block=0.01, sigma_halfspace=0.
     hz = [(cs, npad, -1.5), (cs, ncz), (cs, npad, 1.5)]
     mesh = TensorMesh([hx, hy, hz], "CCC")
     sigma = np.ones(mesh.nC) * sigma_halfspace
-    blk_ind = utils.ModelBuilder.getIndicesBlock(
+    blk_ind = utils.ModelBuilder.get_indices_block(
         np.r_[-40, -40, -160], np.r_[40, 40, -80], mesh.gridCC
     )
     sigma[mesh.gridCC[:, 2] > 0.0] = 1e-8
diff --git a/geoscilabs/inversion/TomographyInversion.py b/geoscilabs/inversion/TomographyInversion.py
index 0be8865..a7b4c75 100644
--- a/geoscilabs/inversion/TomographyInversion.py
+++ b/geoscilabs/inversion/TomographyInversion.py
@@ -130,7 +130,7 @@ def add_noise(self):
     def get_block_index(self, xc=50, yc=50, dx=20, dy=20):
         p0 = np.array([xc - dx / 2.0, yc + dy / 2])
         p1 = np.array([xc + dx / 2.0, yc - dy / 2])
-        index = utils.model_builder.getIndicesBlock(p0, p1, self.mesh_prop.gridCC)
+        index = utils.model_builder.get_indices_block(p0, p1, self.mesh_prop.gridCC)
         return index
 
     def get_block_points(self, xc=50, yc=50, dx=20, dy=20):