choclo.prism.gravity_e

Contents

choclo.prism.gravity_e#

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

Easting component of the gravitational acceleration due to a prism.

Returns the easting component of the gravitational acceleration 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_efloat

Easting component of the gravitational acceleration generated by the rectangular prism on the observation point in m/s2.

Notes

Returns the easting component gx(p) of the gravitational acceleration g 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 ρ:

gx(p)=Gρ|||kx(x,y,z)|X1X2|Y1Y2|Z1Z2

where

kx(x,y,z)=[ysafe-ln(z,r)+zsafe-ln(y,r)xsafe-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-ln and safe-arctan functions are defined as follows:

safe-ln(x,r)={0x=0,r=0ln(x+r)x0ln((r2x2)/(rx))x<0,rxln(2x)x<0,r=x
safe-arctan(y,x)={arctan(yx)x0π2x=0andy>0π2x=0andy<00x=0andy=0

These were defined after [Fukushima2020] and guarantee a good accuracy on any observation point.

References