LPO and SWS plot

Q&A
1

I am trying to plot a comparison between LPO at different depths and shear wave splitting measurements to see variations in a GMT map.
I am not sure how can I plot my LPO set (see format below) to compare with the traditional SWS (Azimuth/Strenght).

Any suggestion.

Cheers,

Lat Long. Depth. TIx TIy TIz
1.151 -10.037 181.9000 1.1142 0.2043 -0.6533
1.764 -10.037 188.3000 1.1676 0.2440 -0.7084
2.464 -10.037 206.8000 1.1855 0.2847 -0.7496
3.724 -10.037 182.0000 1.4140 0.2540 -1.1025
-69.701 -10.043 46.2000 0.6203 -0.5051 1.3565
-63.169 -10.043 175.6000 -0.0328 1.5369 0.6394
-62.384 -10.043 161.8000 0.0032 1.4176 0.5972
-61.696 -10.043 158.5000 0.0247 1.3230 0.5688

2

I’d convert to azimuth and projected strength with something like
gawk ‘BEGIN{pi = 4*atan2(1,1);pif=180/pi} {if(NR>1){amp=sqrt($42+$52+$62);hamp=sqrt(amp2-$6**2);print($2,$3,atan2($4,$5)*pif,hamp)}’

1 Like
3

should have used inline code, the post turned some squared properties into bold font for some reason, but you get it

4

thanks @thwbecker. I just wondered if hamp=sqrt(amp-$62) or hamp=sqrt(amp2-$62) or hamp=sqrt(amp2-$62) ?

same for print($1,$2,atan2($4,$5)*pif,hamp)}’ or print($2,$3,atan2($4,$5)*pif,hamp)}’.

5

Ok to be sure I understood well your comment @thwbecker :
Azimuth= atan2(TIx,TIy)*180/pi
Strength= sqrt( (TIx.TIx+TIy.TIy+TIz.TIz)- TIz.Tiz)

where TIx=$4, TIy=$5, TIz=$6.

6

yes, so as to project the vector

1 Like