[gradsusr] gancgrid error

[pigiottini at libero.it]

Hi all,I have some output data in NetCDF4 format from NEMO model.
I'd like to display my output variables.

I used the "open" command to open the file.nc because "sdfopen" had failed. So I wrote a GrADS descriptor file .ctl 
When I try to display a variable, for example ga-> d vomecrty, I get the following error: 
  gancgrid error: nc_get_vara_double failed; NetCDF: Index exceeds dimension bound
  Data Request Error:  Error for variable 'vomecrty'
  Error ocurred at column 1
  DISPLAY error:  Invalid expression 
  Expression = vomecrty

This is my AMM12_5d_20070101_20070101_grid_V.ctl:
DSET ^AMM12_5d_20070101_20070101_grid_V.nc
DTYPE netcdf
TITLE 4-D Ocean Variables: Meridional Components
UNDEF 9.96921e+36f
OPTIONS zrev
XDEF 198 linear -19.8333 0.16667 
    
YDEF 224 linear 40.16667 0.11111
  
ZDEF 33 levels  3.107505 10.04362 18.13503 27.66259 38.97338 52.49526 
    68.75394 88.39286 112.1955 141.1104 176.2764 219.0491 271.0224 
    334.0439 410.218 501.8885 611.5969 742.007 895.795 1075.506 
    1283.386 1521.216 1790.151 2090.625 2422.302 2784.11 3174.342 
    3590.799 4030.965 4492.172 4971.748 5467.131 5975.949   
    
TDEF 0 linear 01jan2007 5dy
VARS 2
vomecrty    33  t,z,y,x   ocean current along j-axis(m/s)
sometauy     0  t,y,x  Wind Stress along j-axis (N/m^2)
ENDVARS
___
And This is the first part of my AMM12_5d_20070101_20070101_grid_V.cdl (output of ncdump -b):
netcdf AMM12_5d_20070101_20070101_grid_V {
dimensions:
	x = 198 ;
	y = 224 ;
	depthv = 33 ;
	time_counter = UNLIMITED ; // (0 currently)
	tbnds = 2 ;
variables:
	float nav_lon(y, x) ;
		nav_lon:standard_name = "longitude" ;
		nav_lon:units = "degrees_east" ;
		nav_lon:valid_min = 4.002733e-05f ;
		nav_lon:valid_max = 359.8334f ;
		nav_lon:long_name = "Longitude" ;
		nav_lon:nav_model = "Default grid" ;
	float nav_lat(y, x) ;
		nav_lat:standard_name = "latitude" ;
		nav_lat:units = "degrees_north" ;
		nav_lat:valid_min = 40.16667f ;
		nav_lat:valid_max = 64.94442f ;
		nav_lat:long_name = "Latitude" ;
		nav_lat:nav_model = "Default grid" ;
	float depthv(depthv) ;
		depthv:axis = "Z" ;
		depthv:standard_name = "model_level_number" ;
		depthv:units = "m" ;
		depthv:positive = "down" ;
		depthv:valid_min = 3.107505f ;
		depthv:valid_max = 5975.949f ;
		depthv:title = "depthv" ;
		depthv:long_name = "Vertical V levels" ;
	double time_counter(time_counter) ;
		time_counter:axis = "T" ;
		time_counter:standard_name = "time" ;
		time_counter:units = "seconds since 2007-01-01 00:00:00" ;
		time_counter:calendar = "gregorian" ;
		time_counter:title = "Time" ;
		time_counter:long_name = "Time axis" ;
		time_counter:time_origin = " 2007-JAN-01 00:00:00" ;
		time_counter:bounds = "time_counter_bnds" ;
	double time_counter_bnds(time_counter, tbnds) ;
	float vomecrty(time_counter, depthv, y, x) ;
		vomecrty:units = "m/s" ;
		vomecrty:standard_name = "ocean current along j-axis" ;
		vomecrty:_FillValue = 9.96921e+36f ;
		vomecrty:long_name = "ocean current along j-axis" ;
		vomecrty:online_operation = "ave(X)" ;
		vomecrty:interval_operation = 150.f ;
		vomecrty:interval_write = 432000.f ;
		vomecrty:coordinates = "time_counter depthv nav_lat nav_lon" ;
	float vomeeivv(time_counter, depthv, y, x) ;
		vomeeivv:units = "m/s" ;
		vomeeivv:standard_name = "EIV ocean current along j-axis" ;
		vomeeivv:_FillValue = 9.96921e+36f ;
		vomeeivv:long_name = "EIV ocean current along j-axis" ;
		vomeeivv:online_operation = "ave(X)" ;
		vomeeivv:interval_operation = 150.f ;
		vomeeivv:interval_write = 432000.f ;
		vomeeivv:coordinates = "time_counter depthv nav_lat nav_lon" ;
	float sometauy(time_counter, y, x) ;
		sometauy:units = "N/m2" ;
		sometauy:standard_name = "Wind Stress along j-axis" ;
		sometauy:_FillValue = 9.96921e+36f ;
		sometauy:long_name = "Wind Stress along j-axis" ;
		sometauy:online_operation = "ave(X)" ;
		sometauy:interval_operation = 150.f ;
		sometauy:interval_write = 432000.f ;
		sometauy:coordinates = "time_counter nav_lat nav_lon" ;

// global attributes:
		:Conventions = "CF-1.1" ;
		:file_name = "AMM12_5d_20070101_20070101_grid_V.nc" ;
		:production = "An IPSL model" ;
		:TimeStamp = "2012-MAY-25 12:53:46 GMT+0200" ;
 .......

Thanks in advance for your help!
Regards
Giorgia
 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120529/bba6abce/attachment-0003.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: AMM12_5d_20070101_20070101_grid_V.cdl
Type: application/octet-stream
Size: 918816 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20120529/bba6abce/attachment-0003.obj