geodat.math

geodat.math.div(ux, uy, dx, dy, xaxis=-1, yaxis=-2)

Compute the divergence of a vector field ux,uy dx, dy are either 1-d array or a scalar xaxis - integer indicating the location of x axis (default = -1) yaxis - integer indicating the location of y axis (default = -2)

geodat.math.gradient(f, dx=1.0, axis=0, mask_boundary=False)

Compute central difference for a pariticular axis Input: f - numpy ndarray dx - spacing

geodat.math.integrate(data, axes=None, iax=None, **kwargs)

Integrate data along a selected set of axes Dimension is maintained using keepdims.sum data - numpy.ndarray axes - a list of axes iax - a list of integers that select which axes are integrated along