|
From: Eric F. <ef...@ha...> - 2013-12-01 19:20:28
|
On 2013/12/01 5:11 AM, lin...@eu... wrote: > Z = X**2 + Y**2 > Z[(X**2+Y**2 < r**2) | (X**2+Y**2 > R**2)] = np.ma.masked The problem is that Z is not a masked array. Change the first line above to Z = np.ma.array(X**2 + Y**2). Also, use pcolormesh instead of pcolor. Pcolormesh is much faster. Eric |