[ncl-talk] error reading netcdf4 file
Marston Johnston
shejo284 at gmail.com
Wed Nov 30 07:24:50 MST 2016
Hi Karin,
I have tried that and it did not work. I forgot to write that I tested that iteration as well.
fatal:syntax error: line 15 in file process_t511.ncl before or near .128
g = a->p100.128
----------------^
fatal:Syntax Error in block, block not executed
fatal:error at line 17 in file process_t511.ncl
I’m going to try using cdo to rename the variable to remove the “.” And test again.
But the link should be updated to remove this example where the variable name contains a “.”.
/M
From: <ncl-talk-bounces at ucar.edu> on behalf of Karin Meier-Fleischer <meier-fleischer at dkrz.de>
Date: Wednesday, 30 November 2016 at 15:18
To: Marston Johnston <marston.johnston at gu.se>
Cc: "ncl-talk [ncl-talk]" <ncl-talk at ucar.edu>
Subject: Re: [ncl-talk] error reading netcdf4 file
Hi Marston,
the variable name in your netCDF file is p100.128 and not /p100.128. NCLs error message is exactly telling you what is wrong in line 18.
Because the test script you send is not the script producing the error I guess that the line
g = short2flt(a->/p100.128)
is the problem. Change it to
g = short2flt(a->p100.128)
Bye,
Karin
Am 30.11.2016 um 14:59 schrieb Marston Johnston <marston.johnston at gu.se>:
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
_______________________________________________
ncl-talk mailing list
ncl-talk at ucar.edu
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________ ncl-talk mailing list ncl-talk at ucar.edu List instructions, subscriber options, unsubscribe: http://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161130/6eb34469/attachment.html
More information about the ncl-talk
mailing list