geosubset problem with stretched vertical coordinates?

[GoogleCodeExporter]

url=['http://michigan.glin.net:8080/thredds/dodsC/',...
'glos/glcfs/michigan/ncasf_his3d_fc/',...
'Lake_Michigan_Forecast_Collection_3D_best.ncd']
nc = ncgeodataset(url);
nc.variables
tvar = nc.geovariable('temp');
%geosubset takes a structure as input
s.lon = [-88.0244  -85.9673];
s.lat = [41.5265   45.2357];
%  s.time = [now_utc-3/24  now];
s.time = {'10-Nov-2011 09:00:00', '10-Nov-2011 12:00:00'};
s.h_stride = [2 2];
s.z_index=[1 5]; 
s.v_stride=2; 
tsub = tvar.geosubset(s);

Error using ncgeovariable/grid_interop (line 180)
There is a problem applying the vertical coordinate tranform and
subsetting the resuting values.

Error in ncgeovariable/geoij (line 516)
            g = obj.grid_interop(first, last, stride);

Error in ncgeovariable/geosubset (line 384)
            [indstart_r indend_r indstart_c indend_c] =
            obj.geoij(struct);

Error in ncgeovariable/subsref (line 671)
                            sref = builtin('subsref',obj,s);

Error in glos (line 15)
tsub = tvar.geosubset(s);

Original issue reported on code.google.com by rsignell on 8 Nov 2011 at 3:19

• Merged into: #27

[GoogleCodeExporter]

Can other implementations of netcdf-java calculate the vertical transform?

Original comment by crosb...@gmail.com on 8 Nov 2011 at 4:16

[GoogleCodeExporter]

Original comment by crosb...@gmail.com on 8 Nov 2011 at 4:19

• Changed state: Duplicate