John<br>I use MET2.0. I used attached file and rebuild MET. Point_stat works well with U,V,WIND now.<br><br>I looked at FUTURE WORK: ADD SUPPORT FOR THE MERCATOR PROJECTION" in METv2.0 and did not find it anywhere. I guess I was totally overflowed yesterday night and went to METv1.0 to find it. Sorry for this to make you confused.<br>
<br>Thanks again for your helpful support,<br>Ruifang<br><br><div class="gmail_quote">On Wed, Jul 15, 2009 at 7:20 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>
It was in METv2.0 that we added support for Mercator projections into MET. And it looks like we missed a change we needed to make in one spot that's casing the error message you're seeing. I'm<br>
surprised we hadn't heard about this yet, but I don't think that the Mercator projection is used very frequently.<br>
<br>
This bug has an easy fix, just adding a few lines to one file. I'm wondering though, what version of MET are you running? That comment "FUTURE WORK: ADD SUPPORT FOR THE MERCATOR PROJECTION"<br>
shouldn't exist anywhere in the METv2.0 source code. It was in METv1.1.<br>
<br>
In order to use Mercator projection, you'll need to use METv2.0. And to fix this bug, you'll need to replace the file "METv2.0/lib/vx_met_util/read_grid.cc" with the version that's attached. So<br>
please try updating that file, rebuilding MET, and rerunning. And please let me know if that fixes your problem. If so, I'll post that fix to the MET website.<br>
<br>
Thanks for finding this issue.<br>
<div><div></div><div class="h5"><br>
Thanks,<br>
John<br>
<br>
Ruifang Li wrote:<br>
> John,<br>
> I get this error when I run point_stat tool:<br>
><br>
> *** Running POINT_STAT on WRF OUTPUT ***<br>
> *** UGRD ***<br>
> GSL_RNG_TYPE=mt19937<br>
> GSL_RNG_SEED=1666342218<br>
> Forecast File:<br>
> /ptmp/lir/test/data/2009_t8_15km/verify_met2.0/gsi/wpp/2007081512/postprd/wrfprs_d01.000<br>
> Climatology File: none<br>
> Configuration File: ./config/PointStatConfig.UGRD<br>
> Observation File:<br>
> /ptmp/lir/test/data/2009_t8_15km/verify_met2.0/ob_nc/2007081512/<a href="http://ob.nc" target="_blank">ob.nc</a><br>
><br>
> ----------------------------------------<br>
><br>
> Reading records for UGRD/P1000.<br>
> For UGRD/P1000 found 1 forecast levels and 0 climatology levels.<br>
><br>
><br>
> ERROR: is_grid_relative() -> Unsupported grid type value: ^A<br>
><br>
><br>
> But point_stat works well with TMP. I looked at code read_grib.cc and found<br>
> out<br>
><br>
> // FUTURE WORK: ADD SUPPORT FOR THE MERCATOR PROJECTION<br>
><br>
> //<br>
> // Mercator Projection Grid<br>
> //<br>
><br>
> else if(r.gds->type == 1) {<br>
><br>
> cerr << "\n\nERROR: read_gds() -> "<br>
> << "Mercator grid type " << (int) r.gds->type<br>
> << " not currently supported.\n\n"<br>
> << flush;<br>
> exit(1);<br>
> }<br>
><br>
> My wrf forecast use Mercator projection, It looks like that I can not do U,<br>
> V, WIND using MET2.0, right?<br>
> Are there other ways that I can try?<br>
><br>
> Thanks,<br>
> Ruifang<br>
><br>
><br>
> On Tue, Jul 14, 2009 at 10:18 AM, John Halley Gotway <<a href="mailto:johnhg@rap.ucar.edu">johnhg@rap.ucar.edu</a>>wrote:<br>
><br>
>> 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<br>
>> the observations you'd like to retain. But to answer you question, does<br>
>> 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<br>
>> marked "CAR". I'm sure that doesn't line up exactly with NCEP Grid 181,<br>
>> 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>
>> John<br>
>><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><br>
>>> 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<br>
>> the<br>
>>>> changes are listed in red, including the addition of grid #181. This<br>
>> 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<br>
>> of<br>
>>>> definitions. Please try replacing your original version<br>
>>>> of the file with the attached one, and then recompile MET being sure to<br>
>> 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<br>
>> 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>
><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>
</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>