Create new notebook template

assets/.virtual_documents/new_notebook_template.ipynb

@@ -0,0 +1,59 @@
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+# SimPEG functionality
+from simpeg.potential_fields import gravity
+from simpeg.utils import plot2Ddata, model_builder
+from simpeg import maps
+
+# discretize functionality
+from discretize import TensorMesh
+from discretize.utils import mkvc, active_from_xyz
+
+# Common Python functionality
+import numpy as np
+from scipy.interpolate import LinearNDInterpolator
+import matplotlib as mpl
+
+mpl.rcParams.update({"font.size": 14})
+import matplotlib.pyplot as plt
+import os
+
+save_output = False  # Optional
+
+
+
+
+
+[x_topo, y_topo] = np.meshgrid(np.linspace(-200, 200, 41), np.linspace(-200, 200, 41))
+
+rng = np.random.default_rng(seed=737)
+z_topo = (
+    -15 * np.exp(-(x_topo**2 + y_topo**2) / 80**2)
+    + 100.0
+    + rng.uniform(low=0.0, high=0.5, size=x_topo.shape)
+)
+
+
+fig = plt.figure(figsize=(6, 6))
+ax = fig.add_axes([0.1, 0.1, 0.8, 0.8], projection="3d")
+ax.set_zlim([z_topo.max() - 40, z_topo.max()])
+ax.plot_surface(x_topo, y_topo, z_topo, color="r", edgecolor="k", linewidth=0.5)
+ax.set_box_aspect(aspect=None, zoom=0.85)
+ax.set_xlabel("X (m)", labelpad=10)
+ax.set_ylabel("Y (m)", labelpad=10)
+ax.set_zlabel("Z (m)", labelpad=10)
+ax.set_title("Topography (Exaggerated z-axis)", fontsize=16, pad=-20)
+ax.view_init(elev=20.0)