P_wave_gen {ezECM}R Documentation

Generation of noisy p-wave arrival times

Description

Similar utility to time_fn, however multiple seismometer locations can be provided simultaneously and normally distributed noise is added to the arrival time.

Usage

P_wave_gen(
  Si = NULL,
  S0 = NULL,
  Sig = NULL,
  neg.obs = TRUE,
  eps = sqrt(.Machine$double.eps)
)

Arguments

Si

Numeric matrix providing seismometer locations. Must contain 3 columns corresponding to (X,Y) corrdinates and depth.

S0

Numeric 3 element vector stipulating the location of an event, elements correspond to (X, Y, Z)

Sig

Numeric vector, or diagonal matrix, providing the variance in observed arrival times at each seismometer.

neg.obs

Logical indicating whether to allow negative observations of time (eg. the observed time of p-wave arrival is before the true time for the event).

eps

Numeric. If neg.obs = FALSE sets of observations are redrawn until all t_i - t_0 \leq eps.

Value

Numeric vector of observation times that correspond to the rows of Si

Examples


pwave.obs <- P_wave_gen(Si = c(100,200,3), S0 = c(400, 500, 4), Sig = 0.05)


[Package ezECM version 1.0.0 Index]