[Met_help] MET toolkit
John Halley Gotway
johnhg at rap.ucar.edu
Tue Mar 16 16:14:54 MDT 2010
Jackie,
I'm forwarding your message to the "met_help at ucar.edu" email address.
Unfortunately, I'm out of the office at all week. I'm with my son's who's
recovering from a burst appendix. He's recovering nicely, but we'll be
here with him all week.
Perhaps one of the other MET team members could help you figure out why
you're not getting verification output for 10-meter winds. I'd really
suggest sending us the files you're using to run Point-Stat:
(1) Forecast GRIB file
(2) NetCDF point observation file
(3) Point-Stat configuration file
(4) Any files you're using to define the masking verification regions
You can post them to our anonymous ftp site:
ftp ftp.rap.ucar.edu
username=anonymous
password="your email address"
cd incoming/irap/met_help/miller_data
put "your files"
bye
I suspect that there's some issue preventing you from finding any matched
pairs for 10-meter winds. When you run Point-Stat with the "-v 2" option,
it should tell you how many matched pairs it found for 10-meter winds. I
suspect that's 0.
I should be available to help out sometime next week.
Thanks,
John
> Hi John,
>
>
>
> I just wanted to touch base again about the grid_stat output. I went
> ahead and interrogated the input grib file like you recommended using
> the following command.
>
> wgrib in_file.grib | egrep "33|34|11"
>
> And here is a screenshot of the first few lines of the output:
>
>
>
>
>
> As you suspected, the kpds7 value was "2" for temp and "10" for U and V.
> So I went back and changed the config file to:
>
>
>
> fcst_field[] = [ "7/Z0","7/P850","7/P500","7/P250","33/Z10",
>
> "34/Z10","33/P850","34/P850","33/P500",
>
> "34/P500","33/P250","34/P250","11/Z0",
>
> "11/P850","11/P500","11/P250" ];
>
> And the output file that was created did not return results for the UGRD
> and VGRD at 10 meters. Am I missing something? I tried with Z0 and got
> the same result. Just to see what would happen, I tried the following:
>
>
>
> fcst_field[] = [ "7/Z0","7/P850","7/P500","7/P250","33/P1000",
>
> "34/P1000","33/P850","34/P850","33/P500",
>
> "34/P500","33/P250","34/P250","11/Z0",
>
> "11/P850","11/P500","11/P250" ];
>
>
>
> And the output file did include results for the UGRD and VGRD at P1000.
> Is this the best approach?
>
>
>
>
>
> Jackie
>
>
>
> Jackie Miller
> Meteorologist
>
> Vaisala Inc.
>
> Phone +1 303 885 3251
> Mobile +1 303 885 3251
> www.vaisala.com
> <https://exchange.ou.edu/owa/redir.aspx?C=b302bd7f7cb24f099b2eb70ba1b9bf
> 5f&URL=https%3a%2f%2fwebmail.vaisala.com%2fexchweb%2fbin%2fredir.asp%3fU
> RL%3dhttp%3a%2f%2fwww.vaisala.com>
>
>
>
> This electronic message contains a communication which is strictly
> confidential and intended solely for the addressee(s). Any
> non-addressee is prohibited from reading, disseminating, distributing,
> or copying the communication contained herein. If you are in possession
> of the communication in error, please immediately notify the sender via
> an electronic response and destroy the original communication. Thank
> you.
>
>
>
> U.S. Export Restrictions & Disclaimer: Export of any technical
> information contained in this email and/or its attachments is subject to
> the export control laws and regulations of the U.S. Government and may
> require a valid license or written approval prior to export.
>
>
>
>
More information about the Met_help
mailing list