update sphere to block

simpeg · Jan 24, 2024 · 8c2652a · 8c2652a
1 parent ae8d75c
commit 8c2652a
Show file tree

Hide file tree

Showing 7 changed files with 320 additions and 305 deletions.
diff --git a/data/dc_data.obs b/data/dc_data.obs
diff --git a/exercises/01-forward-dc-resisitivity-2d-exercise.ipynb b/exercises/01-forward-dc-resisitivity-2d-exercise.ipynb
@@ -342,9 +342,9 @@
     "\n",
     "* Near surface resistivity of 100 $\\Omega m$\n",
     "* Basement resistivity (z < -16 m) of 10 $\\Omega m$\n",
-    "* An infinite cylinder centered at (10, -5) with a radius of 5 m with a resistivity of 1000 $\\Omega m$\n",
+    "* A block with a resistivity of 1000 $\\Omega m$. The bottom-left corner is at (5, -10) and the top-right corner is at (15, 0)\n",
     "\n",
-    "*For simple geometric shapes, we can use utilities like* [model_builder.get_indices_sphere](https://docs.simpeg.xyz/content/api/generated/SimPEG.utils.model_builder.get_indices_sphere.html).\n",
+    "*For simple geometric shapes, we can use utilities like* [model_builder.get_indices_block](https://docs.simpeg.xyz/content/api/generated/SimPEG.utils.model_builder.get_indices_block.html).\n",
     "\n",
     "### <span style=\"color:darkorange\">Exercise (advanced):</span>\n",
     "\n",

diff --git a/exercises/02-inversion-dc-resistivity-2d-exercise.ipynb b/exercises/02-inversion-dc-resistivity-2d-exercise.ipynb
@@ -167,7 +167,7 @@
    "outputs": [],
    "source": [
     "data_url = \"https://github.com/simpeg/agrogeo24/raw/main/data/dc_data.obs\"\n",
-    "data_hash = \"9BF25B010DD97D3CE4FE14B3C330BC5C7E92D95CB3440CD611F8F8DCD8A2622C\"\n",
+    "data_hash = \"EBA4F058B5189DB974E261F4C0B4C18B7841F29E6304F3138E94FACB436A2E58\"\n",
     "\n",
     "data_filename = pooch.retrieve(data_url, known_hash=data_hash)"
    ]
@@ -1057,7 +1057,7 @@
    "source": [
     "air_resistivity = 1e8\n",
     "background_resistivity = 1e2\n",
-    "sphere_resistivity = 1e3\n",
+    "block_resistivity = 1e3\n",
     "basement_resistivity = 1e1"
    ]
   },
@@ -1073,8 +1073,8 @@
     "ind_basement = mesh.cell_centers[:, -1] < -16.0\n",
     "true_resistivity[ind_basement] = basement_resistivity\n",
     "\n",
-    "ind_sphere = model_builder.get_indices_sphere(np.r_[10.0, -5.0], 5.0, mesh.cell_centers)\n",
-    "true_resistivity[ind_sphere] = sphere_resistivity\n",
+    "ind_block = model_builder.get_indices_block(np.r_[10.0, -5.0], 5.0, mesh.cell_centers)\n",
+    "true_resistivity[ind_block] = block_resistivity\n",
     "\n",
     "true_resistivity = true_resistivity[active_cells]"
    ]

diff --git a/images/dcr_fwd.png b/images/dcr_fwd.png
diff --git a/images/dcr_inv.png b/images/dcr_inv.png
diff --git a/notebooks/01-forward-dc-resisitivity-2d.ipynb b/notebooks/01-forward-dc-resisitivity-2d.ipynb
diff --git a/notebooks/02-inversion-dc-resistivity-2d.ipynb b/notebooks/02-inversion-dc-resistivity-2d.ipynb