Potential Theory in Gravity and Magnetic Applications

Subroutines


The textbook contains an appendix of computer subroutines written in FORTRAN that provide insight into underlying theories discussed in the text. The subroutines are used in some of the problem sets that follow each chapter, and they provide a reference source with which readers can develop their own computer programs. The subroutines are listed in the following table. They can be downloaded individually by selecting the appropriate subroutine name, or they can be downloaded en masse if preferred.

Name Function
contin Analytically continue a gridded potential field from one horizontal level to another
cross Calculate vector products
cylind Calculate the gravitational attraction of an infinitely extended cylinder
dipole Calculate the magnetic induction of a dipole
dircos Calculate direction cosines
expand Add tapered rows and columns to a grid
fac Calculate factorials
facmag Calculate magnetic induction of one polygonal facet of a polyhedron
fork Calculate the one-dimensional Fourier transform and its inverse
fourn Calculate an n-dimensional Fourier transform and its inverse
gbox Calculate the gravitational attraction of a right rectangular prism
gfilt Calculate the earth filter (gravity case) for a horizontal layer
glayer Calculate the gravitational attraction of a flat, horizontal layer
gpoly Calculate the gravitational attraction of an infinitely extended polygonal prism
hgrad Calculate the maximum horizontal gradient of a gridded dataset
kvalue Calculate wavenumber coordinates
line Calculate the intersection of two lines
mbox Calculate the total-field anomaly of a right rectangular prism
mfilt Calculate the earth filter (magnetic case) for a horizontal layer
mlayer Calculate the magnetic induction of a flat, horizontal layer
mtopo Calculate the total-field anomaly of a layer with topographic surface
newvec Transform the direction of magnetization or regional field vectors in gridded magnetic anomalies (e.g., reduction to the pole)
plane Calculate the intersection of a plane and a perpendicular line
pseudo Transform a gridded magnetic anomaly into a pseudogravity anomaly
ribbon Calculate magnetic induction of an infinitely extended ribbon of magnetic charge
rot Find the sense of rotation of one vector with respect to another
schmit Calculate normalized associated Legendre polynomials
sphere Calculate the gravitational attraction of a sphere
verder Calculate the nth order vertical derivative of a gridded potential field
Return to Potential Theory

Cambridge University Press (United Kingdom)

Cambridge University Press (North America)