John,<br>Thank you very much for your suggestion and help. I have also been thinking to turn off MPR output in my next verifaction. <br><br>Thanks again,<br>Ruifang<br><br><div class="gmail_quote">On Fri, Jul 31, 2009 at 5:18 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>
Thanks for sending those files. Let me start by making a suggestion that you turn OFF the matched pair output line (MPR) from Point-Stat. Do this by changing the value for it in the "output_flag"<br>
parameter. You're generating files that are unnecessarily large. The MPR line type is only meant to be used in debugging or testing mode. You probably shouldn't use it for large runs!<br>
<br>
I did look at the MPR line type in this case. I looked at the matched pairs for TMP, Z02, ADPSFC, UW_MEAN, and 1. For the FULL domain, I see 2742 pairs and for the T8 domain, I see 2744 pairs.<br>
Those numbers are very close but they're counter intuitive. It seems like they should be switched. You should have more matched pairs in the FULL domain than in the T8 domain. I looked at the<br>
lat/lon's of the observation points for the 2 that are in T8 but not in full. They both have longitude = -104.8, which I'm guessing is right on the boundary of your domain.<br>
<br>
I have no idea why those 2 are included in T8 but not in FULL. I went hunting through the code to see if I could figure out why, but I didn't turn anything up. I think it might be a rounding issue<br>
along the boundary when converting back and forth between lat/lon and x/y.<br>
<br>
Thanks,<br>
<font color="#888888">John<br>
</font><div><div></div><div class="h5"><br>
Ruifang Li wrote:<br>
> Hi, John,<br>
> I tested only 2007081512 for you to look at this time.<br>
><br>
> You can look at all my data on<br>
> /ptmp/lir/test/data/2009_t8_15km/verify_met2.0/gsi/point_stat<br>
><br>
> Thanks for help,<br>
> Ruifang<br>
><br>
><br>
> On Thu, Jul 30, 2009 at 10:58 AM, Ruifang Li <<a href="mailto:lir@ucar.edu">lir@ucar.edu</a>> wrote:<br>
><br>
>> Hi, John,<br>
>> I am a little bit confused about mask_grid and mask_poly in point_stat<br>
>> config file.<br>
>><br>
>> In pb2nc, I defined mask_grid=" ", mask_poly="MET_BASE/data/poly/T8.poly",<br>
>> so I got obs for T8 domain since my WRF fcst are on this domain.<br>
>> In point_stat, I did several test about mask_grid and mask_poly , here is<br>
>> the result<br>
>><br>
>> * mask_grid[] mask_poly[] domain<br>
>> in .stat .stat file size*<br>
>><br>
>> FULL T8<br>
>> FULL,T8 7,017,368<br>
>> T8<br>
>> T8 3,010,624<br>
>> -- FULL<br>
>> FULL 3,507,520<br>
>><br>
>> My question are:<br>
>><br>
>> - it's better to define both of mask_poly and mask_grid or only one of<br>
>> them is enough?<br>
>> - If I define both of mask_poly and mask_grid, the FULL and T8 records<br>
>> show up in .stat file, is there any difference between FULL and T8 records?<br>
>><br>
>><br>
>> Thanks for your 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>
>><br>
><br>
><br>
><br>
><br>
</div></div><div><div></div><div class="h5">> ------------------------------------------------------------------------<br>
><br>
> _______________________________________________<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>
</div></div></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>