John,<br>Thanks for your code. I used it to filter out some obs(UGRD &gt;= 30.0 m/s) and compare result with WRF-Var verification. The U RMSE is pretty close this time and I also checked these ob/fcst paires in WRF-Var verification and found out that they are treated as missing value. So it looks like that quality control of obs makes big difference between MET and WRF-Var.<br>

<br>Thanks,<br>Ruifang<br><br><div class="gmail_quote">On Fri, Dec 11, 2009 at 10:36 AM, John Halley Gotway <span dir="ltr">&lt;<a href="mailto:johnhg@ucar.edu" target="_blank">johnhg@ucar.edu</a>&gt;</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 and Kathryn,<br>
<br>
Please take a look at the attached file &quot;point_stat.cc&quot;.  Search in there for my initials, &quot;JHG&quot;.  You&#39;ll see that I&#39;ve added a section of code right after we read in the observation values.  In<br>


there, I added a check to skip any observations of UGRD (u-component of wind) that are &gt;= 60.0 m/s.  Feel free to modify this however you&#39;d like - change that value of 60, add other filtering<br>
options... whatever makes sense for your application.<br>
<br>
Also, please let me know how this works out for you.  If this type of simple filtering provides you with the functionality you need, we could consider adding it to the next release of MET.<br>
<br>
So take the attached file, copy it into &quot;METv2.0/src/point_stat/point_stat.cc&quot; and recompile MET.<br>
<br>
Also, if you haven&#39;t already done so, you should grab the latest patches for METv2.0 from:<br>
<div><a href="http://www.dtcenter.org/met/users/support/known_issues/METv2.0/index.php" target="_blank">http://www.dtcenter.org/met/users/support/known_issues/METv2.0/index.php</a><br>
</div>Read and follow the directions at the top of the page.  Please note, if you update METv2.0 with the patches, do so FIRST and then copy over the attached version of point_stat.cc SECOND.<br>
<br>
Hope that does the trick.<br>
<font color="#888888"><br>
John<br>
</font></blockquote></div><br><br clear="all"><br>-- <br>Ruifang <br><br>Mesoscale &amp; Microscale Meteorology Division<br>Phone: 303-497-8938<br>Office:  FL3-3085<br>Email:  <a href="mailto:lir@ucar.edu" target="_blank">lir@ucar.edu</a><br>