[ncl-talk] Problem running ncl_convert2nc (Error unpacking GRIB record section 4)
Dennis Shea
shea at ucar.edu
Wed Dec 28 08:07:54 MST 2016
Hi Ligea
FYI ... a tool unrelated to NCL also fails.
------
There is another tool available from NCEP: wgrib2
http://www.cpc.ncep.noaa.gov/products/wesley/wgrib2/
There is an option to create netCDF
http://www.cpc.ncep.noaa.gov/products/wesley/wgrib2/netcdf.html
===
When I use this tool on the grib2 file you sent, it also fails:
%> wgrib2 gonzalo08l.2014101412.hwrfsat.global.0p25.f024.grb2 -netcdf gon.nc
---
The wgrib2 tool prints many warning lines like:
code_table_4.5a: product definition template #32 not supported
code_table_4.5b: product definition template #32 not supported
....
code_table_4.5b: product definition template #32 not supported
1:0:d=2014101412:var discipline=0 master_table=8 parmcat=5
parm=7:no_level:24 hour fcst:
code_table_4.5a: product definition template #32 not supported
...
code_table_4.5b: product definition template #32 not supported
2:279767:d=2014101412:var discipline=0 master_table=8 parmcat=5
parm=7:no_level:24 hour fcst:
code_table_4.5a: product definition template #32 not supported
...
etc, etc
Like 'ncl_convert2nc', wgrib2 produces just one variable within the output
netCDF. Specifically,
%> ncddump -h gon.nc
netcdf gon {
dimensions:
latitude = 721 ;
longitude = 1440 ;
time = UNLIMITED ; // (1 currently)
variables:
double latitude(latitude) ;
latitude:units = "degrees_north" ;
latitude:long_name = "latitude" ;
double longitude(longitude) ;
longitude:units = "degrees_east" ;
longitude:long_name = "longitude" ;
double time(time) ;
time:units = "seconds since 1970-01-01 00:00:00.0 0:00" ;
time:long_name = "verification time generated by wgrib2 function
verftime()" ;
time:reference_time = 1413288000. ;
time:reference_time_type = 3 ;
time:reference_date = "2014.10.14 12:00:00 UTC" ;
time:reference_time_description = "forecast or accumulated,
reference date is fixed" ;
time:time_step_setting = "auto" ;
time:time_step = 0. ;
float var0_5_7_no_level(time, latitude, longitude) ;
var0_5_7_no_level:_FillValue = 9.999e+20f ;
var0_5_7_no_level:short_name = "var0_5_7_no_level" ;
var0_5_7_no_level:long_name = "desc" ;
var0_5_7_no_level:level = "no_level" ;
var0_5_7_no_level:units = "unit" ;
// global attributes:
:Conventions = "COARDS" ;
:History = "created by wgrib2" ;
:GRIB2_grid_template = 0 ;
==========================================
On Tue, Dec 27, 2016 at 5:01 PM, Dennis Shea <shea at ucar.edu> wrote:
> 6.3.0 will unpack one variable also. The name of the variable will be
> different because (! wild guess !) the needed table was not complete or
> available.
>
> GRIB is a bit of a 'moving target' :-(
>
> I have opened a JIRA ticket: NCL-2533
>
> Unfortunately, this is not a great week to get things done.
>
> Maybe, an IO developer will look?
>
> D
>
>
>
>
>
> On Tue, Dec 27, 2016 at 4:48 PM, Ligia Bernardet - NOAA Affiliate <
> ligia.bernardet at noaa.gov> wrote:
>
>> Hi Dennis,
>>
>> There are 18 fields in the GRIB2 file. In the NetCDF file you sent back
>> to me there is 1.
>> But more importantly, I need to find a way to do this myself, so I can
>> process the multiple files I have.
>>
>> The machine I am using, a NOAA's R&D platform, I was able to find NCL
>> v6.3. Would this version be usable to read this GRIB2 file, until a newer
>> version is available?
>>
>> Thanks,
>> Ligia
>>
>>
>>
>> On Tue, Dec 27, 2016 at 4:11 PM, Dennis Shea <shea at ucar.edu> wrote:
>>
>>> Attached is a gzip'd version of the netCDF file.
>>>
>>> The 'soon-to-be-released' beta 6.4.0 returns a warning message. I will
>>> open a JIRA ticket.
>>> Note: I believe the numbers are correct.
>>>
>>> =========
>>>
>>> NCAR Command Language Version 6.4.0 -27Dec2016_0148
>>>
>>> warning:NclGRIB2: BRTEMP_P32 contains records that NCL cannot currently
>>> differentiate. One or more records will be ignored.
>>>
>>> Variable: f
>>> Type: file
>>> filename: gonzalo08l.2014101412 <(201)%20410-1412>.hwrfsat.
>>> global.0p25.f024
>>> path: /Users/shea/Data/GRIB/gonzalo08l.2014101412 <(201)%20410-1412>
>>> .hwrfsat.global.0p25.f024.grb2
>>> file global attributes:
>>> dimensions:
>>> lat_0 = 721
>>> lon_0 = 1440
>>> variables:
>>> float BRTEMP_P32_GLL0 ( lat_0, lon_0 )
>>> center : US National Weather Service - NCEP (WMC)
>>> production_status : Operational products
>>> long_name : Brightness temperature
>>> units : K
>>> _FillValue : 1e+20
>>> grid_type : Latitude/longitude
>>> parameter_discipline_and_category : Meteorological
>>> products, Long wave radiation
>>> parameter_template_discipline_category_number : ( 32,
>>> 0, 5, 7 )
>>> forecast_time : 0
>>> forecast_time_units : hours
>>> initial_time : 10/14/2014 (12:00)
>>>
>>> float lat_0 ( lat_0 )
>>> long_name : latitude
>>> grid_type : Latitude/Longitude
>>> units : degrees_north
>>> Dj : 0.25
>>> Di : 0.25
>>> Lo2 : 359.75
>>> La2 : -90
>>> Lo1 : 0
>>> La1 : 90
>>>
>>> float lon_0 ( lon_0 )
>>> long_name : longitude
>>> grid_type : Latitude/Longitude
>>> units : degrees_east
>>> Dj : 0.25
>>> Di : 0.25
>>> Lo2 : 359.75
>>> La2 : -90
>>> Lo1 : 0
>>> La1 : 90
>>>
>>> =========
>>> The min and max are:
>>>
>>> (0) Brightness temperature (K) : min=184.67 max=319.49
>>>
>>>
>>>
>>> On Tue, Dec 27, 2016 at 1:13 PM, Ligia Bernardet - NOAA Affiliate <
>>> ligia.bernardet at noaa.gov> wrote:
>>>
>>>> Hello,
>>>>
>>>> I have transferred file gonzalo08l.2014101412 <(201)%20410-1412>
>>>> .hwrfsat.global.0p25.f024.grb2 to /ftp/incoming.
>>>>
>>>> Thanks,
>>>> Ligia
>>>>
>>>> On Tue, Dec 27, 2016 at 1:06 PM, Dennis Shea <shea at ucar.edu> wrote:
>>>>
>>>>> [1]
>>>>> 6.1.2 (Feb. 2013) is almost 4 years old. Much has happened since.
>>>>>
>>>>> [2]
>>>>> To be able to determine the problem or if the current version of NCL
>>>>> has addressed the issue, we would need the file.
>>>>>
>>>>> ftp ftp.cgd.ucar.edu
>>>>> anonymous
>>>>> your_email
>>>>> cd incoming
>>>>> prompt
>>>>> put gonzalo08l.2014101412 <%28201%29%20410-1412>.hwrfsat.
>>>>> global.0p25.f024.grb2
>>>>> quit
>>>>>
>>>>> We can not look at /ftp/incoming, so ....
>>>>>
>>>>> Let ncl-talk know when the file has been transferred.
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> On Tue, Dec 27, 2016 at 12:56 PM, Ligia Bernardet - NOAA Affiliate <
>>>>> ligia.bernardet at noaa.gov> wrote:
>>>>>
>>>>>> Hello,
>>>>>>
>>>>>> I would like to use convert a GRIB2 file generated by NCEP's
>>>>>> Hurricane WRF (HWRF) model to a NetCDF file. However, I am encountering the
>>>>>> error below. Any suggestions?
>>>>>>
>>>>>> Thanks,
>>>>>> Ligia
>>>>>>
>>>>>>
>>>>>> *ncl_convert2nc gonzalo08l.2014101412
>>>>>> <(201)%20410-1412>.hwrfsat.global.0p25.f024.grb2*
>>>>>> Processing file: gonzalo08l.2014101412 <(201)%20410-1412>.hwrfsat.
>>>>>> global.0p25.f024.grb2...
>>>>>> Copyright (C) 1995-2013 - All Rights Reserved
>>>>>> University Corporation for Atmospheric Research
>>>>>> NCAR Command Language Version 6.1.2
>>>>>> The use of this software is governed by a License Agreement.
>>>>>> See http://www.ncl.ucar.edu/ for more details.
>>>>>> getpdstemplate: PDS Template 4.32 not defined.
>>>>>> warning:Error unpacking GRIB record section 4.
>>>>>> fatal:Could not open (./gonzalo08l.2014101412 <(201)%20410-1412>
>>>>>> .hwrfsat.global.0p25.f024.grb2)
>>>>>> Segmentation fault
>>>>>>
>>>>>> _______________________________________________
>>>>>> 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/20161228/927ce580/attachment.html
More information about the ncl-talk
mailing list