Yes, that's the behaviour that I expect setting the 'shrink' keyword to 'False' > Now, just to be clear, you'd want > 'np.ma.masked_values(...,shrink=False) to create a maked array w/ a > full boolean mask by default, right ?