astropy.utils
astropy.utils.masked
Like numpy.nancumsum, skipping masked values as well.
numpy.nancumsum
Masked/NaN values are replaced with 0. The output is not masked.