Get Rake Sense
Get the sense of the focal mechanism rake, from the U, V, P, T vectors
GetRakeSense(uaz, upl, vaz, vpl, paz, ppl, taz, tpl)
uaz |
Azimuth of U vector |
upl |
dip of U vector |
vaz |
Azimuth of V vector |
vpl |
dip of V vector |
paz |
Azimuth of P vector |
ppl |
dip of P vector |
taz |
Azimuth of T vector |
tpl |
dip of T vector |
1, 0 to make sure the region of the T-axis is shaded and the P-axis is blank.
The convention is for the T-axis to be shaded, so this subroutine determines the order of the polygons to be plotted so that the appropriate regins are filled.
Jonathan M. Lees <jonathan.lees@unc.edu>
GetRake
mc =CONVERTSDR(65,25,13) angsense = GetRakeSense(mc$U$az, mc$U$dip, mc$V$az, mc$V$dip,mc$P$az, mc$P$dip,mc$T$az, mc$T$dip)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.