[Met_help] question about point_stat

[John Halley Gotway]

Ruifang,

I took a look at your input file on bluefire in:
/ptmp/lir/test/data/2009_t8_15km/verify_met2.0/gsi/wpp/2007081512/postprd/wrfprs_d01.000

And I ran the following wgrib command on it to see the levels for which
you have UGRD records in your file:
wgrib wrfprs_d01.000 | grep UGRD | cut -d':' -f12

And there appear to be no UGRD records for pressure levels 10, 20, 30, and
70 mb in your GRIB file.  That's why you're getting the error.  If you
think you should have output for these levels, take a look at how you have
WPP configured.

Here's the wgrib output from that command above:
hybrid lev 1
hybrid lev 2
50 mb
100 mb
150 mb
200 mb
250 mb
300 mb
350 mb
400 mb
450 mb
500 mb
550 mb
600 mb
650 mb
700 mb
750 mb
775 mb
800 mb
825 mb
850 mb
875 mb
900 mb
925 mb
950 mb
975 mb
1000 mb
10 m above gnd
30-0 mb above gnd

Thanks,
John

> 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
> _______________________________________________
> Met_help mailing list
> Met_help at mailman.ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/met_help
>