ncgeodataset grid_interop problem bora_feb dataset

[GoogleCodeExporter]

grid_interop does not work with the gridded data in this dataset, even though 
this was a test dataset in the njTBX, and works okay in Tools-UI:

url 
='http://geoport.whoi.edu/thredds/dodsC/usgs/vault0/models/examples/bora_feb.nc'
;
nc=ncgeodataset(url);
tvar=nc.geovariable('temp');
t=tvar.data(2,14,:,:);   % "data" works okay
t=tvar.data(2,14,:,:);   % "grid" works okay
t=tvar.grid_interop(2,14,:,:);  % grid_interop bombs

No method 'getCoordinateArray' with matching signature found for class
'ucar.unidata.geoloc.vertical.VerticalTransformSubset'.

Error in ncgeovariable/grid_interop (line 169)
   array = subgrid.getCoordinateArray();

Error in ncgeovariable/subsref (line 627)
   sref = obj.grid_interop(first, last, stride);

Original issue reported on code.google.com by rsignell on 13 Oct 2011 at 6:30

[GoogleCodeExporter]

Seems to be an issue with the java code itself or some unidentified cdm problem 
within the cdm for this dataset in particular. Dataset appears to be 
conventional. Error occurs using java method on java object that works fine for 
other datasets.

Original comment by crosb...@gmail.com on 13 Oct 2011 at 7:44

[GoogleCodeExporter]

Checking the methodology that njtbx used to do the time dependent vertical 
coordinate transformation on this dataset.

Original comment by crosb...@gmail.com on 18 Oct 2011 at 1:49

• Changed state: Started

[GoogleCodeExporter]

Found a solution sort of. Can you test this out Rich?

I was using a java subsetting method on the vertical transform and then getting 
the values, but njtbx did subsetting in the custom java module after returning 
all of vertical coordinates for a given timestep, which seems to avoid the 
issue on this particular dataset.

Original comment by crosb...@gmail.com on 18 Oct 2011 at 4:59

• Changed state: Fixed

[GoogleCodeExporter]

Verification failed.

url 
='http://geoport.whoi.edu/thredds/dodsC/usgs/vault0/models/examples/bora_feb.nc'
;
nc=ncgeodataset(url);
tvar=nc.geovariable('temp');
t=tvar.grid_interop(2,14,:,:); 
size(t.z)
ans =

     2     1    60   160

This is incorrect.  Only 1 time value should be returned, not 2.  Size should 
be:
1 x 1 x 60 x 160.

Original comment by rsignell on 18 Oct 2011 at 5:49

[GoogleCodeExporter]

Original comment by crosb...@gmail.com on 18 Oct 2011 at 9:24

• Changed state: Verified