This Python Script Arctool plot earthquake focal mechanism beach balls as shapefile. Data needed is earthquakes catalog point featureclass including attributes of lon, lat, M,strike, dip, rake, all in numerical value. The fields name can be any.
Polarity is stored in 'pola' field on which the beach balls and PT axis can be symbolized after balls formed.
The original codes is modified so as no advanced arcmap function of "Feature to Polygon" is needed.
V2.0 modification
a. Repaired a bug: acos domain overflow arised by floating point truncation error.
b. A smoothness option added to make balance between smoothness and speed.
c. Dispersing beach balls with earthquake focals.