def nanprod(a, axis=None, dtype=None, out=None, keepdims=0):
        a, mask = _replace_nan(a, 1)
        return np.prod(a, axis=axis, dtype=dtype, out=out, keepdims=keepdims)
Return the product of array elements over a given axis treating Not a
dataset/ETHPy150Open pydata/xarray/xarray/core/npcompat.py/nanprod