John,<br>I get this error when I run point_stat tool:<br><br><span style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">*** Running POINT_STAT on WRF OUTPUT ***</span><br style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">
<span style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">*** UGRD ***</span><br style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);"><span style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">GSL_RNG_TYPE=mt19937</span><br style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">
<span style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">GSL_RNG_SEED=1666342218</span><br style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);"><span style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">Forecast File: /ptmp/lir/test/data/2009_t8_15km/verify_met2.0/gsi/wpp/2007081512/postprd/wrfprs_d01.000</span><br style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">
<span style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">Climatology File: none</span><br style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);"><span style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">Configuration File: ./config/PointStatConfig.UGRD</span><br style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">
<span style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">Observation File: /ptmp/lir/test/data/2009_t8_15km/verify_met2.0/ob_nc/2007081512/<a href="http://ob.nc">ob.nc</a></span><br style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">
<br style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);"><span style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">----------------------------------------</span><br style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">
<br style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);"><span style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">Reading records for UGRD/P1000.</span><br style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">
<span style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">For UGRD/P1000 found 1 forecast levels and 0 climatology levels.</span><br style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">
<br style="background-color: rgb(153, 153, 153);"><br style="background-color: rgb(153, 153, 153);"><span style="background-color: rgb(153, 153, 153); color: rgb(255, 255, 255);">ERROR: is_grid_relative() -> Unsupported grid type value: ^A</span><br>
<br><br>But point_stat works well with TMP. I looked at code read_grib.cc and found out<br><br><span style="color: rgb(153, 153, 153);">// FUTURE WORK: ADD SUPPORT FOR THE MERCATOR PROJECTION</span><br style="color: rgb(153, 153, 153);">
<br style="color: rgb(153, 153, 153);"><span style="color: rgb(153, 153, 153);"> //</span><br style="color: rgb(153, 153, 153);"><span style="color: rgb(153, 153, 153);"> // Mercator Projection Grid</span><br style="color: rgb(153, 153, 153);">
<span style="color: rgb(153, 153, 153);"> //</span><br style="color: rgb(153, 153, 153);"><br style="color: rgb(153, 153, 153);"><span style="color: rgb(153, 153, 153);">else if(r.gds->type == 1) {</span><br style="color: rgb(153, 153, 153);">
<br style="color: rgb(153, 153, 153);"><span style="color: rgb(153, 153, 153);"> cerr << "\n\nERROR: read_gds() -> "</span><br style="color: rgb(153, 153, 153);"><span style="color: rgb(153, 153, 153);"> << "Mercator grid type " << (int) r.gds->type</span><br style="color: rgb(153, 153, 153);">
<span style="color: rgb(153, 153, 153);"> << " not currently supported.\n\n"</span><br style="color: rgb(153, 153, 153);"><span style="color: rgb(153, 153, 153);"> << flush;</span><br style="color: rgb(153, 153, 153);">
<span style="color: rgb(153, 153, 153);"> exit(1);</span><br style="color: rgb(153, 153, 153);"><span style="color: rgb(153, 153, 153);"> }</span><br><br>My wrf forecast use Mercator projection, It looks like that I can not do U, V, WIND using MET2.0, right?<br>
Are there other ways that I can try?<br><br>Thanks,<br>Ruifang<br><br> <br><div class="gmail_quote">On Tue, Jul 14, 2009 at 10:18 AM, 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>
Sorry about that, I forgot the attachment. It's attached now.<br>
<br>
You're right that you could use a polyline to define the masking region for the observations you'd like to retain. But to answer you question, does G181 = CAR.poly? I'd say the answer is NO.<br>
<br>
Look in the file "METv2.0/data/poly/ncep_vx_regions.pdf", at the region marked "CAR". I'm sure that doesn't line up exactly with NCEP Grid 181, although I can't find a plot of Grid 181 right now.<br>
<br>
You can do the masking with a grid or a polyline, whichever you prefer.<br>
<br>
Thanks,<br>
<font color="#888888">John<br>
</font><div><div></div><div class="h5"><br>
Ruifang Li wrote:<br>
> John,<br>
> Thanks for quick reply.<br>
> One more question, If I define these variables in the config file<br>
><br>
> mask_grid=" "<br>
> maks_poly="MET_BASE/data/poly/CAR.poly"<br>
><br>
> Is it as same as<br>
><br>
> mask_grid="G181 "<br>
> maks_poly=" "<br>
><br>
> By the way, I could not see your attached file.<br>
><br>
> Thanks again,<br>
> Ruifang<br>
><br>
><br>
> On Tue, Jul 14, 2009 at 6:37 AM, John Halley Gotway <<a href="mailto:johnhg@rap.ucar.edu">johnhg@rap.ucar.edu</a>>wrote:<br>
><br>
>> Hello,<br>
>><br>
>> You raise a good point. Here's the webpage from which we originally<br>
>> extracted the grid definitions:<br>
>> <a href="http://www.nco.ncep.noaa.gov/pmb/docs/on388/tableb.html" target="_blank">http://www.nco.ncep.noaa.gov/pmb/docs/on388/tableb.html</a><br>
>><br>
>> I see that the information on there was updated on 9/15/2008, and that the<br>
>> changes are listed in red, including the addition of grid #181. This will<br>
>> require some code changes to grab all of these<br>
>> updated grid definitions. In the meantime, I've modified the file<br>
>> "METv2.0/lib/vx_data_grids/grid_defs.f" to add NCEP Grid 181 to the list of<br>
>> definitions. Please try replacing your original version<br>
>> of the file with the attached one, and then recompile MET being sure to do<br>
>> a "make clean" first. Please let me know if you have any troubles using<br>
>> grid #181 after recompiling.<br>
>><br>
>> When I get a chance sometime over the next month, I'll update that file<br>
>> "METv2.0/lib/vx_data_grids/grid_defs.f" with all of the changes and post it<br>
>> to the MET website.<br>
>><br>
>> Thanks for raising this issue!<br>
>><br>
>> John Halley Gotway<br>
>> <a href="mailto:johnhg@ucar.edu">johnhg@ucar.edu</a><br>
>><br>
>> Ruifang Li wrote:<br>
>>> Hi, John,<br>
>>> I have a question about variable mask_grid = "G212" in the default<br>
>> config<br>
>>> file of PB2NC. I tried to use mask_grid = "G181"<br>
>>> for Caribean domain and it says:<br>
>>> ERROR: process_config_file() -> the mask_grid requested "G181" is not<br>
>>> defined.<br>
>>><br>
>>> But G212 does work.<br>
>>><br>
>>> How to use mask_grid = "G181"<br>
>>><br>
>>> Thanks,<br>
>>><br>
>>><br>
>>> ------------------------------------------------------------------------<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>
><br>
><br>
><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>