[gradsusr] Question on polar GRIB2 data in southern hemisphere

[Frost, Mr. Michael]

Hi:

I am creating GRIB2 data for polar projected data in the southern hemisphere, but something has gone wrong and I can't figure out what.  I can do this fine for the northern hemisphere, but not the southern.

Here is what I am seeing...  Both grb2ctl.pl and the underlying wgrib2 say that it's in the northern hemisphere, but when I use the wgrib2 to go directly to the octet that hold the projection center it shows a value of 1 which means southern.

Attached is a picture of the domain.

>From grb2ctl.pl
===
[frost at typhoon src]$ grb2ctl.pl -verf a.grib2 
dset ^a.grib2
index ^a.grib2.idx
undef 9.999E+20
title a.grib2
dtype grib2
pdef 61 61 nps 1016.71127527321 -779.674789691898 1.439000 54
xdef 8 linear -129.896894 0.537539142857144
ydef 3 linear -71.040194 0.635782999999996
tdef 1 linear 14Z14jun2012 1mo
zdef 1 linear 1 1
vars 1
PRMSLmsl   0,101,0   0,3,1 ** mean sea level Pressure Reduced to MSL [Pa]
ENDVARS
===

>From wgrib2
===
[frost at typhoon src]$ wgrib2 -d 1 -V a.grib2 
1:0:vt=2012061414:mean sea level:120 min fcst:PRMSL Pressure Reduced to MSL [Pa]:
    ndata=3721:undef=0:mean=99840.6:min=96134.2:max=103220
    grid_template=20:winds(N/S):
        polar stereographic grid: (61 x 61) input WE:SN output WE:SN res 48
        North pole lat1 -70.496832 lon1 233.893152 latD -60.000000 lonV 1.439000 dx 54000.000000 m dy 54000.000000 m
===

>From wgrib2, going in directly to the octet that contains the Projection centre flag (http://www.nco.ncep.noaa.gov/pmb/docs/grib2/grib2_temp3-20.shtml)
===
[frost at typhoon src]$ wgrib2 -d 1 -get_byte 3 64 1 a.grib2 
1:0:3-64=1
===


Thanks,
Mike Frost
-------------- next part --------------
A non-text attachment was scrubbed...
Name: southpole.png
Type: image/png
Size: 43577 bytes
Desc: southpole.png
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20120614/61e81a9e/attachment-0003.png