desh: Display Mesh

In geophys: Geophysics, Continuum Mechanics, Gravity Modeling

Description

Display Mesh

Usage

desh(M, add = TRUE, PTS = TRUE, colmesh = grey(0.8), colpts = grey(0.5),
...)

Arguments

M	Mesh output of meshgrid
add	logical, add to plot
PTS	logical, add points
colmesh	color for mesh
colpts	color for points
...	graphical parameters from par

Details

Used for showing strain mesh plots.

Value

graphical side effects

Author(s)

Jonathan M. Lees<[e-mail]>

See Also

meshgrid

Examples

###  warping of the ground from a Mogi source
 P=1e5
 a = 0.01
         f = .1
      nu=0.25
         E = 10e13
     mu = E/(2*(1+nu));
    
EX = seq(from=-3, by=.4, to= 3)
     mm = RPMG::meshgrid(EX, EX)
  rho = sqrt(mm$x^2 +  mm$y^2)

cen = list(x=mean(EX), y=mean(EX))

are = sqrt( (mm$x-cen$x)^2  + (mm$y-cen$y)^2 )


 o1 = mogi1(a, P, rho, mu, f)


shiftx =  o1$ur*(mm$x-cen$x)/are
shifty =  o1$ur*(mm$y-cen$y)/are 

newmesh = list(x=mm$x+shiftx , y=mm$y+shifty)

desh(mm, add=FALSE, PTS=FALSE, colmesh=grey(.8) )


desh(newmesh, add=TRUE, PTS=FALSE, colmesh=rgb(.6,.8, .6) )

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