Klocate: Earthquake Hypocenter Location

In Rquake: Seismic Hypocenter Determination

Description

Earthquake Hypocenter Location

Usage

Klocate(Ldat, sol = c(0, 0, 0, 0), vel=defaultVEL(6),
distwt = 20, errtol = c(0.01, 0.01, 0.01), maxit = 20,
Lambda = 1, guessdepth = 6, APLOT = FALSE,
stas = list(name = "", lat = NA, lon = NA, z = NA))

Arguments

Ldat	swig pick list
sol	vector, initial solution
vel	velocity list
distwt	distance weight parameter
errtol	error tolerance
maxit	Maximum number of iterations
Lambda	damping parameter
guessdepth	initial depth for guess
APLOT	logical, plot intermediate solutions
stas	station list

Details

Inversion is done with SVD.

Value

Event location in Lat-Lon-Z-T.

Note

Damped least squares.

Author(s)

Jonathan M. Lees<[e-mail]>

See Also

swig, defaultVEL

Examples

## Not run: 

LF = list.files(path=pdir, pattern="p$", full.names=TRUE )

GYPSY2 = vector(mode="list")

for(i in 1:length(LF))
  {
    g1 = getpfile( LF[i], sta=staf ) 
    ##  points(g1$H$lon, g1$H$lat, pch=8, col='red')

    w1 = which(!is.na(g1$STAS$lat))
    sec = g1$STAS$sec[w1]

    N = length(sec)
    Ldat =    list(
      name = g1$STAS$name[w1],
      sec = g1$STAS$sec[w1],
      phase = g1$STAS$phase[w1],
      lat=g1$STAS$lat[w1],
      lon = g1$STAS$lon[w1],
      z = g1$STAS$z[w1],
      err= g1$STAS$err[w1],
      yr = rep(g1$LOC$yr , times=N),
      jd = rep(g1$LOC$jd, times=N),
      mo = rep(g1$LOC$mo, times=N),
      dom = rep(g1$LOC$dom, times=N),
      hr =rep( g1$LOC$hr, times=N),
      mi = rep(g1$LOC$mi, times=N) )

 
    NEW = Klocate(Ldat, sol=c(MYLOC$y, MYLOC$x, 6, 0) )

    GYPSY2[[i]] = NEW
    
 
  }



## End(Not run)

Rquake documentation built on Dec. 16, 2020, 5:06 p.m.