irisws: Access the suite of IRIS Web Services

Description Usage Arguments Details Value Author(s) References See Also Examples

From [1]: The distance-azimuth service will calculate the great-circle angular distance, azimuth, and backazimuth between two geographic coordinate pairs. All results are reported in degrees, with azimuth and backazimuth measured clockwise from North. The azimuth is the angle from the station to the event, while the backazimuth is the angle from the event to the station. Latitudes are converted to geocentric latitudes using the WGS84 spheroid to correct for ellipticity.

ws.distaz(station.latlon = c(0, 0), event.latlon = c(0, 0),
  filename = NULL, verbose = TRUE)

distaz.ws(station.latlon = c(0, 0), event.latlon = c(0, 0),
  filename = NULL, verbose = TRUE)

`station.latlon`	numeric; the decimal latitude and longitude of the station
`event.latlon`	numeric; the decimal latitude and longitude of the event (the 'source')
`filename`	`NA` for an auto-generated filename based on the inputs; `NULL` for a temporary filename; or, a character string of the user's choosing. See Details for details about the auto-generated name.
`verbose`	logical; should messages be given?

The query is generated with constructor2, and executed with query.iris, which throws errors based on check.query.

distaz.ws is simply a pointer to ws.distaz

Filename options (`filename=`)

NA: An auto-generated file.
NULL: A temporary file.
<character string>: Any desired name.

Autogenerated files are of the form:

iriswsQ.distaz

Temporary files are generated within query.iris.

A list (invisibly) with the query string, and data from the result

AJ Barbour

[1] http://service.iris.edu/irisws/distaz/1/

Other WebServices: flinnengdahl, irisws-webservices, resp, timeseries, traveltime

## Not run: 

# exact same point -- should give zero
qr <- ws.distaz()
querydata(qr)

# two different locations -- note the coordinates
# must be in the order of latitude, longitude
qr <- ws.distaz(c(0.,0.),c(30.,-100.))
querydata(qr)
# should give
#   azimuth backAzimuth distance
#1 84.98686    300.2136 98.66374


## End(Not run)