# NOT RUN {
# }
# NOT RUN {
library(RSEIS)
data(GH)
data(coso_sta_LLZ)
sta = coso_sta_LLZ
g1 = GH$pickfile
proj = GEOmap::setPROJ(type=2, LAT0 =median(sta$lat) , LON0 = median(sta$lon))
grcol = grey(seq(from=0.3, to=0.95, length=50))
contPFarrivals(g1, sta, proj=proj,cont=TRUE, POINTS=TRUE,
image=TRUE , col=grcol, phase="P",
add=FALSE )
# }
Run the code above in your browser using DataLab