rotation                package:ssEDA                R Documentation

_R_o_t_a_t_e_s _E_v_e_n_t_s _t_o _V_i_e_w _P_l_a_t_e _B_o_u_n_d_a_r_y

_D_e_s_c_r_i_p_t_i_o_n:

     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.

_U_s_a_g_e:

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

_A_r_g_u_m_e_n_t_s:

  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. 

_S_e_e _A_l_s_o:

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

_E_x_a_m_p_l_e_s:

     #    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")

