[Met_help] Zero output

[John Halley Gotway]

Joe,

I see that you're running the Point-Stat tool but are getting zero matched pairs.  There are many reasons why this could be happening, such as...
- Perhaps the observations don't fall into the region you're verifying - maybe you have a sign wrong in the lat/lon location.
- Perhaps the forecast and observation times don't match up.
- Perhaps the observation message type doesn't match what you're trying to verify.

It may be easiest for you to just send me a sample forecast and observation file, and I'll take a look and see what's going on.  Could you post a sample forecast and observation file to our anonymous
ftp site?  Here's how you'd do that:

ftp ftp.rap.ucar.edu
username = anonymous
password = "your email address"
cd incoming/irap/met_help/eastman_data
put "your forecast file"
put "your observation file"
bye

To answer you're second question... the MET stats tools (Grid-Stat, Point-Stat, Wavelet-Stat, and MODE) were designed to run at a single point in time.  So you should run Point-Stat once for each
output time... in your case, every 10 minutes.  To aggregate statistics through time, you'll need to run the STAT-Analysis tool.  Now I see that you're trying to verify over the whole continental US
as well as several individual stations.  When verifying at stations, each individual station is treated as it's own "masking region".  For each run of Point-Stat, you should only be getting 1 matched
pair for each station.  So you really will need to run your output through the STAT-Analysis tool to aggregate through time, which is what I assume you're really interested in.  I realize that this
design is cumbersome for what you're trying to do - but that's how it currently works.

Based on what you appear to be doing, here are some suggestions:
- If you haven't already done so, retrieve and install the latest set of patches for METv2.0 from:
     http://www.dtcenter.org/met/users/support/known_issues/METv2.0/index.php
- Do NOT combine the GRIB files for multiple times into a single file.  MET expects that GRIB output from separate times be in separate files.
- I see that you've set "beg_ds" and "end_ds" to +/- 60 seconds in the config file which will match +/- 1 minute around the forecast valid time.  That looks like a good choice to me for your data.

Let's try to get a single case running for now, and once that's working well, I'll be happy to help you aggregate them through time using STAT-Analysis.

Please let me know if you're able to post those forecast/observation files.  And I'll take a look to see if I can provide better advice.

Thanks,
John Halley Gotway
johnhg at ucar.edu

Joe Eastman wrote:
> Hello,
> 
> I am new to using METv2.0, so this is probably something silly. I run with the command:
> 
> ~/METv2.0-update_bin/point_stat myj_d03_2002-05-05_06:00:00 okmet.nc pointstat -outdir . -v 1000 >&m.out&
> 
> What I end up with is zero output in the txt and stat files. It prints the headers and nothing else. I have converted some Oklahoma mesonet data successfully using the ascii2nc program. I included my station file too. The other question that is not obvious to is how does one do multiple time periods. I have data every 5 minutes. I write the grib WRF output every 10 minutes and would like evaluation over all forecast times in a single file. Is this possible? Do I need to write all WRF output to a single grib file?
> 
> Thanks for your time.
> 
> Take care,
> Joe
> 
> 
> Dr. Joseph L. Eastman
> Senior Atmospheric Scientist
> WindLogics Inc.
> 201 4th St NW
> Grand Rapids, MN
> 55744
> 
> c - 410.279.9702
> w - 651.556.4297
> jeastman at windlogics.com<mailto:jeastman at windlogics.com>
> 
> 
> 
> ________________________________
> CONFIDENTIALITY NOTICE: This e-mail message is for the sole use of the intended recipient(s) and may contain confidential and privileged information. Any unauthorized review, use, disclosure or distribution of any kind is strictly prohibited. If you are not the intended recipient, please contact the sender via reply e-mail and destroy all copies of the original message. Thank you.
> 
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> Met_help mailing list
> Met_help at mailman.ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/met_help