## Not run:library(RSEIS)library(GEOmap)################### ###### set up stationssta=list()sta$'nam'=c("CAL","KAM","DOM","LAV","SMI","CAS")sta$'lat'=c(14.7421759974747,14.7471948493068,14.7422049415205,14.7204249827467,14.7543726234568,14.710961318972)sta$'lon'=c(-91.5659793619529,-91.5698443123368,-91.5775586192333,-91.5716896307798,-91.5518522222222,-91.5702146825397)sta$'el'=c(2.37596727272727,2.29854436407474,2.31819590643275,1.64286335403727,3.65216666666667,1.44584353741497)sta$'das'=c("CAL","KAM","DOM","LAV","SMI","CAS")sta$'sensor1'=c("60T","60T","60T","40T","INF","3T")sta$'comp1'=c("VNE","VNE","VNE","VNE","VNE","VNE")sta$'sensor2'=c("INF","INF","INF","INF","INF","INF")sta$'comp2'=c("IJK","IJK","IJK","IJK","IJK","IJK")sta$'dasSN'=c("9FF2","9FFE","9FFB","9024","A881","9026")sta$'sensorSN'=c("Unknown","Unknown","Unknown","T41034","Unknown","T3A28")sta$'start'=c("2008:366:16:02:59:615","2008:366:20:50:18:615",###### "2008:366:00:58:23:849", "2008:365:23:01:21:315","2008:366:23:57:10:244","2008:365:20:47:51:529")sta$'end'=c("2009:004:18:02:58:615","2009:004:17:50:17:615",###### "2009:004:16:58:22:849", "2009:006:15:01:20:315","2009:004:16:57:09:244","2009:005:22:47:50:529")sta$'name'=c("CAL","KAM","DOM","LAV","SMI","CAS")############## get earthquake epicenterseq1=list()eq1$'yr'=c(2008,2009,2009,2009,2008,2009,2009,2009,2009,2009,2009,2009,2009,2009,2009)eq1$'mo'=c(12,1,1,1,12,1,1,1,1,1,1,1,1,1,1)eq1$'dom'=c(30,1,3,4,30,1,2,3,3,3,3,3,4,4,6)eq1$'lat'=c(14.06,14.73,13.93,15.23,-4.3,-34.84,0.62,-0.41,-0.59,36.42,-0.32,-0.69,-0.4,36.44,-0.66)eq1$'lon'=c(-92.21,-91.39,-91.74,-92.06,101.22,-107.65,-26.66,132.88,133.36,70.74,132.88,133.3,132.76,70.88,133.43)eq1$'mag'=c(4.3,4.7,4,4.7,5.9,5.8,5.6,7.6,5.6,5.8,5.6,7.4,5.9,5.7,6)eq1$'depth'=c(9,169,61,177,20,10,10,17,35,204,29,23,35,186,16)eq1$'hr'=c(23,11,9,19,19,6,19,19,19,20,21,22,7,23,22)eq1$'mi'=c(12,44,16,2,49,27,42,43,53,23,49,33,14,12,48)eq1$'sec'=c(57,51.68,0.8,23,52.61,51.22,27.19,50.65,18.9,20.18,30.88,40.29,0.55,59.29,27.25)eq1$'z'=c(9,169,61,177,20,10,10,17,35,204,29,23,35,186,16)eq1$'jd'=c(365,1,3,4,365,1,2,3,3,3,3,3,4,4,6)########################## use the projection that is derived from the ########################## station file - these are based on the median station locations stinfo = list(mlat=median(sta$lat), mlon=median(sta$lon))proj = setPROJ(6, LAT0=stinfo$mlat, LON0=stinfo$mlon )###### get distances - this is so we can separate regional from teleseismic eventseqdists = distaz(stinfo$mlat , stinfo$mlon, eq1$lat, eq1$lon)mylist = list()for(j in1:length(eq1$sec)){mylist[[j]]= list(yr=eq1$yr[j], jd=eq1$jd[j], mo=eq1$mo[j], dom=eq1$dom[j], hr=eq1$hr[j],mi=eq1$mi[j], sec=eq1$sec[j], lat=eq1$lat[j], lon=eq1$lon[j], z=eq1$z[j], mag=eq1$mag[j])}library(geomapdata) data(worldmap) mapTeleSeis(sta, mylist, worldmap=worldmap)## End(Not run)