[Met_help] [rt.rap.ucar.edu #96443] History for mask file for Metplus
Julie Prestopnik via RT
met_help at ucar.edu
Thu Oct 1 09:04:20 MDT 2020
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hello,
How do I create a msk file to be use in Metplus (point stat)? to limit my
analysis
Thanks
Tatiana
--
*Tatiana D. Gonzalez, PhD.*
*Analysis & Nowcast Branch*
*Analysis and Mission Support Division*
*Analyze, Forecast and Support Office*
*National Weather Service**SSMC2, Room 14115*
*Silver Spring, Maryland*
*Ph. (301) 427-9500 *
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: mask file for Metplus
From: Julie Prestopnik
Time: Tue Aug 25 15:52:11 2020
Hi Tatiana.
I see you have a question about how to create a mask file. I'd
recommend
running the Gen-Vx-Mask tool to create the desired masking region, and
then
reference that in the Point-Stat config file. Here is a link to our
documentation page on Gen-Vx-Mask, but I also go through an example
below:
https://dtcenter.github.io/MET/Users_Guide/masking.html
To illustrate, I pulled a sample GFS file:
> wget
> https://www.ncei.noaa.gov/data/global-forecast-system/access/grid-
004-0.5-degree/forecast/202008/20200820/gfs_4_20200820_0000_012.grb2
I ran wgrib2 and look for a LAND/SEA mask:
> wgrib2 gfs_4_20200820_0000_012.grb2 | egrep -i LAND
>
which returned:
> 571:98325645:d=2020082000:LAND:surface:12 hour fcst:
I ran plot_data_plane on it (see attached image):
> met-9.1/bin/plot_data_plane gfs_4_20200820_0000_012.grb2
gfs_4_LAND.ps
> 'name="LAND"; level="L0";'
>
> DEBUG 1: Opening data file: gfs_4_20200820_0000_012.grb2
> DEBUG 1: Creating postscript file: gfs_4_LAND.ps
>
I ran gen_vx_mask to create a mask for this verification domain. If
you're
actually verifying on the GFS domain, the command looks a little
redundant:
met-9.1/bin/gen_vx_mask met-9.1/bin/gen_vx_mask \
> gfs_4_20200820_0000_012.grb2 \
> gfs_4_20200820_0000_012.grb2 \
> gfs_4_WATER_mask.nc \
> -type data \
> -mask_field 'name="LAND"; level="L0";' \
> -thresh ==0 \
> -name WATER
DEBUG 1: Input Grid: gfs_4_20200820_0000_012.grb2
> DEBUG 1: Mask File: gfs_4_20200820_0000_012.grb2
> DEBUG 2: Parsed Input Grid: LatLon (720 x 361)
> DEBUG 2: Parsed Mask Grid: LatLon (720 x 361)
> DEBUG 1: Output File: gfs_4_WATER_mask.nc
I defined the threshold as "==0" because I want the WATER points, not
the
LAND points. And I used the "-name WATER" command line option to name
the
output NetCDF variable. This will show up in the VX_MASK output
column of
Point-Stat.
And finally, you would then setup your Point-Stat config file like
this:
mask = {
grid = [];
poly = [ "/path/to/gfs_4_WATER_mask.nc" ];
sid = [];
}
I hope that helps!
Julie
On Tue, Aug 25, 2020 at 1:14 PM Tatiana Gonzalez - NOAA Federal via RT
<
met_help at ucar.edu> wrote:
>
> Tue Aug 25 13:14:03 2020: Request 96443 was acted upon.
> Transaction: Ticket created by tatiana.gonzalez at noaa.gov
> Queue: met_help
> Subject: mask file for Metplus
> Owner: Nobody
> Requestors: tatiana.gonzalez at noaa.gov
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
>
>
> Hello,
> How do I create a msk file to be use in Metplus (point stat)? to
limit my
> analysis
>
> Thanks
>
>
> Tatiana
> --
>
>
> *Tatiana D. Gonzalez, PhD.*
> *Analysis & Nowcast Branch*
> *Analysis and Mission Support Division*
>
> *Analyze, Forecast and Support Office*
>
> *National Weather Service**SSMC2, Room 14115*
> *Silver Spring, Maryland*
>
>
> *Ph. (301) 427-9500 *
>
>
--
Julie Prestopnik (she/her/hers)
Software Engineer
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
My working day may not be your working day. Please do not feel
obliged to
reply to this email outside of your normal working hours.
------------------------------------------------
Subject: mask file for Metplus
From: Tatiana Gonzalez - NOAA Federal
Time: Wed Aug 26 12:03:32 2020
Hi Julie,
Thanks for the instructions. I have the RTMA files in grib2 format for
the
whole CONUS. But what I want is only use the data contained in this
box:
South East 45.125 -121.975
North East 46.188 -121.975
North West 46.188 -123.625
South West 45.125 -123.625
I see there is the Polyline (poly) that might work for me. Do you have
an
example for that case or the ASCII file that is needed.
Thank you
Tatiana
On Tue, Aug 25, 2020 at 5:52 PM Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> Hi Tatiana.
>
> I see you have a question about how to create a mask file. I'd
recommend
> running the Gen-Vx-Mask tool to create the desired masking region,
and then
> reference that in the Point-Stat config file. Here is a link to our
> documentation page on Gen-Vx-Mask, but I also go through an example
below:
> https://dtcenter.github.io/MET/Users_Guide/masking.html
>
> To illustrate, I pulled a sample GFS file:
>
> > wget
> >
> https://www.ncei.noaa.gov/data/global-forecast-system/access/grid-
004-0.5-degree/forecast/202008/20200820/gfs_4_20200820_0000_012.grb2
>
>
> I ran wgrib2 and look for a LAND/SEA mask:
>
> > wgrib2 gfs_4_20200820_0000_012.grb2 | egrep -i LAND
> >
> which returned:
>
> > 571:98325645:d=2020082000:LAND:surface:12 hour fcst:
>
>
> I ran plot_data_plane on it (see attached image):
>
> > met-9.1/bin/plot_data_plane gfs_4_20200820_0000_012.grb2
gfs_4_LAND.ps
> > 'name="LAND"; level="L0";'
> >
>
>
> > DEBUG 1: Opening data file: gfs_4_20200820_0000_012.grb2
> > DEBUG 1: Creating postscript file: gfs_4_LAND.ps
> >
>
> I ran gen_vx_mask to create a mask for this verification domain. If
you're
> actually verifying on the GFS domain, the command looks a little
redundant:
>
> met-9.1/bin/gen_vx_mask met-9.1/bin/gen_vx_mask \
> > gfs_4_20200820_0000_012.grb2 \
> > gfs_4_20200820_0000_012.grb2 \
> > gfs_4_WATER_mask.nc \
> > -type data \
> > -mask_field 'name="LAND"; level="L0";' \
> > -thresh ==0 \
> > -name WATER
>
>
> DEBUG 1: Input Grid: gfs_4_20200820_0000_012.grb2
> > DEBUG 1: Mask File: gfs_4_20200820_0000_012.grb2
> > DEBUG 2: Parsed Input Grid: LatLon (720 x 361)
> > DEBUG 2: Parsed Mask Grid: LatLon (720 x 361)
> > DEBUG 1: Output File: gfs_4_WATER_mask.nc
>
>
> I defined the threshold as "==0" because I want the WATER points,
not the
> LAND points. And I used the "-name WATER" command line option to
name the
> output NetCDF variable. This will show up in the VX_MASK output
column of
> Point-Stat.
>
> And finally, you would then setup your Point-Stat config file like
this:
>
> mask = {
> grid = [];
> poly = [ "/path/to/gfs_4_WATER_mask.nc" ];
> sid = [];
> }
>
> I hope that helps!
>
> Julie
>
> On Tue, Aug 25, 2020 at 1:14 PM Tatiana Gonzalez - NOAA Federal via
RT <
> met_help at ucar.edu> wrote:
>
> >
> > Tue Aug 25 13:14:03 2020: Request 96443 was acted upon.
> > Transaction: Ticket created by tatiana.gonzalez at noaa.gov
> > Queue: met_help
> > Subject: mask file for Metplus
> > Owner: Nobody
> > Requestors: tatiana.gonzalez at noaa.gov
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
> >
> >
> > Hello,
> > How do I create a msk file to be use in Metplus (point stat)? to
limit
> my
> > analysis
> >
> > Thanks
> >
> >
> > Tatiana
> > --
> >
> >
> > *Tatiana D. Gonzalez, PhD.*
> > *Analysis & Nowcast Branch*
> > *Analysis and Mission Support Division*
> >
> > *Analyze, Forecast and Support Office*
> >
> > *National Weather Service**SSMC2, Room 14115*
> > *Silver Spring, Maryland*
> >
> >
> > *Ph. (301) 427-9500 *
> >
> >
>
> --
> Julie Prestopnik (she/her/hers)
> Software Engineer
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
> My working day may not be your working day. Please do not feel
obliged to
> reply to this email outside of your normal working hours.
>
>
--
*Tatiana D. Gonzalez, PhD.*
*Analysis & Nowcast Branch*
*Analysis and Mission Support Division*
*Analyze, Forecast and Support Office*
*National Weather Service**SSMC2, Room 14115*
*Silver Spring, Maryland*
*Ph. (301) 427-9500 *
------------------------------------------------
Subject: mask file for Metplus
From: Julie Prestopnik
Time: Wed Aug 26 12:54:35 2020
Hi Tatiana.
To create a polyline mask, you simply create an ASCII file containing
a
list of lat/lons that define the region of interest. For example,
take a
look at met-9.1/data/poly/WEST.poly. The first line of the file
specifies
the name for the polyline. And the remaining lines contain lat/lon
values
for the boundary points. The last point is connected back up to the
first
point to complete the area. You need a minimum of 3 points to define a
polyline, and there is no real maximum. If your polyline contains a
large
number of points (a few hundred) or you have many grid-points, the
process
of determining which grid points are inside/outside the polyline can
be
slow. In that case, I'd suggest running your polyline through the
gen_vx_mask tool.
Julie
On Wed, Aug 26, 2020 at 12:03 PM Tatiana Gonzalez - NOAA Federal via
RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
>
> Hi Julie,
> Thanks for the instructions. I have the RTMA files in grib2 format
for the
> whole CONUS. But what I want is only use the data contained in this
box:
>
> South East 45.125 -121.975
>
> North East 46.188 -121.975
>
> North West 46.188 -123.625
>
> South West 45.125 -123.625
>
> I see there is the Polyline (poly) that might work for me. Do you
have an
> example for that case or the ASCII file that is needed.
>
> Thank you
>
> Tatiana
>
> On Tue, Aug 25, 2020 at 5:52 PM Julie Prestopnik via RT
<met_help at ucar.edu
> >
> wrote:
>
> > Hi Tatiana.
> >
> > I see you have a question about how to create a mask file. I'd
recommend
> > running the Gen-Vx-Mask tool to create the desired masking region,
and
> then
> > reference that in the Point-Stat config file. Here is a link to
our
> > documentation page on Gen-Vx-Mask, but I also go through an
example
> below:
> > https://dtcenter.github.io/MET/Users_Guide/masking.html
> >
> > To illustrate, I pulled a sample GFS file:
> >
> > > wget
> > >
> >
> https://www.ncei.noaa.gov/data/global-forecast-system/access/grid-
004-0.5-degree/forecast/202008/20200820/gfs_4_20200820_0000_012.grb2
> >
> >
> > I ran wgrib2 and look for a LAND/SEA mask:
> >
> > > wgrib2 gfs_4_20200820_0000_012.grb2 | egrep -i LAND
> > >
> > which returned:
> >
> > > 571:98325645:d=2020082000:LAND:surface:12 hour fcst:
> >
> >
> > I ran plot_data_plane on it (see attached image):
> >
> > > met-9.1/bin/plot_data_plane gfs_4_20200820_0000_012.grb2
gfs_4_LAND.ps
> > > 'name="LAND"; level="L0";'
> > >
> >
> >
> > > DEBUG 1: Opening data file: gfs_4_20200820_0000_012.grb2
> > > DEBUG 1: Creating postscript file: gfs_4_LAND.ps
> > >
> >
> > I ran gen_vx_mask to create a mask for this verification domain.
If
> you're
> > actually verifying on the GFS domain, the command looks a little
> redundant:
> >
> > met-9.1/bin/gen_vx_mask met-9.1/bin/gen_vx_mask \
> > > gfs_4_20200820_0000_012.grb2 \
> > > gfs_4_20200820_0000_012.grb2 \
> > > gfs_4_WATER_mask.nc \
> > > -type data \
> > > -mask_field 'name="LAND"; level="L0";' \
> > > -thresh ==0 \
> > > -name WATER
> >
> >
> > DEBUG 1: Input Grid: gfs_4_20200820_0000_012.grb2
> > > DEBUG 1: Mask File: gfs_4_20200820_0000_012.grb2
> > > DEBUG 2: Parsed Input Grid: LatLon (720 x 361)
> > > DEBUG 2: Parsed Mask Grid: LatLon (720 x 361)
> > > DEBUG 1: Output File: gfs_4_WATER_mask.nc
> >
> >
> > I defined the threshold as "==0" because I want the WATER points,
not the
> > LAND points. And I used the "-name WATER" command line option to
name
> the
> > output NetCDF variable. This will show up in the VX_MASK output
column
> of
> > Point-Stat.
> >
> > And finally, you would then setup your Point-Stat config file like
this:
> >
> > mask = {
> > grid = [];
> > poly = [ "/path/to/gfs_4_WATER_mask.nc" ];
> > sid = [];
> > }
> >
> > I hope that helps!
> >
> > Julie
> >
> > On Tue, Aug 25, 2020 at 1:14 PM Tatiana Gonzalez - NOAA Federal
via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > Tue Aug 25 13:14:03 2020: Request 96443 was acted upon.
> > > Transaction: Ticket created by tatiana.gonzalez at noaa.gov
> > > Queue: met_help
> > > Subject: mask file for Metplus
> > > Owner: Nobody
> > > Requestors: tatiana.gonzalez at noaa.gov
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
> >
> > >
> > >
> > > Hello,
> > > How do I create a msk file to be use in Metplus (point stat)?
to limit
> > my
> > > analysis
> > >
> > > Thanks
> > >
> > >
> > > Tatiana
> > > --
> > >
> > >
> > > *Tatiana D. Gonzalez, PhD.*
> > > *Analysis & Nowcast Branch*
> > > *Analysis and Mission Support Division*
> > >
> > > *Analyze, Forecast and Support Office*
> > >
> > > *National Weather Service**SSMC2, Room 14115*
> > > *Silver Spring, Maryland*
> > >
> > >
> > > *Ph. (301) 427-9500 *
> > >
> > >
> >
> > --
> > Julie Prestopnik (she/her/hers)
> > Software Engineer
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
> > My working day may not be your working day. Please do not feel
obliged
> to
> > reply to this email outside of your normal working hours.
> >
> >
>
> --
>
>
> *Tatiana D. Gonzalez, PhD.*
> *Analysis & Nowcast Branch*
> *Analysis and Mission Support Division*
>
> *Analyze, Forecast and Support Office*
>
> *National Weather Service**SSMC2, Room 14115*
> *Silver Spring, Maryland*
>
>
> *Ph. (301) 427-9500 *
>
>
--
Julie Prestopnik (she/her/hers)
Software Engineer
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
My working day may not be your working day. Please do not feel
obliged to
reply to this email outside of your normal working hours.
------------------------------------------------
Subject: mask file for Metplus
From: Tatiana Gonzalez - NOAA Federal
Time: Thu Sep 10 07:10:48 2020
Hello Jullie,
I created a file with 4 coordinates and ran the polyline through the
gen_vx_mask tool. Then on the config file I use the new "file.nc"
But it looks like none of the stations that are within my grid were
included. Is there an order to the coordinates? My file was in the
following order (south east, north east, north west, south west)
MIAMI
25.653243 -79.832631
28.073608 -79.832631
28.073608 -80.938611
25.602028 -80.938611
The grib file I used contained the whole CONUS. Is there something I'm
doing wrong?
Thanks
Tatiana
On Wed, Aug 26, 2020 at 2:54 PM Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> Hi Tatiana.
>
> To create a polyline mask, you simply create an ASCII file
containing a
> list of lat/lons that define the region of interest. For example,
take a
> look at met-9.1/data/poly/WEST.poly. The first line of the file
specifies
> the name for the polyline. And the remaining lines contain lat/lon
values
> for the boundary points. The last point is connected back up to the
first
> point to complete the area. You need a minimum of 3 points to define
a
> polyline, and there is no real maximum. If your polyline contains a
large
> number of points (a few hundred) or you have many grid-points, the
process
> of determining which grid points are inside/outside the polyline can
be
> slow. In that case, I'd suggest running your polyline through the
> gen_vx_mask tool.
>
> Julie
>
> On Wed, Aug 26, 2020 at 12:03 PM Tatiana Gonzalez - NOAA Federal via
RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
> >
> > Hi Julie,
> > Thanks for the instructions. I have the RTMA files in grib2 format
for
> the
> > whole CONUS. But what I want is only use the data contained in
this box:
> >
> > South East 45.125 -121.975
> >
> > North East 46.188 -121.975
> >
> > North West 46.188 -123.625
> >
> > South West 45.125 -123.625
> >
> > I see there is the Polyline (poly) that might work for me. Do you
have an
> > example for that case or the ASCII file that is needed.
> >
> > Thank you
> >
> > Tatiana
> >
> > On Tue, Aug 25, 2020 at 5:52 PM Julie Prestopnik via RT <
> met_help at ucar.edu
> > >
> > wrote:
> >
> > > Hi Tatiana.
> > >
> > > I see you have a question about how to create a mask file. I'd
> recommend
> > > running the Gen-Vx-Mask tool to create the desired masking
region, and
> > then
> > > reference that in the Point-Stat config file. Here is a link to
our
> > > documentation page on Gen-Vx-Mask, but I also go through an
example
> > below:
> > > https://dtcenter.github.io/MET/Users_Guide/masking.html
> > >
> > > To illustrate, I pulled a sample GFS file:
> > >
> > > > wget
> > > >
> > >
> >
> https://www.ncei.noaa.gov/data/global-forecast-system/access/grid-
004-0.5-degree/forecast/202008/20200820/gfs_4_20200820_0000_012.grb2
> > >
> > >
> > > I ran wgrib2 and look for a LAND/SEA mask:
> > >
> > > > wgrib2 gfs_4_20200820_0000_012.grb2 | egrep -i LAND
> > > >
> > > which returned:
> > >
> > > > 571:98325645:d=2020082000:LAND:surface:12 hour fcst:
> > >
> > >
> > > I ran plot_data_plane on it (see attached image):
> > >
> > > > met-9.1/bin/plot_data_plane gfs_4_20200820_0000_012.grb2
> gfs_4_LAND.ps
> > > > 'name="LAND"; level="L0";'
> > > >
> > >
> > >
> > > > DEBUG 1: Opening data file: gfs_4_20200820_0000_012.grb2
> > > > DEBUG 1: Creating postscript file: gfs_4_LAND.ps
> > > >
> > >
> > > I ran gen_vx_mask to create a mask for this verification domain.
If
> > you're
> > > actually verifying on the GFS domain, the command looks a little
> > redundant:
> > >
> > > met-9.1/bin/gen_vx_mask met-9.1/bin/gen_vx_mask \
> > > > gfs_4_20200820_0000_012.grb2 \
> > > > gfs_4_20200820_0000_012.grb2 \
> > > > gfs_4_WATER_mask.nc \
> > > > -type data \
> > > > -mask_field 'name="LAND"; level="L0";' \
> > > > -thresh ==0 \
> > > > -name WATER
> > >
> > >
> > > DEBUG 1: Input Grid: gfs_4_20200820_0000_012.grb2
> > > > DEBUG 1: Mask File: gfs_4_20200820_0000_012.grb2
> > > > DEBUG 2: Parsed Input Grid: LatLon (720 x 361)
> > > > DEBUG 2: Parsed Mask Grid: LatLon (720 x 361)
> > > > DEBUG 1: Output File: gfs_4_WATER_mask.nc
> > >
> > >
> > > I defined the threshold as "==0" because I want the WATER
points, not
> the
> > > LAND points. And I used the "-name WATER" command line option
to name
> > the
> > > output NetCDF variable. This will show up in the VX_MASK output
column
> > of
> > > Point-Stat.
> > >
> > > And finally, you would then setup your Point-Stat config file
like
> this:
> > >
> > > mask = {
> > > grid = [];
> > > poly = [ "/path/to/gfs_4_WATER_mask.nc" ];
> > > sid = [];
> > > }
> > >
> > > I hope that helps!
> > >
> > > Julie
> > >
> > > On Tue, Aug 25, 2020 at 1:14 PM Tatiana Gonzalez - NOAA Federal
via RT
> <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > Tue Aug 25 13:14:03 2020: Request 96443 was acted upon.
> > > > Transaction: Ticket created by tatiana.gonzalez at noaa.gov
> > > > Queue: met_help
> > > > Subject: mask file for Metplus
> > > > Owner: Nobody
> > > > Requestors: tatiana.gonzalez at noaa.gov
> > > > Status: new
> > > > Ticket <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
> > >
> > > >
> > > >
> > > > Hello,
> > > > How do I create a msk file to be use in Metplus (point stat)?
to
> limit
> > > my
> > > > analysis
> > > >
> > > > Thanks
> > > >
> > > >
> > > > Tatiana
> > > > --
> > > >
> > > >
> > > > *Tatiana D. Gonzalez, PhD.*
> > > > *Analysis & Nowcast Branch*
> > > > *Analysis and Mission Support Division*
> > > >
> > > > *Analyze, Forecast and Support Office*
> > > >
> > > > *National Weather Service**SSMC2, Room 14115*
> > > > *Silver Spring, Maryland*
> > > >
> > > >
> > > > *Ph. (301) 427-9500 *
> > > >
> > > >
> > >
> > > --
> > > Julie Prestopnik (she/her/hers)
> > > Software Engineer
> > > National Center for Atmospheric Research
> > > Research Applications Laboratory
> > > Phone: 303.497.8399
> > > Email: jpresto at ucar.edu
> > >
> > > My working day may not be your working day. Please do not feel
obliged
> > to
> > > reply to this email outside of your normal working hours.
> > >
> > >
> >
> > --
> >
> >
> > *Tatiana D. Gonzalez, PhD.*
> > *Analysis & Nowcast Branch*
> > *Analysis and Mission Support Division*
> >
> > *Analyze, Forecast and Support Office*
> >
> > *National Weather Service**SSMC2, Room 14115*
> > *Silver Spring, Maryland*
> >
> >
> > *Ph. (301) 427-9500 *
> >
> >
>
> --
> Julie Prestopnik (she/her/hers)
> Software Engineer
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
> My working day may not be your working day. Please do not feel
obliged to
> reply to this email outside of your normal working hours.
>
>
--
*Tatiana D. Gonzalez, PhD.*
*Analysis & Nowcast Branch*
*Analysis and Mission Support Division*
*Analyze, Forecast and Support Office*
*National Weather Service**SSMC2, Room 14115*
*Silver Spring, Maryland*
*Ph. (301) 427-9500 *
------------------------------------------------
Subject: mask file for Metplus
From: Julie Prestopnik
Time: Thu Sep 10 08:38:46 2020
Hi Tatiana.
Could you please either attach file.nc or upload file.nc to our FTP
server as
described under "How to Send Us Data" on the following page?
https://dtcenter.org/community-code/model-evaluation-tools-met/met-
help-desk
I'll take a look.
Thanks!
Julie
On Thu, Sep 10, 2020, 7:10 AM Tatiana Gonzalez - NOAA Federal via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
>
> Hello Jullie,
> I created a file with 4 coordinates and ran the polyline through the
> gen_vx_mask tool. Then on the config file I use the new "file.nc"
> But it looks like none of the stations that are within my grid were
> included. Is there an order to the coordinates? My file was in the
> following order (south east, north east, north west, south west)
>
> MIAMI
>
> 25.653243 -79.832631
>
> 28.073608 -79.832631
>
> 28.073608 -80.938611
>
> 25.602028 -80.938611
>
>
> The grib file I used contained the whole CONUS. Is there something
I'm
> doing wrong?
>
> Thanks
>
> Tatiana
>
>
>
> On Wed, Aug 26, 2020 at 2:54 PM Julie Prestopnik via RT
<met_help at ucar.edu
> >
> wrote:
>
> > Hi Tatiana.
> >
> > To create a polyline mask, you simply create an ASCII file
containing a
> > list of lat/lons that define the region of interest. For example,
take a
> > look at met-9.1/data/poly/WEST.poly. The first line of the file
> specifies
> > the name for the polyline. And the remaining lines contain
lat/lon
> values
> > for the boundary points. The last point is connected back up to
the
> first
> > point to complete the area. You need a minimum of 3 points to
define a
> > polyline, and there is no real maximum. If your polyline contains
a
> large
> > number of points (a few hundred) or you have many grid-points, the
> process
> > of determining which grid points are inside/outside the polyline
can be
> > slow. In that case, I'd suggest running your polyline through the
> > gen_vx_mask tool.
> >
> > Julie
> >
> > On Wed, Aug 26, 2020 at 12:03 PM Tatiana Gonzalez - NOAA Federal
via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
> > >
> > > Hi Julie,
> > > Thanks for the instructions. I have the RTMA files in grib2
format for
> > the
> > > whole CONUS. But what I want is only use the data contained in
this
> box:
> > >
> > > South East 45.125 -121.975
> > >
> > > North East 46.188 -121.975
> > >
> > > North West 46.188 -123.625
> > >
> > > South West 45.125 -123.625
> > >
> > > I see there is the Polyline (poly) that might work for me. Do
you have
> an
> > > example for that case or the ASCII file that is needed.
> > >
> > > Thank you
> > >
> > > Tatiana
> > >
> > > On Tue, Aug 25, 2020 at 5:52 PM Julie Prestopnik via RT <
> > met_help at ucar.edu
> > > >
> > > wrote:
> > >
> > > > Hi Tatiana.
> > > >
> > > > I see you have a question about how to create a mask file.
I'd
> > recommend
> > > > running the Gen-Vx-Mask tool to create the desired masking
region,
> and
> > > then
> > > > reference that in the Point-Stat config file. Here is a link
to our
> > > > documentation page on Gen-Vx-Mask, but I also go through an
example
> > > below:
> > > > https://dtcenter.github.io/MET/Users_Guide/masking.html
> > > >
> > > > To illustrate, I pulled a sample GFS file:
> > > >
> > > > > wget
> > > > >
> > > >
> > >
> >
> https://www.ncei.noaa.gov/data/global-forecast-system/access/grid-
004-0.5-degree/forecast/202008/20200820/gfs_4_20200820_0000_012.grb2
> > > >
> > > >
> > > > I ran wgrib2 and look for a LAND/SEA mask:
> > > >
> > > > > wgrib2 gfs_4_20200820_0000_012.grb2 | egrep -i LAND
> > > > >
> > > > which returned:
> > > >
> > > > > 571:98325645:d=2020082000:LAND:surface:12 hour fcst:
> > > >
> > > >
> > > > I ran plot_data_plane on it (see attached image):
> > > >
> > > > > met-9.1/bin/plot_data_plane gfs_4_20200820_0000_012.grb2
> > gfs_4_LAND.ps
> > > > > 'name="LAND"; level="L0";'
> > > > >
> > > >
> > > >
> > > > > DEBUG 1: Opening data file: gfs_4_20200820_0000_012.grb2
> > > > > DEBUG 1: Creating postscript file: gfs_4_LAND.ps
> > > > >
> > > >
> > > > I ran gen_vx_mask to create a mask for this verification
domain. If
> > > you're
> > > > actually verifying on the GFS domain, the command looks a
little
> > > redundant:
> > > >
> > > > met-9.1/bin/gen_vx_mask met-9.1/bin/gen_vx_mask \
> > > > > gfs_4_20200820_0000_012.grb2 \
> > > > > gfs_4_20200820_0000_012.grb2 \
> > > > > gfs_4_WATER_mask.nc \
> > > > > -type data \
> > > > > -mask_field 'name="LAND"; level="L0";' \
> > > > > -thresh ==0 \
> > > > > -name WATER
> > > >
> > > >
> > > > DEBUG 1: Input Grid: gfs_4_20200820_0000_012.grb2
> > > > > DEBUG 1: Mask File: gfs_4_20200820_0000_012.grb2
> > > > > DEBUG 2: Parsed Input Grid: LatLon (720 x 361)
> > > > > DEBUG 2: Parsed Mask Grid: LatLon (720 x 361)
> > > > > DEBUG 1: Output File: gfs_4_WATER_mask.nc
> > > >
> > > >
> > > > I defined the threshold as "==0" because I want the WATER
points, not
> > the
> > > > LAND points. And I used the "-name WATER" command line option
to
> name
> > > the
> > > > output NetCDF variable. This will show up in the VX_MASK
output
> column
> > > of
> > > > Point-Stat.
> > > >
> > > > And finally, you would then setup your Point-Stat config file
like
> > this:
> > > >
> > > > mask = {
> > > > grid = [];
> > > > poly = [ "/path/to/gfs_4_WATER_mask.nc" ];
> > > > sid = [];
> > > > }
> > > >
> > > > I hope that helps!
> > > >
> > > > Julie
> > > >
> > > > On Tue, Aug 25, 2020 at 1:14 PM Tatiana Gonzalez - NOAA
Federal via
> RT
> > <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > Tue Aug 25 13:14:03 2020: Request 96443 was acted upon.
> > > > > Transaction: Ticket created by tatiana.gonzalez at noaa.gov
> > > > > Queue: met_help
> > > > > Subject: mask file for Metplus
> > > > > Owner: Nobody
> > > > > Requestors: tatiana.gonzalez at noaa.gov
> > > > > Status: new
> > > > > Ticket <URL:
> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
> > > >
> > > > >
> > > > >
> > > > > Hello,
> > > > > How do I create a msk file to be use in Metplus (point
stat)? to
> > limit
> > > > my
> > > > > analysis
> > > > >
> > > > > Thanks
> > > > >
> > > > >
> > > > > Tatiana
> > > > > --
> > > > >
> > > > >
> > > > > *Tatiana D. Gonzalez, PhD.*
> > > > > *Analysis & Nowcast Branch*
> > > > > *Analysis and Mission Support Division*
> > > > >
> > > > > *Analyze, Forecast and Support Office*
> > > > >
> > > > > *National Weather Service**SSMC2, Room 14115*
> > > > > *Silver Spring, Maryland*
> > > > >
> > > > >
> > > > > *Ph. (301) 427-9500 *
> > > > >
> > > > >
> > > >
> > > > --
> > > > Julie Prestopnik (she/her/hers)
> > > > Software Engineer
> > > > National Center for Atmospheric Research
> > > > Research Applications Laboratory
> > > > Phone: 303.497.8399
> > > > Email: jpresto at ucar.edu
> > > >
> > > > My working day may not be your working day. Please do not
feel
> obliged
> > > to
> > > > reply to this email outside of your normal working hours.
> > > >
> > > >
> > >
> > > --
> > >
> > >
> > > *Tatiana D. Gonzalez, PhD.*
> > > *Analysis & Nowcast Branch*
> > > *Analysis and Mission Support Division*
> > >
> > > *Analyze, Forecast and Support Office*
> > >
> > > *National Weather Service**SSMC2, Room 14115*
> > > *Silver Spring, Maryland*
> > >
> > >
> > > *Ph. (301) 427-9500 *
> > >
> > >
> >
> > --
> > Julie Prestopnik (she/her/hers)
> > Software Engineer
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
> > My working day may not be your working day. Please do not feel
obliged
> to
> > reply to this email outside of your normal working hours.
> >
> >
>
> --
>
>
> *Tatiana D. Gonzalez, PhD.*
> *Analysis & Nowcast Branch*
> *Analysis and Mission Support Division*
>
> *Analyze, Forecast and Support Office*
>
> *National Weather Service**SSMC2, Room 14115*
> *Silver Spring, Maryland*
>
>
> *Ph. (301) 427-9500 *
>
>
------------------------------------------------
Subject: mask file for Metplus
From: Tatiana Gonzalez - NOAA Federal
Time: Thu Sep 10 08:58:49 2020
Jullie,
Thank you for reviewing the file
case_2.nc
<https://drive.google.com/file/d/1tWThlec9seGOruztsodx73wHAD3OU-
gV/view?usp=drive_web>
Tatiana
On Thu, Sep 10, 2020 at 10:38 AM Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> Hi Tatiana.
>
> Could you please either attach file.nc or upload file.nc to our FTP
> server as
> described under "How to Send Us Data" on the following page?
>
> https://dtcenter.org/community-code/model-evaluation-tools-met/met-
help-desk
>
>
> I'll take a look.
>
> Thanks!
>
> Julie
>
>
> On Thu, Sep 10, 2020, 7:10 AM Tatiana Gonzalez - NOAA Federal via RT
<
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
> >
> > Hello Jullie,
> > I created a file with 4 coordinates and ran the polyline through
the
> > gen_vx_mask tool. Then on the config file I use the new "file.nc"
> > But it looks like none of the stations that are within my grid
were
> > included. Is there an order to the coordinates? My file was in the
> > following order (south east, north east, north west, south west)
> >
> > MIAMI
> >
> > 25.653243 -79.832631
> >
> > 28.073608 -79.832631
> >
> > 28.073608 -80.938611
> >
> > 25.602028 -80.938611
> >
> >
> > The grib file I used contained the whole CONUS. Is there something
I'm
> > doing wrong?
> >
> > Thanks
> >
> > Tatiana
> >
> >
> >
> > On Wed, Aug 26, 2020 at 2:54 PM Julie Prestopnik via RT <
> met_help at ucar.edu
> > >
> > wrote:
> >
> > > Hi Tatiana.
> > >
> > > To create a polyline mask, you simply create an ASCII file
containing a
> > > list of lat/lons that define the region of interest. For
example,
> take a
> > > look at met-9.1/data/poly/WEST.poly. The first line of the file
> > specifies
> > > the name for the polyline. And the remaining lines contain
lat/lon
> > values
> > > for the boundary points. The last point is connected back up to
the
> > first
> > > point to complete the area. You need a minimum of 3 points to
define a
> > > polyline, and there is no real maximum. If your polyline
contains a
> > large
> > > number of points (a few hundred) or you have many grid-points,
the
> > process
> > > of determining which grid points are inside/outside the polyline
can be
> > > slow. In that case, I'd suggest running your polyline through
the
> > > gen_vx_mask tool.
> > >
> > > Julie
> > >
> > > On Wed, Aug 26, 2020 at 12:03 PM Tatiana Gonzalez - NOAA Federal
via
> RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
>
> > > >
> > > > Hi Julie,
> > > > Thanks for the instructions. I have the RTMA files in grib2
format
> for
> > > the
> > > > whole CONUS. But what I want is only use the data contained in
this
> > box:
> > > >
> > > > South East 45.125 -121.975
> > > >
> > > > North East 46.188 -121.975
> > > >
> > > > North West 46.188 -123.625
> > > >
> > > > South West 45.125 -123.625
> > > >
> > > > I see there is the Polyline (poly) that might work for me. Do
you
> have
> > an
> > > > example for that case or the ASCII file that is needed.
> > > >
> > > > Thank you
> > > >
> > > > Tatiana
> > > >
> > > > On Tue, Aug 25, 2020 at 5:52 PM Julie Prestopnik via RT <
> > > met_help at ucar.edu
> > > > >
> > > > wrote:
> > > >
> > > > > Hi Tatiana.
> > > > >
> > > > > I see you have a question about how to create a mask file.
I'd
> > > recommend
> > > > > running the Gen-Vx-Mask tool to create the desired masking
region,
> > and
> > > > then
> > > > > reference that in the Point-Stat config file. Here is a
link to
> our
> > > > > documentation page on Gen-Vx-Mask, but I also go through an
example
> > > > below:
> > > > > https://dtcenter.github.io/MET/Users_Guide/masking.html
> > > > >
> > > > > To illustrate, I pulled a sample GFS file:
> > > > >
> > > > > > wget
> > > > > >
> > > > >
> > > >
> > >
> >
> https://www.ncei.noaa.gov/data/global-forecast-system/access/grid-
004-0.5-degree/forecast/202008/20200820/gfs_4_20200820_0000_012.grb2
> > > > >
> > > > >
> > > > > I ran wgrib2 and look for a LAND/SEA mask:
> > > > >
> > > > > > wgrib2 gfs_4_20200820_0000_012.grb2 | egrep -i LAND
> > > > > >
> > > > > which returned:
> > > > >
> > > > > > 571:98325645:d=2020082000:LAND:surface:12 hour fcst:
> > > > >
> > > > >
> > > > > I ran plot_data_plane on it (see attached image):
> > > > >
> > > > > > met-9.1/bin/plot_data_plane gfs_4_20200820_0000_012.grb2
> > > gfs_4_LAND.ps
> > > > > > 'name="LAND"; level="L0";'
> > > > > >
> > > > >
> > > > >
> > > > > > DEBUG 1: Opening data file: gfs_4_20200820_0000_012.grb2
> > > > > > DEBUG 1: Creating postscript file: gfs_4_LAND.ps
> > > > > >
> > > > >
> > > > > I ran gen_vx_mask to create a mask for this verification
domain.
> If
> > > > you're
> > > > > actually verifying on the GFS domain, the command looks a
little
> > > > redundant:
> > > > >
> > > > > met-9.1/bin/gen_vx_mask met-9.1/bin/gen_vx_mask \
> > > > > > gfs_4_20200820_0000_012.grb2 \
> > > > > > gfs_4_20200820_0000_012.grb2 \
> > > > > > gfs_4_WATER_mask.nc \
> > > > > > -type data \
> > > > > > -mask_field 'name="LAND"; level="L0";' \
> > > > > > -thresh ==0 \
> > > > > > -name WATER
> > > > >
> > > > >
> > > > > DEBUG 1: Input Grid: gfs_4_20200820_0000_012.grb2
> > > > > > DEBUG 1: Mask File: gfs_4_20200820_0000_012.grb2
> > > > > > DEBUG 2: Parsed Input Grid: LatLon (720 x 361)
> > > > > > DEBUG 2: Parsed Mask Grid: LatLon (720 x 361)
> > > > > > DEBUG 1: Output File: gfs_4_WATER_mask.nc
> > > > >
> > > > >
> > > > > I defined the threshold as "==0" because I want the WATER
points,
> not
> > > the
> > > > > LAND points. And I used the "-name WATER" command line
option to
> > name
> > > > the
> > > > > output NetCDF variable. This will show up in the VX_MASK
output
> > column
> > > > of
> > > > > Point-Stat.
> > > > >
> > > > > And finally, you would then setup your Point-Stat config
file like
> > > this:
> > > > >
> > > > > mask = {
> > > > > grid = [];
> > > > > poly = [ "/path/to/gfs_4_WATER_mask.nc" ];
> > > > > sid = [];
> > > > > }
> > > > >
> > > > > I hope that helps!
> > > > >
> > > > > Julie
> > > > >
> > > > > On Tue, Aug 25, 2020 at 1:14 PM Tatiana Gonzalez - NOAA
Federal via
> > RT
> > > <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > > >
> > > > > > Tue Aug 25 13:14:03 2020: Request 96443 was acted upon.
> > > > > > Transaction: Ticket created by tatiana.gonzalez at noaa.gov
> > > > > > Queue: met_help
> > > > > > Subject: mask file for Metplus
> > > > > > Owner: Nobody
> > > > > > Requestors: tatiana.gonzalez at noaa.gov
> > > > > > Status: new
> > > > > > Ticket <URL:
> > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
> > > > >
> > > > > >
> > > > > >
> > > > > > Hello,
> > > > > > How do I create a msk file to be use in Metplus (point
stat)? to
> > > limit
> > > > > my
> > > > > > analysis
> > > > > >
> > > > > > Thanks
> > > > > >
> > > > > >
> > > > > > Tatiana
> > > > > > --
> > > > > >
> > > > > >
> > > > > > *Tatiana D. Gonzalez, PhD.*
> > > > > > *Analysis & Nowcast Branch*
> > > > > > *Analysis and Mission Support Division*
> > > > > >
> > > > > > *Analyze, Forecast and Support Office*
> > > > > >
> > > > > > *National Weather Service**SSMC2, Room 14115*
> > > > > > *Silver Spring, Maryland*
> > > > > >
> > > > > >
> > > > > > *Ph. (301) 427-9500 *
> > > > > >
> > > > > >
> > > > >
> > > > > --
> > > > > Julie Prestopnik (she/her/hers)
> > > > > Software Engineer
> > > > > National Center for Atmospheric Research
> > > > > Research Applications Laboratory
> > > > > Phone: 303.497.8399
> > > > > Email: jpresto at ucar.edu
> > > > >
> > > > > My working day may not be your working day. Please do not
feel
> > obliged
> > > > to
> > > > > reply to this email outside of your normal working hours.
> > > > >
> > > > >
> > > >
> > > > --
> > > >
> > > >
> > > > *Tatiana D. Gonzalez, PhD.*
> > > > *Analysis & Nowcast Branch*
> > > > *Analysis and Mission Support Division*
> > > >
> > > > *Analyze, Forecast and Support Office*
> > > >
> > > > *National Weather Service**SSMC2, Room 14115*
> > > > *Silver Spring, Maryland*
> > > >
> > > >
> > > > *Ph. (301) 427-9500 *
> > > >
> > > >
> > >
> > > --
> > > Julie Prestopnik (she/her/hers)
> > > Software Engineer
> > > National Center for Atmospheric Research
> > > Research Applications Laboratory
> > > Phone: 303.497.8399
> > > Email: jpresto at ucar.edu
> > >
> > > My working day may not be your working day. Please do not feel
obliged
> > to
> > > reply to this email outside of your normal working hours.
> > >
> > >
> >
> > --
> >
> >
> > *Tatiana D. Gonzalez, PhD.*
> > *Analysis & Nowcast Branch*
> > *Analysis and Mission Support Division*
> >
> > *Analyze, Forecast and Support Office*
> >
> > *National Weather Service**SSMC2, Room 14115*
> > *Silver Spring, Maryland*
> >
> >
> > *Ph. (301) 427-9500 *
> >
> >
>
>
--
*Tatiana D. Gonzalez, PhD.*
*Analysis & Nowcast Branch*
*Analysis and Mission Support Division*
*Analyze, Forecast and Support Office*
*National Weather Service**SSMC2, Room 14115*
*Silver Spring, Maryland*
*Ph. (301) 427-9500 *
------------------------------------------------
Subject: mask file for Metplus
From: Julie Prestopnik
Time: Thu Sep 10 09:35:49 2020
Hi Tatiana.
I sent a request for access through google drive as I did not have
permissions on the file.
Julie
On Thu, Sep 10, 2020 at 8:58 AM Tatiana Gonzalez - NOAA Federal via RT
<
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
>
> Jullie,
> Thank you for reviewing the file
> case_2.nc
> <
> https://drive.google.com/file/d/1tWThlec9seGOruztsodx73wHAD3OU-
gV/view?usp=drive_web
> >
>
>
> Tatiana
>
> On Thu, Sep 10, 2020 at 10:38 AM Julie Prestopnik via RT <
> met_help at ucar.edu>
> wrote:
>
> > Hi Tatiana.
> >
> > Could you please either attach file.nc or upload file.nc to our
FTP
> > server as
> > described under "How to Send Us Data" on the following page?
> >
> >
> https://dtcenter.org/community-code/model-evaluation-tools-met/met-
help-desk
> >
> >
> > I'll take a look.
> >
> > Thanks!
> >
> > Julie
> >
> >
> > On Thu, Sep 10, 2020, 7:10 AM Tatiana Gonzalez - NOAA Federal via
RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
> > >
> > > Hello Jullie,
> > > I created a file with 4 coordinates and ran the polyline through
the
> > > gen_vx_mask tool. Then on the config file I use the new
"file.nc"
> > > But it looks like none of the stations that are within my grid
were
> > > included. Is there an order to the coordinates? My file was in
the
> > > following order (south east, north east, north west, south west)
> > >
> > > MIAMI
> > >
> > > 25.653243 -79.832631
> > >
> > > 28.073608 -79.832631
> > >
> > > 28.073608 -80.938611
> > >
> > > 25.602028 -80.938611
> > >
> > >
> > > The grib file I used contained the whole CONUS. Is there
something I'm
> > > doing wrong?
> > >
> > > Thanks
> > >
> > > Tatiana
> > >
> > >
> > >
> > > On Wed, Aug 26, 2020 at 2:54 PM Julie Prestopnik via RT <
> > met_help at ucar.edu
> > > >
> > > wrote:
> > >
> > > > Hi Tatiana.
> > > >
> > > > To create a polyline mask, you simply create an ASCII file
> containing a
> > > > list of lat/lons that define the region of interest. For
example,
> > take a
> > > > look at met-9.1/data/poly/WEST.poly. The first line of the
file
> > > specifies
> > > > the name for the polyline. And the remaining lines contain
lat/lon
> > > values
> > > > for the boundary points. The last point is connected back up
to the
> > > first
> > > > point to complete the area. You need a minimum of 3 points to
define
> a
> > > > polyline, and there is no real maximum. If your polyline
contains a
> > > large
> > > > number of points (a few hundred) or you have many grid-points,
the
> > > process
> > > > of determining which grid points are inside/outside the
polyline can
> be
> > > > slow. In that case, I'd suggest running your polyline through
the
> > > > gen_vx_mask tool.
> > > >
> > > > Julie
> > > >
> > > > On Wed, Aug 26, 2020 at 12:03 PM Tatiana Gonzalez - NOAA
Federal via
> > RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
> > > > >
> > > > > Hi Julie,
> > > > > Thanks for the instructions. I have the RTMA files in grib2
format
> > for
> > > > the
> > > > > whole CONUS. But what I want is only use the data contained
in this
> > > box:
> > > > >
> > > > > South East 45.125 -121.975
> > > > >
> > > > > North East 46.188 -121.975
> > > > >
> > > > > North West 46.188 -123.625
> > > > >
> > > > > South West 45.125 -123.625
> > > > >
> > > > > I see there is the Polyline (poly) that might work for me.
Do you
> > have
> > > an
> > > > > example for that case or the ASCII file that is needed.
> > > > >
> > > > > Thank you
> > > > >
> > > > > Tatiana
> > > > >
> > > > > On Tue, Aug 25, 2020 at 5:52 PM Julie Prestopnik via RT <
> > > > met_help at ucar.edu
> > > > > >
> > > > > wrote:
> > > > >
> > > > > > Hi Tatiana.
> > > > > >
> > > > > > I see you have a question about how to create a mask file.
I'd
> > > > recommend
> > > > > > running the Gen-Vx-Mask tool to create the desired masking
> region,
> > > and
> > > > > then
> > > > > > reference that in the Point-Stat config file. Here is a
link to
> > our
> > > > > > documentation page on Gen-Vx-Mask, but I also go through
an
> example
> > > > > below:
> > > > > > https://dtcenter.github.io/MET/Users_Guide/masking.html
> > > > > >
> > > > > > To illustrate, I pulled a sample GFS file:
> > > > > >
> > > > > > > wget
> > > > > > >
> > > > > >
> > > > >
> > > >
> > >
> >
> https://www.ncei.noaa.gov/data/global-forecast-system/access/grid-
004-0.5-degree/forecast/202008/20200820/gfs_4_20200820_0000_012.grb2
> > > > > >
> > > > > >
> > > > > > I ran wgrib2 and look for a LAND/SEA mask:
> > > > > >
> > > > > > > wgrib2 gfs_4_20200820_0000_012.grb2 | egrep -i LAND
> > > > > > >
> > > > > > which returned:
> > > > > >
> > > > > > > 571:98325645:d=2020082000:LAND:surface:12 hour fcst:
> > > > > >
> > > > > >
> > > > > > I ran plot_data_plane on it (see attached image):
> > > > > >
> > > > > > > met-9.1/bin/plot_data_plane gfs_4_20200820_0000_012.grb2
> > > > gfs_4_LAND.ps
> > > > > > > 'name="LAND"; level="L0";'
> > > > > > >
> > > > > >
> > > > > >
> > > > > > > DEBUG 1: Opening data file: gfs_4_20200820_0000_012.grb2
> > > > > > > DEBUG 1: Creating postscript file: gfs_4_LAND.ps
> > > > > > >
> > > > > >
> > > > > > I ran gen_vx_mask to create a mask for this verification
domain.
> > If
> > > > > you're
> > > > > > actually verifying on the GFS domain, the command looks a
little
> > > > > redundant:
> > > > > >
> > > > > > met-9.1/bin/gen_vx_mask met-9.1/bin/gen_vx_mask \
> > > > > > > gfs_4_20200820_0000_012.grb2 \
> > > > > > > gfs_4_20200820_0000_012.grb2 \
> > > > > > > gfs_4_WATER_mask.nc \
> > > > > > > -type data \
> > > > > > > -mask_field 'name="LAND"; level="L0";' \
> > > > > > > -thresh ==0 \
> > > > > > > -name WATER
> > > > > >
> > > > > >
> > > > > > DEBUG 1: Input Grid: gfs_4_20200820_0000_012.grb2
> > > > > > > DEBUG 1: Mask File: gfs_4_20200820_0000_012.grb2
> > > > > > > DEBUG 2: Parsed Input Grid: LatLon (720 x 361)
> > > > > > > DEBUG 2: Parsed Mask Grid: LatLon (720 x 361)
> > > > > > > DEBUG 1: Output File: gfs_4_WATER_mask.nc
> > > > > >
> > > > > >
> > > > > > I defined the threshold as "==0" because I want the WATER
points,
> > not
> > > > the
> > > > > > LAND points. And I used the "-name WATER" command line
option to
> > > name
> > > > > the
> > > > > > output NetCDF variable. This will show up in the VX_MASK
output
> > > column
> > > > > of
> > > > > > Point-Stat.
> > > > > >
> > > > > > And finally, you would then setup your Point-Stat config
file
> like
> > > > this:
> > > > > >
> > > > > > mask = {
> > > > > > grid = [];
> > > > > > poly = [ "/path/to/gfs_4_WATER_mask.nc" ];
> > > > > > sid = [];
> > > > > > }
> > > > > >
> > > > > > I hope that helps!
> > > > > >
> > > > > > Julie
> > > > > >
> > > > > > On Tue, Aug 25, 2020 at 1:14 PM Tatiana Gonzalez - NOAA
Federal
> via
> > > RT
> > > > <
> > > > > > met_help at ucar.edu> wrote:
> > > > > >
> > > > > > >
> > > > > > > Tue Aug 25 13:14:03 2020: Request 96443 was acted upon.
> > > > > > > Transaction: Ticket created by tatiana.gonzalez at noaa.gov
> > > > > > > Queue: met_help
> > > > > > > Subject: mask file for Metplus
> > > > > > > Owner: Nobody
> > > > > > > Requestors: tatiana.gonzalez at noaa.gov
> > > > > > > Status: new
> > > > > > > Ticket <URL:
> > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
> > > > > >
> > > > > > >
> > > > > > >
> > > > > > > Hello,
> > > > > > > How do I create a msk file to be use in Metplus (point
stat)?
> to
> > > > limit
> > > > > > my
> > > > > > > analysis
> > > > > > >
> > > > > > > Thanks
> > > > > > >
> > > > > > >
> > > > > > > Tatiana
> > > > > > > --
> > > > > > >
> > > > > > >
> > > > > > > *Tatiana D. Gonzalez, PhD.*
> > > > > > > *Analysis & Nowcast Branch*
> > > > > > > *Analysis and Mission Support Division*
> > > > > > >
> > > > > > > *Analyze, Forecast and Support Office*
> > > > > > >
> > > > > > > *National Weather Service**SSMC2, Room 14115*
> > > > > > > *Silver Spring, Maryland*
> > > > > > >
> > > > > > >
> > > > > > > *Ph. (301) 427-9500 *
> > > > > > >
> > > > > > >
> > > > > >
> > > > > > --
> > > > > > Julie Prestopnik (she/her/hers)
> > > > > > Software Engineer
> > > > > > National Center for Atmospheric Research
> > > > > > Research Applications Laboratory
> > > > > > Phone: 303.497.8399
> > > > > > Email: jpresto at ucar.edu
> > > > > >
> > > > > > My working day may not be your working day. Please do not
feel
> > > obliged
> > > > > to
> > > > > > reply to this email outside of your normal working hours.
> > > > > >
> > > > > >
> > > > >
> > > > > --
> > > > >
> > > > >
> > > > > *Tatiana D. Gonzalez, PhD.*
> > > > > *Analysis & Nowcast Branch*
> > > > > *Analysis and Mission Support Division*
> > > > >
> > > > > *Analyze, Forecast and Support Office*
> > > > >
> > > > > *National Weather Service**SSMC2, Room 14115*
> > > > > *Silver Spring, Maryland*
> > > > >
> > > > >
> > > > > *Ph. (301) 427-9500 *
> > > > >
> > > > >
> > > >
> > > > --
> > > > Julie Prestopnik (she/her/hers)
> > > > Software Engineer
> > > > National Center for Atmospheric Research
> > > > Research Applications Laboratory
> > > > Phone: 303.497.8399
> > > > Email: jpresto at ucar.edu
> > > >
> > > > My working day may not be your working day. Please do not
feel
> obliged
> > > to
> > > > reply to this email outside of your normal working hours.
> > > >
> > > >
> > >
> > > --
> > >
> > >
> > > *Tatiana D. Gonzalez, PhD.*
> > > *Analysis & Nowcast Branch*
> > > *Analysis and Mission Support Division*
> > >
> > > *Analyze, Forecast and Support Office*
> > >
> > > *National Weather Service**SSMC2, Room 14115*
> > > *Silver Spring, Maryland*
> > >
> > >
> > > *Ph. (301) 427-9500 *
> > >
> > >
> >
> >
>
> --
>
>
> *Tatiana D. Gonzalez, PhD.*
> *Analysis & Nowcast Branch*
> *Analysis and Mission Support Division*
>
> *Analyze, Forecast and Support Office*
>
> *National Weather Service**SSMC2, Room 14115*
> *Silver Spring, Maryland*
>
>
> *Ph. (301) 427-9500 *
>
>
--
Julie Prestopnik (she/her/hers)
Software Engineer
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
My working day may not be your working day. Please do not feel
obliged to
reply to this email outside of your normal working hours.
------------------------------------------------
Subject: mask file for Metplus
From: Tatiana Gonzalez - NOAA Federal
Time: Thu Sep 10 09:46:02 2020
Try again. I just gave you permission
On Thu, Sep 10, 2020 at 11:36 AM Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> Hi Tatiana.
>
> I sent a request for access through google drive as I did not have
> permissions on the file.
>
> Julie
>
> On Thu, Sep 10, 2020 at 8:58 AM Tatiana Gonzalez - NOAA Federal via
RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
> >
> > Jullie,
> > Thank you for reviewing the file
> > case_2.nc
> > <
> >
> https://drive.google.com/file/d/1tWThlec9seGOruztsodx73wHAD3OU-
gV/view?usp=drive_web
> > >
> >
> >
> > Tatiana
> >
> > On Thu, Sep 10, 2020 at 10:38 AM Julie Prestopnik via RT <
> > met_help at ucar.edu>
> > wrote:
> >
> > > Hi Tatiana.
> > >
> > > Could you please either attach file.nc or upload file.nc to our
FTP
> > > server as
> > > described under "How to Send Us Data" on the following page?
> > >
> > >
> >
> https://dtcenter.org/community-code/model-evaluation-tools-met/met-
help-desk
> > >
> > >
> > > I'll take a look.
> > >
> > > Thanks!
> > >
> > > Julie
> > >
> > >
> > > On Thu, Sep 10, 2020, 7:10 AM Tatiana Gonzalez - NOAA Federal
via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
>
> > > >
> > > > Hello Jullie,
> > > > I created a file with 4 coordinates and ran the polyline
through the
> > > > gen_vx_mask tool. Then on the config file I use the new
"file.nc"
> > > > But it looks like none of the stations that are within my grid
were
> > > > included. Is there an order to the coordinates? My file was in
the
> > > > following order (south east, north east, north west, south
west)
> > > >
> > > > MIAMI
> > > >
> > > > 25.653243 -79.832631
> > > >
> > > > 28.073608 -79.832631
> > > >
> > > > 28.073608 -80.938611
> > > >
> > > > 25.602028 -80.938611
> > > >
> > > >
> > > > The grib file I used contained the whole CONUS. Is there
something
> I'm
> > > > doing wrong?
> > > >
> > > > Thanks
> > > >
> > > > Tatiana
> > > >
> > > >
> > > >
> > > > On Wed, Aug 26, 2020 at 2:54 PM Julie Prestopnik via RT <
> > > met_help at ucar.edu
> > > > >
> > > > wrote:
> > > >
> > > > > Hi Tatiana.
> > > > >
> > > > > To create a polyline mask, you simply create an ASCII file
> > containing a
> > > > > list of lat/lons that define the region of interest. For
example,
> > > take a
> > > > > look at met-9.1/data/poly/WEST.poly. The first line of the
file
> > > > specifies
> > > > > the name for the polyline. And the remaining lines contain
lat/lon
> > > > values
> > > > > for the boundary points. The last point is connected back
up to
> the
> > > > first
> > > > > point to complete the area. You need a minimum of 3 points
to
> define
> > a
> > > > > polyline, and there is no real maximum. If your polyline
contains
> a
> > > > large
> > > > > number of points (a few hundred) or you have many grid-
points, the
> > > > process
> > > > > of determining which grid points are inside/outside the
polyline
> can
> > be
> > > > > slow. In that case, I'd suggest running your polyline
through the
> > > > > gen_vx_mask tool.
> > > > >
> > > > > Julie
> > > > >
> > > > > On Wed, Aug 26, 2020 at 12:03 PM Tatiana Gonzalez - NOAA
Federal
> via
> > > RT <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > > >
> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
> > > > > >
> > > > > > Hi Julie,
> > > > > > Thanks for the instructions. I have the RTMA files in
grib2
> format
> > > for
> > > > > the
> > > > > > whole CONUS. But what I want is only use the data
contained in
> this
> > > > box:
> > > > > >
> > > > > > South East 45.125 -121.975
> > > > > >
> > > > > > North East 46.188 -121.975
> > > > > >
> > > > > > North West 46.188 -123.625
> > > > > >
> > > > > > South West 45.125 -123.625
> > > > > >
> > > > > > I see there is the Polyline (poly) that might work for me.
Do you
> > > have
> > > > an
> > > > > > example for that case or the ASCII file that is needed.
> > > > > >
> > > > > > Thank you
> > > > > >
> > > > > > Tatiana
> > > > > >
> > > > > > On Tue, Aug 25, 2020 at 5:52 PM Julie Prestopnik via RT <
> > > > > met_help at ucar.edu
> > > > > > >
> > > > > > wrote:
> > > > > >
> > > > > > > Hi Tatiana.
> > > > > > >
> > > > > > > I see you have a question about how to create a mask
file. I'd
> > > > > recommend
> > > > > > > running the Gen-Vx-Mask tool to create the desired
masking
> > region,
> > > > and
> > > > > > then
> > > > > > > reference that in the Point-Stat config file. Here is a
link
> to
> > > our
> > > > > > > documentation page on Gen-Vx-Mask, but I also go through
an
> > example
> > > > > > below:
> > > > > > > https://dtcenter.github.io/MET/Users_Guide/masking.html
> > > > > > >
> > > > > > > To illustrate, I pulled a sample GFS file:
> > > > > > >
> > > > > > > > wget
> > > > > > > >
> > > > > > >
> > > > > >
> > > > >
> > > >
> > >
> >
> https://www.ncei.noaa.gov/data/global-forecast-system/access/grid-
004-0.5-degree/forecast/202008/20200820/gfs_4_20200820_0000_012.grb2
> > > > > > >
> > > > > > >
> > > > > > > I ran wgrib2 and look for a LAND/SEA mask:
> > > > > > >
> > > > > > > > wgrib2 gfs_4_20200820_0000_012.grb2 | egrep -i LAND
> > > > > > > >
> > > > > > > which returned:
> > > > > > >
> > > > > > > > 571:98325645:d=2020082000:LAND:surface:12 hour fcst:
> > > > > > >
> > > > > > >
> > > > > > > I ran plot_data_plane on it (see attached image):
> > > > > > >
> > > > > > > > met-9.1/bin/plot_data_plane
gfs_4_20200820_0000_012.grb2
> > > > > gfs_4_LAND.ps
> > > > > > > > 'name="LAND"; level="L0";'
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > > DEBUG 1: Opening data file:
gfs_4_20200820_0000_012.grb2
> > > > > > > > DEBUG 1: Creating postscript file: gfs_4_LAND.ps
> > > > > > > >
> > > > > > >
> > > > > > > I ran gen_vx_mask to create a mask for this verification
> domain.
> > > If
> > > > > > you're
> > > > > > > actually verifying on the GFS domain, the command looks
a
> little
> > > > > > redundant:
> > > > > > >
> > > > > > > met-9.1/bin/gen_vx_mask met-9.1/bin/gen_vx_mask \
> > > > > > > > gfs_4_20200820_0000_012.grb2 \
> > > > > > > > gfs_4_20200820_0000_012.grb2 \
> > > > > > > > gfs_4_WATER_mask.nc \
> > > > > > > > -type data \
> > > > > > > > -mask_field 'name="LAND"; level="L0";' \
> > > > > > > > -thresh ==0 \
> > > > > > > > -name WATER
> > > > > > >
> > > > > > >
> > > > > > > DEBUG 1: Input Grid: gfs_4_20200820_0000_012.grb2
> > > > > > > > DEBUG 1: Mask File: gfs_4_20200820_0000_012.grb2
> > > > > > > > DEBUG 2: Parsed Input Grid: LatLon (720 x 361)
> > > > > > > > DEBUG 2: Parsed Mask Grid: LatLon (720 x 361)
> > > > > > > > DEBUG 1: Output File: gfs_4_WATER_mask.nc
> > > > > > >
> > > > > > >
> > > > > > > I defined the threshold as "==0" because I want the
WATER
> points,
> > > not
> > > > > the
> > > > > > > LAND points. And I used the "-name WATER" command line
option
> to
> > > > name
> > > > > > the
> > > > > > > output NetCDF variable. This will show up in the
VX_MASK
> output
> > > > column
> > > > > > of
> > > > > > > Point-Stat.
> > > > > > >
> > > > > > > And finally, you would then setup your Point-Stat config
file
> > like
> > > > > this:
> > > > > > >
> > > > > > > mask = {
> > > > > > > grid = [];
> > > > > > > poly = [ "/path/to/gfs_4_WATER_mask.nc" ];
> > > > > > > sid = [];
> > > > > > > }
> > > > > > >
> > > > > > > I hope that helps!
> > > > > > >
> > > > > > > Julie
> > > > > > >
> > > > > > > On Tue, Aug 25, 2020 at 1:14 PM Tatiana Gonzalez - NOAA
Federal
> > via
> > > > RT
> > > > > <
> > > > > > > met_help at ucar.edu> wrote:
> > > > > > >
> > > > > > > >
> > > > > > > > Tue Aug 25 13:14:03 2020: Request 96443 was acted
upon.
> > > > > > > > Transaction: Ticket created by
tatiana.gonzalez at noaa.gov
> > > > > > > > Queue: met_help
> > > > > > > > Subject: mask file for Metplus
> > > > > > > > Owner: Nobody
> > > > > > > > Requestors: tatiana.gonzalez at noaa.gov
> > > > > > > > Status: new
> > > > > > > > Ticket <URL:
> > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
> > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > Hello,
> > > > > > > > How do I create a msk file to be use in Metplus (point
stat)?
> > to
> > > > > limit
> > > > > > > my
> > > > > > > > analysis
> > > > > > > >
> > > > > > > > Thanks
> > > > > > > >
> > > > > > > >
> > > > > > > > Tatiana
> > > > > > > > --
> > > > > > > >
> > > > > > > >
> > > > > > > > *Tatiana D. Gonzalez, PhD.*
> > > > > > > > *Analysis & Nowcast Branch*
> > > > > > > > *Analysis and Mission Support Division*
> > > > > > > >
> > > > > > > > *Analyze, Forecast and Support Office*
> > > > > > > >
> > > > > > > > *National Weather Service**SSMC2, Room 14115*
> > > > > > > > *Silver Spring, Maryland*
> > > > > > > >
> > > > > > > >
> > > > > > > > *Ph. (301) 427-9500 *
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > > --
> > > > > > > Julie Prestopnik (she/her/hers)
> > > > > > > Software Engineer
> > > > > > > National Center for Atmospheric Research
> > > > > > > Research Applications Laboratory
> > > > > > > Phone: 303.497.8399
> > > > > > > Email: jpresto at ucar.edu
> > > > > > >
> > > > > > > My working day may not be your working day. Please do
not feel
> > > > obliged
> > > > > > to
> > > > > > > reply to this email outside of your normal working
hours.
> > > > > > >
> > > > > > >
> > > > > >
> > > > > > --
> > > > > >
> > > > > >
> > > > > > *Tatiana D. Gonzalez, PhD.*
> > > > > > *Analysis & Nowcast Branch*
> > > > > > *Analysis and Mission Support Division*
> > > > > >
> > > > > > *Analyze, Forecast and Support Office*
> > > > > >
> > > > > > *National Weather Service**SSMC2, Room 14115*
> > > > > > *Silver Spring, Maryland*
> > > > > >
> > > > > >
> > > > > > *Ph. (301) 427-9500 *
> > > > > >
> > > > > >
> > > > >
> > > > > --
> > > > > Julie Prestopnik (she/her/hers)
> > > > > Software Engineer
> > > > > National Center for Atmospheric Research
> > > > > Research Applications Laboratory
> > > > > Phone: 303.497.8399
> > > > > Email: jpresto at ucar.edu
> > > > >
> > > > > My working day may not be your working day. Please do not
feel
> > obliged
> > > > to
> > > > > reply to this email outside of your normal working hours.
> > > > >
> > > > >
> > > >
> > > > --
> > > >
> > > >
> > > > *Tatiana D. Gonzalez, PhD.*
> > > > *Analysis & Nowcast Branch*
> > > > *Analysis and Mission Support Division*
> > > >
> > > > *Analyze, Forecast and Support Office*
> > > >
> > > > *National Weather Service**SSMC2, Room 14115*
> > > > *Silver Spring, Maryland*
> > > >
> > > >
> > > > *Ph. (301) 427-9500 *
> > > >
> > > >
> > >
> > >
> >
> > --
> >
> >
> > *Tatiana D. Gonzalez, PhD.*
> > *Analysis & Nowcast Branch*
> > *Analysis and Mission Support Division*
> >
> > *Analyze, Forecast and Support Office*
> >
> > *National Weather Service**SSMC2, Room 14115*
> > *Silver Spring, Maryland*
> >
> >
> > *Ph. (301) 427-9500 *
> >
> >
>
> --
> Julie Prestopnik (she/her/hers)
> Software Engineer
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
> My working day may not be your working day. Please do not feel
obliged to
> reply to this email outside of your normal working hours.
>
>
--
*Tatiana D. Gonzalez, PhD.*
*Analysis & Nowcast Branch*
*Analysis and Mission Support Division*
*Analyze, Forecast and Support Office*
*National Weather Service**SSMC2, Room 14115*
*Silver Spring, Maryland*
*Ph. (301) 427-9500 *
------------------------------------------------
Subject: mask file for Metplus
From: Tatiana Gonzalez - NOAA Federal
Time: Thu Sep 10 11:19:16 2020
Julie,
I think it was my bad. the coordinates were not correctly converted.
DId
you have a chance to look at the grid? Just want to make sure it is
right
anyways. It should be centered in MIAMI. I don't have a way to plot it
right now not because WCOSS is down.
I apologize for the inconvenience.
Tatiana
On Thu, Sep 10, 2020 at 11:45 AM Tatiana Gonzalez - NOAA Federal <
tatiana.gonzalez at noaa.gov> wrote:
> Try again. I just gave you permission
>
> On Thu, Sep 10, 2020 at 11:36 AM Julie Prestopnik via RT <
> met_help at ucar.edu> wrote:
>
>> Hi Tatiana.
>>
>> I sent a request for access through google drive as I did not have
>> permissions on the file.
>>
>> Julie
>>
>> On Thu, Sep 10, 2020 at 8:58 AM Tatiana Gonzalez - NOAA Federal via
RT <
>> met_help at ucar.edu> wrote:
>>
>> >
>> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
>> >
>> > Jullie,
>> > Thank you for reviewing the file
>> > case_2.nc
>> > <
>> >
>> https://drive.google.com/file/d/1tWThlec9seGOruztsodx73wHAD3OU-
gV/view?usp=drive_web
>> > >
>> >
>> >
>> > Tatiana
>> >
>> > On Thu, Sep 10, 2020 at 10:38 AM Julie Prestopnik via RT <
>> > met_help at ucar.edu>
>> > wrote:
>> >
>> > > Hi Tatiana.
>> > >
>> > > Could you please either attach file.nc or upload file.nc to our
FTP
>> > > server as
>> > > described under "How to Send Us Data" on the following page?
>> > >
>> > >
>> >
>> https://dtcenter.org/community-code/model-evaluation-tools-met/met-
help-desk
>> > >
>> > >
>> > > I'll take a look.
>> > >
>> > > Thanks!
>> > >
>> > > Julie
>> > >
>> > >
>> > > On Thu, Sep 10, 2020, 7:10 AM Tatiana Gonzalez - NOAA Federal
via RT <
>> > > met_help at ucar.edu> wrote:
>> > >
>> > > >
>> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
>
>> > > >
>> > > > Hello Jullie,
>> > > > I created a file with 4 coordinates and ran the polyline
through the
>> > > > gen_vx_mask tool. Then on the config file I use the new
"file.nc"
>> > > > But it looks like none of the stations that are within my
grid were
>> > > > included. Is there an order to the coordinates? My file was
in the
>> > > > following order (south east, north east, north west, south
west)
>> > > >
>> > > > MIAMI
>> > > >
>> > > > 25.653243 -79.832631
>> > > >
>> > > > 28.073608 -79.832631
>> > > >
>> > > > 28.073608 -80.938611
>> > > >
>> > > > 25.602028 -80.938611
>> > > >
>> > > >
>> > > > The grib file I used contained the whole CONUS. Is there
something
>> I'm
>> > > > doing wrong?
>> > > >
>> > > > Thanks
>> > > >
>> > > > Tatiana
>> > > >
>> > > >
>> > > >
>> > > > On Wed, Aug 26, 2020 at 2:54 PM Julie Prestopnik via RT <
>> > > met_help at ucar.edu
>> > > > >
>> > > > wrote:
>> > > >
>> > > > > Hi Tatiana.
>> > > > >
>> > > > > To create a polyline mask, you simply create an ASCII file
>> > containing a
>> > > > > list of lat/lons that define the region of interest. For
example,
>> > > take a
>> > > > > look at met-9.1/data/poly/WEST.poly. The first line of the
file
>> > > > specifies
>> > > > > the name for the polyline. And the remaining lines contain
>> lat/lon
>> > > > values
>> > > > > for the boundary points. The last point is connected back
up to
>> the
>> > > > first
>> > > > > point to complete the area. You need a minimum of 3 points
to
>> define
>> > a
>> > > > > polyline, and there is no real maximum. If your polyline
>> contains a
>> > > > large
>> > > > > number of points (a few hundred) or you have many grid-
points, the
>> > > > process
>> > > > > of determining which grid points are inside/outside the
polyline
>> can
>> > be
>> > > > > slow. In that case, I'd suggest running your polyline
through the
>> > > > > gen_vx_mask tool.
>> > > > >
>> > > > > Julie
>> > > > >
>> > > > > On Wed, Aug 26, 2020 at 12:03 PM Tatiana Gonzalez - NOAA
Federal
>> via
>> > > RT <
>> > > > > met_help at ucar.edu> wrote:
>> > > > >
>> > > > > >
>> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
>> > > > > >
>> > > > > > Hi Julie,
>> > > > > > Thanks for the instructions. I have the RTMA files in
grib2
>> format
>> > > for
>> > > > > the
>> > > > > > whole CONUS. But what I want is only use the data
contained in
>> this
>> > > > box:
>> > > > > >
>> > > > > > South East 45.125 -121.975
>> > > > > >
>> > > > > > North East 46.188 -121.975
>> > > > > >
>> > > > > > North West 46.188 -123.625
>> > > > > >
>> > > > > > South West 45.125 -123.625
>> > > > > >
>> > > > > > I see there is the Polyline (poly) that might work for
me. Do
>> you
>> > > have
>> > > > an
>> > > > > > example for that case or the ASCII file that is needed.
>> > > > > >
>> > > > > > Thank you
>> > > > > >
>> > > > > > Tatiana
>> > > > > >
>> > > > > > On Tue, Aug 25, 2020 at 5:52 PM Julie Prestopnik via RT <
>> > > > > met_help at ucar.edu
>> > > > > > >
>> > > > > > wrote:
>> > > > > >
>> > > > > > > Hi Tatiana.
>> > > > > > >
>> > > > > > > I see you have a question about how to create a mask
file.
>> I'd
>> > > > > recommend
>> > > > > > > running the Gen-Vx-Mask tool to create the desired
masking
>> > region,
>> > > > and
>> > > > > > then
>> > > > > > > reference that in the Point-Stat config file. Here is
a link
>> to
>> > > our
>> > > > > > > documentation page on Gen-Vx-Mask, but I also go
through an
>> > example
>> > > > > > below:
>> > > > > > > https://dtcenter.github.io/MET/Users_Guide/masking.html
>> > > > > > >
>> > > > > > > To illustrate, I pulled a sample GFS file:
>> > > > > > >
>> > > > > > > > wget
>> > > > > > > >
>> > > > > > >
>> > > > > >
>> > > > >
>> > > >
>> > >
>> >
>> https://www.ncei.noaa.gov/data/global-forecast-system/access/grid-
004-0.5-degree/forecast/202008/20200820/gfs_4_20200820_0000_012.grb2
>> > > > > > >
>> > > > > > >
>> > > > > > > I ran wgrib2 and look for a LAND/SEA mask:
>> > > > > > >
>> > > > > > > > wgrib2 gfs_4_20200820_0000_012.grb2 | egrep -i LAND
>> > > > > > > >
>> > > > > > > which returned:
>> > > > > > >
>> > > > > > > > 571:98325645:d=2020082000:LAND:surface:12 hour fcst:
>> > > > > > >
>> > > > > > >
>> > > > > > > I ran plot_data_plane on it (see attached image):
>> > > > > > >
>> > > > > > > > met-9.1/bin/plot_data_plane
gfs_4_20200820_0000_012.grb2
>> > > > > gfs_4_LAND.ps
>> > > > > > > > 'name="LAND"; level="L0";'
>> > > > > > > >
>> > > > > > >
>> > > > > > >
>> > > > > > > > DEBUG 1: Opening data file:
gfs_4_20200820_0000_012.grb2
>> > > > > > > > DEBUG 1: Creating postscript file: gfs_4_LAND.ps
>> > > > > > > >
>> > > > > > >
>> > > > > > > I ran gen_vx_mask to create a mask for this
verification
>> domain.
>> > > If
>> > > > > > you're
>> > > > > > > actually verifying on the GFS domain, the command looks
a
>> little
>> > > > > > redundant:
>> > > > > > >
>> > > > > > > met-9.1/bin/gen_vx_mask met-9.1/bin/gen_vx_mask \
>> > > > > > > > gfs_4_20200820_0000_012.grb2 \
>> > > > > > > > gfs_4_20200820_0000_012.grb2 \
>> > > > > > > > gfs_4_WATER_mask.nc \
>> > > > > > > > -type data \
>> > > > > > > > -mask_field 'name="LAND"; level="L0";' \
>> > > > > > > > -thresh ==0 \
>> > > > > > > > -name WATER
>> > > > > > >
>> > > > > > >
>> > > > > > > DEBUG 1: Input Grid: gfs_4_20200820_0000_012.grb2
>> > > > > > > > DEBUG 1: Mask File: gfs_4_20200820_0000_012.grb2
>> > > > > > > > DEBUG 2: Parsed Input Grid: LatLon (720 x 361)
>> > > > > > > > DEBUG 2: Parsed Mask Grid: LatLon (720 x 361)
>> > > > > > > > DEBUG 1: Output File: gfs_4_WATER_mask.nc
>> > > > > > >
>> > > > > > >
>> > > > > > > I defined the threshold as "==0" because I want the
WATER
>> points,
>> > > not
>> > > > > the
>> > > > > > > LAND points. And I used the "-name WATER" command line
>> option to
>> > > > name
>> > > > > > the
>> > > > > > > output NetCDF variable. This will show up in the
VX_MASK
>> output
>> > > > column
>> > > > > > of
>> > > > > > > Point-Stat.
>> > > > > > >
>> > > > > > > And finally, you would then setup your Point-Stat
config file
>> > like
>> > > > > this:
>> > > > > > >
>> > > > > > > mask = {
>> > > > > > > grid = [];
>> > > > > > > poly = [ "/path/to/gfs_4_WATER_mask.nc" ];
>> > > > > > > sid = [];
>> > > > > > > }
>> > > > > > >
>> > > > > > > I hope that helps!
>> > > > > > >
>> > > > > > > Julie
>> > > > > > >
>> > > > > > > On Tue, Aug 25, 2020 at 1:14 PM Tatiana Gonzalez - NOAA
>> Federal
>> > via
>> > > > RT
>> > > > > <
>> > > > > > > met_help at ucar.edu> wrote:
>> > > > > > >
>> > > > > > > >
>> > > > > > > > Tue Aug 25 13:14:03 2020: Request 96443 was acted
upon.
>> > > > > > > > Transaction: Ticket created by
tatiana.gonzalez at noaa.gov
>> > > > > > > > Queue: met_help
>> > > > > > > > Subject: mask file for Metplus
>> > > > > > > > Owner: Nobody
>> > > > > > > > Requestors: tatiana.gonzalez at noaa.gov
>> > > > > > > > Status: new
>> > > > > > > > Ticket <URL:
>> > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
>> > > > > > >
>> > > > > > > >
>> > > > > > > >
>> > > > > > > > Hello,
>> > > > > > > > How do I create a msk file to be use in Metplus
(point
>> stat)?
>> > to
>> > > > > limit
>> > > > > > > my
>> > > > > > > > analysis
>> > > > > > > >
>> > > > > > > > Thanks
>> > > > > > > >
>> > > > > > > >
>> > > > > > > > Tatiana
>> > > > > > > > --
>> > > > > > > >
>> > > > > > > >
>> > > > > > > > *Tatiana D. Gonzalez, PhD.*
>> > > > > > > > *Analysis & Nowcast Branch*
>> > > > > > > > *Analysis and Mission Support Division*
>> > > > > > > >
>> > > > > > > > *Analyze, Forecast and Support Office*
>> > > > > > > >
>> > > > > > > > *National Weather Service**SSMC2, Room 14115*
>> > > > > > > > *Silver Spring, Maryland*
>> > > > > > > >
>> > > > > > > >
>> > > > > > > > *Ph. (301) 427-9500 *
>> > > > > > > >
>> > > > > > > >
>> > > > > > >
>> > > > > > > --
>> > > > > > > Julie Prestopnik (she/her/hers)
>> > > > > > > Software Engineer
>> > > > > > > National Center for Atmospheric Research
>> > > > > > > Research Applications Laboratory
>> > > > > > > Phone: 303.497.8399
>> > > > > > > Email: jpresto at ucar.edu
>> > > > > > >
>> > > > > > > My working day may not be your working day. Please do
not
>> feel
>> > > > obliged
>> > > > > > to
>> > > > > > > reply to this email outside of your normal working
hours.
>> > > > > > >
>> > > > > > >
>> > > > > >
>> > > > > > --
>> > > > > >
>> > > > > >
>> > > > > > *Tatiana D. Gonzalez, PhD.*
>> > > > > > *Analysis & Nowcast Branch*
>> > > > > > *Analysis and Mission Support Division*
>> > > > > >
>> > > > > > *Analyze, Forecast and Support Office*
>> > > > > >
>> > > > > > *National Weather Service**SSMC2, Room 14115*
>> > > > > > *Silver Spring, Maryland*
>> > > > > >
>> > > > > >
>> > > > > > *Ph. (301) 427-9500 *
>> > > > > >
>> > > > > >
>> > > > >
>> > > > > --
>> > > > > Julie Prestopnik (she/her/hers)
>> > > > > Software Engineer
>> > > > > National Center for Atmospheric Research
>> > > > > Research Applications Laboratory
>> > > > > Phone: 303.497.8399
>> > > > > Email: jpresto at ucar.edu
>> > > > >
>> > > > > My working day may not be your working day. Please do not
feel
>> > obliged
>> > > > to
>> > > > > reply to this email outside of your normal working hours.
>> > > > >
>> > > > >
>> > > >
>> > > > --
>> > > >
>> > > >
>> > > > *Tatiana D. Gonzalez, PhD.*
>> > > > *Analysis & Nowcast Branch*
>> > > > *Analysis and Mission Support Division*
>> > > >
>> > > > *Analyze, Forecast and Support Office*
>> > > >
>> > > > *National Weather Service**SSMC2, Room 14115*
>> > > > *Silver Spring, Maryland*
>> > > >
>> > > >
>> > > > *Ph. (301) 427-9500 *
>> > > >
>> > > >
>> > >
>> > >
>> >
>> > --
>> >
>> >
>> > *Tatiana D. Gonzalez, PhD.*
>> > *Analysis & Nowcast Branch*
>> > *Analysis and Mission Support Division*
>> >
>> > *Analyze, Forecast and Support Office*
>> >
>> > *National Weather Service**SSMC2, Room 14115*
>> > *Silver Spring, Maryland*
>> >
>> >
>> > *Ph. (301) 427-9500 *
>> >
>> >
>>
>> --
>> Julie Prestopnik (she/her/hers)
>> Software Engineer
>> National Center for Atmospheric Research
>> Research Applications Laboratory
>> Phone: 303.497.8399
>> Email: jpresto at ucar.edu
>>
>> My working day may not be your working day. Please do not feel
obliged to
>> reply to this email outside of your normal working hours.
>>
>>
>
> --
>
>
> *Tatiana D. Gonzalez, PhD.*
> *Analysis & Nowcast Branch*
> *Analysis and Mission Support Division*
>
> *Analyze, Forecast and Support Office*
>
> *National Weather Service**SSMC2, Room 14115*
> *Silver Spring, Maryland*
>
>
> *Ph. (301) 427-9500 *
>
>
--
*Tatiana D. Gonzalez, PhD.*
*Analysis & Nowcast Branch*
*Analysis and Mission Support Division*
*Analyze, Forecast and Support Office*
*National Weather Service**SSMC2, Room 14115*
*Silver Spring, Maryland*
*Ph. (301) 427-9500 *
------------------------------------------------
Subject: mask file for Metplus
From: Julie Prestopnik
Time: Thu Sep 10 11:41:17 2020
I had not yet had a chance, because I was in a meeting, but I just
took a
look. I ran it through plot_data_plane:
plot_data_plane case_2.nc case.ps 'name="MIAMI"; level="(*,*)";'
I am attaching the resulting case.ps file. Is that where you expected
it
to be?
Julie
On Thu, Sep 10, 2020 at 11:19 AM Tatiana Gonzalez - NOAA Federal via
RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
>
> Julie,
> I think it was my bad. the coordinates were not correctly converted.
DId
> you have a chance to look at the grid? Just want to make sure it is
right
> anyways. It should be centered in MIAMI. I don't have a way to plot
it
> right now not because WCOSS is down.
> I apologize for the inconvenience.
>
> Tatiana
>
>
> On Thu, Sep 10, 2020 at 11:45 AM Tatiana Gonzalez - NOAA Federal <
> tatiana.gonzalez at noaa.gov> wrote:
>
> > Try again. I just gave you permission
> >
> > On Thu, Sep 10, 2020 at 11:36 AM Julie Prestopnik via RT <
> > met_help at ucar.edu> wrote:
> >
> >> Hi Tatiana.
> >>
> >> I sent a request for access through google drive as I did not
have
> >> permissions on the file.
> >>
> >> Julie
> >>
> >> On Thu, Sep 10, 2020 at 8:58 AM Tatiana Gonzalez - NOAA Federal
via RT <
> >> met_help at ucar.edu> wrote:
> >>
> >> >
> >> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
> >> >
> >> > Jullie,
> >> > Thank you for reviewing the file
> >> > case_2.nc
> >> > <
> >> >
> >>
> https://drive.google.com/file/d/1tWThlec9seGOruztsodx73wHAD3OU-
gV/view?usp=drive_web
> >> > >
> >> >
> >> >
> >> > Tatiana
> >> >
> >> > On Thu, Sep 10, 2020 at 10:38 AM Julie Prestopnik via RT <
> >> > met_help at ucar.edu>
> >> > wrote:
> >> >
> >> > > Hi Tatiana.
> >> > >
> >> > > Could you please either attach file.nc or upload file.nc to
our FTP
> >> > > server as
> >> > > described under "How to Send Us Data" on the following page?
> >> > >
> >> > >
> >> >
> >>
> https://dtcenter.org/community-code/model-evaluation-tools-met/met-
help-desk
> >> > >
> >> > >
> >> > > I'll take a look.
> >> > >
> >> > > Thanks!
> >> > >
> >> > > Julie
> >> > >
> >> > >
> >> > > On Thu, Sep 10, 2020, 7:10 AM Tatiana Gonzalez - NOAA Federal
via
> RT <
> >> > > met_help at ucar.edu> wrote:
> >> > >
> >> > > >
> >> > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
> >> > > >
> >> > > > Hello Jullie,
> >> > > > I created a file with 4 coordinates and ran the polyline
through
> the
> >> > > > gen_vx_mask tool. Then on the config file I use the new
"file.nc"
> >> > > > But it looks like none of the stations that are within my
grid
> were
> >> > > > included. Is there an order to the coordinates? My file was
in the
> >> > > > following order (south east, north east, north west, south
west)
> >> > > >
> >> > > > MIAMI
> >> > > >
> >> > > > 25.653243 -79.832631
> >> > > >
> >> > > > 28.073608 -79.832631
> >> > > >
> >> > > > 28.073608 -80.938611
> >> > > >
> >> > > > 25.602028 -80.938611
> >> > > >
> >> > > >
> >> > > > The grib file I used contained the whole CONUS. Is there
something
> >> I'm
> >> > > > doing wrong?
> >> > > >
> >> > > > Thanks
> >> > > >
> >> > > > Tatiana
> >> > > >
> >> > > >
> >> > > >
> >> > > > On Wed, Aug 26, 2020 at 2:54 PM Julie Prestopnik via RT <
> >> > > met_help at ucar.edu
> >> > > > >
> >> > > > wrote:
> >> > > >
> >> > > > > Hi Tatiana.
> >> > > > >
> >> > > > > To create a polyline mask, you simply create an ASCII
file
> >> > containing a
> >> > > > > list of lat/lons that define the region of interest. For
> example,
> >> > > take a
> >> > > > > look at met-9.1/data/poly/WEST.poly. The first line of
the file
> >> > > > specifies
> >> > > > > the name for the polyline. And the remaining lines
contain
> >> lat/lon
> >> > > > values
> >> > > > > for the boundary points. The last point is connected
back up to
> >> the
> >> > > > first
> >> > > > > point to complete the area. You need a minimum of 3
points to
> >> define
> >> > a
> >> > > > > polyline, and there is no real maximum. If your polyline
> >> contains a
> >> > > > large
> >> > > > > number of points (a few hundred) or you have many grid-
points,
> the
> >> > > > process
> >> > > > > of determining which grid points are inside/outside the
polyline
> >> can
> >> > be
> >> > > > > slow. In that case, I'd suggest running your polyline
through
> the
> >> > > > > gen_vx_mask tool.
> >> > > > >
> >> > > > > Julie
> >> > > > >
> >> > > > > On Wed, Aug 26, 2020 at 12:03 PM Tatiana Gonzalez - NOAA
Federal
> >> via
> >> > > RT <
> >> > > > > met_help at ucar.edu> wrote:
> >> > > > >
> >> > > > > >
> >> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
> >
> >> > > > > >
> >> > > > > > Hi Julie,
> >> > > > > > Thanks for the instructions. I have the RTMA files in
grib2
> >> format
> >> > > for
> >> > > > > the
> >> > > > > > whole CONUS. But what I want is only use the data
contained in
> >> this
> >> > > > box:
> >> > > > > >
> >> > > > > > South East 45.125 -121.975
> >> > > > > >
> >> > > > > > North East 46.188 -121.975
> >> > > > > >
> >> > > > > > North West 46.188 -123.625
> >> > > > > >
> >> > > > > > South West 45.125 -123.625
> >> > > > > >
> >> > > > > > I see there is the Polyline (poly) that might work for
me. Do
> >> you
> >> > > have
> >> > > > an
> >> > > > > > example for that case or the ASCII file that is needed.
> >> > > > > >
> >> > > > > > Thank you
> >> > > > > >
> >> > > > > > Tatiana
> >> > > > > >
> >> > > > > > On Tue, Aug 25, 2020 at 5:52 PM Julie Prestopnik via RT
<
> >> > > > > met_help at ucar.edu
> >> > > > > > >
> >> > > > > > wrote:
> >> > > > > >
> >> > > > > > > Hi Tatiana.
> >> > > > > > >
> >> > > > > > > I see you have a question about how to create a mask
file.
> >> I'd
> >> > > > > recommend
> >> > > > > > > running the Gen-Vx-Mask tool to create the desired
masking
> >> > region,
> >> > > > and
> >> > > > > > then
> >> > > > > > > reference that in the Point-Stat config file. Here
is a
> link
> >> to
> >> > > our
> >> > > > > > > documentation page on Gen-Vx-Mask, but I also go
through an
> >> > example
> >> > > > > > below:
> >> > > > > > >
https://dtcenter.github.io/MET/Users_Guide/masking.html
> >> > > > > > >
> >> > > > > > > To illustrate, I pulled a sample GFS file:
> >> > > > > > >
> >> > > > > > > > wget
> >> > > > > > > >
> >> > > > > > >
> >> > > > > >
> >> > > > >
> >> > > >
> >> > >
> >> >
> >>
> https://www.ncei.noaa.gov/data/global-forecast-system/access/grid-
004-0.5-degree/forecast/202008/20200820/gfs_4_20200820_0000_012.grb2
> >> > > > > > >
> >> > > > > > >
> >> > > > > > > I ran wgrib2 and look for a LAND/SEA mask:
> >> > > > > > >
> >> > > > > > > > wgrib2 gfs_4_20200820_0000_012.grb2 | egrep -i
LAND
> >> > > > > > > >
> >> > > > > > > which returned:
> >> > > > > > >
> >> > > > > > > > 571:98325645:d=2020082000:LAND:surface:12 hour
fcst:
> >> > > > > > >
> >> > > > > > >
> >> > > > > > > I ran plot_data_plane on it (see attached image):
> >> > > > > > >
> >> > > > > > > > met-9.1/bin/plot_data_plane
gfs_4_20200820_0000_012.grb2
> >> > > > > gfs_4_LAND.ps
> >> > > > > > > > 'name="LAND"; level="L0";'
> >> > > > > > > >
> >> > > > > > >
> >> > > > > > >
> >> > > > > > > > DEBUG 1: Opening data file:
gfs_4_20200820_0000_012.grb2
> >> > > > > > > > DEBUG 1: Creating postscript file: gfs_4_LAND.ps
> >> > > > > > > >
> >> > > > > > >
> >> > > > > > > I ran gen_vx_mask to create a mask for this
verification
> >> domain.
> >> > > If
> >> > > > > > you're
> >> > > > > > > actually verifying on the GFS domain, the command
looks a
> >> little
> >> > > > > > redundant:
> >> > > > > > >
> >> > > > > > > met-9.1/bin/gen_vx_mask met-9.1/bin/gen_vx_mask \
> >> > > > > > > > gfs_4_20200820_0000_012.grb2 \
> >> > > > > > > > gfs_4_20200820_0000_012.grb2 \
> >> > > > > > > > gfs_4_WATER_mask.nc \
> >> > > > > > > > -type data \
> >> > > > > > > > -mask_field 'name="LAND"; level="L0";' \
> >> > > > > > > > -thresh ==0 \
> >> > > > > > > > -name WATER
> >> > > > > > >
> >> > > > > > >
> >> > > > > > > DEBUG 1: Input Grid: gfs_4_20200820_0000_012.grb2
> >> > > > > > > > DEBUG 1: Mask File: gfs_4_20200820_0000_012.grb2
> >> > > > > > > > DEBUG 2: Parsed Input Grid: LatLon (720 x 361)
> >> > > > > > > > DEBUG 2: Parsed Mask Grid: LatLon (720 x 361)
> >> > > > > > > > DEBUG 1: Output File: gfs_4_WATER_mask.nc
> >> > > > > > >
> >> > > > > > >
> >> > > > > > > I defined the threshold as "==0" because I want the
WATER
> >> points,
> >> > > not
> >> > > > > the
> >> > > > > > > LAND points. And I used the "-name WATER" command
line
> >> option to
> >> > > > name
> >> > > > > > the
> >> > > > > > > output NetCDF variable. This will show up in the
VX_MASK
> >> output
> >> > > > column
> >> > > > > > of
> >> > > > > > > Point-Stat.
> >> > > > > > >
> >> > > > > > > And finally, you would then setup your Point-Stat
config
> file
> >> > like
> >> > > > > this:
> >> > > > > > >
> >> > > > > > > mask = {
> >> > > > > > > grid = [];
> >> > > > > > > poly = [ "/path/to/gfs_4_WATER_mask.nc" ];
> >> > > > > > > sid = [];
> >> > > > > > > }
> >> > > > > > >
> >> > > > > > > I hope that helps!
> >> > > > > > >
> >> > > > > > > Julie
> >> > > > > > >
> >> > > > > > > On Tue, Aug 25, 2020 at 1:14 PM Tatiana Gonzalez -
NOAA
> >> Federal
> >> > via
> >> > > > RT
> >> > > > > <
> >> > > > > > > met_help at ucar.edu> wrote:
> >> > > > > > >
> >> > > > > > > >
> >> > > > > > > > Tue Aug 25 13:14:03 2020: Request 96443 was acted
upon.
> >> > > > > > > > Transaction: Ticket created by
tatiana.gonzalez at noaa.gov
> >> > > > > > > > Queue: met_help
> >> > > > > > > > Subject: mask file for Metplus
> >> > > > > > > > Owner: Nobody
> >> > > > > > > > Requestors: tatiana.gonzalez at noaa.gov
> >> > > > > > > > Status: new
> >> > > > > > > > Ticket <URL:
> >> > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
> >> > > > > > >
> >> > > > > > > >
> >> > > > > > > >
> >> > > > > > > > Hello,
> >> > > > > > > > How do I create a msk file to be use in Metplus
(point
> >> stat)?
> >> > to
> >> > > > > limit
> >> > > > > > > my
> >> > > > > > > > analysis
> >> > > > > > > >
> >> > > > > > > > Thanks
> >> > > > > > > >
> >> > > > > > > >
> >> > > > > > > > Tatiana
> >> > > > > > > > --
> >> > > > > > > >
> >> > > > > > > >
> >> > > > > > > > *Tatiana D. Gonzalez, PhD.*
> >> > > > > > > > *Analysis & Nowcast Branch*
> >> > > > > > > > *Analysis and Mission Support Division*
> >> > > > > > > >
> >> > > > > > > > *Analyze, Forecast and Support Office*
> >> > > > > > > >
> >> > > > > > > > *National Weather Service**SSMC2, Room 14115*
> >> > > > > > > > *Silver Spring, Maryland*
> >> > > > > > > >
> >> > > > > > > >
> >> > > > > > > > *Ph. (301) 427-9500 *
> >> > > > > > > >
> >> > > > > > > >
> >> > > > > > >
> >> > > > > > > --
> >> > > > > > > Julie Prestopnik (she/her/hers)
> >> > > > > > > Software Engineer
> >> > > > > > > National Center for Atmospheric Research
> >> > > > > > > Research Applications Laboratory
> >> > > > > > > Phone: 303.497.8399
> >> > > > > > > Email: jpresto at ucar.edu
> >> > > > > > >
> >> > > > > > > My working day may not be your working day. Please
do not
> >> feel
> >> > > > obliged
> >> > > > > > to
> >> > > > > > > reply to this email outside of your normal working
hours.
> >> > > > > > >
> >> > > > > > >
> >> > > > > >
> >> > > > > > --
> >> > > > > >
> >> > > > > >
> >> > > > > > *Tatiana D. Gonzalez, PhD.*
> >> > > > > > *Analysis & Nowcast Branch*
> >> > > > > > *Analysis and Mission Support Division*
> >> > > > > >
> >> > > > > > *Analyze, Forecast and Support Office*
> >> > > > > >
> >> > > > > > *National Weather Service**SSMC2, Room 14115*
> >> > > > > > *Silver Spring, Maryland*
> >> > > > > >
> >> > > > > >
> >> > > > > > *Ph. (301) 427-9500 *
> >> > > > > >
> >> > > > > >
> >> > > > >
> >> > > > > --
> >> > > > > Julie Prestopnik (she/her/hers)
> >> > > > > Software Engineer
> >> > > > > National Center for Atmospheric Research
> >> > > > > Research Applications Laboratory
> >> > > > > Phone: 303.497.8399
> >> > > > > Email: jpresto at ucar.edu
> >> > > > >
> >> > > > > My working day may not be your working day. Please do
not feel
> >> > obliged
> >> > > > to
> >> > > > > reply to this email outside of your normal working hours.
> >> > > > >
> >> > > > >
> >> > > >
> >> > > > --
> >> > > >
> >> > > >
> >> > > > *Tatiana D. Gonzalez, PhD.*
> >> > > > *Analysis & Nowcast Branch*
> >> > > > *Analysis and Mission Support Division*
> >> > > >
> >> > > > *Analyze, Forecast and Support Office*
> >> > > >
> >> > > > *National Weather Service**SSMC2, Room 14115*
> >> > > > *Silver Spring, Maryland*
> >> > > >
> >> > > >
> >> > > > *Ph. (301) 427-9500 *
> >> > > >
> >> > > >
> >> > >
> >> > >
> >> >
> >> > --
> >> >
> >> >
> >> > *Tatiana D. Gonzalez, PhD.*
> >> > *Analysis & Nowcast Branch*
> >> > *Analysis and Mission Support Division*
> >> >
> >> > *Analyze, Forecast and Support Office*
> >> >
> >> > *National Weather Service**SSMC2, Room 14115*
> >> > *Silver Spring, Maryland*
> >> >
> >> >
> >> > *Ph. (301) 427-9500 *
> >> >
> >> >
> >>
> >> --
> >> Julie Prestopnik (she/her/hers)
> >> Software Engineer
> >> National Center for Atmospheric Research
> >> Research Applications Laboratory
> >> Phone: 303.497.8399
> >> Email: jpresto at ucar.edu
> >>
> >> My working day may not be your working day. Please do not feel
obliged
> to
> >> reply to this email outside of your normal working hours.
> >>
> >>
> >
> > --
> >
> >
> > *Tatiana D. Gonzalez, PhD.*
> > *Analysis & Nowcast Branch*
> > *Analysis and Mission Support Division*
> >
> > *Analyze, Forecast and Support Office*
> >
> > *National Weather Service**SSMC2, Room 14115*
> > *Silver Spring, Maryland*
> >
> >
> > *Ph. (301) 427-9500 *
> >
> >
>
> --
>
>
> *Tatiana D. Gonzalez, PhD.*
> *Analysis & Nowcast Branch*
> *Analysis and Mission Support Division*
>
> *Analyze, Forecast and Support Office*
>
> *National Weather Service**SSMC2, Room 14115*
> *Silver Spring, Maryland*
>
>
> *Ph. (301) 427-9500 *
>
>
--
Julie Prestopnik (she/her/hers)
Software Engineer
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
My working day may not be your working day. Please do not feel
obliged to
reply to this email outside of your normal working hours.
------------------------------------------------
Subject: mask file for Metplus
From: Tatiana Gonzalez - NOAA Federal
Time: Thu Sep 10 12:27:03 2020
Yes, thank you! That seems right to me. Thanks a lot for your help
On Thu, Sep 10, 2020 at 1:41 PM Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> I had not yet had a chance, because I was in a meeting, but I just
took a
> look. I ran it through plot_data_plane:
>
> plot_data_plane case_2.nc case.ps 'name="MIAMI"; level="(*,*)";'
>
>
> I am attaching the resulting case.ps file. Is that where you
expected it
> to be?
>
> Julie
>
>
>
> On Thu, Sep 10, 2020 at 11:19 AM Tatiana Gonzalez - NOAA Federal via
RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
> >
> > Julie,
> > I think it was my bad. the coordinates were not correctly
converted. DId
> > you have a chance to look at the grid? Just want to make sure it
is right
> > anyways. It should be centered in MIAMI. I don't have a way to
plot it
> > right now not because WCOSS is down.
> > I apologize for the inconvenience.
> >
> > Tatiana
> >
> >
> > On Thu, Sep 10, 2020 at 11:45 AM Tatiana Gonzalez - NOAA Federal <
> > tatiana.gonzalez at noaa.gov> wrote:
> >
> > > Try again. I just gave you permission
> > >
> > > On Thu, Sep 10, 2020 at 11:36 AM Julie Prestopnik via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > >> Hi Tatiana.
> > >>
> > >> I sent a request for access through google drive as I did not
have
> > >> permissions on the file.
> > >>
> > >> Julie
> > >>
> > >> On Thu, Sep 10, 2020 at 8:58 AM Tatiana Gonzalez - NOAA Federal
via
> RT <
> > >> met_help at ucar.edu> wrote:
> > >>
> > >> >
> > >> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
>
> > >> >
> > >> > Jullie,
> > >> > Thank you for reviewing the file
> > >> > case_2.nc
> > >> > <
> > >> >
> > >>
> >
> https://drive.google.com/file/d/1tWThlec9seGOruztsodx73wHAD3OU-
gV/view?usp=drive_web
> > >> > >
> > >> >
> > >> >
> > >> > Tatiana
> > >> >
> > >> > On Thu, Sep 10, 2020 at 10:38 AM Julie Prestopnik via RT <
> > >> > met_help at ucar.edu>
> > >> > wrote:
> > >> >
> > >> > > Hi Tatiana.
> > >> > >
> > >> > > Could you please either attach file.nc or upload file.nc to
our
> FTP
> > >> > > server as
> > >> > > described under "How to Send Us Data" on the following
page?
> > >> > >
> > >> > >
> > >> >
> > >>
> >
> https://dtcenter.org/community-code/model-evaluation-tools-met/met-
help-desk
> > >> > >
> > >> > >
> > >> > > I'll take a look.
> > >> > >
> > >> > > Thanks!
> > >> > >
> > >> > > Julie
> > >> > >
> > >> > >
> > >> > > On Thu, Sep 10, 2020, 7:10 AM Tatiana Gonzalez - NOAA
Federal via
> > RT <
> > >> > > met_help at ucar.edu> wrote:
> > >> > >
> > >> > > >
> > >> > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
> > >> > > >
> > >> > > > Hello Jullie,
> > >> > > > I created a file with 4 coordinates and ran the polyline
through
> > the
> > >> > > > gen_vx_mask tool. Then on the config file I use the new "
> file.nc"
> > >> > > > But it looks like none of the stations that are within my
grid
> > were
> > >> > > > included. Is there an order to the coordinates? My file
was in
> the
> > >> > > > following order (south east, north east, north west,
south west)
> > >> > > >
> > >> > > > MIAMI
> > >> > > >
> > >> > > > 25.653243 -79.832631
> > >> > > >
> > >> > > > 28.073608 -79.832631
> > >> > > >
> > >> > > > 28.073608 -80.938611
> > >> > > >
> > >> > > > 25.602028 -80.938611
> > >> > > >
> > >> > > >
> > >> > > > The grib file I used contained the whole CONUS. Is there
> something
> > >> I'm
> > >> > > > doing wrong?
> > >> > > >
> > >> > > > Thanks
> > >> > > >
> > >> > > > Tatiana
> > >> > > >
> > >> > > >
> > >> > > >
> > >> > > > On Wed, Aug 26, 2020 at 2:54 PM Julie Prestopnik via RT <
> > >> > > met_help at ucar.edu
> > >> > > > >
> > >> > > > wrote:
> > >> > > >
> > >> > > > > Hi Tatiana.
> > >> > > > >
> > >> > > > > To create a polyline mask, you simply create an ASCII
file
> > >> > containing a
> > >> > > > > list of lat/lons that define the region of interest.
For
> > example,
> > >> > > take a
> > >> > > > > look at met-9.1/data/poly/WEST.poly. The first line of
the
> file
> > >> > > > specifies
> > >> > > > > the name for the polyline. And the remaining lines
contain
> > >> lat/lon
> > >> > > > values
> > >> > > > > for the boundary points. The last point is connected
back up
> to
> > >> the
> > >> > > > first
> > >> > > > > point to complete the area. You need a minimum of 3
points to
> > >> define
> > >> > a
> > >> > > > > polyline, and there is no real maximum. If your
polyline
> > >> contains a
> > >> > > > large
> > >> > > > > number of points (a few hundred) or you have many grid-
points,
> > the
> > >> > > > process
> > >> > > > > of determining which grid points are inside/outside the
> polyline
> > >> can
> > >> > be
> > >> > > > > slow. In that case, I'd suggest running your polyline
through
> > the
> > >> > > > > gen_vx_mask tool.
> > >> > > > >
> > >> > > > > Julie
> > >> > > > >
> > >> > > > > On Wed, Aug 26, 2020 at 12:03 PM Tatiana Gonzalez -
NOAA
> Federal
> > >> via
> > >> > > RT <
> > >> > > > > met_help at ucar.edu> wrote:
> > >> > > > >
> > >> > > > > >
> > >> > > > > > <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
> > >
> > >> > > > > >
> > >> > > > > > Hi Julie,
> > >> > > > > > Thanks for the instructions. I have the RTMA files in
grib2
> > >> format
> > >> > > for
> > >> > > > > the
> > >> > > > > > whole CONUS. But what I want is only use the data
contained
> in
> > >> this
> > >> > > > box:
> > >> > > > > >
> > >> > > > > > South East 45.125 -121.975
> > >> > > > > >
> > >> > > > > > North East 46.188 -121.975
> > >> > > > > >
> > >> > > > > > North West 46.188 -123.625
> > >> > > > > >
> > >> > > > > > South West 45.125 -123.625
> > >> > > > > >
> > >> > > > > > I see there is the Polyline (poly) that might work
for me.
> Do
> > >> you
> > >> > > have
> > >> > > > an
> > >> > > > > > example for that case or the ASCII file that is
needed.
> > >> > > > > >
> > >> > > > > > Thank you
> > >> > > > > >
> > >> > > > > > Tatiana
> > >> > > > > >
> > >> > > > > > On Tue, Aug 25, 2020 at 5:52 PM Julie Prestopnik via
RT <
> > >> > > > > met_help at ucar.edu
> > >> > > > > > >
> > >> > > > > > wrote:
> > >> > > > > >
> > >> > > > > > > Hi Tatiana.
> > >> > > > > > >
> > >> > > > > > > I see you have a question about how to create a
mask file.
> > >> I'd
> > >> > > > > recommend
> > >> > > > > > > running the Gen-Vx-Mask tool to create the desired
masking
> > >> > region,
> > >> > > > and
> > >> > > > > > then
> > >> > > > > > > reference that in the Point-Stat config file. Here
is a
> > link
> > >> to
> > >> > > our
> > >> > > > > > > documentation page on Gen-Vx-Mask, but I also go
through
> an
> > >> > example
> > >> > > > > > below:
> > >> > > > > > >
https://dtcenter.github.io/MET/Users_Guide/masking.html
> > >> > > > > > >
> > >> > > > > > > To illustrate, I pulled a sample GFS file:
> > >> > > > > > >
> > >> > > > > > > > wget
> > >> > > > > > > >
> > >> > > > > > >
> > >> > > > > >
> > >> > > > >
> > >> > > >
> > >> > >
> > >> >
> > >>
> >
> https://www.ncei.noaa.gov/data/global-forecast-system/access/grid-
004-0.5-degree/forecast/202008/20200820/gfs_4_20200820_0000_012.grb2
> > >> > > > > > >
> > >> > > > > > >
> > >> > > > > > > I ran wgrib2 and look for a LAND/SEA mask:
> > >> > > > > > >
> > >> > > > > > > > wgrib2 gfs_4_20200820_0000_012.grb2 | egrep -i
LAND
> > >> > > > > > > >
> > >> > > > > > > which returned:
> > >> > > > > > >
> > >> > > > > > > > 571:98325645:d=2020082000:LAND:surface:12 hour
fcst:
> > >> > > > > > >
> > >> > > > > > >
> > >> > > > > > > I ran plot_data_plane on it (see attached image):
> > >> > > > > > >
> > >> > > > > > > > met-9.1/bin/plot_data_plane
gfs_4_20200820_0000_012.grb2
> > >> > > > > gfs_4_LAND.ps
> > >> > > > > > > > 'name="LAND"; level="L0";'
> > >> > > > > > > >
> > >> > > > > > >
> > >> > > > > > >
> > >> > > > > > > > DEBUG 1: Opening data file:
gfs_4_20200820_0000_012.grb2
> > >> > > > > > > > DEBUG 1: Creating postscript file: gfs_4_LAND.ps
> > >> > > > > > > >
> > >> > > > > > >
> > >> > > > > > > I ran gen_vx_mask to create a mask for this
verification
> > >> domain.
> > >> > > If
> > >> > > > > > you're
> > >> > > > > > > actually verifying on the GFS domain, the command
looks a
> > >> little
> > >> > > > > > redundant:
> > >> > > > > > >
> > >> > > > > > > met-9.1/bin/gen_vx_mask met-9.1/bin/gen_vx_mask \
> > >> > > > > > > > gfs_4_20200820_0000_012.grb2 \
> > >> > > > > > > > gfs_4_20200820_0000_012.grb2 \
> > >> > > > > > > > gfs_4_WATER_mask.nc \
> > >> > > > > > > > -type data \
> > >> > > > > > > > -mask_field 'name="LAND"; level="L0";' \
> > >> > > > > > > > -thresh ==0 \
> > >> > > > > > > > -name WATER
> > >> > > > > > >
> > >> > > > > > >
> > >> > > > > > > DEBUG 1: Input Grid: gfs_4_20200820_0000_012.grb2
> > >> > > > > > > > DEBUG 1: Mask File: gfs_4_20200820_0000_012.grb2
> > >> > > > > > > > DEBUG 2: Parsed Input Grid: LatLon (720 x 361)
> > >> > > > > > > > DEBUG 2: Parsed Mask Grid: LatLon (720 x 361)
> > >> > > > > > > > DEBUG 1: Output File: gfs_4_WATER_mask.nc
> > >> > > > > > >
> > >> > > > > > >
> > >> > > > > > > I defined the threshold as "==0" because I want the
WATER
> > >> points,
> > >> > > not
> > >> > > > > the
> > >> > > > > > > LAND points. And I used the "-name WATER" command
line
> > >> option to
> > >> > > > name
> > >> > > > > > the
> > >> > > > > > > output NetCDF variable. This will show up in the
VX_MASK
> > >> output
> > >> > > > column
> > >> > > > > > of
> > >> > > > > > > Point-Stat.
> > >> > > > > > >
> > >> > > > > > > And finally, you would then setup your Point-Stat
config
> > file
> > >> > like
> > >> > > > > this:
> > >> > > > > > >
> > >> > > > > > > mask = {
> > >> > > > > > > grid = [];
> > >> > > > > > > poly = [ "/path/to/gfs_4_WATER_mask.nc" ];
> > >> > > > > > > sid = [];
> > >> > > > > > > }
> > >> > > > > > >
> > >> > > > > > > I hope that helps!
> > >> > > > > > >
> > >> > > > > > > Julie
> > >> > > > > > >
> > >> > > > > > > On Tue, Aug 25, 2020 at 1:14 PM Tatiana Gonzalez -
NOAA
> > >> Federal
> > >> > via
> > >> > > > RT
> > >> > > > > <
> > >> > > > > > > met_help at ucar.edu> wrote:
> > >> > > > > > >
> > >> > > > > > > >
> > >> > > > > > > > Tue Aug 25 13:14:03 2020: Request 96443 was acted
upon.
> > >> > > > > > > > Transaction: Ticket created by
> tatiana.gonzalez at noaa.gov
> > >> > > > > > > > Queue: met_help
> > >> > > > > > > > Subject: mask file for Metplus
> > >> > > > > > > > Owner: Nobody
> > >> > > > > > > > Requestors: tatiana.gonzalez at noaa.gov
> > >> > > > > > > > Status: new
> > >> > > > > > > > Ticket <URL:
> > >> > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
> > >> > > > > > >
> > >> > > > > > > >
> > >> > > > > > > >
> > >> > > > > > > > Hello,
> > >> > > > > > > > How do I create a msk file to be use in Metplus
(point
> > >> stat)?
> > >> > to
> > >> > > > > limit
> > >> > > > > > > my
> > >> > > > > > > > analysis
> > >> > > > > > > >
> > >> > > > > > > > Thanks
> > >> > > > > > > >
> > >> > > > > > > >
> > >> > > > > > > > Tatiana
> > >> > > > > > > > --
> > >> > > > > > > >
> > >> > > > > > > >
> > >> > > > > > > > *Tatiana D. Gonzalez, PhD.*
> > >> > > > > > > > *Analysis & Nowcast Branch*
> > >> > > > > > > > *Analysis and Mission Support Division*
> > >> > > > > > > >
> > >> > > > > > > > *Analyze, Forecast and Support Office*
> > >> > > > > > > >
> > >> > > > > > > > *National Weather Service**SSMC2, Room 14115*
> > >> > > > > > > > *Silver Spring, Maryland*
> > >> > > > > > > >
> > >> > > > > > > >
> > >> > > > > > > > *Ph. (301) 427-9500 *
> > >> > > > > > > >
> > >> > > > > > > >
> > >> > > > > > >
> > >> > > > > > > --
> > >> > > > > > > Julie Prestopnik (she/her/hers)
> > >> > > > > > > Software Engineer
> > >> > > > > > > National Center for Atmospheric Research
> > >> > > > > > > Research Applications Laboratory
> > >> > > > > > > Phone: 303.497.8399
> > >> > > > > > > Email: jpresto at ucar.edu
> > >> > > > > > >
> > >> > > > > > > My working day may not be your working day. Please
do not
> > >> feel
> > >> > > > obliged
> > >> > > > > > to
> > >> > > > > > > reply to this email outside of your normal working
hours.
> > >> > > > > > >
> > >> > > > > > >
> > >> > > > > >
> > >> > > > > > --
> > >> > > > > >
> > >> > > > > >
> > >> > > > > > *Tatiana D. Gonzalez, PhD.*
> > >> > > > > > *Analysis & Nowcast Branch*
> > >> > > > > > *Analysis and Mission Support Division*
> > >> > > > > >
> > >> > > > > > *Analyze, Forecast and Support Office*
> > >> > > > > >
> > >> > > > > > *National Weather Service**SSMC2, Room 14115*
> > >> > > > > > *Silver Spring, Maryland*
> > >> > > > > >
> > >> > > > > >
> > >> > > > > > *Ph. (301) 427-9500 *
> > >> > > > > >
> > >> > > > > >
> > >> > > > >
> > >> > > > > --
> > >> > > > > Julie Prestopnik (she/her/hers)
> > >> > > > > Software Engineer
> > >> > > > > National Center for Atmospheric Research
> > >> > > > > Research Applications Laboratory
> > >> > > > > Phone: 303.497.8399
> > >> > > > > Email: jpresto at ucar.edu
> > >> > > > >
> > >> > > > > My working day may not be your working day. Please do
not
> feel
> > >> > obliged
> > >> > > > to
> > >> > > > > reply to this email outside of your normal working
hours.
> > >> > > > >
> > >> > > > >
> > >> > > >
> > >> > > > --
> > >> > > >
> > >> > > >
> > >> > > > *Tatiana D. Gonzalez, PhD.*
> > >> > > > *Analysis & Nowcast Branch*
> > >> > > > *Analysis and Mission Support Division*
> > >> > > >
> > >> > > > *Analyze, Forecast and Support Office*
> > >> > > >
> > >> > > > *National Weather Service**SSMC2, Room 14115*
> > >> > > > *Silver Spring, Maryland*
> > >> > > >
> > >> > > >
> > >> > > > *Ph. (301) 427-9500 *
> > >> > > >
> > >> > > >
> > >> > >
> > >> > >
> > >> >
> > >> > --
> > >> >
> > >> >
> > >> > *Tatiana D. Gonzalez, PhD.*
> > >> > *Analysis & Nowcast Branch*
> > >> > *Analysis and Mission Support Division*
> > >> >
> > >> > *Analyze, Forecast and Support Office*
> > >> >
> > >> > *National Weather Service**SSMC2, Room 14115*
> > >> > *Silver Spring, Maryland*
> > >> >
> > >> >
> > >> > *Ph. (301) 427-9500 *
> > >> >
> > >> >
> > >>
> > >> --
> > >> Julie Prestopnik (she/her/hers)
> > >> Software Engineer
> > >> National Center for Atmospheric Research
> > >> Research Applications Laboratory
> > >> Phone: 303.497.8399
> > >> Email: jpresto at ucar.edu
> > >>
> > >> My working day may not be your working day. Please do not feel
> obliged
> > to
> > >> reply to this email outside of your normal working hours.
> > >>
> > >>
> > >
> > > --
> > >
> > >
> > > *Tatiana D. Gonzalez, PhD.*
> > > *Analysis & Nowcast Branch*
> > > *Analysis and Mission Support Division*
> > >
> > > *Analyze, Forecast and Support Office*
> > >
> > > *National Weather Service**SSMC2, Room 14115*
> > > *Silver Spring, Maryland*
> > >
> > >
> > > *Ph. (301) 427-9500 *
> > >
> > >
> >
> > --
> >
> >
> > *Tatiana D. Gonzalez, PhD.*
> > *Analysis & Nowcast Branch*
> > *Analysis and Mission Support Division*
> >
> > *Analyze, Forecast and Support Office*
> >
> > *National Weather Service**SSMC2, Room 14115*
> > *Silver Spring, Maryland*
> >
> >
> > *Ph. (301) 427-9500 *
> >
> >
>
> --
> Julie Prestopnik (she/her/hers)
> Software Engineer
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
> My working day may not be your working day. Please do not feel
obliged to
> reply to this email outside of your normal working hours.
>
>
--
*Tatiana D. Gonzalez, PhD.*
*Analysis & Nowcast Branch*
*Analysis and Mission Support Division*
*Analyze, Forecast and Support Office*
*National Weather Service**SSMC2, Room 14115*
*Silver Spring, Maryland*
*Ph. (301) 427-9500 *
------------------------------------------------
Subject: mask file for Metplus
From: John Halley Gotway
Time: Thu Sep 10 12:30:52 2020
Tatiana,
Hello, this is John Halley Gotway. I work with Julie on the
development and
support of MET. I see that you have some questions about creating a
relatively small masking region around Miami.
If you're running the Point-Stat tool in MET, there are multiple
options
for defining masks. You're doing it spatially, but another option is
specifying a list of station id's to be used. Below is an excerpt from
the
user's guide about this:
https://dtcenter.github.io/MET/Users_Guide/data_io.html
// - The "sid" entry is an array of strings which define groups of
// observation station ID's over which to compute statistics.
Each entry
// in the array is either a filename of a comma-separated list.
// - For a filename, the strings are whitespace-separated. The
first
// string is the mask "name" and the remaining strings are the
station
// ID's to be used.
// - For a comma-separated list, optionally use a colon to
specify a name.
// For "MY_LIST:SID1,SID2", name = MY_LIST and values = SID1
and SID2.
// - For a comma-separated list of length one with no name
specified, the
// mask "name" and value are both set to the single station ID
string.
// For "SID1", name = SID1 and value = SID1.
// - For a comma-separated list of length greater than one with
no name
// specified, the name is set to MASK_SID and the values are
the station
// ID's to be used.
// For "SID1,SID2", name = MASK_SID and values = SID1 and
SID2.
// - The "name" of the station ID mask is written to the VX_MASK
column
// of the MET output files.
mask = {
...
sid = [ "CONUS.stations" ];
...
}
So if you have a specific set of stations in mind, it may be more
convenient to specify them as a list of strings instead of a spatial
area.
Thanks,
John
On Thu, Sep 10, 2020 at 11:41 AM Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
>
> I had not yet had a chance, because I was in a meeting, but I just
took a
> look. I ran it through plot_data_plane:
>
> plot_data_plane case_2.nc case.ps 'name="MIAMI"; level="(*,*)";'
>
>
> I am attaching the resulting case.ps file. Is that where you
expected it
> to be?
>
> Julie
>
>
>
> On Thu, Sep 10, 2020 at 11:19 AM Tatiana Gonzalez - NOAA Federal via
RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
> >
> > Julie,
> > I think it was my bad. the coordinates were not correctly
converted. DId
> > you have a chance to look at the grid? Just want to make sure it
is right
> > anyways. It should be centered in MIAMI. I don't have a way to
plot it
> > right now not because WCOSS is down.
> > I apologize for the inconvenience.
> >
> > Tatiana
> >
> >
> > On Thu, Sep 10, 2020 at 11:45 AM Tatiana Gonzalez - NOAA Federal <
> > tatiana.gonzalez at noaa.gov> wrote:
> >
> > > Try again. I just gave you permission
> > >
> > > On Thu, Sep 10, 2020 at 11:36 AM Julie Prestopnik via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > >> Hi Tatiana.
> > >>
> > >> I sent a request for access through google drive as I did not
have
> > >> permissions on the file.
> > >>
> > >> Julie
> > >>
> > >> On Thu, Sep 10, 2020 at 8:58 AM Tatiana Gonzalez - NOAA Federal
via
> RT <
> > >> met_help at ucar.edu> wrote:
> > >>
> > >> >
> > >> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
>
> > >> >
> > >> > Jullie,
> > >> > Thank you for reviewing the file
> > >> > case_2.nc
> > >> > <
> > >> >
> > >>
> >
> https://drive.google.com/file/d/1tWThlec9seGOruztsodx73wHAD3OU-
gV/view?usp=drive_web
> > >> > >
> > >> >
> > >> >
> > >> > Tatiana
> > >> >
> > >> > On Thu, Sep 10, 2020 at 10:38 AM Julie Prestopnik via RT <
> > >> > met_help at ucar.edu>
> > >> > wrote:
> > >> >
> > >> > > Hi Tatiana.
> > >> > >
> > >> > > Could you please either attach file.nc or upload file.nc to
our
> FTP
> > >> > > server as
> > >> > > described under "How to Send Us Data" on the following
page?
> > >> > >
> > >> > >
> > >> >
> > >>
> >
> https://dtcenter.org/community-code/model-evaluation-tools-met/met-
help-desk
> > >> > >
> > >> > >
> > >> > > I'll take a look.
> > >> > >
> > >> > > Thanks!
> > >> > >
> > >> > > Julie
> > >> > >
> > >> > >
> > >> > > On Thu, Sep 10, 2020, 7:10 AM Tatiana Gonzalez - NOAA
Federal via
> > RT <
> > >> > > met_help at ucar.edu> wrote:
> > >> > >
> > >> > > >
> > >> > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
> > >> > > >
> > >> > > > Hello Jullie,
> > >> > > > I created a file with 4 coordinates and ran the polyline
through
> > the
> > >> > > > gen_vx_mask tool. Then on the config file I use the new "
> file.nc"
> > >> > > > But it looks like none of the stations that are within my
grid
> > were
> > >> > > > included. Is there an order to the coordinates? My file
was in
> the
> > >> > > > following order (south east, north east, north west,
south west)
> > >> > > >
> > >> > > > MIAMI
> > >> > > >
> > >> > > > 25.653243 -79.832631
> > >> > > >
> > >> > > > 28.073608 -79.832631
> > >> > > >
> > >> > > > 28.073608 -80.938611
> > >> > > >
> > >> > > > 25.602028 -80.938611
> > >> > > >
> > >> > > >
> > >> > > > The grib file I used contained the whole CONUS. Is there
> something
> > >> I'm
> > >> > > > doing wrong?
> > >> > > >
> > >> > > > Thanks
> > >> > > >
> > >> > > > Tatiana
> > >> > > >
> > >> > > >
> > >> > > >
> > >> > > > On Wed, Aug 26, 2020 at 2:54 PM Julie Prestopnik via RT <
> > >> > > met_help at ucar.edu
> > >> > > > >
> > >> > > > wrote:
> > >> > > >
> > >> > > > > Hi Tatiana.
> > >> > > > >
> > >> > > > > To create a polyline mask, you simply create an ASCII
file
> > >> > containing a
> > >> > > > > list of lat/lons that define the region of interest.
For
> > example,
> > >> > > take a
> > >> > > > > look at met-9.1/data/poly/WEST.poly. The first line of
the
> file
> > >> > > > specifies
> > >> > > > > the name for the polyline. And the remaining lines
contain
> > >> lat/lon
> > >> > > > values
> > >> > > > > for the boundary points. The last point is connected
back up
> to
> > >> the
> > >> > > > first
> > >> > > > > point to complete the area. You need a minimum of 3
points to
> > >> define
> > >> > a
> > >> > > > > polyline, and there is no real maximum. If your
polyline
> > >> contains a
> > >> > > > large
> > >> > > > > number of points (a few hundred) or you have many grid-
points,
> > the
> > >> > > > process
> > >> > > > > of determining which grid points are inside/outside the
> polyline
> > >> can
> > >> > be
> > >> > > > > slow. In that case, I'd suggest running your polyline
through
> > the
> > >> > > > > gen_vx_mask tool.
> > >> > > > >
> > >> > > > > Julie
> > >> > > > >
> > >> > > > > On Wed, Aug 26, 2020 at 12:03 PM Tatiana Gonzalez -
NOAA
> Federal
> > >> via
> > >> > > RT <
> > >> > > > > met_help at ucar.edu> wrote:
> > >> > > > >
> > >> > > > > >
> > >> > > > > > <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
> > >
> > >> > > > > >
> > >> > > > > > Hi Julie,
> > >> > > > > > Thanks for the instructions. I have the RTMA files in
grib2
> > >> format
> > >> > > for
> > >> > > > > the
> > >> > > > > > whole CONUS. But what I want is only use the data
contained
> in
> > >> this
> > >> > > > box:
> > >> > > > > >
> > >> > > > > > South East 45.125 -121.975
> > >> > > > > >
> > >> > > > > > North East 46.188 -121.975
> > >> > > > > >
> > >> > > > > > North West 46.188 -123.625
> > >> > > > > >
> > >> > > > > > South West 45.125 -123.625
> > >> > > > > >
> > >> > > > > > I see there is the Polyline (poly) that might work
for me.
> Do
> > >> you
> > >> > > have
> > >> > > > an
> > >> > > > > > example for that case or the ASCII file that is
needed.
> > >> > > > > >
> > >> > > > > > Thank you
> > >> > > > > >
> > >> > > > > > Tatiana
> > >> > > > > >
> > >> > > > > > On Tue, Aug 25, 2020 at 5:52 PM Julie Prestopnik via
RT <
> > >> > > > > met_help at ucar.edu
> > >> > > > > > >
> > >> > > > > > wrote:
> > >> > > > > >
> > >> > > > > > > Hi Tatiana.
> > >> > > > > > >
> > >> > > > > > > I see you have a question about how to create a
mask file.
> > >> I'd
> > >> > > > > recommend
> > >> > > > > > > running the Gen-Vx-Mask tool to create the desired
masking
> > >> > region,
> > >> > > > and
> > >> > > > > > then
> > >> > > > > > > reference that in the Point-Stat config file. Here
is a
> > link
> > >> to
> > >> > > our
> > >> > > > > > > documentation page on Gen-Vx-Mask, but I also go
through
> an
> > >> > example
> > >> > > > > > below:
> > >> > > > > > >
https://dtcenter.github.io/MET/Users_Guide/masking.html
> > >> > > > > > >
> > >> > > > > > > To illustrate, I pulled a sample GFS file:
> > >> > > > > > >
> > >> > > > > > > > wget
> > >> > > > > > > >
> > >> > > > > > >
> > >> > > > > >
> > >> > > > >
> > >> > > >
> > >> > >
> > >> >
> > >>
> >
> https://www.ncei.noaa.gov/data/global-forecast-system/access/grid-
004-0.5-degree/forecast/202008/20200820/gfs_4_20200820_0000_012.grb2
> > >> > > > > > >
> > >> > > > > > >
> > >> > > > > > > I ran wgrib2 and look for a LAND/SEA mask:
> > >> > > > > > >
> > >> > > > > > > > wgrib2 gfs_4_20200820_0000_012.grb2 | egrep -i
LAND
> > >> > > > > > > >
> > >> > > > > > > which returned:
> > >> > > > > > >
> > >> > > > > > > > 571:98325645:d=2020082000:LAND:surface:12 hour
fcst:
> > >> > > > > > >
> > >> > > > > > >
> > >> > > > > > > I ran plot_data_plane on it (see attached image):
> > >> > > > > > >
> > >> > > > > > > > met-9.1/bin/plot_data_plane
gfs_4_20200820_0000_012.grb2
> > >> > > > > gfs_4_LAND.ps
> > >> > > > > > > > 'name="LAND"; level="L0";'
> > >> > > > > > > >
> > >> > > > > > >
> > >> > > > > > >
> > >> > > > > > > > DEBUG 1: Opening data file:
gfs_4_20200820_0000_012.grb2
> > >> > > > > > > > DEBUG 1: Creating postscript file: gfs_4_LAND.ps
> > >> > > > > > > >
> > >> > > > > > >
> > >> > > > > > > I ran gen_vx_mask to create a mask for this
verification
> > >> domain.
> > >> > > If
> > >> > > > > > you're
> > >> > > > > > > actually verifying on the GFS domain, the command
looks a
> > >> little
> > >> > > > > > redundant:
> > >> > > > > > >
> > >> > > > > > > met-9.1/bin/gen_vx_mask met-9.1/bin/gen_vx_mask \
> > >> > > > > > > > gfs_4_20200820_0000_012.grb2 \
> > >> > > > > > > > gfs_4_20200820_0000_012.grb2 \
> > >> > > > > > > > gfs_4_WATER_mask.nc \
> > >> > > > > > > > -type data \
> > >> > > > > > > > -mask_field 'name="LAND"; level="L0";' \
> > >> > > > > > > > -thresh ==0 \
> > >> > > > > > > > -name WATER
> > >> > > > > > >
> > >> > > > > > >
> > >> > > > > > > DEBUG 1: Input Grid: gfs_4_20200820_0000_012.grb2
> > >> > > > > > > > DEBUG 1: Mask File: gfs_4_20200820_0000_012.grb2
> > >> > > > > > > > DEBUG 2: Parsed Input Grid: LatLon (720 x 361)
> > >> > > > > > > > DEBUG 2: Parsed Mask Grid: LatLon (720 x 361)
> > >> > > > > > > > DEBUG 1: Output File: gfs_4_WATER_mask.nc
> > >> > > > > > >
> > >> > > > > > >
> > >> > > > > > > I defined the threshold as "==0" because I want the
WATER
> > >> points,
> > >> > > not
> > >> > > > > the
> > >> > > > > > > LAND points. And I used the "-name WATER" command
line
> > >> option to
> > >> > > > name
> > >> > > > > > the
> > >> > > > > > > output NetCDF variable. This will show up in the
VX_MASK
> > >> output
> > >> > > > column
> > >> > > > > > of
> > >> > > > > > > Point-Stat.
> > >> > > > > > >
> > >> > > > > > > And finally, you would then setup your Point-Stat
config
> > file
> > >> > like
> > >> > > > > this:
> > >> > > > > > >
> > >> > > > > > > mask = {
> > >> > > > > > > grid = [];
> > >> > > > > > > poly = [ "/path/to/gfs_4_WATER_mask.nc" ];
> > >> > > > > > > sid = [];
> > >> > > > > > > }
> > >> > > > > > >
> > >> > > > > > > I hope that helps!
> > >> > > > > > >
> > >> > > > > > > Julie
> > >> > > > > > >
> > >> > > > > > > On Tue, Aug 25, 2020 at 1:14 PM Tatiana Gonzalez -
NOAA
> > >> Federal
> > >> > via
> > >> > > > RT
> > >> > > > > <
> > >> > > > > > > met_help at ucar.edu> wrote:
> > >> > > > > > >
> > >> > > > > > > >
> > >> > > > > > > > Tue Aug 25 13:14:03 2020: Request 96443 was acted
upon.
> > >> > > > > > > > Transaction: Ticket created by
> tatiana.gonzalez at noaa.gov
> > >> > > > > > > > Queue: met_help
> > >> > > > > > > > Subject: mask file for Metplus
> > >> > > > > > > > Owner: Nobody
> > >> > > > > > > > Requestors: tatiana.gonzalez at noaa.gov
> > >> > > > > > > > Status: new
> > >> > > > > > > > Ticket <URL:
> > >> > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
> > >> > > > > > >
> > >> > > > > > > >
> > >> > > > > > > >
> > >> > > > > > > > Hello,
> > >> > > > > > > > How do I create a msk file to be use in Metplus
(point
> > >> stat)?
> > >> > to
> > >> > > > > limit
> > >> > > > > > > my
> > >> > > > > > > > analysis
> > >> > > > > > > >
> > >> > > > > > > > Thanks
> > >> > > > > > > >
> > >> > > > > > > >
> > >> > > > > > > > Tatiana
> > >> > > > > > > > --
> > >> > > > > > > >
> > >> > > > > > > >
> > >> > > > > > > > *Tatiana D. Gonzalez, PhD.*
> > >> > > > > > > > *Analysis & Nowcast Branch*
> > >> > > > > > > > *Analysis and Mission Support Division*
> > >> > > > > > > >
> > >> > > > > > > > *Analyze, Forecast and Support Office*
> > >> > > > > > > >
> > >> > > > > > > > *National Weather Service**SSMC2, Room 14115*
> > >> > > > > > > > *Silver Spring, Maryland*
> > >> > > > > > > >
> > >> > > > > > > >
> > >> > > > > > > > *Ph. (301) 427-9500 *
> > >> > > > > > > >
> > >> > > > > > > >
> > >> > > > > > >
> > >> > > > > > > --
> > >> > > > > > > Julie Prestopnik (she/her/hers)
> > >> > > > > > > Software Engineer
> > >> > > > > > > National Center for Atmospheric Research
> > >> > > > > > > Research Applications Laboratory
> > >> > > > > > > Phone: 303.497.8399
> > >> > > > > > > Email: jpresto at ucar.edu
> > >> > > > > > >
> > >> > > > > > > My working day may not be your working day. Please
do not
> > >> feel
> > >> > > > obliged
> > >> > > > > > to
> > >> > > > > > > reply to this email outside of your normal working
hours.
> > >> > > > > > >
> > >> > > > > > >
> > >> > > > > >
> > >> > > > > > --
> > >> > > > > >
> > >> > > > > >
> > >> > > > > > *Tatiana D. Gonzalez, PhD.*
> > >> > > > > > *Analysis & Nowcast Branch*
> > >> > > > > > *Analysis and Mission Support Division*
> > >> > > > > >
> > >> > > > > > *Analyze, Forecast and Support Office*
> > >> > > > > >
> > >> > > > > > *National Weather Service**SSMC2, Room 14115*
> > >> > > > > > *Silver Spring, Maryland*
> > >> > > > > >
> > >> > > > > >
> > >> > > > > > *Ph. (301) 427-9500 *
> > >> > > > > >
> > >> > > > > >
> > >> > > > >
> > >> > > > > --
> > >> > > > > Julie Prestopnik (she/her/hers)
> > >> > > > > Software Engineer
> > >> > > > > National Center for Atmospheric Research
> > >> > > > > Research Applications Laboratory
> > >> > > > > Phone: 303.497.8399
> > >> > > > > Email: jpresto at ucar.edu
> > >> > > > >
> > >> > > > > My working day may not be your working day. Please do
not
> feel
> > >> > obliged
> > >> > > > to
> > >> > > > > reply to this email outside of your normal working
hours.
> > >> > > > >
> > >> > > > >
> > >> > > >
> > >> > > > --
> > >> > > >
> > >> > > >
> > >> > > > *Tatiana D. Gonzalez, PhD.*
> > >> > > > *Analysis & Nowcast Branch*
> > >> > > > *Analysis and Mission Support Division*
> > >> > > >
> > >> > > > *Analyze, Forecast and Support Office*
> > >> > > >
> > >> > > > *National Weather Service**SSMC2, Room 14115*
> > >> > > > *Silver Spring, Maryland*
> > >> > > >
> > >> > > >
> > >> > > > *Ph. (301) 427-9500 *
> > >> > > >
> > >> > > >
> > >> > >
> > >> > >
> > >> >
> > >> > --
> > >> >
> > >> >
> > >> > *Tatiana D. Gonzalez, PhD.*
> > >> > *Analysis & Nowcast Branch*
> > >> > *Analysis and Mission Support Division*
> > >> >
> > >> > *Analyze, Forecast and Support Office*
> > >> >
> > >> > *National Weather Service**SSMC2, Room 14115*
> > >> > *Silver Spring, Maryland*
> > >> >
> > >> >
> > >> > *Ph. (301) 427-9500 *
> > >> >
> > >> >
> > >>
> > >> --
> > >> Julie Prestopnik (she/her/hers)
> > >> Software Engineer
> > >> National Center for Atmospheric Research
> > >> Research Applications Laboratory
> > >> Phone: 303.497.8399
> > >> Email: jpresto at ucar.edu
> > >>
> > >> My working day may not be your working day. Please do not feel
> obliged
> > to
> > >> reply to this email outside of your normal working hours.
> > >>
> > >>
> > >
> > > --
> > >
> > >
> > > *Tatiana D. Gonzalez, PhD.*
> > > *Analysis & Nowcast Branch*
> > > *Analysis and Mission Support Division*
> > >
> > > *Analyze, Forecast and Support Office*
> > >
> > > *National Weather Service**SSMC2, Room 14115*
> > > *Silver Spring, Maryland*
> > >
> > >
> > > *Ph. (301) 427-9500 *
> > >
> > >
> >
> > --
> >
> >
> > *Tatiana D. Gonzalez, PhD.*
> > *Analysis & Nowcast Branch*
> > *Analysis and Mission Support Division*
> >
> > *Analyze, Forecast and Support Office*
> >
> > *National Weather Service**SSMC2, Room 14115*
> > *Silver Spring, Maryland*
> >
> >
> > *Ph. (301) 427-9500 *
> >
> >
>
> --
> Julie Prestopnik (she/her/hers)
> Software Engineer
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
> My working day may not be your working day. Please do not feel
obliged to
> reply to this email outside of your normal working hours.
>
>
------------------------------------------------
Subject: mask file for Metplus
From: Tatiana Gonzalez - NOAA Federal
Time: Thu Sep 10 13:54:15 2020
Hi John,
Thank you. I do see how this might help. l will try this methods for
some
cases that I have the stations
regards
Tatiana
On Thu, Sep 10, 2020 at 2:40 PM John Halley Gotway via RT
<met_help at ucar.edu>
wrote:
> Tatiana,
>
> Hello, this is John Halley Gotway. I work with Julie on the
development and
> support of MET. I see that you have some questions about creating a
> relatively small masking region around Miami.
>
> If you're running the Point-Stat tool in MET, there are multiple
options
> for defining masks. You're doing it spatially, but another option is
> specifying a list of station id's to be used. Below is an excerpt
from the
> user's guide about this:
> https://dtcenter.github.io/MET/Users_Guide/data_io.html
>
> // - The "sid" entry is an array of strings which define groups
of
> // observation station ID's over which to compute statistics.
Each
> entry
> // in the array is either a filename of a comma-separated list.
> // - For a filename, the strings are whitespace-separated. The
first
> // string is the mask "name" and the remaining strings are
the
> station
> // ID's to be used.
> // - For a comma-separated list, optionally use a colon to
specify a
> name.
> // For "MY_LIST:SID1,SID2", name = MY_LIST and values = SID1
and
> SID2.
> // - For a comma-separated list of length one with no name
> specified, the
> // mask "name" and value are both set to the single station
ID
> string.
> // For "SID1", name = SID1 and value = SID1.
> // - For a comma-separated list of length greater than one
with no
> name
> // specified, the name is set to MASK_SID and the values are
the
> station
> // ID's to be used.
> // For "SID1,SID2", name = MASK_SID and values = SID1 and
SID2.
> // - The "name" of the station ID mask is written to the
VX_MASK
> column
> // of the MET output files.
>
> mask = {
> ...
> sid = [ "CONUS.stations" ];
>
> ...
> }
>
> So if you have a specific set of stations in mind, it may be more
> convenient to specify them as a list of strings instead of a spatial
area.
>
> Thanks,
> John
>
> On Thu, Sep 10, 2020 at 11:41 AM Julie Prestopnik via RT <
> met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
> >
> > I had not yet had a chance, because I was in a meeting, but I just
took a
> > look. I ran it through plot_data_plane:
> >
> > plot_data_plane case_2.nc case.ps 'name="MIAMI"; level="(*,*)";'
> >
> >
> > I am attaching the resulting case.ps file. Is that where you
expected
> it
> > to be?
> >
> > Julie
> >
> >
> >
> > On Thu, Sep 10, 2020 at 11:19 AM Tatiana Gonzalez - NOAA Federal
via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
> > >
> > > Julie,
> > > I think it was my bad. the coordinates were not correctly
converted.
> DId
> > > you have a chance to look at the grid? Just want to make sure it
is
> right
> > > anyways. It should be centered in MIAMI. I don't have a way to
plot it
> > > right now not because WCOSS is down.
> > > I apologize for the inconvenience.
> > >
> > > Tatiana
> > >
> > >
> > > On Thu, Sep 10, 2020 at 11:45 AM Tatiana Gonzalez - NOAA Federal
<
> > > tatiana.gonzalez at noaa.gov> wrote:
> > >
> > > > Try again. I just gave you permission
> > > >
> > > > On Thu, Sep 10, 2020 at 11:36 AM Julie Prestopnik via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > >> Hi Tatiana.
> > > >>
> > > >> I sent a request for access through google drive as I did not
have
> > > >> permissions on the file.
> > > >>
> > > >> Julie
> > > >>
> > > >> On Thu, Sep 10, 2020 at 8:58 AM Tatiana Gonzalez - NOAA
Federal via
> > RT <
> > > >> met_help at ucar.edu> wrote:
> > > >>
> > > >> >
> > > >> > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443 >
> > > >> >
> > > >> > Jullie,
> > > >> > Thank you for reviewing the file
> > > >> > case_2.nc
> > > >> > <
> > > >> >
> > > >>
> > >
> >
> https://drive.google.com/file/d/1tWThlec9seGOruztsodx73wHAD3OU-
gV/view?usp=drive_web
> > > >> > >
> > > >> >
> > > >> >
> > > >> > Tatiana
> > > >> >
> > > >> > On Thu, Sep 10, 2020 at 10:38 AM Julie Prestopnik via RT <
> > > >> > met_help at ucar.edu>
> > > >> > wrote:
> > > >> >
> > > >> > > Hi Tatiana.
> > > >> > >
> > > >> > > Could you please either attach file.nc or upload file.nc
to our
> > FTP
> > > >> > > server as
> > > >> > > described under "How to Send Us Data" on the following
page?
> > > >> > >
> > > >> > >
> > > >> >
> > > >>
> > >
> >
> https://dtcenter.org/community-code/model-evaluation-tools-met/met-
help-desk
> > > >> > >
> > > >> > >
> > > >> > > I'll take a look.
> > > >> > >
> > > >> > > Thanks!
> > > >> > >
> > > >> > > Julie
> > > >> > >
> > > >> > >
> > > >> > > On Thu, Sep 10, 2020, 7:10 AM Tatiana Gonzalez - NOAA
Federal
> via
> > > RT <
> > > >> > > met_help at ucar.edu> wrote:
> > > >> > >
> > > >> > > >
> > > >> > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
> >
> > > >> > > >
> > > >> > > > Hello Jullie,
> > > >> > > > I created a file with 4 coordinates and ran the
polyline
> through
> > > the
> > > >> > > > gen_vx_mask tool. Then on the config file I use the new
"
> > file.nc"
> > > >> > > > But it looks like none of the stations that are within
my grid
> > > were
> > > >> > > > included. Is there an order to the coordinates? My file
was in
> > the
> > > >> > > > following order (south east, north east, north west,
south
> west)
> > > >> > > >
> > > >> > > > MIAMI
> > > >> > > >
> > > >> > > > 25.653243 -79.832631
> > > >> > > >
> > > >> > > > 28.073608 -79.832631
> > > >> > > >
> > > >> > > > 28.073608 -80.938611
> > > >> > > >
> > > >> > > > 25.602028 -80.938611
> > > >> > > >
> > > >> > > >
> > > >> > > > The grib file I used contained the whole CONUS. Is
there
> > something
> > > >> I'm
> > > >> > > > doing wrong?
> > > >> > > >
> > > >> > > > Thanks
> > > >> > > >
> > > >> > > > Tatiana
> > > >> > > >
> > > >> > > >
> > > >> > > >
> > > >> > > > On Wed, Aug 26, 2020 at 2:54 PM Julie Prestopnik via RT
<
> > > >> > > met_help at ucar.edu
> > > >> > > > >
> > > >> > > > wrote:
> > > >> > > >
> > > >> > > > > Hi Tatiana.
> > > >> > > > >
> > > >> > > > > To create a polyline mask, you simply create an ASCII
file
> > > >> > containing a
> > > >> > > > > list of lat/lons that define the region of interest.
For
> > > example,
> > > >> > > take a
> > > >> > > > > look at met-9.1/data/poly/WEST.poly. The first line
of the
> > file
> > > >> > > > specifies
> > > >> > > > > the name for the polyline. And the remaining lines
contain
> > > >> lat/lon
> > > >> > > > values
> > > >> > > > > for the boundary points. The last point is connected
back
> up
> > to
> > > >> the
> > > >> > > > first
> > > >> > > > > point to complete the area. You need a minimum of 3
points
> to
> > > >> define
> > > >> > a
> > > >> > > > > polyline, and there is no real maximum. If your
polyline
> > > >> contains a
> > > >> > > > large
> > > >> > > > > number of points (a few hundred) or you have many
> grid-points,
> > > the
> > > >> > > > process
> > > >> > > > > of determining which grid points are inside/outside
the
> > polyline
> > > >> can
> > > >> > be
> > > >> > > > > slow. In that case, I'd suggest running your
polyline
> through
> > > the
> > > >> > > > > gen_vx_mask tool.
> > > >> > > > >
> > > >> > > > > Julie
> > > >> > > > >
> > > >> > > > > On Wed, Aug 26, 2020 at 12:03 PM Tatiana Gonzalez -
NOAA
> > Federal
> > > >> via
> > > >> > > RT <
> > > >> > > > > met_help at ucar.edu> wrote:
> > > >> > > > >
> > > >> > > > > >
> > > >> > > > > > <URL:
> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
> > > >
> > > >> > > > > >
> > > >> > > > > > Hi Julie,
> > > >> > > > > > Thanks for the instructions. I have the RTMA files
in
> grib2
> > > >> format
> > > >> > > for
> > > >> > > > > the
> > > >> > > > > > whole CONUS. But what I want is only use the data
> contained
> > in
> > > >> this
> > > >> > > > box:
> > > >> > > > > >
> > > >> > > > > > South East 45.125 -121.975
> > > >> > > > > >
> > > >> > > > > > North East 46.188 -121.975
> > > >> > > > > >
> > > >> > > > > > North West 46.188 -123.625
> > > >> > > > > >
> > > >> > > > > > South West 45.125 -123.625
> > > >> > > > > >
> > > >> > > > > > I see there is the Polyline (poly) that might work
for me.
> > Do
> > > >> you
> > > >> > > have
> > > >> > > > an
> > > >> > > > > > example for that case or the ASCII file that is
needed.
> > > >> > > > > >
> > > >> > > > > > Thank you
> > > >> > > > > >
> > > >> > > > > > Tatiana
> > > >> > > > > >
> > > >> > > > > > On Tue, Aug 25, 2020 at 5:52 PM Julie Prestopnik
via RT <
> > > >> > > > > met_help at ucar.edu
> > > >> > > > > > >
> > > >> > > > > > wrote:
> > > >> > > > > >
> > > >> > > > > > > Hi Tatiana.
> > > >> > > > > > >
> > > >> > > > > > > I see you have a question about how to create a
mask
> file.
> > > >> I'd
> > > >> > > > > recommend
> > > >> > > > > > > running the Gen-Vx-Mask tool to create the
desired
> masking
> > > >> > region,
> > > >> > > > and
> > > >> > > > > > then
> > > >> > > > > > > reference that in the Point-Stat config file.
Here is a
> > > link
> > > >> to
> > > >> > > our
> > > >> > > > > > > documentation page on Gen-Vx-Mask, but I also go
through
> > an
> > > >> > example
> > > >> > > > > > below:
> > > >> > > > > > >
https://dtcenter.github.io/MET/Users_Guide/masking.html
> > > >> > > > > > >
> > > >> > > > > > > To illustrate, I pulled a sample GFS file:
> > > >> > > > > > >
> > > >> > > > > > > > wget
> > > >> > > > > > > >
> > > >> > > > > > >
> > > >> > > > > >
> > > >> > > > >
> > > >> > > >
> > > >> > >
> > > >> >
> > > >>
> > >
> >
> https://www.ncei.noaa.gov/data/global-forecast-system/access/grid-
004-0.5-degree/forecast/202008/20200820/gfs_4_20200820_0000_012.grb2
> > > >> > > > > > >
> > > >> > > > > > >
> > > >> > > > > > > I ran wgrib2 and look for a LAND/SEA mask:
> > > >> > > > > > >
> > > >> > > > > > > > wgrib2 gfs_4_20200820_0000_012.grb2 | egrep -i
LAND
> > > >> > > > > > > >
> > > >> > > > > > > which returned:
> > > >> > > > > > >
> > > >> > > > > > > > 571:98325645:d=2020082000:LAND:surface:12 hour
fcst:
> > > >> > > > > > >
> > > >> > > > > > >
> > > >> > > > > > > I ran plot_data_plane on it (see attached image):
> > > >> > > > > > >
> > > >> > > > > > > > met-9.1/bin/plot_data_plane
> gfs_4_20200820_0000_012.grb2
> > > >> > > > > gfs_4_LAND.ps
> > > >> > > > > > > > 'name="LAND"; level="L0";'
> > > >> > > > > > > >
> > > >> > > > > > >
> > > >> > > > > > >
> > > >> > > > > > > > DEBUG 1: Opening data file:
> gfs_4_20200820_0000_012.grb2
> > > >> > > > > > > > DEBUG 1: Creating postscript file:
gfs_4_LAND.ps
> > > >> > > > > > > >
> > > >> > > > > > >
> > > >> > > > > > > I ran gen_vx_mask to create a mask for this
verification
> > > >> domain.
> > > >> > > If
> > > >> > > > > > you're
> > > >> > > > > > > actually verifying on the GFS domain, the command
looks
> a
> > > >> little
> > > >> > > > > > redundant:
> > > >> > > > > > >
> > > >> > > > > > > met-9.1/bin/gen_vx_mask met-9.1/bin/gen_vx_mask \
> > > >> > > > > > > > gfs_4_20200820_0000_012.grb2 \
> > > >> > > > > > > > gfs_4_20200820_0000_012.grb2 \
> > > >> > > > > > > > gfs_4_WATER_mask.nc \
> > > >> > > > > > > > -type data \
> > > >> > > > > > > > -mask_field 'name="LAND"; level="L0";' \
> > > >> > > > > > > > -thresh ==0 \
> > > >> > > > > > > > -name WATER
> > > >> > > > > > >
> > > >> > > > > > >
> > > >> > > > > > > DEBUG 1: Input Grid: gfs_4_20200820_0000_012.grb2
> > > >> > > > > > > > DEBUG 1: Mask File:
gfs_4_20200820_0000_012.grb2
> > > >> > > > > > > > DEBUG 2: Parsed Input Grid: LatLon (720 x 361)
> > > >> > > > > > > > DEBUG 2: Parsed Mask Grid: LatLon (720 x 361)
> > > >> > > > > > > > DEBUG 1: Output File: gfs_4_WATER_mask.nc
> > > >> > > > > > >
> > > >> > > > > > >
> > > >> > > > > > > I defined the threshold as "==0" because I want
the
> WATER
> > > >> points,
> > > >> > > not
> > > >> > > > > the
> > > >> > > > > > > LAND points. And I used the "-name WATER"
command line
> > > >> option to
> > > >> > > > name
> > > >> > > > > > the
> > > >> > > > > > > output NetCDF variable. This will show up in the
> VX_MASK
> > > >> output
> > > >> > > > column
> > > >> > > > > > of
> > > >> > > > > > > Point-Stat.
> > > >> > > > > > >
> > > >> > > > > > > And finally, you would then setup your Point-Stat
config
> > > file
> > > >> > like
> > > >> > > > > this:
> > > >> > > > > > >
> > > >> > > > > > > mask = {
> > > >> > > > > > > grid = [];
> > > >> > > > > > > poly = [ "/path/to/gfs_4_WATER_mask.nc" ];
> > > >> > > > > > > sid = [];
> > > >> > > > > > > }
> > > >> > > > > > >
> > > >> > > > > > > I hope that helps!
> > > >> > > > > > >
> > > >> > > > > > > Julie
> > > >> > > > > > >
> > > >> > > > > > > On Tue, Aug 25, 2020 at 1:14 PM Tatiana Gonzalez
- NOAA
> > > >> Federal
> > > >> > via
> > > >> > > > RT
> > > >> > > > > <
> > > >> > > > > > > met_help at ucar.edu> wrote:
> > > >> > > > > > >
> > > >> > > > > > > >
> > > >> > > > > > > > Tue Aug 25 13:14:03 2020: Request 96443 was
acted
> upon.
> > > >> > > > > > > > Transaction: Ticket created by
> > tatiana.gonzalez at noaa.gov
> > > >> > > > > > > > Queue: met_help
> > > >> > > > > > > > Subject: mask file for Metplus
> > > >> > > > > > > > Owner: Nobody
> > > >> > > > > > > > Requestors: tatiana.gonzalez at noaa.gov
> > > >> > > > > > > > Status: new
> > > >> > > > > > > > Ticket <URL:
> > > >> > > > >
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96443
> > > >> > > > > > >
> > > >> > > > > > > >
> > > >> > > > > > > >
> > > >> > > > > > > > Hello,
> > > >> > > > > > > > How do I create a msk file to be use in Metplus
(point
> > > >> stat)?
> > > >> > to
> > > >> > > > > limit
> > > >> > > > > > > my
> > > >> > > > > > > > analysis
> > > >> > > > > > > >
> > > >> > > > > > > > Thanks
> > > >> > > > > > > >
> > > >> > > > > > > >
> > > >> > > > > > > > Tatiana
> > > >> > > > > > > > --
> > > >> > > > > > > >
> > > >> > > > > > > >
> > > >> > > > > > > > *Tatiana D. Gonzalez, PhD.*
> > > >> > > > > > > > *Analysis & Nowcast Branch*
> > > >> > > > > > > > *Analysis and Mission Support Division*
> > > >> > > > > > > >
> > > >> > > > > > > > *Analyze, Forecast and Support Office*
> > > >> > > > > > > >
> > > >> > > > > > > > *National Weather Service**SSMC2, Room 14115*
> > > >> > > > > > > > *Silver Spring, Maryland*
> > > >> > > > > > > >
> > > >> > > > > > > >
> > > >> > > > > > > > *Ph. (301) 427-9500 *
> > > >> > > > > > > >
> > > >> > > > > > > >
> > > >> > > > > > >
> > > >> > > > > > > --
> > > >> > > > > > > Julie Prestopnik (she/her/hers)
> > > >> > > > > > > Software Engineer
> > > >> > > > > > > National Center for Atmospheric Research
> > > >> > > > > > > Research Applications Laboratory
> > > >> > > > > > > Phone: 303.497.8399
> > > >> > > > > > > Email: jpresto at ucar.edu
> > > >> > > > > > >
> > > >> > > > > > > My working day may not be your working day.
Please do
> not
> > > >> feel
> > > >> > > > obliged
> > > >> > > > > > to
> > > >> > > > > > > reply to this email outside of your normal
working
> hours.
> > > >> > > > > > >
> > > >> > > > > > >
> > > >> > > > > >
> > > >> > > > > > --
> > > >> > > > > >
> > > >> > > > > >
> > > >> > > > > > *Tatiana D. Gonzalez, PhD.*
> > > >> > > > > > *Analysis & Nowcast Branch*
> > > >> > > > > > *Analysis and Mission Support Division*
> > > >> > > > > >
> > > >> > > > > > *Analyze, Forecast and Support Office*
> > > >> > > > > >
> > > >> > > > > > *National Weather Service**SSMC2, Room 14115*
> > > >> > > > > > *Silver Spring, Maryland*
> > > >> > > > > >
> > > >> > > > > >
> > > >> > > > > > *Ph. (301) 427-9500 *
> > > >> > > > > >
> > > >> > > > > >
> > > >> > > > >
> > > >> > > > > --
> > > >> > > > > Julie Prestopnik (she/her/hers)
> > > >> > > > > Software Engineer
> > > >> > > > > National Center for Atmospheric Research
> > > >> > > > > Research Applications Laboratory
> > > >> > > > > Phone: 303.497.8399
> > > >> > > > > Email: jpresto at ucar.edu
> > > >> > > > >
> > > >> > > > > My working day may not be your working day. Please
do not
> > feel
> > > >> > obliged
> > > >> > > > to
> > > >> > > > > reply to this email outside of your normal working
hours.
> > > >> > > > >
> > > >> > > > >
> > > >> > > >
> > > >> > > > --
> > > >> > > >
> > > >> > > >
> > > >> > > > *Tatiana D. Gonzalez, PhD.*
> > > >> > > > *Analysis & Nowcast Branch*
> > > >> > > > *Analysis and Mission Support Division*
> > > >> > > >
> > > >> > > > *Analyze, Forecast and Support Office*
> > > >> > > >
> > > >> > > > *National Weather Service**SSMC2, Room 14115*
> > > >> > > > *Silver Spring, Maryland*
> > > >> > > >
> > > >> > > >
> > > >> > > > *Ph. (301) 427-9500 *
> > > >> > > >
> > > >> > > >
> > > >> > >
> > > >> > >
> > > >> >
> > > >> > --
> > > >> >
> > > >> >
> > > >> > *Tatiana D. Gonzalez, PhD.*
> > > >> > *Analysis & Nowcast Branch*
> > > >> > *Analysis and Mission Support Division*
> > > >> >
> > > >> > *Analyze, Forecast and Support Office*
> > > >> >
> > > >> > *National Weather Service**SSMC2, Room 14115*
> > > >> > *Silver Spring, Maryland*
> > > >> >
> > > >> >
> > > >> > *Ph. (301) 427-9500 *
> > > >> >
> > > >> >
> > > >>
> > > >> --
> > > >> Julie Prestopnik (she/her/hers)
> > > >> Software Engineer
> > > >> National Center for Atmospheric Research
> > > >> Research Applications Laboratory
> > > >> Phone: 303.497.8399
> > > >> Email: jpresto at ucar.edu
> > > >>
> > > >> My working day may not be your working day. Please do not
feel
> > obliged
> > > to
> > > >> reply to this email outside of your normal working hours.
> > > >>
> > > >>
> > > >
> > > > --
> > > >
> > > >
> > > > *Tatiana D. Gonzalez, PhD.*
> > > > *Analysis & Nowcast Branch*
> > > > *Analysis and Mission Support Division*
> > > >
> > > > *Analyze, Forecast and Support Office*
> > > >
> > > > *National Weather Service**SSMC2, Room 14115*
> > > > *Silver Spring, Maryland*
> > > >
> > > >
> > > > *Ph. (301) 427-9500 *
> > > >
> > > >
> > >
> > > --
> > >
> > >
> > > *Tatiana D. Gonzalez, PhD.*
> > > *Analysis & Nowcast Branch*
> > > *Analysis and Mission Support Division*
> > >
> > > *Analyze, Forecast and Support Office*
> > >
> > > *National Weather Service**SSMC2, Room 14115*
> > > *Silver Spring, Maryland*
> > >
> > >
> > > *Ph. (301) 427-9500 *
> > >
> > >
> >
> > --
> > Julie Prestopnik (she/her/hers)
> > Software Engineer
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
> > My working day may not be your working day. Please do not feel
obliged
> to
> > reply to this email outside of your normal working hours.
> >
> >
>
>
--
*Tatiana D. Gonzalez, PhD.*
*Analysis & Nowcast Branch*
*Analysis and Mission Support Division*
*Analyze, Forecast and Support Office*
*National Weather Service**SSMC2, Room 14115*
*Silver Spring, Maryland*
*Ph. (301) 427-9500 *
------------------------------------------------
More information about the Met_help
mailing list