[Met_help] question about point_stat

Ruifang Li lir at ucar.edu
Tue Aug 4 15:56:12 MDT 2009


Hi, John,
I run point_stat for 2007081512 fcst/obs and get this error:

Reading records for UGRD/P30.

ERROR: process_grib_codes() -> no records matching GRIB code 33 with level
indicator of P30 found in GRIB file: /ptmp/lir/test/data/2009_t8_15k
m/verify_met2.0/gsi/wpp/2007081512/postprd/wrfprs_d01.000

The ERROR happens for P70, P30, P20, P10. When I remove these levels,
point_stat run well.
So how do you handle single levels inside the code?  It should be better to
use these levels in my verification.

Here is fcst_field in my config:

fcst_field[] =
["UGRD/Z10","UGRD/P1000","UGRD/P925","UGRD/P850","UGRD/P700","UGRD/P500","UGRD/P400","UGRD/P300","UGRD/P250","UGRD/P200","UGRD/P
150","UGRD/P100","UGRD/P70","UGRD/P50","UGRD/P30","UGRD/P20","UGRD/P10" ];
obs_field[]  = [];


Thanks for help,
Ruifang

Mesoscale & Microscale Meteorology Division
Phone: 303-497-8938
Office:  FL3-3085
Email:  lir at ucar.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/met_help/attachments/20090804/57751582/attachment.html 


More information about the Met_help mailing list