[Met_help] PB2NC config file
John Halley Gotway
johnhg at rap.ucar.edu
Wed Jul 15 07:20:30 MDT 2009
Ruifang,
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
surprised we hadn't heard about this yet, but I don't think that the Mercator projection is used very frequently.
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"
shouldn't exist anywhere in the METv2.0 source code. It was in METv1.1.
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
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.
Thanks for finding this issue.
Thanks,
John
Ruifang Li wrote:
> John,
> I get this error when I run point_stat tool:
>
> *** Running POINT_STAT on WRF OUTPUT ***
> *** UGRD ***
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=1666342218
> Forecast File:
> /ptmp/lir/test/data/2009_t8_15km/verify_met2.0/gsi/wpp/2007081512/postprd/wrfprs_d01.000
> Climatology File: none
> Configuration File: ./config/PointStatConfig.UGRD
> Observation File:
> /ptmp/lir/test/data/2009_t8_15km/verify_met2.0/ob_nc/2007081512/ob.nc
>
> ----------------------------------------
>
> Reading records for UGRD/P1000.
> For UGRD/P1000 found 1 forecast levels and 0 climatology levels.
>
>
> ERROR: is_grid_relative() -> Unsupported grid type value: ^A
>
>
> But point_stat works well with TMP. I looked at code read_grib.cc and found
> out
>
> // FUTURE WORK: ADD SUPPORT FOR THE MERCATOR PROJECTION
>
> //
> // Mercator Projection Grid
> //
>
> else if(r.gds->type == 1) {
>
> cerr << "\n\nERROR: read_gds() -> "
> << "Mercator grid type " << (int) r.gds->type
> << " not currently supported.\n\n"
> << flush;
> exit(1);
> }
>
> My wrf forecast use Mercator projection, It looks like that I can not do U,
> V, WIND using MET2.0, right?
> Are there other ways that I can try?
>
> Thanks,
> Ruifang
>
>
> On Tue, Jul 14, 2009 at 10:18 AM, John Halley Gotway <johnhg at rap.ucar.edu>wrote:
>
>> Ruifang,
>>
>> Sorry about that, I forgot the attachment. It's attached now.
>>
>> 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.
>>
>> 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.
>>
>> You can do the masking with a grid or a polyline, whichever you prefer.
>>
>> Thanks,
>> John
>>
>> Ruifang Li wrote:
>>> John,
>>> Thanks for quick reply.
>>> One more question, If I define these variables in the config file
>>>
>>> mask_grid=" "
>>> maks_poly="MET_BASE/data/poly/CAR.poly"
>>>
>>> Is it as same as
>>>
>>> mask_grid="G181 "
>>> maks_poly=" "
>>>
>>> By the way, I could not see your attached file.
>>>
>>> Thanks again,
>>> Ruifang
>>>
>>>
>>> On Tue, Jul 14, 2009 at 6:37 AM, John Halley Gotway <johnhg at rap.ucar.edu
>>> wrote:
>>>
>>>> Hello,
>>>>
>>>> You raise a good point. Here's the webpage from which we originally
>>>> extracted the grid definitions:
>>>> http://www.nco.ncep.noaa.gov/pmb/docs/on388/tableb.html
>>>>
>>>> I see that the information on there was updated on 9/15/2008, and that
>> the
>>>> changes are listed in red, including the addition of grid #181. This
>> will
>>>> require some code changes to grab all of these
>>>> updated grid definitions. In the meantime, I've modified the file
>>>> "METv2.0/lib/vx_data_grids/grid_defs.f" to add NCEP Grid 181 to the list
>> of
>>>> definitions. Please try replacing your original version
>>>> of the file with the attached one, and then recompile MET being sure to
>> do
>>>> a "make clean" first. Please let me know if you have any troubles using
>>>> grid #181 after recompiling.
>>>>
>>>> When I get a chance sometime over the next month, I'll update that file
>>>> "METv2.0/lib/vx_data_grids/grid_defs.f" with all of the changes and post
>> it
>>>> to the MET website.
>>>>
>>>> Thanks for raising this issue!
>>>>
>>>> John Halley Gotway
>>>> johnhg at ucar.edu
>>>>
>>>> Ruifang Li wrote:
>>>>> Hi, John,
>>>>> I have a question about variable mask_grid = "G212" in the default
>>>> config
>>>>> file of PB2NC. I tried to use mask_grid = "G181"
>>>>> for Caribean domain and it says:
>>>>> ERROR: process_config_file() -> the mask_grid requested "G181" is not
>>>>> defined.
>>>>>
>>>>> But G212 does work.
>>>>>
>>>>> How to use mask_grid = "G181"
>>>>>
>>>>> Thanks,
>>>>>
>>>>>
>>>>>
>> ------------------------------------------------------------------------
>>>>> _______________________________________________
>>>>> Met_help mailing list
>>>>> Met_help at mailman.ucar.edu
>>>>> http://mailman.ucar.edu/mailman/listinfo/met_help
>>>
>>>
>
>
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> Met_help mailing list
> Met_help at mailman.ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/met_help
-------------- next part --------------
A non-text attachment was scrubbed...
Name: read_grib.cc
Type: text/x-c++src
Size: 50125 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/met_help/attachments/20090715/b9fc7aed/attachment-0001.bin
More information about the Met_help
mailing list