NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
NOTICE: This version of the NSF Unidata web site (archive.unidata.ucar.edu) is no longer being updated.
Current content can be found at unidata.ucar.edu.
To learn about what's going on, see About the Archive Site.
Hi Chris,I think it fails because the grid.makeSubset() creates a Variable that is a "section" view of the underlying grid variable. In other words, it isn't a fully contained new Variable, any data access on the new Variable reads from the underlying Variable. Our code doesn't allow a slice to be made on an already sectioned Variable. Instead it throws the exception you are seeing.
Your Array stuff works because you are actually reading the data into memory and creating a new Variable.
There may be a better way to deal with this. But I'm going to let John Caron, the author of most of this code, answer that question. I think he is back in the office next week.
Sorry I don't have a better answer for you, Ethan Christopher Mueller wrote:
Hi Ethan,
Thanks for the response. I actually tried Variable.slice() but it throws an
invalid range exception (below):
ucar.ma2.InvalidRangeException: Variable.slice: cannot slice a section
at ucar.nc2.Variable.makeSlice(Variable.java:475)
at ucar.nc2.dataset.VariableDS.slice(VariableDS.java:98)
I should probably mention that I'm obtaining the variable from a
GridDatatype ("grid" below) using GridDatatype.getVariable(). I know that
the index I'm trying to use in the slice is correct (levelRange is 1 value -
so there's only 1 value for the level dimension).
I have come up with a way of doing what I needed to do, but I'm not sure if
it's really a "good" way of doing it. It runs fine (and dumps the 1-length
dimension) perfectly using the code below. The "Variable sliceV..." line is
where it crashes (I comment it out when it works properly).
if((llbb != null) || (timeRange != null) || (stride_h >= 1))
grid = grid.makeSubset(timeRange, levelRange, llbb, 1, stride_h, stride_h);
Variable gridV = (Variable) grid.getVariable();
Variable sliceV = gridV.slice(gridV.findDimensionIndex(levelName), 0);
Array arrV = gridV.read().reduce(gridV.findDimensionIndex(levelName));
Variable newGridV = new Variable(gds.getNetcdfFile(), null, null,
gridV.getShortName());
newGridV.setDataType(gridV.getDataType());
for(Attribute a : (List<Attribute>)gridV.getAttributes()){
newGridV.addAttribute(a); }
Dimension levelD = gridV.getDimension(gridV.findDimensionIndex(levelName));
List<Dimension> dims = gridV.getDimensions();
if(levelD != null){ dims.remove(levelD); }
newGridV.setDimensions(dims);
newGridV.setCachedData(arrV, false);
Does this look like a reasonable way of doing it? Any idea why I'm getting
the error when I try to slice? That would be the preferred way in my mind -
it seems more flexible (I can do it to any of the dimensions easily)...
Thanks,
Best,
Chris
-- Ethan R. Davis Telephone: (303) 497-8155 Software Engineer Fax: (303) 497-8690 UCAR Unidata Program Center E-mail: edavis@xxxxxxxx P.O. Box 3000 Boulder, CO 80307-3000 http://www.unidata.ucar.edu/ ---------------------------------------------------------------------------
netcdf-java archives: