verticalfault class

class csi.verticalfault(name, utmzone=None, ellps='WGS84', lon0=None, lat0=None)
BuildPatchesVarResolution(depths, Depthpoints, Resolpoints, interpolation='linear', minpatchsize=0.1, extrap=None)

Patchizes the fault with a variable patch size at depth. The variable patch size is given by the respoints table. Depthpoints = [depth1, depth2, depth3, …., depthN] Resolpoints = [Resol1, Resol2, Resol3, …., ResolN] The final resolution is interpolated given the ‘interpolation’ method. Interpolation can be ‘linear’, ‘cubic’.

associatePatch2PDFs(directory='.', prefix='step_001_param')

Associates a patch with a pdf called directory/prefix_{#}.dat. import AltarExplore….

build_patches()

Builds rectangular patches from the discretized fault. A patch is a list of 4 corners.

computetotalslip()

Computes the total slip.

cutPatchesVertically(iP, cuttingDepth)

Cut a patche into 2 patches at depth given by cuttingDepth. Args:

  • iP : patch index or list of patch indexes.

  • cuttingDepth : Depth where patch is going to be split in 2.

differentiateGFs(datas)

Uses the Delaunay triangulation to prepare a differential Green’s function matrix, data vector and data covariance matrix. Args:

  • datas : List of dataset concerned

extrap1d(interpolator)

Linear extrapolation routine. Found on StackOverflow by sastanin.

extrapolate(length_added=50, tol=2.0, fracstep=5.0, extrap='ud')

Extrapolates the surface trace. This is usefull when building deep patches for interseismic loading. Args:

  • length_added : Length to add when extrapolating.

  • tol : Tolerance to find the good length.

  • fracstep : control each jump size.

  • extrapif u in extrap -> extrapolates at the end

    if d in extrap -> extrapolates at the beginning default is ‘ud’

rotationHoriz(center, angle)

Rotates the geometry of the fault around center, of an angle. Args:

  • center : [lon,lat]

  • angle : degrees

setDepth(depth, top=0, num=5)

Set the maximum depth of the fault patches.

Args:

  • depth : Depth of the fault patches.

  • num : Number of fault patches at depth.

translationHoriz(dx, dy)

Translates the patches. Args:

  • dx : Translation along x (km)

  • dy : Translation along y (km)