grib2 fie read

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
Received on Fri Sep 21 14:43:01 2012