Get Seismic Gap

Get Pand S travel times and derivatives

1 | ```
GETpsTT(phase, eqz = 6, staz = 0, delx = 1, dely = 1, deltadis = 6, vel)
``` |

`phase` |
character vector, phase |

`eqz` |
event depth |

`staz` |
station elevation |

`delx` |
km, delta X |

`dely` |
km, delta Y |

`deltadis` |
km, distance |

`vel` |
velocity models (P and S) |

Creates a vector of travel times, and a matrix and derivatives used for inversion.

list:

`TT` |
travel time vector |

`Derivs` |
matrix of derivatives, dtdx, dtdy, dtdz |

Jonathan M. Lees<jonathan.lees@unc.edu>

many.time1D

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 | ```
## Not run:
library(RSEIS)
library(GEOmap)
data(GH)
data(VELMOD1D)
vel = VELMOD1D
p1 = GH$pickfile$STAS
loc = GH$pickfile$LOC
proj = GEOmap::setPROJ(type = 2, LAT0 =loc$lat, LON0 = loc$lon)
XYsta = GEOmap::GLOB.XY(p1$lat, p1$lon, proj)
XYq = GEOmap::GLOB.XY(loc$lat, loc$lon, proj)
delx = XYq$x-XYsta$x
dely = XYq$y-XYsta$y
dists = sqrt(delx^2+dely^2)
G1 = GETpsTT(p1$phase, eqz=loc$z, staz=0, delx=delx, dely=dely, deltadis=dists , vel)
## End(Not run)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.