ModelG: Interactive 2.5D Gravity Modeling

In geophys: Geophysics, Continuum Mechanics, Gravity Modeling

Description

Interactive 2.5D Gravity Modeling

Usage

ModelG(Data, labs = c("Done"), obs=NULL, ZCOLS=RPMG::pastel.colors(24, seed=2 ))

Arguments

Data	list: xmin, xmax, zmin, zmax
labs	Labels
obs	Observed Profile Gravity (x,g)
ZCOLS	Color Palette for plotting blocks

Details

Forward modeling using the code of Won and Bevis.

For input, at a minimum, one must provide the boundary of the model region, xmin, xmax, ymin, ymax (or depth). more comprehensive input, i.e. with a starting model, will provide a model in addition: xmin, xmax, zmin, zmax, mod, cens, n. The mod is a list of polygons, each with xy coordinates and a delta-rho value.

Value

Model

Note

You type in the polygons interactively and execute the gravity forward modeling.

Author(s)

Jonathan M. Lees<[e-mail]>

References

Won and Bevis (1987) Computing the gravitational and magnetic anomalies due to a polygon: Algorithms and Fortran subroutines <doi:https://doi.org/10.1190/1.1442298>

See Also

DGzx

Examples

## Not run: 
blum = list(xmin=0, xmax=26000, zmin=0, zmax=-1000)
gill = ModelG(blum)



## End(Not run)

geophys documentation built on May 1, 2019, 9:26 p.m.