Klocate: Earthquake Hypocenter Location

In Rquake: Seismic Hypocenter Determination

Earthquake Hypocenter Location

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

######  could read from a list of files on disk
###  LF = list.files(path=pdir, pattern="p$", full.names=TRUE )

   data(GH, package='RSEIS')

g1 = GH$pickfile




    ##  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) )

 ######  let the code determine the initial guess
    NEW = Klocate(Ldat  )

  

Rquake documentation built on June 22, 2024, 11:54 a.m.