Isofill Blank Spots on Solid with Indices/Levels

[chaosphere2112]

@aashish24 @sankhesh @doutriaux1 Just ran into this one while testing my performance branch, then went back and found it on master.

import vcs, cdms2
x = vcs.init()
f = cdms2.open(vcs.sample_data + "/clt.nc")
s = f("clt")

iso = x.createisofill()
iso.levels = range(0, 100, 10)
iso.fillareaindices = range(1, 11)
iso.ext_2 = True

x.plot(s, iso)

will give you this:

Obviously, fillareaindices shouldn't have any effect on a solid isofill, so that's a pretty weird behavior.

[chaosphere2112]

The ext_2 is insignificant, but was in the script that I used to generate the image.

[chaosphere2112]

Also with patterns:

import vcs, cdms2
x = vcs.init()
f = cdms2.open(vcs.sample_data + "/clt.nc")
s = f("clt")

iso = x.createisofill()
iso.levels = range(0, 110, 10)
iso.fillareastyle = "pattern"
iso.fillareaindices = range(1, 11)
iso.fillareacolors = range(20, 220, 20)
iso.fillareaopacity = [100] * 10

x.plot(s, iso)

[sankhesh]

@danlipsa Could you please look into this next?

[aashish24]

@danlipsa I am hoping this would a tiny bug somewhere. @doutriaux1 we may have to patch our release.

[danlipsa]

Seems the problem is how vtkBandedPolyDataContourFilter fills a band. Attached are two pictures: first isocontours at 40 and 50 then a band for 40-50. The triangle below the point labeled 53.9 never gets filled.