rotation {ssEDA}R Documentation

Rotates Events to View Plate Boundary

Description

Plots a cross section of depth for the selected events. Depth is on the vertical axis, and the view is theta degrees west of north.

Usage

rotation(events, theta=0, km=TRUE, criteria=TRUE)

Arguments

events defines the events to be plotted. It is an object of class "subset", generally created by subsetcircle, subsetpolygon, subsetrect or subsetsphere.
theta number of degrees of the direction of view from north, positive to the west. For example, theta=0 is viewing to the north (default) and for theta=-45 one would be viewing towards the NE.
km if TRUE units of kilometres are used on the horizontal scale, if FALSE units of degrees are used.
criteria boolean variable. If TRUE (default) the subsetting parameters are printed below the graph.

See Also

depth.hist, epicentres, freq.cusum, freq.magnitude, magnitude.cusum, magnitude.time, threeD, timeplot, multigraph

Examples

#    These examples require the NZ catalogue

require(ssNZ)
b <- subsetrect(NZ, minlong=170, maxlong=180, minlat=-43,
                 maxlat=-35, minmag=3, minday=julian(1,1,1970),
                 maxday=julian(1,1,1993))

rotation(b, theta=-47)
title(main="NZ Plate Boundary")

#--------------------------------------------------------------------

#   Make the Wellington Catalogue
require(ssNZ)
as.catalogue(subsetrect(NZ, minlat=-42.2, maxlat=-40.5,
             minlong=173.6, maxlong=176.0, minday=julian(1,1,1978)),
             catname="Wellington")
    
b <- subsetrect(Wellington, minlong=173.6, maxlong=176.1, minlat=-42.0,
                maxlat=-40.3, minmag=2, minday=julian(1,1,1978),
                maxday=julian(1,1,1992), maxdepth=200)

rotation(b, theta=-40, km=FALSE)
title(main="Wellington Plate Boundary")
[Package ssEDA version 2.2-1 Index]