Computation of 3D gravity anomalies for arbitrary-shaped 3D bodies, which are approximated by rectangular prisms (blocks). It is considered that a parabolic density-contrast function can reasonably approximate the arbitrary depth-dependent density contrast of geologic bodies (i.e. sedimentary basins)
Computes the theoretical gravity anomaly of 3d rectangular/square block
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=dp) | :: | x | ||||
| real(kind=dp) | :: | t | ||||
| real(kind=dp) | :: | y | ||||
| real(kind=dp) | :: | z | ||||
| real(kind=dp) | :: | sd | ||||
| real(kind=dp) | :: | alpha |
Compute the field (f array) for all ('m') or some (ista:iend) points in the xy plane, as the sum of the contributions of 'n' rectangular prisms
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer | :: | ista | ||||
| integer | :: | iend | ||||
| integer | :: | n | ||||
| integer | :: | m | ||||
| real(kind=dp) | :: | sd | ||||
| real(kind=dp) | :: | alpha | ||||
| real(kind=dp) | :: | dx | ||||
| real(kind=dp) | :: | dy | ||||
| real(kind=dp), | dimension (n) | :: | z | |||
| real(kind=dp), | dimension (n) | :: | xprm | |||
| real(kind=dp), | dimension (n) | :: | yprm | |||
| real(kind=dp), | dimension (m) | :: | xrec | |||
| real(kind=dp), | dimension (m) | :: | yrec | |||
| real(kind=dp), | dimension (m) | :: | f | |||
| character(len=5), | optional | :: | pl_opt |