Re: empty header in reading grib2

From: David Brown <dbrown_at_nyahnyahspammersnyahnyah>
Date: Tue, 24 Feb 2009 12:27:09 -0700

Hi YuGuang He,

I got your sample file, and NCL seems to read it pretty well, except
that there is one code table entry that it does not recognize.

ncl 0> f = addfile("gfs.t00z.pgrbf240.grib2","r")
warning: Entry (203) not found in code table file /usr/local/ncarg/
lib/ncarg/grib2_codetables/ncep/4/4.2.2.3.table
ncl 1> print(f)
Variable: f (file variable)

filename: gfs.t00z.pgrbf240
path: gfs.t00z.pgrbf240.grib2
    file global attributes:
    dimensions:
       lat_0 = 73
       lon_0 = 144
       lv_ISBL0 = 26
       lv_AMSL1 = 3
       lv_DBLL2 = 4
       lv_PVL3 = 2
       lv_ISBL4 = 21
       lv_SIGL5 = 4
       lv_ISBL6 = 2
       lv_ISBL7 = 6
    variables:
       float TMP_P0_L1_GLL0 ( lat_0, lon_0 )
          center : US National Weather Service - NCEP (WMC)
          production_status : Operational products
          long_name : Temperature
          units : K
          _FillValue : 1e+20
          grid_type : Latitude/longitude
          parameter_discipline_and_category : Meteorological
products, Temperature
          parameter_template_discipline_category_number : ( 0,
0, 0, 0 )
          level_type : Ground or water surface
          level : 0
          forecast_time : 240
          forecast_time_units : hours
          initial_time : 02/20/2009 (00:00)
... etc.

Are you saying that when you read it, there do not seem to be any
variables or dimensions?
That is quite strange.
My best suggestion would be to wait for 5.1.0, which is in final test
phase right now and will be released next week.
However, I truly do not know why you would have this sort of problem
with 5.0.0.
  -dave

On Feb 24, 2009, at 12:03 PM, YuGuang He wrote:

> Dear NCL Supports,
>
> I am having an issue of reading the grib2 format files from NCEP
> in the updated NCL version 5.0.0 on Fedora 9 LINUX. I was wondering
> if I could get some help and shed some light on it.
> Basically, when I opened a file by
>
> fh = addfile("gfs.t00z.pgrbf240.grib2","r")
> print(fh)
>
> it shows,
>
> Variable: fh (file variable)
>
> filename: gfs.t00z.pgrbf240
> path: gfs.t00z.pgrbf240.grib2
> file global attributes:
> dimensions:
> variables:
>
> which does not contain the header information I was looking for. I
> have been using version 4.3.0 all along for the past few years and
> there has not been such an issue for me. Also, there is no problem
> for me to read grib1 and netCDF files with this updated version
> 5.0.0 of NCL.
>
> Thank you very much for your help.
>
> - Yuguang He
>
> --
> contact : Yuguang He
> company : AER, Inc.
> email : yhe_at_aer.com
> subject : Technical Support
> phone : 781-761-2288 (ext. 217)
> facsimile : 781-761-2299
> address : 131 Hartwell Ave
> city : Lexington
> state : MA
> postalcode : 02421
> country : USA
> website : http://www.aer.com
> referrer : http://weather.aer.com/contact.php
>
> <yhe.vcf>_______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Feb 24 2009 - 12:27:09 MST

This archive was generated by hypermail 2.2.0 : Tue Mar 03 2009 - 09:53:57 MST