[Met_help] question about point_stat

[Ruifang Li]

John,
You are right. I made a change in WPP configure file and now get levels
output for P70,P30,P20,P10.

Thanks for your help,
Ruifang

On Tue, Aug 4, 2009 at 4:16 PM, John Halley Gotway <johnhg at rap.ucar.edu>wrote:

> 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
> >
>
>
>


-- 
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/20090805/5eaedf41/attachment.html