[gradsusr] Gribmap not MATCHING with SREF data

Dan Leins theedge981 at gmail.com
Sun Sep 9 08:46:22 EDT 2012


All,

I am working with SREF grib2 data and am extracting a subset of
probabilities (i.e. Probability of visibilities below 1 mile) from the
grib2 files found on NCEP's server (sref.pgrb212.prob_3hrly.grib2).
 Without getting into all the gorey details, I extract the probabilities
I'm interested in using wgrib2 and write out new grib2 files for each
forecast hour.  ("wgrib2 sref.pgrb212.prob_3hrly.grib2 -match ":${fcst}
hour" -grib SREF.f${hr}"). This gives me 30 individual grib2 files
containing probabilities of visibilities below 1 mile
(VIS.1609.f00, VIS.1609.f03, VIS.1609.f06, etc...).

I would like to use a control file template so I don't have to open 30
different files.  Running *g2ctl VIS.1609.f%f2 *gives me a control file
which looks like this:

dset ^VIS.1609.f%f2
index ^VIS.1609.f00.idx
undef 9.999E+20
title VIS.1609.f00
* produced by g2ctl v0.0.8.4
* command line options: VIS.1609.f%f2
* griddef=1:0:(185 x 129):grid_template=30:winds(grid): Lambert Conformal:
(185 x 129) input WE:SN output WE:SN res 8 Lat1 12.190000 Lon1 226.541000
LoV 265.000000 LatD 25.000000 Latin1 25.000000 Latin2 25.000000 LatSP
0.000000 LonSP 0.000000 North P
dtype grib2
options template
pdef 185 129 lccr 12.190000 -133.459 1 1 25.000000 25.000000 -95
40635.000000 40635.000000
xdef 274 linear -152.855107 0.377930376864113
ydef 133 linear 12.289391 0.369409090909091
tdef 30 linear 03Z09sep2012 3hr
zdef 1 linear 1 1
vars 1
VISsfc   0,1,0   0,19,0 ** surface Visibility [prob]
ENDVARS

This is what I expected. However when I run gribmap on this control file, I
don't get anything to MATCH.

[dan at leins devel]$ gribmap -v -i control
grib2map: scanning GRIB2 file: VIS.1609.f00
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f03
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f06
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f09
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f12
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f15
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f18
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f21
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f24
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f27
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f30
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f33
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f36
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f39
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f42
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f45
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f48
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f51
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f54
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f57
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f60
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f63
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f66
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f69
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f72
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f75
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f78
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f81
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f84
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: scanning GRIB2 file: VIS.1609.f87
1.1: lev1=1,0 var=0,19,0 pdt=5
grib2map: reached end of files
grib2map: Writing out the index file

Now, if I were to create individual control files for each file and run
gribmap, that seems to work fine. I've also verified that each grib file
contains data and it plots correctly in GrADS.  What did I do wrong that I
can't seem to template these data properly?

Thanks for any help!
Dan Leins
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120909/2410071d/attachment-0003.html 


More information about the gradsusr mailing list