[ncl-talk] error reading netcdf4 file

[Marston Johnston]

Hi,

I writing a small script to read a netcdf4 file that was created by grib_api (version 1.17.0). But despite following the examples given here: http://www.ncl.ucar.edu/Applications/netcdf4.shtml. At first I thought the “.” in the variable name was causing the problem but in the examples given in the link above, such a file naming type was given. I’ve tried using $”/p100.128”$, treating it as a group: a=>/p100.128, but to no avail.

I can upload the file for testing.

****************** script *********************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
  idir="/home/ "
  fil="ICMGGglxc+2007var100.nc"
  setfileoption("nc", "FileStructure", "Advanced")
  setfileoption("nc", "Format",  "NetCDF4")
  a = addfile(idir+fil,"r")
  printVarSummary(a)
  vNames = getfilevarnames (a)
  print(vNames)
  g = short2flt(a->/p100.128)
  printVarInfo(g,"g")
end
*****************  results  ********************
[nclcodes]$ ncdump -k /home/ICMGGglxc+2007var100.nc
netCDF-4
[nclcodes]$
ncdump -g /p100.128 /home/ICMGGglxc+2007var100.nc
ncdump: /p100.128: No such group

[nclcodes]$ ncdump -h /home/ICMGGglxc+2007var100.nc
netcdf ICMGGglxc+2007var100 {
dimensions:
                longitude = 1024 ;
                latitude = 512 ;
                level = 11 ;
                time = UNLIMITED ; // (2920 currently)
variables:
                float longitude(longitude) ;
                                longitude:standard_name = "longitude" ;
                                longitude:long_name = "longitude" ;
                                longitude:units = "degrees_east" ;
                                longitude:axis = "X" ;
                float latitude(latitude) ;
                                latitude:standard_name = "latitude" ;
                                latitude:long_name = "latitude" ;
                                latitude:units = "degrees_north" ;
                                latitude:axis = "Y" ;
                double level(level) ;
                                level:long_name = "model_level_number" ;
                                level:axis = "Z" ;
                double time(time) ;
                                time:standard_name = "time" ;
                                time:long_name = "time" ;
                                time:units = "hours since 1900-1-1 00:00:00" ;
                                time:calendar = "standard" ;
                                time:axis = "T" ;
                short p100.128(time, level, latitude, longitude) ;
                                p100.128:long_name = "Experimental product" ;
                                p100.128:units = "~" ;
                                p100.128:grid_type = "gaussian" ;
                                p100.128:add_offset = 2826.8713701289 ;
                                p100.128:scale_factor = 0.0914882578243023 ;
                                p100.128:_FillValue = -32767s ;
                                p100.128:missing_value = -32767s ;

// global attributes:
                                :CDI = "Climate Data Interface version 1.7.0 (http://mpimet.mpg.de/cdi)" ;
                                :history = "Mon Nov 28 23:20:13 2016: cdo -O -v mergetime ICMGGglxc+200701_var100.nc ICMGGglxc+200702_var100.nc ICMGGglxc+200703_var100.nc ICMGGglxc+200704_var100.nc ICMGGglxc+200705_var100.nc ICMGGglxc+200706_var100.nc ICMGGglxc+200707_var100.nc ICMGGglxc+200708_var100.nc ICMGGglxc+200709_var100.nc ICMGGglxc+200710_var100.nc ICMGGglxc+200711_var100.nc ICMGGglxc+200712_var100.nc ICMGGglxc+2007var100.nc\n2016-11-28 22:01:41 GMT by grib_to_netcdf-1.17.0: grib_to_netcdf -k 3 -o ICMGGglxc+200712_var100.nc ICMGGglxc+200712_var100.grb" ;
                                :Conventions = "CF-1.6" ;
                                :CDO = "Climate Data Operators version 1.7.0 (http://mpimet.mpg.de/cdo)" ;
}
[nclcodes]$

[nclcodes]$ ncl process_t511.ncl
Copyright (C) 1995-2015 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 6.3.0
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.

Variable: a
Type: file
File path                :               /home/ICMGGglxc+2007var100.nc

Number of global attributes            :               4
Number of dimensions      :               4
Number of chunk_dimensions         :               1
Number of variables          :               5

Variable: vNames
Type: string
Total Size: 40 bytes
            5 values
Number of Dimensions: 1
Dimensions and sizes:        [5]
Coordinates:
(0)           /longitude
(1)           /latitude
(2)           /level
(3)           /time
(4)           /p100.128
fatal:["Execute.c":132]:variable (/p100.128) is not in file (/home/ICMGGglxc+2007var100.nc)
fatal:["Execute.c":6316]:variable (/p100.128) is not in file (a)
fatal:["Execute.c":8575]:Execute: Error occurred at or near line 18 in file process_t511.ncl

Cheers,
/M

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161130/a20f6198/attachment.html