[Met_help] [rt.rap.ucar.edu #74090] History for masking region by latitude
John Halley Gotway via RT
met_help at ucar.edu
Wed Nov 25 08:14:34 MST 2015
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hello,
I am trying to compute the rmse and acc using only the regions 20-80
degrees N and 20-80 degrees S. Is there a simple way to do this with MET
5.1? I see there are various .poly files that make it easy to mask these
regions but I do not think any of them are the regions that I want.
Best,
--
Shannon Rees
Programming Scientist
Engility
Geophysical Fluid Dynamics Lab
201 Forrestal Rd Princeton, NJ
609-452-5384
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: masking region by latitude
From: Shannon Rees - NOAA Affiliate
Time: Thu Nov 12 13:36:43 2015
Here is an update:
I created my own simple .poly file:
20_90N
20.00000 0.00000
20.00000 180.00000
20.00000 359.90000
89.00000 0.00000
89.00000 180.00000
89.00000 359.90000
I got this error when trying to run gen_vx_mask:
DEBUG 1: Input File:
/home/Shannon.McElhinney/met-
5.1/data/sample_obs/ST2ml/ST2ml2005080700.Grb_G212
DEBUG 1: Mask File: test.poly
DEBUG 2: Parsed Data Grid: Lambert Conformal (185 x 129)
ERROR :
ERROR : MaskPoly::load() -> read error in mask poly file "test.poly"
ERROR :
On Thu, Nov 12, 2015 at 10:43 AM, met_help at ucar.edu via RT <
met_help at ucar.edu> wrote:
> Greetings,
>
> This message has been automatically generated in response to the
> creation of a trouble ticket regarding:
> "masking region by latitude",
> a summary of which appears below.
>
> There is no need to reply to this message right now. Your ticket
has been
> assigned an ID of [rt.rap.ucar.edu #74090].
>
> Please include the string:
>
> [rt.rap.ucar.edu #74090]
>
> in the subject line of all future correspondence about this issue.
To do
> so,
> you may reply to this message.
>
> Thank you,
> met_help at ucar.edu
>
>
-------------------------------------------------------------------------
> Hello,
>
> I am trying to compute the rmse and acc using only the regions 20-80
> degrees N and 20-80 degrees S. Is there a simple way to do this
with MET
> 5.1? I see there are various .poly files that make it easy to mask
these
> regions but I do not think any of them are the regions that I want.
>
> Best,
> --
> Shannon Rees
> Programming Scientist
> Engility
> Geophysical Fluid Dynamics Lab
> 201 Forrestal Rd Princeton, NJ
> 609-452-5384
>
>
--
Shannon Rees
Programming Scientist
Engility
Geophysical Fluid Dynamics Lab
201 Forrestal Rd Princeton, NJ
609-452-5384
------------------------------------------------
Subject: masking region by latitude
From: John Halley Gotway
Time: Thu Nov 12 14:28:51 2015
Hi Shannon,
I am out on travel today and can't get VPN running to access my
machine
remotely.
I see that you're trying to generate a mask of all grid points that
are
20-80 degrees north and a second mask for 20-80 degrees south.
Do you actually want those as a single mask or separate masks? I'll
assume
separate for now.
Here's the mask I'd try for the North Hemisphere (but I haven't tested
with
this yet)...
20_80N
20 -180
80 -180
80 180
20 180
For the Southern Hemisphere, you'd just switch the 20 and 80 latitude
to
-20 and -80.
Another cool option is this... suppose you want all grid points other
than
-20 to 20 degrees. You could use the following mask and then use the
"-complement" option for gen_vx_mask. This will define the -20 to 20
degree band, but using "-complement" will essentially make it reverse-
video:
20N_20S
-20 -180
20 -180
20 180
-20 180
Please let me know how it goes.
Thanks,
John
On Thu, Nov 12, 2015 at 2:36 PM, Shannon Rees - NOAA Affiliate via RT
<
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74090 >
>
> Here is an update:
>
> I created my own simple .poly file:
> 20_90N
> 20.00000 0.00000
> 20.00000 180.00000
> 20.00000 359.90000
> 89.00000 0.00000
> 89.00000 180.00000
> 89.00000 359.90000
>
> I got this error when trying to run gen_vx_mask:
> DEBUG 1: Input File:
>
> /home/Shannon.McElhinney/met-
5.1/data/sample_obs/ST2ml/ST2ml2005080700.Grb_G212
> DEBUG 1: Mask File: test.poly
> DEBUG 2: Parsed Data Grid: Lambert Conformal (185 x 129)
> ERROR :
> ERROR : MaskPoly::load() -> read error in mask poly file
"test.poly"
> ERROR :
>
>
>
> On Thu, Nov 12, 2015 at 10:43 AM, met_help at ucar.edu via RT <
> met_help at ucar.edu> wrote:
>
> > Greetings,
> >
> > This message has been automatically generated in response to the
> > creation of a trouble ticket regarding:
> > "masking region by latitude",
> > a summary of which appears below.
> >
> > There is no need to reply to this message right now. Your ticket
has
> been
> > assigned an ID of [rt.rap.ucar.edu #74090].
> >
> > Please include the string:
> >
> > [rt.rap.ucar.edu #74090]
> >
> > in the subject line of all future correspondence about this issue.
To do
> > so,
> > you may reply to this message.
> >
> > Thank you,
> > met_help at ucar.edu
> >
> >
-------------------------------------------------------------------------
> > Hello,
> >
> > I am trying to compute the rmse and acc using only the regions 20-
80
> > degrees N and 20-80 degrees S. Is there a simple way to do this
with MET
> > 5.1? I see there are various .poly files that make it easy to
mask these
> > regions but I do not think any of them are the regions that I
want.
> >
> > Best,
> > --
> > Shannon Rees
> > Programming Scientist
> > Engility
> > Geophysical Fluid Dynamics Lab
> > 201 Forrestal Rd Princeton, NJ
> > 609-452-5384
> >
> >
>
>
> --
> Shannon Rees
> Programming Scientist
> Engility
> Geophysical Fluid Dynamics Lab
> 201 Forrestal Rd Princeton, NJ
> 609-452-5384
>
>
------------------------------------------------
Subject: masking region by latitude
From: Shannon Rees - NOAA Affiliate
Time: Thu Nov 12 16:46:47 2015
Hi John,
Sorry to disturb your travels. I do want two separate masks for the
northern and southern hemispheres. I tried your suggestion and got an
error.
met-5.1/bin/gen_vx_mask \
met-5.1/data/sample_obs/ST2ml/ST2ml2005080700.Grb_G212 \
20_80N.poly \
out/gen_vx_mask/20_80N_mask.nc
ERROR : MaskPoly::load() -> read error in mask poly file
"20_80N.poly"
I get the same error using different input files as well. Do I need
to use
a special option?
Thanks
On Thu, Nov 12, 2015 at 4:28 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Hi Shannon,
>
> I am out on travel today and can't get VPN running to access my
machine
> remotely.
>
> I see that you're trying to generate a mask of all grid points that
are
> 20-80 degrees north and a second mask for 20-80 degrees south.
>
> Do you actually want those as a single mask or separate masks? I'll
assume
> separate for now.
>
> Here's the mask I'd try for the North Hemisphere (but I haven't
tested with
> this yet)...
>
> 20_80N
> 20 -180
> 80 -180
> 80 180
> 20 180
>
> For the Southern Hemisphere, you'd just switch the 20 and 80
latitude to
> -20 and -80.
>
> Another cool option is this... suppose you want all grid points
other than
> -20 to 20 degrees. You could use the following mask and then use
the
> "-complement" option for gen_vx_mask. This will define the -20 to
20
> degree band, but using "-complement" will essentially make it
> reverse-video:
>
> 20N_20S
> -20 -180
> 20 -180
> 20 180
> -20 180
>
> Please let me know how it goes.
>
> Thanks,
> John
>
> On Thu, Nov 12, 2015 at 2:36 PM, Shannon Rees - NOAA Affiliate via
RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74090 >
> >
> > Here is an update:
> >
> > I created my own simple .poly file:
> > 20_90N
> > 20.00000 0.00000
> > 20.00000 180.00000
> > 20.00000 359.90000
> > 89.00000 0.00000
> > 89.00000 180.00000
> > 89.00000 359.90000
> >
> > I got this error when trying to run gen_vx_mask:
> > DEBUG 1: Input File:
> >
> >
> /home/Shannon.McElhinney/met-
5.1/data/sample_obs/ST2ml/ST2ml2005080700.Grb_G212
> > DEBUG 1: Mask File: test.poly
> > DEBUG 2: Parsed Data Grid: Lambert Conformal (185 x 129)
> > ERROR :
> > ERROR : MaskPoly::load() -> read error in mask poly file
"test.poly"
> > ERROR :
> >
> >
> >
> > On Thu, Nov 12, 2015 at 10:43 AM, met_help at ucar.edu via RT <
> > met_help at ucar.edu> wrote:
> >
> > > Greetings,
> > >
> > > This message has been automatically generated in response to the
> > > creation of a trouble ticket regarding:
> > > "masking region by latitude",
> > > a summary of which appears below.
> > >
> > > There is no need to reply to this message right now. Your
ticket has
> > been
> > > assigned an ID of [rt.rap.ucar.edu #74090].
> > >
> > > Please include the string:
> > >
> > > [rt.rap.ucar.edu #74090]
> > >
> > > in the subject line of all future correspondence about this
issue. To
> do
> > > so,
> > > you may reply to this message.
> > >
> > > Thank you,
> > > met_help at ucar.edu
> > >
> > >
>
-------------------------------------------------------------------------
> > > Hello,
> > >
> > > I am trying to compute the rmse and acc using only the regions
20-80
> > > degrees N and 20-80 degrees S. Is there a simple way to do this
with
> MET
> > > 5.1? I see there are various .poly files that make it easy to
mask
> these
> > > regions but I do not think any of them are the regions that I
want.
> > >
> > > Best,
> > > --
> > > Shannon Rees
> > > Programming Scientist
> > > Engility
> > > Geophysical Fluid Dynamics Lab
> > > 201 Forrestal Rd Princeton, NJ
> > > 609-452-5384
> > >
> > >
> >
> >
> > --
> > Shannon Rees
> > Programming Scientist
> > Engility
> > Geophysical Fluid Dynamics Lab
> > 201 Forrestal Rd Princeton, NJ
> > 609-452-5384
> >
> >
>
>
--
Shannon Rees
Programming Scientist
Engility
Geophysical Fluid Dynamics Lab
201 Forrestal Rd Princeton, NJ
609-452-5384
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #74090] AutoReply: masking region by latitude
From: John Halley Gotway
Time: Thu Nov 12 18:44:35 2015
Shannon,
I don't see any obvious problems in the command line. I'll replicate
the behavior tomorrow and let you know what I find.
Thanks
John
Sent from my iPhone
On Nov 12, 2015, at 5:46 PM, "Shannon Rees - NOAA Affiliate via RT"
<met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74090 >
>
> Hi John,
>
> Sorry to disturb your travels. I do want two separate masks for the
> northern and southern hemispheres. I tried your suggestion and got
an
> error.
>
> met-5.1/bin/gen_vx_mask \
> met-5.1/data/sample_obs/ST2ml/ST2ml2005080700.Grb_G212 \
> 20_80N.poly \
> out/gen_vx_mask/20_80N_mask.nc
>
> ERROR : MaskPoly::load() -> read error in mask poly file
"20_80N.poly"
>
> I get the same error using different input files as well. Do I need
to use
> a special option?
>
> Thanks
>
>
> On Thu, Nov 12, 2015 at 4:28 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
>> Hi Shannon,
>>
>> I am out on travel today and can't get VPN running to access my
machine
>> remotely.
>>
>> I see that you're trying to generate a mask of all grid points that
are
>> 20-80 degrees north and a second mask for 20-80 degrees south.
>>
>> Do you actually want those as a single mask or separate masks?
I'll assume
>> separate for now.
>>
>> Here's the mask I'd try for the North Hemisphere (but I haven't
tested with
>> this yet)...
>>
>> 20_80N
>> 20 -180
>> 80 -180
>> 80 180
>> 20 180
>>
>> For the Southern Hemisphere, you'd just switch the 20 and 80
latitude to
>> -20 and -80.
>>
>> Another cool option is this... suppose you want all grid points
other than
>> -20 to 20 degrees. You could use the following mask and then use
the
>> "-complement" option for gen_vx_mask. This will define the -20 to
20
>> degree band, but using "-complement" will essentially make it
>> reverse-video:
>>
>> 20N_20S
>> -20 -180
>> 20 -180
>> 20 180
>> -20 180
>>
>> Please let me know how it goes.
>>
>> Thanks,
>> John
>>
>> On Thu, Nov 12, 2015 at 2:36 PM, Shannon Rees - NOAA Affiliate via
RT <
>> met_help at ucar.edu> wrote:
>>
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74090 >
>>>
>>> Here is an update:
>>>
>>> I created my own simple .poly file:
>>> 20_90N
>>> 20.00000 0.00000
>>> 20.00000 180.00000
>>> 20.00000 359.90000
>>> 89.00000 0.00000
>>> 89.00000 180.00000
>>> 89.00000 359.90000
>>>
>>> I got this error when trying to run gen_vx_mask:
>>> DEBUG 1: Input File:
>>>
>>>
>> /home/Shannon.McElhinney/met-
5.1/data/sample_obs/ST2ml/ST2ml2005080700.Grb_G212
>>> DEBUG 1: Mask File: test.poly
>>> DEBUG 2: Parsed Data Grid: Lambert Conformal (185 x 129)
>>> ERROR :
>>> ERROR : MaskPoly::load() -> read error in mask poly file
"test.poly"
>>> ERROR :
>>>
>>>
>>>
>>> On Thu, Nov 12, 2015 at 10:43 AM, met_help at ucar.edu via RT <
>>> met_help at ucar.edu> wrote:
>>>
>>>> Greetings,
>>>>
>>>> This message has been automatically generated in response to the
>>>> creation of a trouble ticket regarding:
>>>> "masking region by latitude",
>>>> a summary of which appears below.
>>>>
>>>> There is no need to reply to this message right now. Your ticket
has
>>> been
>>>> assigned an ID of [rt.rap.ucar.edu #74090].
>>>>
>>>> Please include the string:
>>>>
>>>> [rt.rap.ucar.edu #74090]
>>>>
>>>> in the subject line of all future correspondence about this
issue. To
>> do
>>>> so,
>>>> you may reply to this message.
>>>>
>>>> Thank you,
>>>> met_help at ucar.edu
>>>>
>>>>
>>
-------------------------------------------------------------------------
>>>> Hello,
>>>>
>>>> I am trying to compute the rmse and acc using only the regions
20-80
>>>> degrees N and 20-80 degrees S. Is there a simple way to do this
with
>> MET
>>>> 5.1? I see there are various .poly files that make it easy to
mask
>> these
>>>> regions but I do not think any of them are the regions that I
want.
>>>>
>>>> Best,
>>>> --
>>>> Shannon Rees
>>>> Programming Scientist
>>>> Engility
>>>> Geophysical Fluid Dynamics Lab
>>>> 201 Forrestal Rd Princeton, NJ
>>>> 609-452-5384
>>>>
>>>>
>>>
>>>
>>> --
>>> Shannon Rees
>>> Programming Scientist
>>> Engility
>>> Geophysical Fluid Dynamics Lab
>>> 201 Forrestal Rd Princeton, NJ
>>> 609-452-5384
>>>
>>>
>>
>>
>
>
> --
> Shannon Rees
> Programming Scientist
> Engility
> Geophysical Fluid Dynamics Lab
> 201 Forrestal Rd Princeton, NJ
> 609-452-5384
>
------------------------------------------------
Subject: masking region by latitude
From: John Halley Gotway
Time: Fri Nov 13 14:36:11 2015
Shannon,
In my first test replicating your issue, the gen_vx_mask ran fine
without
error.
For my second test, I added extra empty blank lines at the end of the
ASCII
polyline file, and that did the trick. I was able to reproduce your
error:
ERROR : MaskPoly::load() -> read error in mask poly file
"20_80N.poly"
For now, please just delete all trailing empty lines from your
polyline
file. Hopefully that'll get it to run.
And I'll look into updating the MET code to handle this situation more
gracefully.
Thanks,
John
On Thu, Nov 12, 2015 at 6:44 PM, John NCAR <johnhg at ucar.edu> wrote:
> Shannon,
>
> I don't see any obvious problems in the command line. I'll
replicate the
> behavior tomorrow and let you know what I find.
>
> Thanks
> John
>
> Sent from my iPhone
>
> On Nov 12, 2015, at 5:46 PM, "Shannon Rees - NOAA Affiliate via RT"
<
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74090 >
> >
> > Hi John,
> >
> > Sorry to disturb your travels. I do want two separate masks for
the
> > northern and southern hemispheres. I tried your suggestion and
got an
> > error.
> >
> > met-5.1/bin/gen_vx_mask \
> > met-5.1/data/sample_obs/ST2ml/ST2ml2005080700.Grb_G212 \
> > 20_80N.poly \
> > out/gen_vx_mask/20_80N_mask.nc
> >
> > ERROR : MaskPoly::load() -> read error in mask poly file
"20_80N.poly"
> >
> > I get the same error using different input files as well. Do I
need to
> use
> > a special option?
> >
> > Thanks
> >
> >
> > On Thu, Nov 12, 2015 at 4:28 PM, John Halley Gotway via RT <
> > met_help at ucar.edu> wrote:
> >
> >> Hi Shannon,
> >>
> >> I am out on travel today and can't get VPN running to access my
machine
> >> remotely.
> >>
> >> I see that you're trying to generate a mask of all grid points
that are
> >> 20-80 degrees north and a second mask for 20-80 degrees south.
> >>
> >> Do you actually want those as a single mask or separate masks?
I'll
> assume
> >> separate for now.
> >>
> >> Here's the mask I'd try for the North Hemisphere (but I haven't
tested
> with
> >> this yet)...
> >>
> >> 20_80N
> >> 20 -180
> >> 80 -180
> >> 80 180
> >> 20 180
> >>
> >> For the Southern Hemisphere, you'd just switch the 20 and 80
latitude to
> >> -20 and -80.
> >>
> >> Another cool option is this... suppose you want all grid points
other
> than
> >> -20 to 20 degrees. You could use the following mask and then use
the
> >> "-complement" option for gen_vx_mask. This will define the -20
to 20
> >> degree band, but using "-complement" will essentially make it
> >> reverse-video:
> >>
> >> 20N_20S
> >> -20 -180
> >> 20 -180
> >> 20 180
> >> -20 180
> >>
> >> Please let me know how it goes.
> >>
> >> Thanks,
> >> John
> >>
> >> On Thu, Nov 12, 2015 at 2:36 PM, Shannon Rees - NOAA Affiliate
via RT <
> >> met_help at ucar.edu> wrote:
> >>
> >>>
> >>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74090 >
> >>>
> >>> Here is an update:
> >>>
> >>> I created my own simple .poly file:
> >>> 20_90N
> >>> 20.00000 0.00000
> >>> 20.00000 180.00000
> >>> 20.00000 359.90000
> >>> 89.00000 0.00000
> >>> 89.00000 180.00000
> >>> 89.00000 359.90000
> >>>
> >>> I got this error when trying to run gen_vx_mask:
> >>> DEBUG 1: Input File:
> >>>
> >>>
> >>
> /home/Shannon.McElhinney/met-
5.1/data/sample_obs/ST2ml/ST2ml2005080700.Grb_G212
> >>> DEBUG 1: Mask File: test.poly
> >>> DEBUG 2: Parsed Data Grid: Lambert Conformal (185 x 129)
> >>> ERROR :
> >>> ERROR : MaskPoly::load() -> read error in mask poly file
"test.poly"
> >>> ERROR :
> >>>
> >>>
> >>>
> >>> On Thu, Nov 12, 2015 at 10:43 AM, met_help at ucar.edu via RT <
> >>> met_help at ucar.edu> wrote:
> >>>
> >>>> Greetings,
> >>>>
> >>>> This message has been automatically generated in response to
the
> >>>> creation of a trouble ticket regarding:
> >>>> "masking region by latitude",
> >>>> a summary of which appears below.
> >>>>
> >>>> There is no need to reply to this message right now. Your
ticket has
> >>> been
> >>>> assigned an ID of [rt.rap.ucar.edu #74090].
> >>>>
> >>>> Please include the string:
> >>>>
> >>>> [rt.rap.ucar.edu #74090]
> >>>>
> >>>> in the subject line of all future correspondence about this
issue. To
> >> do
> >>>> so,
> >>>> you may reply to this message.
> >>>>
> >>>> Thank you,
> >>>> met_help at ucar.edu
> >>>>
> >>>>
> >>
>
-------------------------------------------------------------------------
> >>>> Hello,
> >>>>
> >>>> I am trying to compute the rmse and acc using only the regions
20-80
> >>>> degrees N and 20-80 degrees S. Is there a simple way to do
this with
> >> MET
> >>>> 5.1? I see there are various .poly files that make it easy to
mask
> >> these
> >>>> regions but I do not think any of them are the regions that I
want.
> >>>>
> >>>> Best,
> >>>> --
> >>>> Shannon Rees
> >>>> Programming Scientist
> >>>> Engility
> >>>> Geophysical Fluid Dynamics Lab
> >>>> 201 Forrestal Rd Princeton, NJ
> >>>> 609-452-5384
> >>>>
> >>>>
> >>>
> >>>
> >>> --
> >>> Shannon Rees
> >>> Programming Scientist
> >>> Engility
> >>> Geophysical Fluid Dynamics Lab
> >>> 201 Forrestal Rd Princeton, NJ
> >>> 609-452-5384
> >>>
> >>>
> >>
> >>
> >
> >
> > --
> > Shannon Rees
> > Programming Scientist
> > Engility
> > Geophysical Fluid Dynamics Lab
> > 201 Forrestal Rd Princeton, NJ
> > 609-452-5384
> >
>
------------------------------------------------
Subject: masking region by latitude
From: Shannon Rees - NOAA Affiliate
Time: Fri Nov 13 14:54:12 2015
John,
That fixed it. Thank you!
-Shannon
On Fri, Nov 13, 2015 at 4:36 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Shannon,
>
> In my first test replicating your issue, the gen_vx_mask ran fine
without
> error.
>
> For my second test, I added extra empty blank lines at the end of
the ASCII
> polyline file, and that did the trick. I was able to reproduce your
error:
> ERROR : MaskPoly::load() -> read error in mask poly file
"20_80N.poly"
>
> For now, please just delete all trailing empty lines from your
polyline
> file. Hopefully that'll get it to run.
>
> And I'll look into updating the MET code to handle this situation
more
> gracefully.
>
> Thanks,
> John
>
>
> On Thu, Nov 12, 2015 at 6:44 PM, John NCAR <johnhg at ucar.edu> wrote:
>
> > Shannon,
> >
> > I don't see any obvious problems in the command line. I'll
replicate the
> > behavior tomorrow and let you know what I find.
> >
> > Thanks
> > John
> >
> > Sent from my iPhone
> >
> > On Nov 12, 2015, at 5:46 PM, "Shannon Rees - NOAA Affiliate via
RT" <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74090 >
> > >
> > > Hi John,
> > >
> > > Sorry to disturb your travels. I do want two separate masks for
the
> > > northern and southern hemispheres. I tried your suggestion and
got an
> > > error.
> > >
> > > met-5.1/bin/gen_vx_mask \
> > > met-5.1/data/sample_obs/ST2ml/ST2ml2005080700.Grb_G212 \
> > > 20_80N.poly \
> > > out/gen_vx_mask/20_80N_mask.nc
> > >
> > > ERROR : MaskPoly::load() -> read error in mask poly file
"20_80N.poly"
> > >
> > > I get the same error using different input files as well. Do I
need to
> > use
> > > a special option?
> > >
> > > Thanks
> > >
> > >
> > > On Thu, Nov 12, 2015 at 4:28 PM, John Halley Gotway via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > >> Hi Shannon,
> > >>
> > >> I am out on travel today and can't get VPN running to access my
> machine
> > >> remotely.
> > >>
> > >> I see that you're trying to generate a mask of all grid points
that
> are
> > >> 20-80 degrees north and a second mask for 20-80 degrees south.
> > >>
> > >> Do you actually want those as a single mask or separate masks?
I'll
> > assume
> > >> separate for now.
> > >>
> > >> Here's the mask I'd try for the North Hemisphere (but I haven't
tested
> > with
> > >> this yet)...
> > >>
> > >> 20_80N
> > >> 20 -180
> > >> 80 -180
> > >> 80 180
> > >> 20 180
> > >>
> > >> For the Southern Hemisphere, you'd just switch the 20 and 80
latitude
> to
> > >> -20 and -80.
> > >>
> > >> Another cool option is this... suppose you want all grid points
other
> > than
> > >> -20 to 20 degrees. You could use the following mask and then
use the
> > >> "-complement" option for gen_vx_mask. This will define the -20
to 20
> > >> degree band, but using "-complement" will essentially make it
> > >> reverse-video:
> > >>
> > >> 20N_20S
> > >> -20 -180
> > >> 20 -180
> > >> 20 180
> > >> -20 180
> > >>
> > >> Please let me know how it goes.
> > >>
> > >> Thanks,
> > >> John
> > >>
> > >> On Thu, Nov 12, 2015 at 2:36 PM, Shannon Rees - NOAA Affiliate
via RT
> <
> > >> met_help at ucar.edu> wrote:
> > >>
> > >>>
> > >>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74090
>
> > >>>
> > >>> Here is an update:
> > >>>
> > >>> I created my own simple .poly file:
> > >>> 20_90N
> > >>> 20.00000 0.00000
> > >>> 20.00000 180.00000
> > >>> 20.00000 359.90000
> > >>> 89.00000 0.00000
> > >>> 89.00000 180.00000
> > >>> 89.00000 359.90000
> > >>>
> > >>> I got this error when trying to run gen_vx_mask:
> > >>> DEBUG 1: Input File:
> > >>>
> > >>>
> > >>
> >
> /home/Shannon.McElhinney/met-
5.1/data/sample_obs/ST2ml/ST2ml2005080700.Grb_G212
> > >>> DEBUG 1: Mask File: test.poly
> > >>> DEBUG 2: Parsed Data Grid: Lambert Conformal (185 x 129)
> > >>> ERROR :
> > >>> ERROR : MaskPoly::load() -> read error in mask poly file
"test.poly"
> > >>> ERROR :
> > >>>
> > >>>
> > >>>
> > >>> On Thu, Nov 12, 2015 at 10:43 AM, met_help at ucar.edu via RT <
> > >>> met_help at ucar.edu> wrote:
> > >>>
> > >>>> Greetings,
> > >>>>
> > >>>> This message has been automatically generated in response to
the
> > >>>> creation of a trouble ticket regarding:
> > >>>> "masking region by latitude",
> > >>>> a summary of which appears below.
> > >>>>
> > >>>> There is no need to reply to this message right now. Your
ticket
> has
> > >>> been
> > >>>> assigned an ID of [rt.rap.ucar.edu #74090].
> > >>>>
> > >>>> Please include the string:
> > >>>>
> > >>>> [rt.rap.ucar.edu #74090]
> > >>>>
> > >>>> in the subject line of all future correspondence about this
issue.
> To
> > >> do
> > >>>> so,
> > >>>> you may reply to this message.
> > >>>>
> > >>>> Thank you,
> > >>>> met_help at ucar.edu
> > >>>>
> > >>>>
> > >>
> >
-------------------------------------------------------------------------
> > >>>> Hello,
> > >>>>
> > >>>> I am trying to compute the rmse and acc using only the
regions 20-80
> > >>>> degrees N and 20-80 degrees S. Is there a simple way to do
this
> with
> > >> MET
> > >>>> 5.1? I see there are various .poly files that make it easy
to mask
> > >> these
> > >>>> regions but I do not think any of them are the regions that I
want.
> > >>>>
> > >>>> Best,
> > >>>> --
> > >>>> Shannon Rees
> > >>>> Programming Scientist
> > >>>> Engility
> > >>>> Geophysical Fluid Dynamics Lab
> > >>>> 201 Forrestal Rd Princeton, NJ
> > >>>> 609-452-5384
> > >>>>
> > >>>>
> > >>>
> > >>>
> > >>> --
> > >>> Shannon Rees
> > >>> Programming Scientist
> > >>> Engility
> > >>> Geophysical Fluid Dynamics Lab
> > >>> 201 Forrestal Rd Princeton, NJ
> > >>> 609-452-5384
> > >>>
> > >>>
> > >>
> > >>
> > >
> > >
> > > --
> > > Shannon Rees
> > > Programming Scientist
> > > Engility
> > > Geophysical Fluid Dynamics Lab
> > > 201 Forrestal Rd Princeton, NJ
> > > 609-452-5384
> > >
> >
>
>
--
Shannon Rees
Programming Scientist
Engility
Geophysical Fluid Dynamics Lab
201 Forrestal Rd Princeton, NJ
609-452-5384
------------------------------------------------
Subject: masking region by latitude
From: John Halley Gotway
Time: Fri Nov 13 14:56:29 2015
Shannon,
Thanks for pointing this out.
In future versions of MET, we could either update the error message to
instruct the user to delete trailing whitespace, or assume we're at
the end
of the polyline file rather than erroring out. I've asked the MET
developers for their opinions.
Thanks and have a good weekend.
John
On Fri, Nov 13, 2015 at 2:54 PM, Shannon Rees - NOAA Affiliate via RT
<
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74090 >
>
> John,
>
> That fixed it. Thank you!
>
> -Shannon
>
> On Fri, Nov 13, 2015 at 4:36 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
> > Shannon,
> >
> > In my first test replicating your issue, the gen_vx_mask ran fine
without
> > error.
> >
> > For my second test, I added extra empty blank lines at the end of
the
> ASCII
> > polyline file, and that did the trick. I was able to reproduce
your
> error:
> > ERROR : MaskPoly::load() -> read error in mask poly file
> "20_80N.poly"
> >
> > For now, please just delete all trailing empty lines from your
polyline
> > file. Hopefully that'll get it to run.
> >
> > And I'll look into updating the MET code to handle this situation
more
> > gracefully.
> >
> > Thanks,
> > John
> >
> >
> > On Thu, Nov 12, 2015 at 6:44 PM, John NCAR <johnhg at ucar.edu>
wrote:
> >
> > > Shannon,
> > >
> > > I don't see any obvious problems in the command line. I'll
replicate
> the
> > > behavior tomorrow and let you know what I find.
> > >
> > > Thanks
> > > John
> > >
> > > Sent from my iPhone
> > >
> > > On Nov 12, 2015, at 5:46 PM, "Shannon Rees - NOAA Affiliate via
RT" <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74090
>
> > > >
> > > > Hi John,
> > > >
> > > > Sorry to disturb your travels. I do want two separate masks
for the
> > > > northern and southern hemispheres. I tried your suggestion
and got
> an
> > > > error.
> > > >
> > > > met-5.1/bin/gen_vx_mask \
> > > > met-5.1/data/sample_obs/ST2ml/ST2ml2005080700.Grb_G212 \
> > > > 20_80N.poly \
> > > > out/gen_vx_mask/20_80N_mask.nc
> > > >
> > > > ERROR : MaskPoly::load() -> read error in mask poly file
> "20_80N.poly"
> > > >
> > > > I get the same error using different input files as well. Do
I need
> to
> > > use
> > > > a special option?
> > > >
> > > > Thanks
> > > >
> > > >
> > > > On Thu, Nov 12, 2015 at 4:28 PM, John Halley Gotway via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > >> Hi Shannon,
> > > >>
> > > >> I am out on travel today and can't get VPN running to access
my
> > machine
> > > >> remotely.
> > > >>
> > > >> I see that you're trying to generate a mask of all grid
points that
> > are
> > > >> 20-80 degrees north and a second mask for 20-80 degrees
south.
> > > >>
> > > >> Do you actually want those as a single mask or separate
masks? I'll
> > > assume
> > > >> separate for now.
> > > >>
> > > >> Here's the mask I'd try for the North Hemisphere (but I
haven't
> tested
> > > with
> > > >> this yet)...
> > > >>
> > > >> 20_80N
> > > >> 20 -180
> > > >> 80 -180
> > > >> 80 180
> > > >> 20 180
> > > >>
> > > >> For the Southern Hemisphere, you'd just switch the 20 and 80
> latitude
> > to
> > > >> -20 and -80.
> > > >>
> > > >> Another cool option is this... suppose you want all grid
points
> other
> > > than
> > > >> -20 to 20 degrees. You could use the following mask and then
use
> the
> > > >> "-complement" option for gen_vx_mask. This will define the
-20 to
> 20
> > > >> degree band, but using "-complement" will essentially make it
> > > >> reverse-video:
> > > >>
> > > >> 20N_20S
> > > >> -20 -180
> > > >> 20 -180
> > > >> 20 180
> > > >> -20 180
> > > >>
> > > >> Please let me know how it goes.
> > > >>
> > > >> Thanks,
> > > >> John
> > > >>
> > > >> On Thu, Nov 12, 2015 at 2:36 PM, Shannon Rees - NOAA
Affiliate via
> RT
> > <
> > > >> met_help at ucar.edu> wrote:
> > > >>
> > > >>>
> > > >>> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74090 >
> > > >>>
> > > >>> Here is an update:
> > > >>>
> > > >>> I created my own simple .poly file:
> > > >>> 20_90N
> > > >>> 20.00000 0.00000
> > > >>> 20.00000 180.00000
> > > >>> 20.00000 359.90000
> > > >>> 89.00000 0.00000
> > > >>> 89.00000 180.00000
> > > >>> 89.00000 359.90000
> > > >>>
> > > >>> I got this error when trying to run gen_vx_mask:
> > > >>> DEBUG 1: Input File:
> > > >>>
> > > >>>
> > > >>
> > >
> >
> /home/Shannon.McElhinney/met-
5.1/data/sample_obs/ST2ml/ST2ml2005080700.Grb_G212
> > > >>> DEBUG 1: Mask File: test.poly
> > > >>> DEBUG 2: Parsed Data Grid: Lambert Conformal (185 x
129)
> > > >>> ERROR :
> > > >>> ERROR : MaskPoly::load() -> read error in mask poly file
> "test.poly"
> > > >>> ERROR :
> > > >>>
> > > >>>
> > > >>>
> > > >>> On Thu, Nov 12, 2015 at 10:43 AM, met_help at ucar.edu via RT <
> > > >>> met_help at ucar.edu> wrote:
> > > >>>
> > > >>>> Greetings,
> > > >>>>
> > > >>>> This message has been automatically generated in response
to the
> > > >>>> creation of a trouble ticket regarding:
> > > >>>> "masking region by latitude",
> > > >>>> a summary of which appears below.
> > > >>>>
> > > >>>> There is no need to reply to this message right now. Your
ticket
> > has
> > > >>> been
> > > >>>> assigned an ID of [rt.rap.ucar.edu #74090].
> > > >>>>
> > > >>>> Please include the string:
> > > >>>>
> > > >>>> [rt.rap.ucar.edu #74090]
> > > >>>>
> > > >>>> in the subject line of all future correspondence about this
issue.
> > To
> > > >> do
> > > >>>> so,
> > > >>>> you may reply to this message.
> > > >>>>
> > > >>>> Thank you,
> > > >>>> met_help at ucar.edu
> > > >>>>
> > > >>>>
> > > >>
> > >
>
-------------------------------------------------------------------------
> > > >>>> Hello,
> > > >>>>
> > > >>>> I am trying to compute the rmse and acc using only the
regions
> 20-80
> > > >>>> degrees N and 20-80 degrees S. Is there a simple way to do
this
> > with
> > > >> MET
> > > >>>> 5.1? I see there are various .poly files that make it easy
to
> mask
> > > >> these
> > > >>>> regions but I do not think any of them are the regions that
I
> want.
> > > >>>>
> > > >>>> Best,
> > > >>>> --
> > > >>>> Shannon Rees
> > > >>>> Programming Scientist
> > > >>>> Engility
> > > >>>> Geophysical Fluid Dynamics Lab
> > > >>>> 201 Forrestal Rd Princeton, NJ
> > > >>>> 609-452-5384
> > > >>>>
> > > >>>>
> > > >>>
> > > >>>
> > > >>> --
> > > >>> Shannon Rees
> > > >>> Programming Scientist
> > > >>> Engility
> > > >>> Geophysical Fluid Dynamics Lab
> > > >>> 201 Forrestal Rd Princeton, NJ
> > > >>> 609-452-5384
> > > >>>
> > > >>>
> > > >>
> > > >>
> > > >
> > > >
> > > > --
> > > > Shannon Rees
> > > > Programming Scientist
> > > > Engility
> > > > Geophysical Fluid Dynamics Lab
> > > > 201 Forrestal Rd Princeton, NJ
> > > > 609-452-5384
> > > >
> > >
> >
> >
>
>
> --
> Shannon Rees
> Programming Scientist
> Engility
> Geophysical Fluid Dynamics Lab
> 201 Forrestal Rd Princeton, NJ
> 609-452-5384
>
>
------------------------------------------------
More information about the Met_help
mailing list