pyacs.lib.euler.pole_matrix_fault

Design matrix relating fault slip components to rotation rate vector.

pyacs.lib.euler.pole_matrix_fault.pole_matrix_fault(coor, strike, order=None)

Design matrix relating fault slip components to a rotation rate vector.

For n sites with [lon, lat] (decimal degrees) and strike (counter-clockwise from north), returns matrix W such that np.dot(W, w) is the flattened [ss1, ns1, ss2, ns2, …] in m/yr for rotation rate vector w in rad/yr.

Parameters:

• coor (array_like, shape (n, 2)) – Site coordinates as [longitude, latitude] in decimal degrees.

• strike (array_like, shape (n,) or (n, 1)) – Fault strike at each site, counter-clockwise from north, decimal degrees.

• order ({None, 'SS_DS'}, optional) – If provided, return matrix ordered as strike-slip then dip-slip rows.

Returns:

Pole matrix. np.dot(W, w).reshape(-1, 2) gives [along-strike, normal] components in two columns, in m/yr. Values refer to the hanging-wall block (right-hand rule polygon).

Return type:

ndarray, shape (2*n, 3)