New fwi notebook overview

.github/workflows/build_firedrake_main.yml

@@ -112,4 +112,4 @@ jobs:
 
       - name: Post-run cleanup
         if: always()
-        run: find . -delete
+        run: find . -delete

.gitignore

@@ -188,3 +188,4 @@ __pycache__/
 fw_outputs*.*
 output/
 forgotten_files/
+PacMeshIgnore/

README.md

@@ -160,11 +160,11 @@ dictionary["parallelism"] = {
 # Define the domain size without the PML. Here we'll assume a 0.75 x 1.50 km
 dictionary["mesh"] = {
     # depth in km - always positive
-    "Lz": 0.75,
+    "length_z": 0.75,
     # width in km - always positive
-    "Lx": 1.50,
+    "length_x": 1.50,
     # thickness in km - always positive
-    "Ly": 0.0,
+    "length_y": 0.0,
     # If we are loading and external .msh mesh file
     "mesh_file": None,
     # options: None (default), firedrake_mesh, user_mesh, or SeismicMesh

benchmarks/benchmark_2d.py

@@ -25,9 +25,9 @@
 }
 
 model["mesh"] = {
-    "Lz": 114.4,  # depth in km - always positive
-    "Lx": 85.8,  # width in km - always positive
-    "Ly": 0.0,  # thickness in km - always positive
+    "length_z": 114.4,  # depth in km - always positive
+    "length_x": 85.8,  # width in km - always positive
+    "length_y": 0.0,  # thickness in km - always positive
     "meshfile": "not_used.msh",
     "initmodel": "not_used.hdf5",
     "truemodel": "not_used.hdf5",
@@ -86,8 +86,8 @@
     M = 1.69
 
 edge_length = 0.286/M
-Real_Lz = model["mesh"]["Lz"] + model["BCs"]["lz"]
-Lx = model["mesh"]["Lx"]
+Real_Lz = model["mesh"]["length_z"] + model["BCs"]["lz"]
+Lx = model["mesh"]["length_x"]
 pad = model["BCs"]["lz"]
 
 bbox = (-Real_Lz, 0.0, -pad, Lx+pad)

benchmarks/benchmark_3d.py

@@ -25,9 +25,9 @@
 }
 
 model["mesh"] = {
-    "Lz": 8.0,  # depth in km - always positive
-    "Lx": 8.0,  # width in km - always positive
-    "Ly": 8.0,  # thickness in km - always positive
+    "length_z": 8.0,  # depth in km - always positive
+    "length_x": 8.0,  # width in km - always positive
+    "length_y": 8.0,  # thickness in km - always positive
     "meshfile": "not_used.msh",
     "initmodel": "not_used.hdf5",
     "truemodel": "not_used.hdf5",
@@ -82,9 +82,9 @@
     M = 3.1
 
 edge_length = 0.286/M
-Real_Lz = model["mesh"]["Lz"] + model["BCs"]["lz"]
-Lx = model["mesh"]["Lx"]
-Ly = model["mesh"]["Ly"]
+Real_Lz = model["mesh"]["length_z"] + model["BCs"]["lz"]
+Lx = model["mesh"]["length_x"]
+Ly = model["mesh"]["length_y"]
 pad = model["BCs"]["lz"]
 
 bbox = (-Real_Lz, 0.0, -pad, Lx+pad, -pad, Ly+pad)

benchmarks/benchmark_cuboid.py

@@ -22,9 +22,9 @@
 }
 
 model["mesh"] = {
-    "Lz": 1.0,  # depth in km - always positive
-    "Lx": 1.0,  # width in km - always positive
-    "Ly": 1.0,  # thickness in km - always positive
+    "length_z": 1.0,  # depth in km - always positive
+    "length_x": 1.0,  # width in km - always positive
+    "length_y": 1.0,  # thickness in km - always positive
     "meshfile": "not_used.msh",
     "initmodel": "not_used.hdf5",
     "truemodel": "not_used.hdf5",

demos/with_automatic_differentiation/utils.py

@@ -35,9 +35,9 @@ def model_settings():
 
     # Define the domain size without the ABL.
     model["mesh"] = {
-        "Lz": 1.0,  # depth in km - always positive
-        "Lx": 1.0,  # width in km - always positive
-        "Ly": 0.0,  # thickness in km - always positive
+        "length_z": 1.0,  # depth in km - always positive
+        "length_x": 1.0,  # width in km - always positive
+        "length_y": 0.0,  # thickness in km - always positive
         "meshfile": "not_used.msh",
         "initmodel": "not_used.hdf5",
         "truemodel": "not_used.hdf5",