[gradsusr] Gribmap not MATCHING with SREF data

[Wesley Ebisuzaki]

Dan,

   The SREF problem is a becoming a more common problem with grib2 files.
For example, the probability that the temperature is below freezing has a
variable name "TMP" with a modifier (Prob < 273.15K).  As you know this is
the same as the degrees K temperature according to GrADS.

   Jennifer, Brian and I have discussed this problem.  Jennifer and Brian
have
decided to add more numbers to the grib2 variable definition in the ctl
file.  No
idea on their progress.

   I am working on an alternative approach.  Here, the grib2 variable
definitions will
be wgrib2 extended names.  The alt_gmp (alt gribmap) is working except for
vertical profiles.  The alt_g2ctl is in rougher shape. Nevertheless I was
able to
display some sref fields.

*g2ctl VIS.1609.f%f2* problem:

Not enough information was provide to figure out your problem.   Try
using a -b option to g2ctl and gribmap.


             Wesley


On Sun, Sep 9, 2012 at 8:46 AM, Dan Leins <theedge981 at gmail.com> wrote:

> 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
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20120910/cc566106/attachment-0003.html