NetCDF reader bounds check incorrect
Affects | Status | Importance | Assigned to | Milestone | |
---|---|---|---|---|---|
Fluidity |
Confirmed
|
Wishlist
|
Stephan Kramer |
Bug Description
Fluidity will error with a 'Point ... not found ...' error when requesting data at a point on the boundary of a NetCDF file.
For example, when requesting z for x=-1.0, y=0.0 from a NetCDF file containing data z(x,y) for x,y in [-1.0,1.0].
To reproduce this behaviour, see test case 'netcdf_
Apply the following patch to tests/netcdf_
10,11c10,11
< x = arange(-1.2, 1.21, 0.2)
< y = arange(-1.2, 1.21, 0.2)
---
> x = arange(-1.0, 1.01, 0.2)
> y = arange(-1.0, 1.01, 0.2)
Related branches
Changed in fluidity: | |
status: | New → Confirmed |
Stephen has agreed to arrange a meeting for interested parties to consolidate and refactor the netcdf code. There are at least three different implementations of the reader and all will show this bug.