[Met_help] Re: Problem with point-stat procedure

[John Halley Gotway]

Valerio,

I've updated 2 library files in MET, and they are attached.

Please copy the attached files "read_grib.cc" and "read_grib.h" into the directory METv1.1/lib/vx_met_util and rebuild MET.

With these changes, I'm still not getting exactly the same forecast value you're seeing.
Instead, for lat/lon = 43.72111, 10.39556, the nearest x/y is (44, 88) and the value at that point is 96.89951.

Does that match what you're seeing, or is there still a problem?

Thanks,
John

capecchi wrote:
> John, I'm writing (again) to you to overtake a problem with the point-stat
> procedure.
> Syntax of the command:
> point_stat meteosalute_2007012100-2007012100_d01.grb obs.ncPointStatConfigMETEOS
> (in attach you can find input files)
> It produce, among the others, the attached continuos statistics (cnt file)
> but as you can see
> the forecast value 72.2 for variable pressure tendency is not the same as
> you can dump from
> the grib file 96.6998....
> I guess there's some problem with the definition of masking region (I tried
> proper(?) definition for mask_grids[] and mask_polys[] and mask_stations)
> but I had the same result.
> Where am I wrong??
> Thanks, Valerio
> 
> 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: read_grib.cc
Type: text/x-c++src
Size: 47465 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/met_help/attachments/20081013/f54c9f9e/read_grib-0002.bin
-------------- next part --------------
A non-text attachment was scrubbed...
Name: read_grib.h
Type: text/x-chdr
Size: 3102 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/met_help/attachments/20081013/f54c9f9e/read_grib-0003.bin