[Met_help] [rt.rap.ucar.edu #100332] History for question on grid_stat masking

[John Halley Gotway via RT]

----------------------------------------------------------------
  Initial Request
----------------------------------------------------------------

Hi,

I'm trying to run the MET grid_stat on some RTOFS model data (regridded 
to cylindrical projection), and I'm having trouble specifying the 
verification masking in the user config file.  I've attached the user 
config file I'm using.

I specified the masking like so (copied out of the attached file)

mask = {
    llpnt   = [ { name       = "LAT30TO40";
                  lat_thresh = >=30&&<=40;
                  lon_thresh = NA; },
                { name       = "BOX";
                  lat_thresh = >=20&&<=40;
                  lon_thresh = >=-110&&<=-90; } ];
}

this is from one of the examples in the docs.

However, when I run grid_stat with this config file, it only uses the 
default full mask (run log attached)

DEBUG 2: Processing masking regions.
DEBUG 3: Processing grid mask: FULL
DEBUG 4: parse_grid_mask() -> parsing grid mask "FULL"

I'd like to be able to set up a set of regions for the major ocean 
basins.  Can you help?

-- 
Dr. Todd Spindler
IMSG at NOAA/NWS/NCEP/EMC
5830 University Research Ct., #2118
College Park, MD 20740
Todd.Spindler at noaa.gov
301-683-3757



----------------------------------------------------------------
  Complete Ticket History
----------------------------------------------------------------

Subject: question on grid_stat masking
From: Minna Win
Time: Fri Jun 25 09:39:29 2021

Hi Todd,

>From the documentation, it looks like llpnt and sid are applicable to
point-stat and ensemble-stat.

*The masking logic for processing point observations in Point-Stat and
Ensemble-Stat fall into two cateogries. The “sid” and “llpnt” options
apply
directly to the point observations. Only those observations for the
specified station id’s are included in the “sid” masks. Only those
observations meeting the latitude and longitude threshold criteria are
included in the “llpnt” masks.*


 Have you tried creating your own poly files and using this syntax:

mask = { grid=[];
               poly=[$METBASE/"MET_BASE/poly/LMV.poly"]
}

Leaving the grid array empty to not use any grid.

You can create your own poly files to delineate the major ocean
basins:

 *The "poly" entry contains a comma-separated list of files that
define
verification masking regions. These masking regions may be specified
in two
ways: as a lat/lon polygon or using a gridded data file such as the
NetCDF
output of the Gen-Vx-Mask tool.  An ASCII file containing a lat/lon
polygon. Latitude in degrees north and longitude in degrees east. The
first
and last polygon points are connected. e.g. "MET_BASE/poly/EAST.poly"
which
consists of n points:  "poly_name lat1 lon1 lat2 lon2... latn lonn"*

Regards,
Minna


---------------
Minna Win
Pronouns: she/her
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax:   303-497-8401
---------------



On Thu, Jun 24, 2021 at 3:55 PM Todd Spindler via RT
<met_help at ucar.edu>
wrote:

>
> Thu Jun 24 15:55:08 2021: Request 100332 was acted upon.
> Transaction: Ticket created by todd.spindler at noaa.gov
>        Queue: met_help
>      Subject: question on grid_stat masking
>        Owner: Nobody
>   Requestors: todd.spindler at noaa.gov
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100332 >
>
>
> Hi,
>
> I'm trying to run the MET grid_stat on some RTOFS model data
(regridded
> to cylindrical projection), and I'm having trouble specifying the
> verification masking in the user config file.  I've attached the
user
> config file I'm using.
>
> I specified the masking like so (copied out of the attached file)
>
> mask = {
>     llpnt   = [ { name       = "LAT30TO40";
>                   lat_thresh = >=30&&<=40;
>                   lon_thresh = NA; },
>                 { name       = "BOX";
>                   lat_thresh = >=20&&<=40;
>                   lon_thresh = >=-110&&<=-90; } ];
> }
>
> this is from one of the examples in the docs.
>
> However, when I run grid_stat with this config file, it only uses
the
> default full mask (run log attached)
>
> DEBUG 2: Processing masking regions.
> DEBUG 3: Processing grid mask: FULL
> DEBUG 4: parse_grid_mask() -> parsing grid mask "FULL"
>
> I'd like to be able to set up a set of regions for the major ocean
> basins.  Can you help?
>
> --
> Dr. Todd Spindler
> IMSG at NOAA/NWS/NCEP/EMC
> 5830 University Research Ct., #2118
> College Park, MD 20740
> Todd.Spindler at noaa.gov
> 301-683-3757
>
>
>

------------------------------------------------