From c07ce1bca97d76f0fbbcc77aa0cef86bb56e81a4 Mon Sep 17 00:00:00 2001 From: "github-actions[bot]" Date: Mon, 11 Mar 2024 16:15:09 +0000 Subject: [PATCH] Deploy dev from ff34dde --- dev/.buildinfo | 2 +- dev/_modules/boule/_ellipsoid.html | 12 +- dev/_modules/boule/_sphere.html | 15 +- dev/_modules/boule/_triaxialellipsoid.html | 117 ++++++++++- dev/_modules/index.html | 4 +- .../api/generated/boule.Sphere.rst.txt | 4 + .../generated/boule.TriaxialEllipsoid.rst.txt | 19 ++ dev/_sources/references.rst.txt | 1 + dev/api/generated/boule.Ellipsoid.html | 30 +-- dev/api/generated/boule.Sphere.html | 36 ++-- .../generated/boule.TriaxialEllipsoid.html | 98 +++++++++- dev/api/index.html | 4 +- dev/changes.html | 4 +- dev/citing.html | 4 +- dev/compatibility.html | 4 +- dev/ellipsoids.html | 4 +- dev/genindex.html | 14 +- dev/glossary.html | 4 +- dev/index.html | 4 +- dev/install.html | 16 +- dev/objects.inv | Bin 1291 -> 1326 bytes dev/overview.html | 53 +++-- dev/py-modindex.html | 4 +- dev/references.html | 7 +- dev/search.html | 4 +- dev/searchindex.js | 2 +- dev/user_guide/coordinates.html | 12 +- dev/user_guide/defining_ellipsoids.html | 4 +- dev/user_guide/gravity_disturbance.html | 183 +++++++++++++----- dev/user_guide/normal_gravity.html | 82 ++++++-- dev/versions.html | 4 +- 31 files changed, 576 insertions(+), 175 deletions(-) diff --git a/dev/.buildinfo b/dev/.buildinfo index 439f06a4..8be51281 100644 --- a/dev/.buildinfo +++ b/dev/.buildinfo @@ -1,4 +1,4 @@ # Sphinx build info version 1 # This file hashes the configuration used when building these files. When it is not found, a full rebuild will be done. -config: bd2cb66756f386d7eaf3b3bcc21e6ee5 +config: 6f05172153ebbe2dfe06684739675159 tags: 645f666f9bcd5a90fca523b33c5a78b7 diff --git a/dev/_modules/boule/_ellipsoid.html b/dev/_modules/boule/_ellipsoid.html index 46d624a3..9d3fe912 100644 --- a/dev/_modules/boule/_ellipsoid.html +++ b/dev/_modules/boule/_ellipsoid.html @@ -649,9 +649,7 @@

Source code for boule._ellipsoid

             / (1.5 * ((1 + 3 * ratio**2) * arctan - 3 * ratio))
         )
         result = (
-            self.geocentric_grav_const
-            * (1 + self._emm * aux)
-            / self.semimajor_axis**2
+            self.geocentric_grav_const * (1 + self._emm * aux) / self.semimajor_axis**2
         )
         return result
 
@@ -791,9 +789,7 @@ 
Source code for boule._ellipsoid
         .. [2] See https://en.wikipedia.org/wiki/Earth_radius#Prime_vertical
 
         """
-        return self.semimajor_axis / np.sqrt(
-            1 - self.first_eccentricity**2 * sinlat**2
-        )

+        return self.semimajor_axis / np.sqrt(1 - self.first_eccentricity**2 * sinlat**2)
[docs] def geodetic_to_spherical(self, longitude, latitude, height): """ @@ -1015,8 +1011,8 @@

Source code for boule._ellipsoid

     

         

   
-      © Copyright 2022, The Boule Developers.

-    Last updated on Oct 27, 2022.

+      © Copyright 2024, The Boule Developers.

+    Last updated on Mar 11, 2024.

 

         

     

diff --git a/dev/_modules/boule/_sphere.html b/dev/_modules/boule/_sphere.html
index c9e014ba..37ea019b 100644
--- a/dev/_modules/boule/_sphere.html
+++ b/dev/_modules/boule/_sphere.html
@@ -473,6 +473,8 @@ 
Source code for boule._sphere
     0
     >>> print(sphere.thirdflattening)
     0
+    >>> print(f"{sphere.volume:.10f} m³")
+    4.1887902048 m³
 
     """
 
@@ -550,6 +552,15 @@ 
Source code for boule._sphere
         """
         return 0
 
