def GetAnomaly(x, axis):
	"""Computes the anomaly of array x along dimension axis.
	INPUTS:
	  x    - array to compute anomalies from
	  axis - axis along dimension for anomalies
	OUTPUTS:
	  x    - anomalous array
	""" #bring axis to the front
	xt = x - x.mean(axis)
	return xt