[gradsusr] Problem with Alaska RTMA (CONUS 2.5km RTMA works fine)
Joe Moore
joe at wxjoe.com
Sat Nov 5 20:01:25 EDT 2011
I'm having an issue that I believe lies in g2ctl, but I'm not absolutely
sure. I've been making images from the CONUS 2.5 km RTMA data from NCEP
just fine for the past month or so. This weekend I decided to tackle the
Alaska RTMA data, but I've hit an issue.
Here's the wgrib2 inventory of the CONUS 2.5 km RTMA, which works:
CONUS 2.5 km RTMA data:
http://nomads.ncep.noaa.gov/pub/data/nccf/com/rtma2p5/prod/ - using the
"anal" files
~# wgrib2 rtma_anal_20111105_22z.gr2
1:0:d=2011110522:HGT:surface:anl:
2:5907518:d=2011110522:PRES:surface:anl:
3:11815036:d=2011110522:TMP:2 m above ground:anl:
4:16614930:d=2011110522:DPT:2 m above ground:anl:
5:19937992:d=2011110522:UGRD:10 m above ground:anl:
6:24737886:d=2011110522:VGRD:10 m above ground:anl:
7:29537780:d=2011110522:SPFH:2 m above ground:anl:
8:35076090:d=2011110522:WDIworksR:10 m above ground:anl:
9:39506776:d=2011110522:WIND:10 m above ground:anl:
Using g2ctl results in the following control file:
dset ^rtma_anal_20111105_22z.gr2
index ^rtma_anal_20111105_22z.gr2.idx
undef 9.999E+20
title rtma_anal_20111105_22z.gr2
* produced by g2ctl v0.0.8.2
* command line options: rtma_anal_20111105_22z.gr2
* griddef=1:0:(2145 x 1377):grid_template=30:winds(grid): Lambert
Conformal: (2145 x 1377) input WE:SN output WE:SN res 8 Lat1 20.191999 Lon1
238.445999 LoV 265.000000 LatD 25.000000 Latin1 25.000000 Latin2 25.000000
LatSP -90.000000 LonSP 0.000000 N
dtype grib2
pdef 2145 1377 lccr 20.191999 -121.554001 1 1 25.000000 25.000000 -95
2539.703000 2539.703000
xdef 2814 linear -130.103387 0.0246001009814572
ydef 1413 linear 20.196426 0.0230882090909091
tdef 1 linear 22Z05nov2011 1mo
zdef 1 linear 1 1
vars 9
DPT2m 0,103,2 0,0,6 ** 2 m above ground Dew Point Temperature [K]
HGTsfc 0,1,0 0,3,5 ** surface Geopotential Height [gpm]
PRESsfc 0,1,0 0,3,0 ** surface Pressure [Pa]
SPFH2m 0,103,2 0,1,0 ** 2 m above ground Specific Humidity [kg/kg]
TMP2m 0,103,2 0,0,0 ** 2 m above ground Temperature [K]
UGRD10m 0,103,10 0,2,2 ** 10 m above ground U-Component of Wind [m/s]
VGRD10m 0,103,10 0,2,3 ** 10 m above ground V-Component of Wind [m/s]
WDIR10m 0,103,10 0,2,0 ** 10 m above ground Wind Direction (from which
blowing) [deg]
WIND10m 0,103,10 0,2,1 ** 10 m above ground Wind Speed [m/s]
ENDVARS
Again, the CONUS data above works as expected. I can plot temp/dew point
just fine.
The problem is with the Alaska data:
Alaska RTMA data:
http://nomads.ncep.noaa.gov/pub/data/nccf/com/akrtma/prod/- using
"anl" files
~# wgrib2 rtma_anal_ak_20111105_22z.gr2
1:0:d=2011110522:HGT:surface:anl:
2:912622:d=2011110522:PRES:surface:anl:
3:1825244:d=2011110522:TMP:2 m above ground:anl:
4:2566782:d=2011110522:DPT:2 m above ground:anl:
5:3137236:d=2011110522:UGRD:10 m above ground:anl:
6:3821746:d=2011110522:VGRD:10 m above ground:anl:
7:4506256:d=2011110522:SPFH:2 m above ground:anl:
8:5247794:d=2011110522:WDIR:10 m above ground:anl:
9:5932304:d=2011110522:WIND:10 m above ground:anl:
10:6616814:d=2011110522:PRES:surface:anl:
11:7472408:d=2011110522:TMP:2 m above ground:anl:
12:7985834:d=2011110522:DPT:2 m above ground:anl:
13:8499260:d=2011110522:UGRD:10 m above ground:anl:
14:9012686:d=2011110522:VGRD:10 m above ground:anl:
15:9526112:d=2011110522:SPFH:2 m above ground:anl:
16:10210622:d=2011110522:WDIR:10 m above ground:anl:
17:10552963:d=2011110522:WIND:10 m above ground:anl:
Using g2ctl results in:
dset ^rtma_anal_ak_20111105_22z.gr2
index ^rtma_anal_ak_20111105_22z.gr2.idx
undef 9.999E+20
title rtma_anal_ak_20111105_22z.gr2
* produced by g2ctl v0.0.8.2
* command line options: rtma_anal_ak_20111105_22z.gr2
* griddef=1:0:(825 x 553):grid_template=20:winds(grid): polar stereographic
grid: (825 x 553) input WE:SN output WE:SN res 8 North pole lat1 40.530101
lon1 181.429000 latD 60.000000 lonV 210.000000 dx 5953.125000 m dy
5953.125000 m:winds(grid)
dtype grib2
pdef 825 553 nps 441.000499063067 808.990794912376 210.000000 5.953125
xdef 2168 linear 150.190944 0.0535853193354869
ydef 651 linear 40.552720 0.0535970923076923
tdef 1 linear 22Z05nov2011 1mo
zdef 1 linear 1 1
vars 9
DPT2m 0,103,2 0,0,6 ** 2 m above ground Dew Point Temperature [K]
HGTsfc 0,1,0 0,3,5 ** surface Geopotential Height [gpm]
PRESsfc 0,1,0 0,3,0 ** surface Pressure [Pa]
SPFH2m 0,103,2 0,1,0 ** 2 m above ground Specific Humidity [kg/kg]
TMP2m 0,103,2 0,0,0 ** 2 m above ground Temperature [K]
UGRD10m 0,103,10 0,2,2 ** 10 m above ground U-Component of Wind [m/s]
VGRD10m 0,103,10 0,2,3 ** 10 m above ground V-Component of Wind [m/s]
WDIR10m 0,103,10 0,2,0 ** 10 m above ground Wind Direction (from which
blowing) [deg]
WIND10m 0,103,10 0,2,1 ** 10 m above ground Wind Speed [m/s]
ENDVARS
While I had some trouble getting the lat lon correct (which I figured out
by trial-and-error), I can't seem to plot the temperature and dew point.
While I can plot the tmp2m, the data is NOT the temperature; it's actually
the temperature error. (Record 11 of the grib2 - not sure if my terminology
is correct, but it's the "11" line of the wgrib2 output above.) I want to
plot the tmp2m (the actually temperature), but there's only one temp in the
control file.
In GrADS, when I open up the control file, querying the file gives the
following:
ga-> q file
File 1 : rtma_anal_ak_20111105_02z.gr2
Descriptor: rtma_anal_ak_20111105_02z.ctl
Binary: rtma_anal_ak_20111105_02z.gr2
Type = Gridded
Xsize = 2168 Ysize = 651 Zsize = 1 Tsize = 1 Esize = 1
Number of Variables = 9
dpt2m 0 0,0,6 ** 2 m above ground Dew Point Temperature [K]
hgtsfc 0 0,3,5 ** surface Geopotential Height [gpm]
pressfc 0 0,3,0 ** surface Pressure [Pa]
spfh2m 0 0,1,0 ** 2 m above ground Specific Humidity [kg/kg]
tmp2m 0 0,0,0 ** 2 m above ground Temperature [K]
ugrd10m 0 0,2,2 ** 10 m above ground U-Component of Wind [m/s]
vgrd10m 0 0,2,3 ** 10 m above ground V-Component of Wind [m/s]
wdir10m 0 0,2,0 ** 10 m above ground Wind Direction (from which
blowing) [deg]
wind10m 0 0,2,1 ** 10 m above ground Wind Speed [m/s]
I know that the tmp2m in the control file is incorrect because when I do a
plot with "set stat on," I see the values only range from ~1-4 Kelvin!
ga-> d tmp2m
Notice: Automatic Grid Interpolation Taking Place
Data Type = grid
Dimensions = 0 1
I Dimension = 1 to 2168 Linear 150.191 0.0535853
J Dimension = 1 to 651 Linear 40.5527 0.0535971
Sizes = 2168 651 1411368
Undef value = -9.99e+08
Undef count = 518479 Valid count = 892889
Min, Max = 1.41458 4
Cmin, cmax, cint = 1.5 3.9 0.3
Stats[sum,sumsqr,root(sumsqr),n]: 3.49332e+06 1.37142e+07 3703.27 892889
Stats[(sum,sumsqr,root(sumsqr))/n]: 3.91238 15.3593 3.9191
Stats[(sum,sumsqr,root(sumsqr))/(n-1)]: 3.91239 15.3594 3.9191
Stats[(sigma,var)(n)]: 0.229326 0.0525902
Stats[(sigma,var)(n-1)]: 0.229326 0.0525903
Contouring: 1.5 to 3.9 interval 0.3
Is this a problem that I should solve myself? I can think of a quick and
dirty fix using a shell script with sed, but *it appears to me* that g2ctl
is not functioning as expected.
Thank you,
Joe Moore
PS: I am using openGrADS with the following version of GrADS: 2.0.a9.oga.1
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20111105/087eaf47/attachment-0003.html
More information about the gradsusr
mailing list