getCMT {RFOC}

Read CMT

Description

Read and reformat CMT solutions downloaded from the web.

Usage

getCMT(fn, skip=1)

Arguments

fn

character file name

skip

number of lines to skip (e.g. header)

Details

Data can be extracted from web site: http://www.globalcmt.org/CMTsearch.html

The file must be cleaned prior to scanning - on download from the web site there are extra lines on top and bottom of file. Delete these. Leave one line on the top that describesthe columns. Data is separated by blanks. The files have a mixture of dates - some with 7 component dates (YYMMDD and others with 14 components YYYYMODDHHMM these are read in separately. Missing hours and minutes areset to zero.

Value

list of CMT solution data:

lon

lon of epicenter

lat

lat of epicenter

str1

strike of fault plane

dip1

dip of fault plane

rake1

rake of fault plane

str2

strike of auxilliary plane

dip2

dip of auxilliary plane

rake2

rake of auxilliary plane

sc

scale?

iexp

exponent?

name

name, includes the date

Elat

exploding latitude, set to lat initially

Elon

exploding longitude, set to lon initially

jd

julian day

yr

year

mo

month

dom

day of month

Note

Use ExplodeSymbols or explode to get new locations for expanding the plotting points.

Author(s)

Jonathan M. Lees<[e-mail]>

References

http://www.globalcmt.org/CMTsearch.html

G. Ekstrom. Rapid earthquake analysis utilizes the internet. Computers in Physics, 8:632-638, 1994.

See Also

ExplodeSymbols, spherefocgeo, ternfocgeo

Examples

## Not run: 


g = getCMT("/home/lees/aleut.cmt")

pg = prepFOCS(g)


plot(range(pg$LONS), range(pg$LATS), type = "n", xlab = "LON",
    ylab = "LAT", asp = 1)


 for (i in 1:length(pg$LATS)) {
    mc = CONVERTSDR(g$str1[i], g$dip1[i], g$rake1[i])
     MEC <- MRake(mc$M)
MEC$UP = FALSE
     Fcol <- foc.color(foc.icolor(MEC$rake1), pal = 1)
     justfocXY(MEC, x = pg$LONS[i], y = pg$LATS[i], focsiz = 0.4,
     fcol = Fcol, xpd = FALSE)
 }




## End(Not run)