examples GridSearch.DoubleCouple.py - beachball.png station azimuths are rotated

[mikehagerty]

I'm trying to follow along with last year's virtual workshop on youtube. On day 1 at [1:13:52] Felix shows the beachball.png produced by the GridSearch.DoubleCouple example.

My results - waveform fits, misfit contours, solution s,d,r - are identical to what he shows, except the station triangles/names plotted on the beachball look rotated by 180 deg. e.g., the beachball is the same but the station azimuths are totally different. For instance, you have TRF plotting at 353 deg (towards north) and mine is rotated to the south.

I thought maybe this could be an upper/lower hemisphere focal mechanism thing, since the upper hemisphere plot is rotated 180deg from the lower hemisphere plot, but as I say, only the station locations (azimuths) are rotated, not the fault planes, so that doesn't make sense.

After some debugging ... I think you might have an error in your example. For instance, here is the tmp station file that polar1 plots:

cat tmp.20090407201255351DC_beachball.png.sta BIGB 345.527769 152.458534 PMR 64.454566 128.320995 KASH 338.666196 118.526754 TUPA 157.288772 105.902894 LSUM 170.675982 103.422575 MPEN 206.794148 103.171812 DEVL 175.365655 101.057717 BLAK 136.074440 100.572742 LSKI 200.106721 99.145813 NSKI 223.868000 98.472389 AVAL 169.692917 98.364355 PERI 129.875238 97.862826 SOLD 213.907196 97.457392 SWD 173.872713 96.070030 HEAD 173.413687 95.478760 DIV 97.918007 93.427289 TRF 353.014642 93.206074 EYAK 113.279988 92.876460 PAX 50.896740 92.049359 BMR 98.847989 92.010388

TRF has an azimuth of 353 but a takeoff angle of 93 which is > 90. so in fact, in the lower hemisphere plot, it is an upgoing ray and should plot 180 deg. away from the azimuth(or towards the south) as it does in my plot (see attached).

I don't think I have a way to know the takeoff angles you used, so maybe in a different velocity model it was < 90. for TRF ?

Thanks! -Mike

Versions: python: 3.11.0 mtuq: 0.2.0 pygmt: 0.8.0