A three-dimensional p x q x n array is referred to as a field, with first two space dimensions p and q, respectively, (e.g. longitude and latitude) and third time dimension n. A matrix p x q is referred to as a map, with first and second dimensions p and q, respectively, (e.g. longitude and latitude). A three-dimensional p x q x n array is referred to as a field, with first two space dimensions p and q, respectively, (e.g. longitude and latitude) and third time dimension n. A matrix p x q is referred to as a map, with first and second dimensions p and q, respectively, (e.g. longitude and latitude). Compute average cluster size for a given three-dimensional p x q x n array of excesses. First two dimensions p and q are space dimensions (e.g. longitude and latitude). Third dimension n is time. Calculate anomalies by subtracting mean annual cycle or long-term mean of a three-dimensional p x q x n array of monthly mean values. Compute basic statistics (e.g. mean, variance, min, max, skewness, quantile) for a given p x q x n three-dimensional array. Same as above, but also allows specification of fraction of non-missing values for the computation of the statistics. Grid points with larger fraction of missing values than specified are excluded. Compute upper bound of excesses for a given 2 x p x q array of Generalized Pareto distribution parameters. First index of the first dimension of the array represents the scale parameter. Second index of the first dimension of the array represents the shape parameter. Extract a subset of data from either a p x q two-dimensional map (matrix) or a p x q x n three-dimensional array (field). Flip a p x q two-dimensional map (matrix) or a p x q x n three-dimensional array up side down. Flip either latitude or longitude or both dimensions of a map or a three-dimensional array. Identify points by clicking on an image plot of the first time slice of a given p x q x n three-dimensional array, or an image plot of a p x q two-dimensional matrix. Also apply a user specified function to each selected points of the three-dimensional array. Write two-dimensional map or three-dimensional array of data, longitude vector, latitude vector and time vector into a netcdf file. Map two p x q dimensional data matrix in either cylindrical equidistant latitude and longitude projection or stereographic projection. Reshape a n x p*q matrix in a three-dimensional n x p x q array or reshape three-dimensional n x p x q array into a n x p*q matrix. Compute return period for a given p x q matrix of excesses and a given 2 x p x q array of Generalized Pareto distribution parameters. Compute extreme dependence measures between a given p x q x n three-dimensional array and a given time series of length n. Same as above, but also allows specification of fraction of non-missing values for the computation of the statistics. Grid points with larger fraction of missing values than specified are excluded. Compute location, shape and scale parameters of a Generalized Extreme Value Distribution for block annual maxima or minima of a given p x q x n three-dimensional array. Compute the intervals estimator for the extremal index, an index for time clusters, for a given time series and threshold. Fit Generalized Pareto distribution with time-varying threshold at each grid point for a given p x q x n three-dimensional array of montly data. Fit Generalized Pareto distribution with time-varying threshold at each grid point for a given p x q x n three-dimensional array of montly data. Allows linear modelling of the paramters. Zoom in on either a p x q two-dimensional map (matrix) or a p x q x n three-dimensional array (field), by selecting two corners that define a rectangular region of the space and a subset of time slices. Source.