Analyze your faults (compute the traction components and the CFF along them). Also compute the angular misfit between the observed or assumed slip direction and the computed shear traction orientation. Use the specified stress magnitudes for a start:
%assume that the principal stresses are appropriate for normal faulting
%conditions so maximum stress is the vertical stress
sv = -26.5.*12; %assume 26.5 MPa per km and 12 km depth
shmin = sv.*0.1; %assume the 1 direction is the minimum horizontal stress and is 10% sv
shmax = sv.*0.5; %assume the 2 direction is intermediate
S = [shmin 0 0;
0 shmax 0;
0 0 sv]