choclo.prism.gravity_ee

Contents

choclo.prism.gravity_ee#

choclo.prism.gravity_ee(easting, northing, upward, prism_west, prism_east, prism_south, prism_north, prism_bottom, prism_top, density)[source]#

Easting-easting component of the gravitational tensor due to a prism.

Returns the easting-easting component of the gravitational tensor produced by a single rectangular prism on a single computation point.

Parameters:
easting, northing, upwardfloat

Easting, northing and upward coordinates of the observation point. Must be in meters.

prism_west, prism_east, prism_south, prism_north, prism_bottom, prism_topfloat

The boundaries of the prism. Must be in meters.

densityfloat

Density of the rectangular prism in kilograms per cubic meter.

Returns:
g_eefloat

Easting-easting component of the gravitational tensor generated by the rectangular prism on the observation point in m/s2. Return numpy.nan if the observation point falls in a singular point: prism vertices or prism edges perpendicular to the easting direction.

Notes

Returns the easting-easting component gxx(p) of the gravitational tensor T on the observation point p=(xp,yp,zp) generated by a single rectangular prism defined by its boundaries x1,x2,y1,y2,z1,z2 and with a density ρ:

gxx(p)=Gρ|||kxx(x,y,z)|X1X2|Y1Y2|Z1Z2

where

kxx(x,y,z)=safe-arctan(yz,xr),
r=x2+y2+z2,

and

X1=x1xpX2=x2xpY1=y1ypY2=y2ypZ1=z1zpZ2=z2zp

are the shifted coordinates of the prism boundaries and G is the Universal Gravitational Constant.

The safe-arctan function is defined as follows:

safe-arctan(y,x)={arctan(yx)x0π2x=0andy>0π2x=0andy<00x=0andy=0

It was defined after [Fukushima2020] and guarantee a good accuracy on any observation point.

References