RFOC (version 2.0-02)

getCMT: 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)

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

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.

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

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], size = c(0.5,
         0.5), fcol = Fcol, xpd = FALSE)
 }