John,<br>You are right. I made a change in WPP configure file and now get levels output for P70,P30,P20,P10.<br><br>Thanks for your help,<br>Ruifang<br><br><div class="gmail_quote">On Tue, Aug 4, 2009 at 4:16 PM, John Halley Gotway <span dir="ltr"><<a href="mailto:johnhg@rap.ucar.edu">johnhg@rap.ucar.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Ruifang,<br>
<br>
I took a look at your input file on bluefire in:<br>
/ptmp/lir/test/data/2009_t8_15km/verify_met2.0/gsi/wpp/2007081512/postprd/wrfprs_d01.000<br>
<br>
And I ran the following wgrib command on it to see the levels for which<br>
you have UGRD records in your file:<br>
wgrib wrfprs_d01.000 | grep UGRD | cut -d':' -f12<br>
<br>
And there appear to be no UGRD records for pressure levels 10, 20, 30, and<br>
70 mb in your GRIB file. That's why you're getting the error. If you<br>
think you should have output for these levels, take a look at how you have<br>
WPP configured.<br>
<br>
Here's the wgrib output from that command above:<br>
hybrid lev 1<br>
hybrid lev 2<br>
50 mb<br>
100 mb<br>
150 mb<br>
200 mb<br>
250 mb<br>
300 mb<br>
350 mb<br>
400 mb<br>
450 mb<br>
500 mb<br>
550 mb<br>
600 mb<br>
650 mb<br>
700 mb<br>
750 mb<br>
775 mb<br>
800 mb<br>
825 mb<br>
850 mb<br>
875 mb<br>
900 mb<br>
925 mb<br>
950 mb<br>
975 mb<br>
1000 mb<br>
10 m above gnd<br>
30-0 mb above gnd<br>
<br>
Thanks,<br>
John<br>
<div><div></div><div class="h5"><br>
> Hi, John,<br>
> I run point_stat for 2007081512 fcst/obs and get this error:<br>
><br>
> Reading records for UGRD/P30.<br>
><br>
> ERROR: process_grib_codes() -> no records matching GRIB code 33 with level<br>
> indicator of P30 found in GRIB file: /ptmp/lir/test/data/2009_t8_15k<br>
> m/verify_met2.0/gsi/wpp/2007081512/postprd/wrfprs_d01.000<br>
><br>
> The ERROR happens for P70, P30, P20, P10. When I remove these levels,<br>
> point_stat run well.<br>
> So how do you handle single levels inside the code? It should be better<br>
> to<br>
> use these levels in my verification.<br>
><br>
> Here is fcst_field in my config:<br>
><br>
> fcst_field[] =<br>
> ["UGRD/Z10","UGRD/P1000","UGRD/P925","UGRD/P850","UGRD/P700","UGRD/P500","UGRD/P400","UGRD/P300","UGRD/P250","UGRD/P200","UGRD/P<br>
> 150","UGRD/P100","UGRD/P70","UGRD/P50","UGRD/P30","UGRD/P20","UGRD/P10" ];<br>
> obs_field[] = [];<br>
><br>
><br>
> Thanks for help,<br>
> Ruifang<br>
><br>
> Mesoscale & Microscale Meteorology Division<br>
> Phone: 303-497-8938<br>
> Office: FL3-3085<br>
> Email: <a href="mailto:lir@ucar.edu">lir@ucar.edu</a><br>
</div></div>> _______________________________________________<br>
> Met_help mailing list<br>
> <a href="mailto:Met_help@mailman.ucar.edu">Met_help@mailman.ucar.edu</a><br>
> <a href="http://mailman.ucar.edu/mailman/listinfo/met_help" target="_blank">http://mailman.ucar.edu/mailman/listinfo/met_help</a><br>
><br>
<br>
<br>
</blockquote></div><br><br clear="all"><br>-- <br>Ruifang <br><br>Mesoscale & Microscale Meteorology Division<br>Phone: 303-497-8938<br>Office: FL3-3085<br>Email: <a href="mailto:lir@ucar.edu">lir@ucar.edu</a><br>