+    @property
+    def volume(self):
+        r"""
+        The volume of the sphere.
+        Definition: :math:`V = \dfrac{4}{3} \pi r^3`.
+        Units: :math:`m^3`.
+        """
+        return (4 / 3 * np.pi) * self.radius**3
+
 
[docs]    def normal_gravity(self, latitude, height, si_units=False):
         r"""
         Normal gravity of the sphere at the given latitude and height.
@@ -760,8 +771,8 @@ 
Source code for boule._sphere
     

         

   
-      © Copyright 2022, The Boule Developers.

-    Last updated on Oct 27, 2022.

+      © Copyright 2024, The Boule Developers.

+    Last updated on Mar 11, 2024.

 

         

     

diff --git a/dev/_modules/boule/_triaxialellipsoid.html b/dev/_modules/boule/_triaxialellipsoid.html
index 18dc6cf5..4e461ef2 100644
--- a/dev/_modules/boule/_triaxialellipsoid.html
+++ b/dev/_modules/boule/_triaxialellipsoid.html
@@ -559,7 +559,118 @@ 
Source code for boule._triaxialellipsoid
             * self.semimajor_axis
             * self.semimedium_axis
             * self.semiminor_axis
-        )

+        )
+
+    @property
+    def equatorial_flattening(self):
+        r"""
+        The equatorial flattening of the ellipsoid.
+        Definition: :math:`f_b = \frac{a - b}{a}`.
+        Units: adimensional.
+        """
+        return (self.semimajor_axis - self.semimedium_axis) / self.semimajor_axis
+
+    @property
+    def meridional_flattening(self):
+        r"""
+        The meridional flattening of the ellipsoid in the meridian plane
+        containing the semi-major axis.
+        Definition: :math:`f_c = \frac{a - c}{a}`.
+        Units: adimensional.
+        """
+        return (self.semimajor_axis - self.semiminor_axis) / self.semimajor_axis
+
+
[docs]    def geocentric_radius(self, longitude, latitude, longitude_semimajor_axis=0.0):
+        r"""
+        Radial distance from the center of the ellipsoid to its surface.
+
+        Assumes geocentric spherical latitude and geocentric spherical
+        longitudes. The geocentric radius is calculated following [Pec1983]_.
+
+        Parameters
+        ----------
+        longitude : float or array
+            Longitude coordinates on spherical coordinate system in degrees.
+        latitude : float or array
+            Latitude coordinates on spherical coordinate system in degrees.
+        longitude_semimajor_axis : float (optional)
+            Longitude coordinate of the meridian containing the semi-major axis
+            on spherical coordinate system in degrees. Optional, default value
+            is 0.0.
+
+        Returns
+        -------
+        geocentric_radius : float or array
+            The geocentric radius for the given spherical latitude(s) and
+            spherical longitude(s) in the same units as the axes of the
+            ellipsoid.
+
+
+        .. tip::
+
+            No elevation is taken into account.
+
+        Notes
+        -----
+
+        Given geocentric spherical latitude :math:`\phi` and geocentric
+        spherical longitude :math:`\lambda`, the geocentric surface radius
+        :math:`R` is computed as (see Eq. 1 of [Pec1983]_)
+
+        .. math::
+
+            R(\phi, \lambda) =
+            \frac{
+                a \, (1 - f_c) \, (1 - f_b)
+            }{
+                \sqrt{
+                    1
+                    - (2 f_c - f_c^2) \cos^2 \phi
+                    - (2 f_b - f_b^2) \sin^2 \phi
+                    - (1 - f_c)^2 (2 f_b - f_b^2)
+                        \cos^2 \phi \cos^2 (\lambda - \lambda_a)
+                }
+            },
+
+        where :math:`f_c` is the meridional flattening
+
+        .. math::
+
+            f_c = \frac{a - c}{a},
+
+        :math:`f_b` is the equatorial flattening
+
+        .. math::
+
+            f_b = \frac{a - b}{a},
+
+        with :math:`a`, :math:`b` and :math:`c` being the semi-major,
+        semi-medium and semi-minor axes of the ellipsoid, and :math:`\lambda_a`
+        being the geocentric spherical longitude of the meridian containing the
+        semi-major axis.
+
+        Note that [Pec1983]_ use geocentric spherical co-latitude, while here
+        we used geocentric spherical latitude.
+        """
+        latitude_rad = np.radians(latitude)
+        longitude_rad = np.radians(longitude)
+        longitude_semimajor_axis_rad = np.radians(longitude_semimajor_axis)
+
+        coslat, sinlat = np.cos(latitude_rad), np.sin(latitude_rad)
+
+        fc = self.meridional_flattening
+        fb = self.equatorial_flattening
+
+        radius = (self.semimajor_axis * (1.0 - fc) * (1.0 - fb)) / np.sqrt(
+            1.0
+            - (2.0 * fc - fc**2) * coslat**2
+            - (2.0 * fb - fb**2) * sinlat**2
+            - (1.0 - fc) ** 2
+            * (2.0 * fb - fb**2)
+            * coslat**2
+            * np.cos(longitude_rad - longitude_semimajor_axis_rad) ** 2
+        )
+        return radius

 

 
               

@@ -576,8 +687,8 @@ 
Source code for boule._triaxialellipsoid
     

         

   
-      © Copyright 2022, The Boule Developers.

-    Last updated on Oct 27, 2022.

+      © Copyright 2024, The Boule Developers.

+    Last updated on Mar 11, 2024.

 

         

     

diff --git a/dev/_modules/index.html b/dev/_modules/index.html
index c05521eb..36776296 100644
--- a/dev/_modules/index.html
+++ b/dev/_modules/index.html
@@ -393,8 +393,8 @@ 
All modules for which code is available

     

         

   
-      © Copyright 2022, The Boule Developers.

-    Last updated on Oct 27, 2022.

+      © Copyright 2024, The Boule Developers.

+    Last updated on Mar 11, 2024.

 

         

     

diff --git a/dev/_sources/api/generated/boule.Sphere.rst.txt b/dev/_sources/api/generated/boule.Sphere.rst.txt
index b23b819c..ff5d936a 100644
--- a/dev/_sources/api/generated/boule.Sphere.rst.txt
+++ b/dev/_sources/api/generated/boule.Sphere.rst.txt
@@ -36,6 +36,10 @@ Attributes
 
 
 
+.. autoattribute:: Sphere.volume
+
+
+
 
 
 Methods
diff --git a/dev/_sources/api/generated/boule.TriaxialEllipsoid.rst.txt b/dev/_sources/api/generated/boule.TriaxialEllipsoid.rst.txt
index 924c0158..b637a429 100644
--- a/dev/_sources/api/generated/boule.TriaxialEllipsoid.rst.txt
+++ b/dev/_sources/api/generated/boule.TriaxialEllipsoid.rst.txt
@@ -12,10 +12,18 @@ Attributes
 
 
 
+.. autoattribute:: TriaxialEllipsoid.equatorial_flattening
+
+
+
 .. autoattribute:: TriaxialEllipsoid.mean_radius
 
 
 
+.. autoattribute:: TriaxialEllipsoid.meridional_flattening
+
+
+
 .. autoattribute:: TriaxialEllipsoid.volume
 
 
@@ -31,6 +39,10 @@ Methods
     
     
     
+    
+        TriaxialEllipsoid.geocentric_radius
+    
+    
 
 .. rubric:: Methods documentation
 
@@ -38,6 +50,13 @@ Methods
 
 
 
+.. automethod:: TriaxialEllipsoid.geocentric_radius
+
+----
+
+
+
+
 .. raw:: html
 
      

diff --git a/dev/_sources/references.rst.txt b/dev/_sources/references.rst.txt
index 8d9dc7ed..840f01c7 100644
--- a/dev/_sources/references.rst.txt
+++ b/dev/_sources/references.rst.txt
@@ -11,3 +11,4 @@ References
 .. [Russell2012] Russell, C. T., Raymond, C. A., Coradini, A., McSween, H. Y., Zuber, M. T., Nathues, A., et al. (2012). Dawn at Vesta: Testing the Protoplanetary Paradigm. Science. doi:`10.1126/science.1219381 `__
 .. [Vermeille2002] Vermeille, H., 2002. Direct transformation from geocentric coordinates to geodetic coordinates. Journal of Geodesy. 76. 451-454. doi:`10.1007/s00190-002-0273-6 `__
 .. [Wieczorek2015] Wieczorek, M. A. (2015). Gravity and Topography of the Terrestrial Planets. In Treatise on Geophysics (pp. 153–193). Elsevier. doi:`10.1016/b978-0-444-53802-4.00169-x `__
+.. [Pec1983] Pěč, K. & Martinec, Z. (1983). Expansion of geoid heights over a triaxial Earth's ellipsoid into a spherical harmonic series.  Studia Geophysica et Geodaetica, 27, 217-232. doi: `10.1007/BF01592791 `__
diff --git a/dev/api/generated/boule.Ellipsoid.html b/dev/api/generated/boule.Ellipsoid.html
index b2049fb8..e3d6f415 100644
--- a/dev/api/generated/boule.Ellipsoid.html
+++ b/dev/api/generated/boule.Ellipsoid.html
@@ -504,23 +504,23 @@ 
boule.Ellipsoid
 
Parameters

 
    
-
  • name (str) – A short name for the ellipsoid, for example "WGS84".
    • 
-
  • semimajor_axis (float) – The semimajor axis of the ellipsoid. The equatorial (large) radius.
+
  • name (str) – A short name for the ellipsoid, for example "WGS84".
    • 
+
  • semimajor_axis (float) – The semimajor axis of the ellipsoid. The equatorial (large) radius.
 Definition: \(a\).
 Units: \(m\).
    • 
-
  • flattening (float) – The (first) flattening of the ellipsoid.
+
  • flattening (float) – The (first) flattening of the ellipsoid.
 Definition: \(f = (a - b)/a\).
 Units: adimensional.
    • 
-
  • geocentric_grav_const (float) – The geocentric gravitational constant. The product of the mass of the
+
  • geocentric_grav_const (float) – The geocentric gravitational constant. The product of the mass of the
 ellipsoid \(M\) and the gravitational constant \(G\).
 Definition: \(GM\). Units:
 \(m^3.s^{-2}\).
    • 
-
  • angular_velocity (float) – The angular velocity of the rotating ellipsoid.
+
  • angular_velocity (float) – The angular velocity of the rotating ellipsoid.
 Definition: \(\omega\).
 Units: \(\\rad.s^{-1}\).
    • 
-
  • long_name (str or None) – A long name for the ellipsoid, for example "World Geodetic System
+
  • long_name (str or None) – A long name for the ellipsoid, for example "World Geodetic System
 1984" (optional).
    • 
-
  • reference (str or None) – Citation for the ellipsoid parameter values (optional).
    • 
+
  • reference (str or None) – Citation for the ellipsoid parameter values (optional).
    • 
 

 

 
@@ -696,8 +696,8 @@ 
Methods
 
Parameters

 
    
-
  • latitude (float or array) – Latitude coordinates on geodetic coordinate system in degrees.
    • 
-
  • geodetic (bool) – If True (default), will assume that latitudes are geodetic
+
  • latitude (float or array) – Latitude coordinates on geodetic coordinate system in degrees.
    • 
+
  • geodetic (bool) – If True (default), will assume that latitudes are geodetic
 latitudes. Otherwise, will assume that they are geocentric
 spherical latitudes.
    • 
 

@@ -804,11 +804,11 @@ 
Methods
 
Parameters

 
    
-
  • latitude (float or array) – The geodetic latitude where the normal gravity will be computed (in
+
  • latitude (float or array) – The geodetic latitude where the normal gravity will be computed (in
 degrees).
    • 
-
  • height (float or array) – The ellipsoidal (geometric) height of computation the point (in
+
  • height (float or array) – The ellipsoidal (geometric) height of computation the point (in
 meters).
    • 
-
  • si_units (bool) – Return the value in mGal (False, default) or m/s² (True)
    • 
+
  • si_units (bool) – Return the value in mGal (False, default) or m/s² (True)
    • 
 

 

 
Returns

@@ -830,7 +830,7 @@ 
Methods
 
Parameters

-
sinlat (float or array-like) – Sine of the geocentric geodetic latitude.

+
sinlat (float or array-like) – Sine of the geocentric geodetic latitude.

 

 
Returns

 
prime_vertical_radius (float or array-like) – Prime vertical radius given in the same units as the semi-major
@@ -914,8 +914,8 @@ 
Methods
         

   
-      © Copyright 2022, The Boule Developers.

-    Last updated on Oct 27, 2022.

+      © Copyright 2024, The Boule Developers.

+    Last updated on Mar 11, 2024.

 

         

     

diff --git a/dev/api/generated/boule.Sphere.html b/dev/api/generated/boule.Sphere.html
index 1192605f..2c898aba 100644
--- a/dev/api/generated/boule.Sphere.html
+++ b/dev/api/generated/boule.Sphere.html
@@ -505,20 +505,20 @@ 
boule.Sphere
 
Parameters

 
    
-
  • name (str) – A short name for the sphere, for example "Moon".
    • 
-
  • radius (float) – The radius of the sphere.
+
  • name (str) – A short name for the sphere, for example "Moon".
    • 
+
  • radius (float) – The radius of the sphere.
 Definition: \(R\).
 Units: \(m\).
    • 
-
  • geocentric_grav_const (float) – The geocentric gravitational constant. The product of the mass of the
+
  • geocentric_grav_const (float) – The geocentric gravitational constant. The product of the mass of the
 sphere \(M\) and the gravitational constant \(G\).
 Definition: \(GM\). Units:
 \(m^3.s^{-2}\).
    • 
-
  • angular_velocity (float) – The angular velocity of the rotating sphere.
+
  • angular_velocity (float) – The angular velocity of the rotating sphere.
 Definition: \(\omega\).
 Units: \(\\rad.s^{-1}\).
    • 
-
  • long_name (str or None) – A long name for the sphere, for example "Moon Reference System"
+
  • long_name (str or None) – A long name for the sphere, for example "Moon Reference System"
 (optional).
    • 
-
  • reference (str or None) – Citation for the sphere parameter values (optional).
    • 
+
  • reference (str or None) – Citation for the sphere parameter values (optional).
    • 
 

 

 
@@ -560,6 +560,8 @@ 
boule.Sphere0
 >>> print(sphere.thirdflattening)
 0
+>>> print(f"{sphere.volume:.10f} m³")
+4.1887902048 m³
@@ -617,6 +619,14 @@

Attributes +
+Sphere.volume#
+

The volume of the sphere. +Definition: \(V = \dfrac{4}{3} \pi r^3\). +Units: \(m^3\).

+
+

Methods#

@@ -654,9 +664,9 @@

Methods
Parameters
    -

  • height (float or array) – The height above the surface of the sphere of the computation point +

  • height (float or array) – The height above the surface of the sphere of the computation point (in meters).

    • -

  • si_units (bool) – Return the value in mGal (False, default) or m/s² (True)

    • +

  • si_units (bool) – Return the value in mGal (False, default) or m/s² (True)

Returns
@@ -702,11 +712,11 @@

Methods
Parameters
    -

  • latitude (float or array) – The spherical latitude where the normal gravity will be computed +

  • latitude (float or array) – The spherical latitude where the normal gravity will be computed (in degrees).

    • -

  • height (float or array) – The height above the surface of the sphere of the computation point +

  • height (float or array) – The height above the surface of the sphere of the computation point (in meters).

    • -

  • si_units (bool) – Return the value in mGal (False, default) or m/s² (True)

    • +

  • si_units (bool) – Return the value in mGal (False, default) or m/s² (True)

Returns
@@ -800,8 +810,8 @@

Methods

- © Copyright 2022, The Boule Developers.
- Last updated on Oct 27, 2022.
+ © Copyright 2024, The Boule Developers.
+ Last updated on Mar 11, 2024.

diff --git a/dev/api/generated/boule.TriaxialEllipsoid.html b/dev/api/generated/boule.TriaxialEllipsoid.html index 8f4b2603..41edaddd 100644 --- a/dev/api/generated/boule.TriaxialEllipsoid.html +++ b/dev/api/generated/boule.TriaxialEllipsoid.html @@ -510,26 +510,26 @@

boule.TriaxialEllipsoid
Parameters
    -

  • name (str) – A short name for the ellipsoid, for example "WGS84".

    • -

  • semimajor_axis (float) – The semimajor (largest) axis of the ellipsoid. +

  • name (str) – A short name for the ellipsoid, for example "WGS84".

    • +

  • semimajor_axis (float) – The semimajor (largest) axis of the ellipsoid. Definition: \(a\). Units: \(m\).

    • -

  • semimedium_axis (float) – The semimedium (middle) axis of the ellipsoid. +

  • semimedium_axis (float) – The semimedium (middle) axis of the ellipsoid. Definition: \(b\). Units: \(m\).

    • -

  • semiminor_axis (float) – The semiminor (smallest) axis of the ellipsoid. +

  • semiminor_axis (float) – The semiminor (smallest) axis of the ellipsoid. Definition: \(c\). Units: \(m\).

    • -

  • geocentric_grav_const (float) – The geocentric gravitational constant. The product of the mass of the +

  • geocentric_grav_const (float) – The geocentric gravitational constant. The product of the mass of the ellipsoid \(M\) and the gravitational constant \(G\). Definition: \(GM\). Units: \(m^3.s^{-2}\).

    • -

  • angular_velocity (float) – The angular velocity of the rotating ellipsoid. +

  • angular_velocity (float) – The angular velocity of the rotating ellipsoid. Definition: \(\omega\). Units: \(\\rad.s^{-1}\).

    • -

  • long_name (str or None) – A long name for the ellipsoid, for example "World Geodetic System +

  • long_name (str or None) – A long name for the ellipsoid, for example "World Geodetic System 1984" (optional).

    • -

  • reference (str or None) – Citation for the ellipsoid parameter values (optional).

    • +

  • reference (str or None) – Citation for the ellipsoid parameter values (optional).

@@ -568,6 +568,14 @@

boule.TriaxialEllipsoid

Attributes#

+
+
+TriaxialEllipsoid.equatorial_flattening#
+

The equatorial flattening of the ellipsoid. +Definition: \(f_b = \frac{a - b}{a}\). +Units: adimensional.

+
+
TriaxialEllipsoid.mean_radius#
@@ -576,6 +584,15 @@

Attributes\(m\).

+
+
+TriaxialEllipsoid.meridional_flattening#
+

The meridional flattening of the ellipsoid in the meridian plane +containing the semi-major axis. +Definition: \(f_c = \frac{a - c}{a}\). +Units: adimensional.

+
+
TriaxialEllipsoid.volume#
@@ -594,9 +611,70 @@

Methods +

TriaxialEllipsoid.geocentric_radius(...[, ...])

+

Radial distance from the center of the ellipsoid to its surface.

+

Methods documentation

+
+
+TriaxialEllipsoid.geocentric_radius(longitude, latitude, longitude_semimajor_axis=0.0)[source]#
+

Radial distance from the center of the ellipsoid to its surface.

+

Assumes geocentric spherical latitude and geocentric spherical +longitudes. The geocentric radius is calculated following [Pec1983].

+
+
Parameters
+
    +

  • longitude (float or array) – Longitude coordinates on spherical coordinate system in degrees.

    • +

  • latitude (float or array) – Latitude coordinates on spherical coordinate system in degrees.

    • +

  • longitude_semimajor_axis (float (optional)) – Longitude coordinate of the meridian containing the semi-major axis +on spherical coordinate system in degrees. Optional, default value +is 0.0.

    • +
+
+
Returns
+

geocentric_radius (float or array) – The geocentric radius for the given spherical latitude(s) and +spherical longitude(s) in the same units as the axes of the +ellipsoid.

+
+
+
+

Tip

+

No elevation is taken into account.

+
+

Notes

+

Given geocentric spherical latitude \(\phi\) and geocentric +spherical longitude \(\lambda\), the geocentric surface radius +\(R\) is computed as (see Eq. 1 of [Pec1983])

+
+\[R(\phi, \lambda) = +\frac{ + a \, (1 - f_c) \, (1 - f_b) +}{ + \sqrt{ + 1 + - (2 f_c - f_c^2) \cos^2 \phi + - (2 f_b - f_b^2) \sin^2 \phi + - (1 - f_c)^2 (2 f_b - f_b^2) + \cos^2 \phi \cos^2 (\lambda - \lambda_a) + } +},\]
+

where \(f_c\) is the meridional flattening

+
+\[f_c = \frac{a - c}{a},\]
+

\(f_b\) is the equatorial flattening

+
+\[f_b = \frac{a - b}{a},\]
+

with \(a\), \(b\) and \(c\) being the semi-major, +semi-medium and semi-minor axes of the ellipsoid, and \(\lambda_a\) +being the geocentric spherical longitude of the meridian containing the +semi-major axis.

+

Note that [Pec1983] use geocentric spherical co-latitude, while here +we used geocentric spherical latitude.

+
+ +

@@ -629,8 +707,8 @@

Methods

- © Copyright 2022, The Boule Developers.
- Last updated on Oct 27, 2022.
+ © Copyright 2024, The Boule Developers.
+ Last updated on Mar 11, 2024.

diff --git a/dev/api/index.html b/dev/api/index.html index 0898995a..ebf408cc 100644 --- a/dev/api/index.html +++ b/dev/api/index.html @@ -492,8 +492,8 @@

List of functions and classes (API)

- © Copyright 2022, The Boule Developers.
- Last updated on Oct 27, 2022.
+ © Copyright 2024, The Boule Developers.
+ Last updated on Mar 11, 2024.

diff --git a/dev/changes.html b/dev/changes.html index 4a622721..f8c43b99 100644 --- a/dev/changes.html +++ b/dev/changes.html @@ -769,8 +769,8 @@

Version 0.0.1

- © Copyright 2022, The Boule Developers.
- Last updated on Oct 27, 2022.
+ © Copyright 2024, The Boule Developers.
+ Last updated on Mar 11, 2024.

diff --git a/dev/citing.html b/dev/citing.html index affea009..9b858024 100644 --- a/dev/citing.html +++ b/dev/citing.html @@ -476,8 +476,8 @@

Citing Boule

- © Copyright 2022, The Boule Developers.
- Last updated on Oct 27, 2022.
+ © Copyright 2024, The Boule Developers.
+ Last updated on Mar 11, 2024.

diff --git a/dev/compatibility.html b/dev/compatibility.html index 1135c2df..346d829e 100644 --- a/dev/compatibility.html +++ b/dev/compatibility.html @@ -586,8 +586,8 @@

Boule version compatibility

- © Copyright 2022, The Boule Developers.
- Last updated on Oct 27, 2022.
+ © Copyright 2024, The Boule Developers.
+ Last updated on Mar 11, 2024.

diff --git a/dev/ellipsoids.html b/dev/ellipsoids.html index 5bbb7596..386215cf 100644 --- a/dev/ellipsoids.html +++ b/dev/ellipsoids.html @@ -685,8 +685,8 @@

Vesta

- © Copyright 2022, The Boule Developers.
- Last updated on Oct 27, 2022.
+ © Copyright 2024, The Boule Developers.
+ Last updated on Mar 11, 2024.

diff --git a/dev/genindex.html b/dev/genindex.html index 689b39c2..751ca0bb 100644 --- a/dev/genindex.html +++ b/dev/genindex.html @@ -425,6 +425,8 @@

E

 @@ -451,7 +453,11 @@

G

  • Geocentric gravitational constant
  • geocentric_radius() (boule.Ellipsoid method) + +
  • geodetic_to_spherical() (boule.Ellipsoid method)
    • @@ -486,6 +492,8 @@

    M

    • latitude
      PandasIndex
      PandasIndex(Float64Index([             -90.0, -89.83333333333333, -89.66666666666667,
+                           -89.5, -89.33333333333333, -89.16666666666667,
+                           -89.0, -88.83333333333333, -88.66666666666667,
+                           -88.5,
+              ...
+                            88.5,  88.66666666666666,  88.83333333333331,
+                            89.0,  89.16666666666666,  89.33333333333331,
+                            89.5,  89.66666666666666,  89.83333333333331,
+                            90.0],
+             dtype='float64', name='latitude', length=1081))
    • longitude
      PandasIndex
      PandasIndex(Float64Index([             -180.0, -179.83333333333334, -179.66666666666666,
+                           -179.5, -179.33333333333334, -179.16666666666666,
+                           -179.0, -178.83333333333334, -178.66666666666666,
+                           -178.5,
+              ...
+                            178.5,  178.66666666666663,  178.83333333333331,
+                            179.0,  179.16666666666663,  179.33333333333331,
+                            179.5,  179.66666666666663,  179.83333333333331,
+                            180.0],
+             dtype='float64', name='longitude', length=2161))
    • @@ -1632,6 +1690,11 @@

    Gravity disturbances

    grid-column: 4; } +.xr-index-preview { + grid-column: 2 / 5; + color: var(--xr-font-color2); +} + .xr-var-name, .xr-var-dims, .xr-var-dtype, @@ -1653,14 +1716,16 @@

    Gravity disturbances

    } .xr-var-attrs, -.xr-var-data { +.xr-var-data, +.xr-index-data { display: none; background-color: var(--xr-background-color) !important; padding-bottom: 5px !important; } .xr-var-attrs-in:checked ~ .xr-var-attrs, -.xr-var-data-in:checked ~ .xr-var-data { +.xr-var-data-in:checked ~ .xr-var-data, +.xr-index-data-in:checked ~ .xr-index-data { display: block; } @@ -1670,13 +1735,16 @@

    Gravity disturbances

    .xr-var-name span, .xr-var-data, +.xr-index-name div, +.xr-index-data, .xr-attrs { padding-left: 25px !important; } .xr-attrs, .xr-var-attrs, -.xr-var-data { +.xr-var-data, +.xr-index-data { grid-column: 1 / -1; } @@ -1714,7 +1782,8 @@

    Gravity disturbances

    } .xr-icon-database, -.xr-icon-file-text2 { +.xr-icon-file-text2, +.xr-no-icon { display: inline-block; vertical-align: middle; width: 1em; @@ -1728,11 +1797,11 @@

    Gravity disturbances

    -35.83509235, -35.83509235], [-34.04097894, -34.04097894, -34.04097894, ..., -34.04097894, -34.04097894, -34.04097894], - [-33.34614029, -33.34614029, -33.34614029, ..., -33.40864029, - -33.40864029, -33.34614029], + [-33.34614042, -33.34614042, -33.34614042, ..., -33.40864042, + -33.40864042, -33.34614042], ..., - [ 11.65385965, 11.59135965, 11.46635965, ..., 11.77885965, - 11.65385965, 11.65385965], + [ 11.65385958, 11.59135958, 11.46635958, ..., 11.77885958, + 11.65385958, 11.65385958], [ 18.14652106, 18.14652106, 18.14652106, ..., 18.14652106, 18.14652106, 18.14652106], [ 15.16490765, 15.16490765, 15.16490765, ..., 15.16490765, @@ -1740,28 +1809,46 @@

    Gravity disturbances

    Coordinates: * longitude (longitude) float64 -180.0 -179.8 -179.7 ... 179.7 179.8 180.0 * latitude (latitude) float64 -90.0 -89.83 -89.67 -89.5 ... 89.67 89.83 90.0 - height (latitude, longitude) float32 1e+04 1e+04 1e+04 ... 1e+04 1e+04 + dtype=float32)
    • longitude
      PandasIndex
      PandasIndex(Float64Index([             -180.0, -179.83333333333334, -179.66666666666666,
+                           -179.5, -179.33333333333334, -179.16666666666666,
+                           -179.0, -178.83333333333334, -178.66666666666666,
+                           -178.5,
+              ...
+                            178.5,  178.66666666666663,  178.83333333333331,
+                            179.0,  179.16666666666663,  179.33333333333331,
+                            179.5,  179.66666666666663,  179.83333333333331,
+                            180.0],
+             dtype='float64', name='longitude', length=2161))
    • latitude
      PandasIndex
      PandasIndex(Float64Index([             -90.0, -89.83333333333333, -89.66666666666667,
+                           -89.5, -89.33333333333333, -89.16666666666667,
+                           -89.0, -88.83333333333333, -88.66666666666667,
+                           -88.5,
+              ...
+                            88.5,  88.66666666666666,  88.83333333333331,
+                            89.0,  89.16666666666666,  89.33333333333331,
+                            89.5,  89.66666666666666,  89.83333333333331,
+                            90.0],
+             dtype='float64', name='latitude', length=1081))

    • Finally, we can display the disturbance in a nice global map:

    @@ -1825,8 +1912,8 @@

    Gravity disturbances

    - © Copyright 2022, The Boule Developers.
    - Last updated on Oct 27, 2022.
    + © Copyright 2024, The Boule Developers.
    + Last updated on Mar 11, 2024.

    diff --git a/dev/user_guide/normal_gravity.html b/dev/user_guide/normal_gravity.html index 2b62b4da..fc478944 100644 --- a/dev/user_guide/normal_gravity.html +++ b/dev/user_guide/normal_gravity.html @@ -540,14 +540,14 @@

    Oblate ellipsoids -
    [980619.77693773 980585.4934137  980551.21168079 980516.93173896
- 980482.65358805 980448.3772279  980414.10265843 980379.82987947
- 980345.55889092 980311.28969268]
+
    [980619.77693773 980585.4934137  980551.21168081 980516.93173899
+ 980482.65358804 980448.3772279  980414.10265843 980379.82987947
+ 980345.55889093 980311.28969268]
    -

    The arrays can be multi-dimensional so we can use verde to generate a +

    The arrays can be multi-dimensional so we can use verde to generate a grid of normal gravity:

    @@ -579,7 +579,7 @@

    Oblate ellipsoidsxarray.Dataset:

    +

    Which can be put in a xarray.Dataset:

    grid = vd.make_xarray_grid(
@@ -851,6 +851,11 @@ 
    Oblate ellipsoidspygmt:
    
+
    And plotted with pygmt:
    
 
    
 
    
 
    import pygmt
@@ -1367,6 +1388,11 @@ 
    Gravity versus gravitation
         
    
   
-      © Copyright 2022, The Boule Developers.
    
-    Last updated on Oct 27, 2022.
    
+      © Copyright 2024, The Boule Developers.
    
+    Last updated on Mar 11, 2024.
    
 
    
         
    
     
    
diff --git a/dev/versions.html b/dev/versions.html
index 5af45854..01cb3774 100644
--- a/dev/versions.html
+++ b/dev/versions.html
@@ -479,8 +479,8 @@ 
    Documentation for other versions
         
    
   
-      © Copyright 2022, The Boule Developers.
    
-    Last updated on Oct 27, 2022.
    
+      © Copyright 2024, The Boule Developers.
    
+    Last updated on Mar 11, 2024.