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. <br>
<br>Here's the wgrib2 inventory of the CONUS 2.5 km RTMA, which works: <br>CONUS 2.5 km RTMA data: <a href="http://nomads.ncep.noaa.gov/pub/data/nccf/com/rtma2p5/prod/">http://nomads.ncep.noaa.gov/pub/data/nccf/com/rtma2p5/prod/</a> - using the "anal" files<br>
<span style="font-family: courier new,monospace;">~# wgrib2 rtma_anal_20111105_22z.gr2</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">1:0:d=2011110522:HGT:surface:anl:</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">2:5907518:d=2011110522:PRES:surface:anl:</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">3:11815036:d=2011110522:TMP:2 m above ground:anl:</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">4:16614930:d=2011110522:DPT:2 m above ground:anl:</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">5:19937992:d=2011110522:UGRD:10 m above ground:anl:</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">6:24737886:d=2011110522:VGRD:10 m above ground:anl:</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">7:29537780:d=2011110522:SPFH:2 m above ground:anl:</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">8:35076090:d=2011110522:WDI</span>works<span style="font-family: courier new,monospace;">R:10 m above ground:anl:</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">9:39506776:d=2011110522:WIND:10 m above ground:anl:</span><br>
<br>Using g2ctl results in the following control file:<br><span style="font-family: courier new,monospace;">dset ^rtma_anal_20111105_22z.gr2</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">index ^rtma_anal_20111105_22z.gr2.idx</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">undef 9.999E+20</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">title rtma_anal_20111105_22z.gr2</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">* produced by g2ctl v0.0.8.2</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">* command line options: rtma_anal_20111105_22z.gr2</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">* 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</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">dtype grib2</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">pdef 2145 1377 lccr 20.191999 -121.554001 1 1 25.000000 25.000000 -95 2539.703000 2539.703000</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">xdef 2814 linear -130.103387 0.0246001009814572</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">ydef 1413 linear 20.196426 0.0230882090909091</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">tdef 1 linear 22Z05nov2011 1mo</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">zdef 1 linear 1 1</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">vars 9</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">DPT2m 0,103,2 0,0,6 ** 2 m above ground Dew Point Temperature [K]</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">HGTsfc 0,1,0 0,3,5 ** surface Geopotential Height [gpm]</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">PRESsfc 0,1,0 0,3,0 ** surface Pressure [Pa]</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">SPFH2m 0,103,2 0,1,0 ** 2 m above ground Specific Humidity [kg/kg]</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">TMP2m 0,103,2 0,0,0 ** 2 m above ground Temperature [K]</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">UGRD10m 0,103,10 0,2,2 ** 10 m above ground U-Component of Wind [m/s]</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">VGRD10m 0,103,10 0,2,3 ** 10 m above ground V-Component of Wind [m/s]</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">WDIR10m 0,103,10 0,2,0 ** 10 m above ground Wind Direction (from which blowing) [deg]</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">WIND10m 0,103,10 0,2,1 ** 10 m above ground Wind Speed [m/s]</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">ENDVARS</span><br><br>Again, the CONUS data above works as expected. I can plot temp/dew point just fine.<br><br>The problem is with the Alaska data:<br><br>Alaska RTMA data: <a href="http://nomads.ncep.noaa.gov/pub/data/nccf/com/akrtma/prod/">http://nomads.ncep.noaa.gov/pub/data/nccf/com/akrtma/prod/</a> - using "anl" files<br>
<span style="font-family: courier new,monospace;">~# wgrib2 rtma_anal_ak_20111105_22z.gr2</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">1:0:d=2011110522:HGT:surface:anl:</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">2:912622:d=2011110522:PRES:surface:anl:</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">3:1825244:d=2011110522:TMP:2 m above ground:anl:</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">4:2566782:d=2011110522:DPT:2 m above ground:anl:</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">5:3137236:d=2011110522:UGRD:10 m above ground:anl:</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">6:3821746:d=2011110522:VGRD:10 m above ground:anl:</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">7:4506256:d=2011110522:SPFH:2 m above ground:anl:</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">8:5247794:d=2011110522:WDIR:10 m above ground:anl:</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">9:5932304:d=2011110522:WIND:10 m above ground:anl:</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">10:6616814:d=2011110522:PRES:surface:anl:</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">11:7472408:d=2011110522:TMP:2 m above ground:anl:</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">12:7985834:d=2011110522:DPT:2 m above ground:anl:</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">13:8499260:d=2011110522:UGRD:10 m above ground:anl:</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">14:9012686:d=2011110522:VGRD:10 m above ground:anl:</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">15:9526112:d=2011110522:SPFH:2 m above ground:anl:</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">16:10210622:d=2011110522:WDIR:10 m above ground:anl:</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">17:10552963:d=2011110522:WIND:10 m above ground:anl:</span><br>
<br>Using g2ctl results in:<br><span style="font-family: courier new,monospace;">dset ^rtma_anal_ak_20111105_22z.gr2</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">index ^rtma_anal_ak_20111105_22z.gr2.idx</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">undef 9.999E+20</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">title rtma_anal_ak_20111105_22z.gr2</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">* produced by g2ctl v0.0.8.2</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">* command line options: rtma_anal_ak_20111105_22z.gr2</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">* 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)</span><br style="font-family: courier new,monospace;">
<br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">dtype grib2</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">pdef 825 553 nps 441.000499063067 808.990794912376 210.000000 5.953125</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">xdef 2168 linear 150.190944 0.0535853193354869</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">ydef 651 linear 40.552720 0.0535970923076923</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">tdef 1 linear 22Z05nov2011 1mo</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">zdef 1 linear 1 1</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">vars 9</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">DPT2m 0,103,2 0,0,6 ** 2 m above ground Dew Point Temperature [K]</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">HGTsfc 0,1,0 0,3,5 ** surface Geopotential Height [gpm]</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">PRESsfc 0,1,0 0,3,0 ** surface Pressure [Pa]</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">SPFH2m 0,103,2 0,1,0 ** 2 m above ground Specific Humidity [kg/kg]</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">TMP2m 0,103,2 0,0,0 ** 2 m above ground Temperature [K]</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">UGRD10m 0,103,10 0,2,2 ** 10 m above ground U-Component of Wind [m/s]</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">VGRD10m 0,103,10 0,2,3 ** 10 m above ground V-Component of Wind [m/s]</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">WDIR10m 0,103,10 0,2,0 ** 10 m above ground Wind Direction (from which blowing) [deg]</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">WIND10m 0,103,10 0,2,1 ** 10 m above ground Wind Speed [m/s]</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">ENDVARS</span><br><br>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.<br>
<br>In GrADS, when I open up the control file, querying the file gives the following:<br><span style="font-family: courier new,monospace;">ga-> q file</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">File 1 : rtma_anal_ak_20111105_02z.gr2</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"> Descriptor: rtma_anal_ak_20111105_02z.ctl</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"> Binary: rtma_anal_ak_20111105_02z.gr2</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"> Type = Gridded</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"> Xsize = 2168 Ysize = 651 Zsize = 1 Tsize = 1 Esize = 1</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"> Number of Variables = 9</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"> dpt2m 0 0,0,6 ** 2 m above ground Dew Point Temperature [K]</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"> hgtsfc 0 0,3,5 ** surface Geopotential Height [gpm]</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"> pressfc 0 0,3,0 ** surface Pressure [Pa]</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"> spfh2m 0 0,1,0 ** 2 m above ground Specific Humidity [kg/kg]</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"> tmp2m 0 0,0,0 ** 2 m above ground Temperature [K]</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"> ugrd10m 0 0,2,2 ** 10 m above ground U-Component of Wind [m/s]</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"> vgrd10m 0 0,2,3 ** 10 m above ground V-Component of Wind [m/s]</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;"> wdir10m 0 0,2,0 ** 10 m above ground Wind Direction (from which blowing) [deg]</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;"> wind10m 0 0,2,1 ** 10 m above ground Wind Speed [m/s]</span><br style="font-family: courier new,monospace;">
<br>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!<br><br><span style="font-family: courier new,monospace;">ga-> d tmp2m</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">Notice: Automatic Grid Interpolation Taking Place</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">Data Type = grid</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">Dimensions = 0 1</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">I Dimension = 1 to 2168 Linear 150.191 0.0535853</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">J Dimension = 1 to 651 Linear 40.5527 0.0535971</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">Sizes = 2168 651 1411368</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">Undef value = -9.99e+08</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">Undef count = 518479 Valid count = 892889</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">Min, Max = 1.41458 4</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">Cmin, cmax, cint = 1.5 3.9 0.3</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">Stats[sum,sumsqr,root(sumsqr),n]: 3.49332e+06 1.37142e+07 3703.27 892889</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">Stats[(sum,sumsqr,root(sumsqr))/n]: 3.91238 15.3593 3.9191</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">Stats[(sum,sumsqr,root(sumsqr))/(n-1)]: 3.91239 15.3594 3.9191</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">Stats[(sigma,var)(n)]: 0.229326 0.0525902</span><br style="font-family: courier new,monospace;">
<span style="font-family: courier new,monospace;">Stats[(sigma,var)(n-1)]: 0.229326 0.0525903</span><br style="font-family: courier new,monospace;"><span style="font-family: courier new,monospace;">Contouring: 1.5 to 3.9 interval 0.3 </span><br>
<br>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 <i>it appears to me</i> that g2ctl is not functioning as expected. <br><br>Thank you,<br>Joe Moore<br>
<br>PS: I am using openGrADS with the following version of GrADS: 2.0.a9.oga.1<br><br>