PDEF and HIRLAM
Kees Kling
ckling at PLANET.NL
Thu Apr 19 05:32:23 EDT 2007
Dear gradsusers,
like many of you I have some problems displaying my data due to rotated pole
coordinates
If I run wgrib -V mygribfile I get the next result
+++++++++++++
rec 1:4:date 2006082700 PRES kpds5=1 kpds6=105 kpds7=0 levels=(0,0) grid=255
0 m above gnd 12hr fcst:
PRES=Pressure [Pa]
timerange 0 P1 12 P2 0 TimeU 1 nx 406 ny 324 GDS grid 10 num_in_ave 0
missing 0
center 150 subcenter 0 process 8 Table 1
rotated LatLon grid lat -35.600000 to 29.000000 lon -34.400000 to 46.600000
nxny 131544 (406 x 324) dx 200 dy 200 scan 64 mode 136
transform: south pole lat -30.000000 lon -15.000000 rot angle 0.000000
min/max data 676.308 1043.46 num bits 17 BDS_Ref 676.308 DecScale 0
BinScale -8
++++++++++++++++
from which I can see that it is a rotetade grid
Running grib2ctl gives the following results
++++++++++++++++
dset ^./grib/H6E_18_2006082700_512_100
index ^./grib/H6E_18_2006082700_512_100.idx
undef 9.999E+20
title ./grib/H6E_18_2006082700_512_100
* produced by grib2ctl v0.9.12.5p39d
dtype grib 255
* Rotated LatLon grid: South pole lat -30.000000 lon -15.000000 rot angle
0.000000
ydef 324 linear -35.600000 0.2
xdef 406 linear -34.400000 0.2
tdef 1 linear 00Z28aug2006 1mo
zdef 1 linear 1 1
vars 3
DPT2m 0 17,105,2 ** 2 m above ground Dew point temp. [K]
PRES0m 0 1,105,0 ** 0 m above ground Pressure [Pa]
TMP2m 0 11,105,2 ** 2 m above ground Temp. [K]
ENDVARS
+++++++++++++++++++
Till now everything is correct and I can display it
Now I need the PDEF information for this grid. I'v read the manuals about it
but it is still not clear for me.
Is there anyone out there who can help me???
Thanks in advance
Kees Kling
More information about the gradsusr
mailing list