doAmap {Rquake}

Plot a map of station locations

Description

Plot a map of station locations

Usage

doAmap(stas, doproj = TRUE)

Arguments

stas

station list

doproj

logical, if TRUE, project (UTM) the data so plot is in units of km with the median lat-lon as the center. If FALSE, use the lat-lon coordinates.

Details

The range of the plot is expanded by 10 percent prior to plotting.

Value

list, GEOmap projection

Author(s)

Jonathan M. Lees<[e-mail]>

See Also

gMAP,expandbound,GLOB.XY

Examples

data(coso_sta_LLZ)
###  or read in from file:
##   fsta = "staLLZ.txt"
##  stas = scan(file=fsta,what=list(name="", lat=0, lon=0, z=0))
##  stas$z = stas$z/1000

stas = coso_sta_LLZ

STA = doAmap(stas, doproj = TRUE)