choclo.prism.kernel_n#

choclo.prism.kernel_n(easting, northing, upward, radius)[source]#

Kernel for northing component of the gradient due to a rectangular prism.

Evaluates the integration kernel for the northing component of the gradient of the potential field generated by a prism [Nagy2000] on a single vertex of the prism. The coordinates that must be passed are shifted coordinates: the coordinates of the vertex from a Cartesian coordinate system whose origin is located in the observation point.

This function makes use of a safe natural logarithmic function and a safe arctangent function [Fukushima2020] that guarantee a good accuracy on every observation point.

Parameters:

easting, northing, upwardfloat

Shifted easting, northing and upward coordinates of the vertex of the prism. Must be in meters.

radiusfloat

Square root of the sum of the squares of the easting, northing and upward shifted coordinates.

Returns:

kernelfloat

Value of the numerical kernel function for the northing component of the gradient of the potential field due to a rectangular prism evaluated on a single vertex.

Notes

Computes the following numerical kernel on the passed shifted coordinates:

$$k_y(x,y,z)=-[z\mathrm{safe}\mathrm{\_}\ln (x,r)+x\mathrm{safe}\mathrm{\_}\ln (z,r)-y\mathrm{safe-arctan}(zx,yr)]$$

where

$$\begin{array}{r}\mathrm{safe}\mathrm{\_}\ln (x,r)=\{\begin{array}{ll}0& r=0\\ \ln (x+r)& x\ge 0\\ \ln ((y^2+z^2)/(r-x))& x<0,r\ne |x|\\ -\ln (-2x)& x<0,r=|x|\end{array}\end{array}$$

and

$$\begin{array}{r}\mathrm{safe-arctan}(y,x)=\{\begin{array}{ll}\text{arctan}\left(\frac{y}{x}\right)& x\ne 0\\ \frac{\pi }{2}& x=0\text{and}y>0\\ -\frac{\pi }{2}& x=0\text{and}y<0\\ 0& x=0\text{and}y=0\end{array}\end{array}$$

Important

In the first equation a minus sign has been added to the one obtained by [Nagy2000] in order to compute the numerical kernel for the northward component instead for the southward one.

References

• [Nagy2000]

• [Nagy2002]

• [Fukushima2020]