Merge branch 'main' into eqs-default-depth

doc/user_guide/equivalent_sources/block-averaged-eqs.rst

@@ -106,7 +106,7 @@ And project the geographic coordinates to plain Cartesian ones:
             frame=[f"WSne+t{title}", "xa10000", "ya10000"],
             x=easting,
             y=northing,
-            color=data.total_field_anomaly_nt,
+            fill=data.total_field_anomaly_nt,
             style="c0.1c",
             cmap=True,
         )
@@ -193,7 +193,7 @@ we are efectivelly upward continuing the data.
             frame=[f"WSne+t{title}", "xa10000", "ya10000"],
             x=easting,
             y=northing,
-            color=data.total_field_anomaly_nt,
+            fill=data.total_field_anomaly_nt,
             style="c0.1c",
             cmap=True,
         )

doc/user_guide/equivalent_sources/eq-sources-spherical.rst

@@ -176,7 +176,7 @@ Lets plot it:
         frame=[f"WSne+t{title}", "xa5", "ya4"],
         x=longitude,
         y=latitude,
-        color=gravity_disturbance,
+        fill=gravity_disturbance,
         style="c0.1c",
         cmap=True,
     )

doc/user_guide/equivalent_sources/gradient-boosted-eqs.rst

@@ -183,7 +183,7 @@ And plot it:
             frame=[f"WSne+t{title}", "xa500000+a15", "ya400000"],
             x=easting,
             y=northing,
-            color=disturbance,
+            fill=disturbance,
             style="c0.1c",
             cmap=True,
         )

doc/user_guide/equivalent_sources/index.rst

@@ -149,7 +149,7 @@ And plot it:
       fig.plot(
          x=easting,
          y=northing,
-         color=disturbance,
+         fill=disturbance,
          cmap=True,
          style="c3p",
          projection=fig_proj,
@@ -165,7 +165,7 @@ And plot it:
       fig.plot(
          x=easting,
          y=northing,
-         color=data.gravity_disturbance_mgal,
+         fill=data.gravity_disturbance_mgal,
          cmap=True,
          style="c3p",
          frame=['ag', f"+t{title}"],

doc/user_guide/forward_modelling/prism.rst

@@ -343,10 +343,10 @@ bottom boundaries.
 It returns a :class:`xarray.Dataset` with the coordinates of the centers of the
 prisms and their corresponding physical properties.
 
-The :class:`harmonica.DatasetAccessorPrismsLayer` Dataset accessor can be used
+The :class:`harmonica.DatasetAccessorPrismLayer` Dataset accessor can be used
 to obtain some properties of the layer like its shape and size or the
 boundaries of any prism in the layer.
-Moreover, we can use the :meth:`harmonica.DatasetAcessorPrismsLayer.gravity`
+Moreover, we can use the :meth:`harmonica.DatasetAcessorPrismLayer.gravity`
 method to compute the gravitational field of the prism layer on any set of
 computation points.
 

doc/user_guide/topographic_correction.rst

@@ -76,7 +76,7 @@ And plot it:
    fig.plot(
       x=data.longitude,
       y=data.latitude,
-      color=data.gravity_disturbance_mgal,
+      fill=data.gravity_disturbance_mgal,
       cmap=True,
       style="c3p",
       projection="M15c",
@@ -129,7 +129,7 @@ We can now compute the Bouguer disturbance and plot it:
    fig.plot(
       x=data.longitude,
       y=data.latitude,
-      color=bouguer_disturbance,
+      fill=bouguer_disturbance,
       cmap=True,
       style="c3p",
       projection="M15c",
@@ -249,7 +249,7 @@ And plot it:
    fig.plot(
       x=data.longitude,
       y=data.latitude,
-      color=topo_free_disturbance,
+      fill=topo_free_disturbance,
       cmap=True,
       style="c3p",
       projection="M15c",

examples/equivalent_sources/block_averaged_sources.py

@@ -123,7 +123,7 @@
         frame=[f"WSne+t{title}", "xa10000", "ya10000"],
         x=easting,
         y=northing,
-        color=data.total_field_anomaly_nt,
+        fill=data.total_field_anomaly_nt,
         style="c0.1c",
         cmap=True,
     )

examples/equivalent_sources/cartesian.py

@@ -112,7 +112,7 @@
         frame=[f"WSne+t{title}", "xa10000", "ya10000"],
         x=easting,
         y=northing,
-        color=data.total_field_anomaly_nt,
+        fill=data.total_field_anomaly_nt,
         style="c0.1c",
         cmap=True,
     )

examples/equivalent_sources/gradient_boosted.py

@@ -131,7 +131,7 @@
         frame=[f"WSne+t{title}", "xa200000+a15", "ya100000"],
         x=easting,
         y=northing,
-        color=data.gravity_disturbance,
+        fill=data.gravity_disturbance,
         style="c0.1c",
         cmap=True,
     )

examples/equivalent_sources/spherical.py

@@ -96,7 +96,7 @@
     frame=["WSne", "xa5", "ya4"],
     x=longitude,
     y=latitude,
-    color=gravity_disturbance,
+    fill=gravity_disturbance,
     style="c0.1c",
     cmap=True,
 )

examples/forward/point_gravity.py

@@ -67,7 +67,7 @@
         cmap=True,
     )
 
-fig.plot(x=easting, y=northing, style="c0.2c", color="grey")
+fig.plot(x=easting, y=northing, style="c0.2c", fill="grey")
 
 fig.colorbar(cmap=True, position="JMR", frame=["a.6f.2", "x+lmGal"])
 

examples/forward/prism_layer.py

@@ -15,7 +15,7 @@
 horizontal directions with variable top and bottom boundaries. It returns
 a :class:`xarray.Dataset` with the coordinates of the centers of the prisms and
 their corresponding physical properties.
-The :class:`harmonica.DatasetAccessorPrismsLayer` Dataset accessor can be used
+The :class:`harmonica.DatasetAccessorPrismLayer` Dataset accessor can be used
 to obtain some properties of the layer like its shape and size or the
 boundaries of any prism in the layer. The methods of this Dataset accessor can
 be used together with the :func:`harmonica.prism_gravity` to compute the

harmonica/tests/test_prism_layer.py

@@ -170,7 +170,7 @@ def test_prism_layer_no_regular_grid(
 
 def test_prism_layer_attributes():
     """
-    Check attributes of the DatasetAccessorPrismsLayer class
+    Check attributes of the DatasetAccessorPrismLayer class
     """
     easting = np.linspace(1, 3, 5)
     northing = np.linspace(7, 10, 4)