xan/awips2
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Issue #189 add more comments to grib decoder.

Browse source 
Former-commit-id: 397ebf02f6 [formerly 67db4ae4e4] [formerly 237a7ea992 [formerly fdc27da79f1da6c2bda0f3061096227d95745c46]]
Former-commit-id: 237a7ea992
Former-commit-id: d57e097e3d
This commit is contained in:
Ben Steffensmeier 2012-11-06 18:29:07 -06:00
parent 76102aa826
commit 929df097c0
1 changed files with 10 additions and 2 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

12
edexOsgi/com.raytheon.edex.plugin.grib/GribDecoder.py
View file

@ -810,9 +810,17 @@ class GribDecoder():
dx = self._divideBy10e6(gdsTemplate[16])
dy = self._divideBy10e6(gdsTemplate[17])
if dx == 65.535:
# According to the grib2 spec 65.535 is completely valid, however it
# is impossible to define anything larger than a 5x2 grid with this
# spacing so we assume it is invalid and try to calculate a better
# value. 65.535 was chosen because it is the value encoded in the
# GFS161 model and it is completely wrong. This value is probably
# an artifact of converting from grib1 to grib2 since in grib1 this
# value would be encoded as an unsigned short with all bits as 1
# which is a special value in grib1, but in grib2 its just wrong
if dx >= 65.535:
dx = abs(lo1-lo2)/nx
if dy == 65.535:
if dy >= 65.535:
dy = abs(la1-la2)/ny
coverage.setSpacingUnit(DEFAULT_SPACING_UNIT2)
coverage.setNx(Integer(nx))