Re: grib2 fie read

Hi CL,

I believe this is related to a problem that another user reported to me offline today. I've filed a ticket on this (NCL-1550 for your reference).

Hopefully we'll have a fix for you early next week.

--Mary

On Sep 21, 2012, at 2:42 PM, Luo, Chao wrote:

> Hi,
>
> I am reading biomass burning area of NOAA satellite grib2 data by using ncl-6.1.0. I print out the latitude of the grib2, and found maximum of latitude is 1.e+12, which is wrong. The error message when I reading the file is:
>
> warning: Entry (3) not found in code table file /data11/cluo/GAfire/ncl-
> 6.1.0/lib/ncarg/grib2_codetables/ncep/4/4.2.2.4.table
>
> I was wondering if this the problem caused by this warning. The NCL script I used is very simple:
>
> begin
> f = addfile("../NOAA_HMS/march_april/grib2/burned_area_20120316_18_23_12km.grib2","r")
> lat = f->gridlat_0(:,:)
> lon = f->gridlon_0(:,:)
> area = f->VAR_2_4_3_P31_GLC0(:,:)
> printMinMax (lat, True)
> printVarSummary (lat)
> end
> ~
> and output is :
>
> Variable: lat
> Type: float
> Total Size: 2987436 bytes
> 746859 values
> Number of Dimensions: 2
> Dimensions and sizes: [ygrid_0 | 723] x [xgrid_0 | 1033]
> Coordinates:
> Number Of Attributes: 11
> corners : ( 45.7454, 40.18576, 1e+12, 1e+12 )
> long_name : latitude
> grid_type : Lambert Conformal (secant, tangent, conical or bipolar)
> units : degrees_north
> Latin2 : 60
> Latin1 : 20
> Dy : 12
> Dx : 12
> Lov : 264
> Lo1 : 161.7229
> La1 : 45.7454
> (0)
> (0) latitude: min=36.144 max=1e+12
>
>
> Thanks very much,
>
> CL
> _______________________________________________
> 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 Fri Sep 21 16:29:04 2012