Update to VisAD library [SEC=UNCLASSIFIED]

[Tom Rink]

Hi Rob,

I missed this little problem because I normally use the
IDV to test changes/enhancements to the Contour classes
and most of the time there's an alpha component used.
In your case, there wasn't, and hence the Exception. 

If your application can be distributed, I'd be happy to
test changes against it as well.

Tom


Rob Hackett wrote:

> Hi
>
> I've had a look at the new release of the source code and I noticed 
> what seems like a minor bug relating to the message from James Kelly 
> which I have quoted below.
>
> The new version of the source code solves James' problem, but shortly 
> afterwards fails with  the following error ...
>
> java.lang.ArrayIndexOutOfBoundsException: 3
>         at visad.Contour2D.contour(Contour2D.java: )
>         at visad.Gridded3DSet.makeIsoLines(Gridded3DSet.java:2180)
>         at visad.ShadowType.makeContour(ShadowType.java:3846)
>         at visad.java3d.ShadowTypeJ3D.makeContour(ShadowTypeJ3D.java:473)
>         at 
> visad.ShadowFunctionOrSetType.doTransform(ShadowFunctionOrSetType.java:2141)
>         at 
> visad.java3d.ShadowFunctionOrSetTypeJ3D.doTransform(ShadowFunctionOrSetTypeJ3D.java:159)
>         at 
> visad.java3d.DefaultRendererJ3D.doTransform(DefaultRendererJ3D.java:134)
>         at visad.java3d.RendererJ3D.doAction(RendererJ3D.java:181)
>         at visad.DisplayImpl.doAction(DisplayImpl.java:1711)
>         at visad.ActionImpl.run(ActionImpl.java:364)
>         at visad.util.ThreadPool$ThreadMinnow.run(ThreadPool.java:95)
>
>
> Looking at the source for Contour2D.java I can see that this 
> corresponds to the following line ...
>
> auxLevels1[3] = new byte[numv1[0]+cnt];
>
> In this case, auxLevels is defined as having a length of 3, hence the 
> array exception.
>
> I notice that in the preceding lines of Contour2D.java, references to 
> auxLevels[3] are all protected by an "if statement" which ensures that 
> the array is big enough. And when I add a similar "if statement" to 
> line 1340 the problem disappears.
>
> Can you confirm that this is a safe solution to this problem?
>
> thanks,
>
> Rob
>
>
>
>
>
> From James Kelly ...
>
>
> Hi,
>
> Under some circumstances, I'm having trouble contouring a simple 2D
> FlatField, using the latest visad release. It works ok with a release
> from last year ( size 4444712 bytes,  Apr  2006).
>
> Here are the details in brief:
> I want to contour an atmospheric pressure field (P: mean sea level
> pressure) at 2hPa intervals. If I load a field that has a size of 80x104
> rows and columns it works ok. If I load the same field, but with a
> larger domain (112x120) it fails with the exception below:
>
> java.lang.ArrayIndexOutOfBoundsException: 21
>         at visad.ContourQuadSet.add(Contour2D.java:2968)
>         at visad.ContourStripSet.add(Contour2D.java:3456)
>         at visad.Contour2D.contour(Contour2D.java:1239)
>         at visad.Gridded3DSet.makeIsoLines(Gridded3DSet.java:2180)
>         at visad.ShadowType.makeContour(ShadowType.java:3805)
>
> It does work if I don't specify the 2hPa contour interval, and use a
> default interval. Also, as mentioned, it works with an earlier version
> of visad.
>
>
> Here's how to reproduce the problem:
>
> Download the following 2 files to your working directory:
> ftp://ftp.ssec.wisc.edu/pub/incoming/usmodel.nc
> ftp://ftp.ssec.wisc.edu/pub/incoming/usmodel_works.nc
>
>
> Save the jython program at the end of this email to your working
> directory (eg as test_contour.py).
> Run it using "./jy test_contour.py"
>
> [comment out the line "ci.setContourInterval(2.0, 0.0, 1050.0, 0.0)" and
> it works]
>
> Thanks for any help,
>
> James
>
> Jython program:
>
> from visad.python.JPythonMethods import *
> from visad import ScalarMap, Display
> from visad.data.netcdf import Plain
>
> import subs, graph, sys
>
> width = 800
> height = 600
>
> filename_works = "usmodel_works.nc"
> filename = "usmodel.nc"
> print "reading from netCDF..."
> netCDFReader = Plain()
> model = netCDFReader.open(filename)
> model_works = netCDFReader.open(filename_works)
>
> dom = domainType(model)
> dom_1 = dom.getComponent(0)
> dom_2 = dom.getComponent(1)
> rng = rangeType(model)
> rngMap = ScalarMap(rng, Display.IsoContour)
> xMap = ScalarMap(dom_1, Display.XAxis)
> yMap = ScalarMap(dom_2, Display.YAxis)
> maps = (xMap, yMap, rngMap)
>
> disp = subs.makeDisplay(maps)
>
> subs.addData("contours", model_works, disp)
> ci = rngMap.getControl()
> ci.enableLabels(1)
> ci.enableContours(1)
> ci.setContourInterval(2.0, 0.0, 1050.0, 0.0)
>
> subs.setBoxSize(disp, .80)
> subs.setAspectRatio(disp, float(width)/float(height))
> subs.showDisplay(disp,width,height,"contour test")
>
> # now show second model field which doesn't work:
> subs.addData("contours don't work", model, disp)
>
>

==============================================================================
To unsubscribe visad, visit:
http://www.unidata.ucar.edu/mailing-list-delete-form.html
==============================================================================