[gradsusr] gancgrid error
Jennifer Adams
jma at cola.iges.org
Tue May 29 13:45:31 EDT 2012
Hi, Georgia --
There are two problems:
1. The time_counter coordinate has size 0, but the variables vomecrty and sometauy list this as one of their dimensions. I don't think I've ever seen a netcdf file like this; I would try removing the "t," from your variable declarations and see if that fixes anything.
2. The variables nav_lon and nav_lat reveal that the data are not on a rectilinear grid -- GrADS requires 1-dimensional grid coordinates.
If you want to look at the data in GrADS on an abstract grid that won't be georegistered, you will have to change your XDEF and YDEF statements to:
XDEF 198 linear 1 1
YDEF 224 linear 1 1
Then use 'set mproj off'. This is obviously not ideal.
I have been working on a solution to this problem, using PDEF and a regridding tool from ESMF called RegridWeightGen that will interpolate this kind of ocean model output on a curvilinear grid to a rectilinear grid that GrADS can handle. It is still in the testing/development stage. If you want to know more about how to make this work, please write to me offline.
--Jennifer
On May 29, 2012, at 12:26 PM, pigiottini at libero.it wrote:
> 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
> <AMM12_5d_20070101_20070101_grid_V.cdl>_______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120529/4e020d74/attachment-0003.html
More information about the gradsusr
mailing list