Return the cumulative sum of the elements along the given axis. The cumulative sum is calculated over the flattened array by default, otherwise over the specified axis.
Masked values are set to 0 internally during the computation. However, their position is saved, and the result will be masked at the same locations.
Parameters:  axis : {None, 1, int}, optional
dtype : {None, dtype}, optional
out : ndarray, optional


Returns:  cumsum : ndarray.

Notes
The mask is lost if out is not a valid MaskedArray !
Arithmetic is modular when using integer types, and no error is raised on overflow.
Examples
>>> marr = np.ma.array(np.arange(10), mask=[0,0,0,1,1,1,0,0,0,0])
>>> print marr.cumsum()
[0 1 3    9 16 24 33]