[Met_help] [rt.rap.ucar.edu #88200] History for Grid Stat Climo File

John Halley Gotway via RT met_help at ucar.edu
Thu Jan 3 14:13:48 MST 2019 

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

John Gotway recently sent me a reply to Grid Stat question ( #88155).  I have on additional question.   I am using NCEP climo data to verify wind probability forecasts.  When I run MET, everything works except that it found 0 climatology files to match PROB(WIND>15)/z10/PROB.  The listing is below:

GSL_RNG_SEED=18446744072601605980
DEBUG 1: Forecast File: /h/data/global/ENSEMBLE/data/post/GEPS/global/2018112600/grib2.2018112600.0006.EPOSTOUT
DEBUG 1: Observation File: /h/data/global/WXQC/data/galwem/1811260600
DEBUG 3: Reading forecast data for PROB(WIND>15.000)/Z10/PROB.
DEBUG 3: MetGrib2DataFile::data_plane() - Found exact match for 'PROB(WIND>15.000)/Z10/PROB' in GRIB2 record 97 field 1 of GRIB2 file '/h/data/global/ENSEMBLE/data/post/GEPS/global/2018112600/grib2.2018112600.0006.EPOSTOUT'
DEBUG 3: Rescaling probabilistic field from [0,100] to [0,1].
DEBUG 3: Reading observation data for WIND/Z10.
DEBUG 3: MetGrib2DataFile::data_plane_array() - Found exact match for "WIND/Z10" in GRIB2 record 99 field 1 of GRIB2 file "/h/data/global/WXQC/data/galwem/1811260600"
DEBUG 3: MetGrib2DataFile::data_plane_array() - Found exact match for "WIND/Z10" in GRIB2 record 100 field 1 of GRIB2 file "/h/data/global/WXQC/data/galwem/1811260600"
DEBUG 3: Deriving wind speed from U and V wind components.
DEBUG 3: MetGrib2DataFile::data_plane() - Found derived match for 'WIND/Z10' in GRIB2 file '/h/data/global/WXQC/data/galwem/1811260600'
DEBUG 1: Regridding field WIND/Z10 to the verification grid.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Attempt to derive winds from U and V components.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading U-wind records.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match for VarInfo "WIND/Z10" in GRIB record 69 of GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match for VarInfo "WIND/Z10" in GRIB record 70 of GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match for VarInfo "WIND/Z10" in GRIB record 71 of GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match for VarInfo "WIND/Z10" in GRIB record 72 of GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB records matching VarInfo "WIND/Z10" in GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading V-wind records.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match for VarInfo "WIND/Z10" in GRIB record 73 of GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match for VarInfo "WIND/Z10" in GRIB record 74 of GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match for VarInfo "WIND/Z10" in GRIB record 75 of GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match for VarInfo "WIND/Z10" in GRIB record 76 of GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB records matching VarInfo "WIND/Z10" in GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
DEBUG 3: Deriving wind speed from U and V wind components.
DEBUG 3: Deriving wind speed from U and V wind components.
DEBUG 3: Deriving wind speed from U and V wind components.
DEBUG 3: Deriving wind speed from U and V wind components.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB records matching VarInfo "WIND/Z10" in GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
DEBUG 3: Found 0 climatology fields.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Attempt to derive winds from U and V components.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading U-wind records.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match for VarInfo "WIND/Z10" in GRIB record 69 of GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match for VarInfo "WIND/Z10" in GRIB record 70 of GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match for VarInfo "WIND/Z10" in GRIB record 71 of GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match for VarInfo "WIND/Z10" in GRIB record 72 of GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB records matching VarInfo "WIND/Z10" in GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading V-wind records.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match for VarInfo "WIND/Z10" in GRIB record 73 of GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match for VarInfo "WIND/Z10" in GRIB record 74 of GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match for VarInfo "WIND/Z10" in GRIB record 75 of GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match for VarInfo "WIND/Z10" in GRIB record 76 of GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB records matching VarInfo "WIND/Z10" in GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
DEBUG 3: Deriving wind speed from U and V wind components.
DEBUG 3: Deriving wind speed from U and V wind components.
DEBUG 3: Deriving wind speed from U and V wind components.
DEBUG 3: Deriving wind speed from U and V wind components.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB records matching VarInfo "WIND/Z10" in GRIB file "/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
DEBUG 3: Found 0 climatology fields.
DEBUG 3: Found 0 climatology mean and 0 climatology standard deviation field(s) for forecast PROB(WIND>15.000)/Z10/PROB.
DEBUG 2: Processing masking regions.
DEBUG 3: Processing grid mask: FULL

So, MET is reading the climo file just fine but does not find a match with the FCST Probability (WIND>15).   Do I have to convert the NCEP climo file in a probabilistic format where winds > 15 is 1 and winds < 15 is 0?  I need the climo field to get a Brier Skill Score.

Thanks
Bob Craig



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

Subject: Grid Stat Climo File
From: Minna Win
Time: Mon Dec 17 16:39:09 2018

Hi Bob,

John is out of the office for the next two weeks.  I will try and find
another content-expert who can handle your questions.

Thanks for your patience.

Regards,
Minna
---------------
Minna Win
NCAR
Research Applications Lab
Phone: 303-497-8423
Fax:   303-497-8401



On Mon, Dec 17, 2018 at 8:14 PM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:

>
> Mon Dec 17 13:13:44 2018: Request 88200 was acted upon.
> Transaction: Ticket created by robert.craig.2 at us.af.mil
>        Queue: met_help
>      Subject: Grid Stat Climo File
>        Owner: Nobody
>   Requestors: robert.craig.2 at us.af.mil
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88200 >
>
>
> John Gotway recently sent me a reply to Grid Stat question (
#88155).  I
> have on additional question.   I am using NCEP climo data to verify
wind
> probability forecasts.  When I run MET, everything works except that
it
> found 0 climatology files to match PROB(WIND>15)/z10/PROB.  The
listing is
> below:
>
> GSL_RNG_SEED=18446744072601605980
> DEBUG 1: Forecast File:
>
/h/data/global/ENSEMBLE/data/post/GEPS/global/2018112600/grib2.2018112600.0006.EPOSTOUT
> DEBUG 1: Observation File:
/h/data/global/WXQC/data/galwem/1811260600
> DEBUG 3: Reading forecast data for PROB(WIND>15.000)/Z10/PROB.
> DEBUG 3: MetGrib2DataFile::data_plane() - Found exact match for
> 'PROB(WIND>15.000)/Z10/PROB' in GRIB2 record 97 field 1 of GRIB2
file
>
'/h/data/global/ENSEMBLE/data/post/GEPS/global/2018112600/grib2.2018112600.0006.EPOSTOUT'
> DEBUG 3: Rescaling probabilistic field from [0,100] to [0,1].
> DEBUG 3: Reading observation data for WIND/Z10.
> DEBUG 3: MetGrib2DataFile::data_plane_array() - Found exact match
for
> "WIND/Z10" in GRIB2 record 99 field 1 of GRIB2 file
> "/h/data/global/WXQC/data/galwem/1811260600"
> DEBUG 3: MetGrib2DataFile::data_plane_array() - Found exact match
for
> "WIND/Z10" in GRIB2 record 100 field 1 of GRIB2 file
> "/h/data/global/WXQC/data/galwem/1811260600"
> DEBUG 3: Deriving wind speed from U and V wind components.
> DEBUG 3: MetGrib2DataFile::data_plane() - Found derived match for
> 'WIND/Z10' in GRIB2 file
'/h/data/global/WXQC/data/galwem/1811260600'
> DEBUG 1: Regridding field WIND/Z10 to the verification grid.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Attempt to derive
winds
> from U and V components.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading U-wind
records.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 69 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 70 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 71 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 72 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
records
> matching VarInfo "WIND/Z10" in GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading V-wind
records.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 73 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 74 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 75 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 76 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
records
> matching VarInfo "WIND/Z10" in GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: Deriving wind speed from U and V wind components.
> DEBUG 3: Deriving wind speed from U and V wind components.
> DEBUG 3: Deriving wind speed from U and V wind components.
> DEBUG 3: Deriving wind speed from U and V wind components.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
records
> matching VarInfo "WIND/Z10" in GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: Found 0 climatology fields.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Attempt to derive
winds
> from U and V components.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading U-wind
records.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 69 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 70 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 71 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 72 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
records
> matching VarInfo "WIND/Z10" in GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading V-wind
records.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 73 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 74 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 75 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 76 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
records
> matching VarInfo "WIND/Z10" in GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: Deriving wind speed from U and V wind components.
> DEBUG 3: Deriving wind speed from U and V wind components.
> DEBUG 3: Deriving wind speed from U and V wind components.
> DEBUG 3: Deriving wind speed from U and V wind components.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
records
> matching VarInfo "WIND/Z10" in GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: Found 0 climatology fields.
> DEBUG 3: Found 0 climatology mean and 0 climatology standard
deviation
> field(s) for forecast PROB(WIND>15.000)/Z10/PROB.
> DEBUG 2: Processing masking regions.
> DEBUG 3: Processing grid mask: FULL
>
> So, MET is reading the climo file just fine but does not find a
match with
> the FCST Probability (WIND>15).   Do I have to convert the NCEP
climo file
> in a probabilistic format where winds > 15 is 1 and winds < 15 is 0?
I
> need the climo field to get a Brier Skill Score.
>
> Thanks
> Bob Craig
>
>
>

------------------------------------------------
Subject: RE: [Non-DoD Source] Re: [rt.rap.ucar.edu #88200] Grid Stat Climo File
From: robert.craig.2 at us.af.mil
Time: Tue Dec 18 15:36:31 2018

Minna, That’s fine - I can wait for John if there is not an easy
answer to this.

Thanks
Bob

-----Original Message-----
From: Minna Win via RT <met_help at ucar.edu>
Sent: Monday, December 17, 2018 5:39 PM
To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #88200] Grid Stat Climo
File

Hi Bob,

John is out of the office for the next two weeks.  I will try and find
another content-expert who can handle your questions.

Thanks for your patience.

Regards,
Minna
---------------
Minna Win
NCAR
Research Applications Lab
Phone: 303-497-8423
Fax:   303-497-8401



On Mon, Dec 17, 2018 at 8:14 PM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:

>
> Mon Dec 17 13:13:44 2018: Request 88200 was acted upon.
> Transaction: Ticket created by robert.craig.2 at us.af.mil
>        Queue: met_help
>      Subject: Grid Stat Climo File
>        Owner: Nobody
>   Requestors: robert.craig.2 at us.af.mil
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88200
> >
>
>
> John Gotway recently sent me a reply to Grid Stat question (
#88155).  I
> have on additional question.   I am using NCEP climo data to verify
wind
> probability forecasts.  When I run MET, everything works except that
> it found 0 climatology files to match PROB(WIND>15)/z10/PROB.  The
> listing is
> below:
>
> GSL_RNG_SEED=18446744072601605980
> DEBUG 1: Forecast File:
>
/h/data/global/ENSEMBLE/data/post/GEPS/global/2018112600/grib2.2018112
> 600.0006.EPOSTOUT DEBUG 1: Observation File:
> /h/data/global/WXQC/data/galwem/1811260600
> DEBUG 3: Reading forecast data for PROB(WIND>15.000)/Z10/PROB.
> DEBUG 3: MetGrib2DataFile::data_plane() - Found exact match for
> 'PROB(WIND>15.000)/Z10/PROB' in GRIB2 record 97 field 1 of GRIB2
file
>
'/h/data/global/ENSEMBLE/data/post/GEPS/global/2018112600/grib2.2018112600.0006.EPOSTOUT'
> DEBUG 3: Rescaling probabilistic field from [0,100] to [0,1].
> DEBUG 3: Reading observation data for WIND/Z10.
> DEBUG 3: MetGrib2DataFile::data_plane_array() - Found exact match
for
> "WIND/Z10" in GRIB2 record 99 field 1 of GRIB2 file
> "/h/data/global/WXQC/data/galwem/1811260600"
> DEBUG 3: MetGrib2DataFile::data_plane_array() - Found exact match
for
> "WIND/Z10" in GRIB2 record 100 field 1 of GRIB2 file
> "/h/data/global/WXQC/data/galwem/1811260600"
> DEBUG 3: Deriving wind speed from U and V wind components.
> DEBUG 3: MetGrib2DataFile::data_plane() - Found derived match for
> 'WIND/Z10' in GRIB2 file
'/h/data/global/WXQC/data/galwem/1811260600'
> DEBUG 1: Regridding field WIND/Z10 to the verification grid.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Attempt to derive
> winds from U and V components.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading U-wind
records.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 69 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 70 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 71 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 72 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
records
> matching VarInfo "WIND/Z10" in GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading V-wind
records.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 73 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 74 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 75 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 76 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
records
> matching VarInfo "WIND/Z10" in GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: Deriving wind speed from U and V wind components.
> DEBUG 3: Deriving wind speed from U and V wind components.
> DEBUG 3: Deriving wind speed from U and V wind components.
> DEBUG 3: Deriving wind speed from U and V wind components.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
records
> matching VarInfo "WIND/Z10" in GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> DEBUG 3: Found 0 climatology fields.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Attempt to derive
> winds from U and V components.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading U-wind
records.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 69 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 70 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 71 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 72 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
records
> matching VarInfo "WIND/Z10" in GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading V-wind
records.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 73 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 74 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 75 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> VarInfo "WIND/Z10" in GRIB record 76 of GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
records
> matching VarInfo "WIND/Z10" in GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: Deriving wind speed from U and V wind components.
> DEBUG 3: Deriving wind speed from U and V wind components.
> DEBUG 3: Deriving wind speed from U and V wind components.
> DEBUG 3: Deriving wind speed from U and V wind components.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
records
> matching VarInfo "WIND/Z10" in GRIB file
>
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> DEBUG 3: Found 0 climatology fields.
> DEBUG 3: Found 0 climatology mean and 0 climatology standard
deviation
> field(s) for forecast PROB(WIND>15.000)/Z10/PROB.
> DEBUG 2: Processing masking regions.
> DEBUG 3: Processing grid mask: FULL
>
> So, MET is reading the climo file just fine but does not find a
match with
> the FCST Probability (WIND>15).   Do I have to convert the NCEP
climo file
> in a probabilistic format where winds > 15 is 1 and winds < 15 is 0?
> I need the climo field to get a Brier Skill Score.
>
> Thanks
> Bob Craig
>
>
>



------------------------------------------------
Subject: Grid Stat Climo File
From: Minna Win
Time: Wed Dec 19 11:51:21 2018

Hi Bob,

Thanks for your patience.  I checked with a statistician and it looks
like
this is a question for John.  I'm assigning this ticket to him so he
can
take a look at it when he returns.

Thanks,
Minna
---------------
Minna Win
NCAR
Research Applications Lab
Phone: 303-497-8423
Fax:   303-497-8401



On Tue, Dec 18, 2018 at 10:37 PM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88200 >
>
> Minna, That’s fine - I can wait for John if there is not an easy
answer to
> this.
>
> Thanks
> Bob
>
> -----Original Message-----
> From: Minna Win via RT <met_help at ucar.edu>
> Sent: Monday, December 17, 2018 5:39 PM
> To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
> Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #88200] Grid Stat
Climo
> File
>
> Hi Bob,
>
> John is out of the office for the next two weeks.  I will try and
find
> another content-expert who can handle your questions.
>
> Thanks for your patience.
>
> Regards,
> Minna
> ---------------
> Minna Win
> NCAR
> Research Applications Lab
> Phone: 303-497-8423
> Fax:   303-497-8401
>
>
>
> On Mon, Dec 17, 2018 at 8:14 PM robert.craig.2 at us.af.mil via RT <
> met_help at ucar.edu> wrote:
>
> >
> > Mon Dec 17 13:13:44 2018: Request 88200 was acted upon.
> > Transaction: Ticket created by robert.craig.2 at us.af.mil
> >        Queue: met_help
> >      Subject: Grid Stat Climo File
> >        Owner: Nobody
> >   Requestors: robert.craig.2 at us.af.mil
> >       Status: new
> >  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88200
> > >
> >
> >
> > John Gotway recently sent me a reply to Grid Stat question (
#88155).  I
> > have on additional question.   I am using NCEP climo data to
verify wind
> > probability forecasts.  When I run MET, everything works except
that
> > it found 0 climatology files to match PROB(WIND>15)/z10/PROB.  The
> > listing is
> > below:
> >
> > GSL_RNG_SEED=18446744072601605980
> > DEBUG 1: Forecast File:
> >
/h/data/global/ENSEMBLE/data/post/GEPS/global/2018112600/grib2.2018112
> > 600.0006.EPOSTOUT DEBUG 1: Observation File:
> > /h/data/global/WXQC/data/galwem/1811260600
> > DEBUG 3: Reading forecast data for PROB(WIND>15.000)/Z10/PROB.
> > DEBUG 3: MetGrib2DataFile::data_plane() - Found exact match for
> > 'PROB(WIND>15.000)/Z10/PROB' in GRIB2 record 97 field 1 of GRIB2
file
> >
>
'/h/data/global/ENSEMBLE/data/post/GEPS/global/2018112600/grib2.2018112600.0006.EPOSTOUT'
> > DEBUG 3: Rescaling probabilistic field from [0,100] to [0,1].
> > DEBUG 3: Reading observation data for WIND/Z10.
> > DEBUG 3: MetGrib2DataFile::data_plane_array() - Found exact match
for
> > "WIND/Z10" in GRIB2 record 99 field 1 of GRIB2 file
> > "/h/data/global/WXQC/data/galwem/1811260600"
> > DEBUG 3: MetGrib2DataFile::data_plane_array() - Found exact match
for
> > "WIND/Z10" in GRIB2 record 100 field 1 of GRIB2 file
> > "/h/data/global/WXQC/data/galwem/1811260600"
> > DEBUG 3: Deriving wind speed from U and V wind components.
> > DEBUG 3: MetGrib2DataFile::data_plane() - Found derived match for
> > 'WIND/Z10' in GRIB2 file
'/h/data/global/WXQC/data/galwem/1811260600'
> > DEBUG 1: Regridding field WIND/Z10 to the verification grid.
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Attempt to derive
> > winds from U and V components.
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading U-wind
records.
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> > VarInfo "WIND/Z10" in GRIB record 69 of GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> > VarInfo "WIND/Z10" in GRIB record 70 of GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> > VarInfo "WIND/Z10" in GRIB record 71 of GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> > VarInfo "WIND/Z10" in GRIB record 72 of GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
records
> > matching VarInfo "WIND/Z10" in GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading V-wind
records.
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> > VarInfo "WIND/Z10" in GRIB record 73 of GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> > VarInfo "WIND/Z10" in GRIB record 74 of GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> > VarInfo "WIND/Z10" in GRIB record 75 of GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> > VarInfo "WIND/Z10" in GRIB record 76 of GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
records
> > matching VarInfo "WIND/Z10" in GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> > DEBUG 3: Deriving wind speed from U and V wind components.
> > DEBUG 3: Deriving wind speed from U and V wind components.
> > DEBUG 3: Deriving wind speed from U and V wind components.
> > DEBUG 3: Deriving wind speed from U and V wind components.
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
records
> > matching VarInfo "WIND/Z10" in GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
> > DEBUG 3: Found 0 climatology fields.
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Attempt to derive
> > winds from U and V components.
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading U-wind
records.
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> > VarInfo "WIND/Z10" in GRIB record 69 of GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> > VarInfo "WIND/Z10" in GRIB record 70 of GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> > VarInfo "WIND/Z10" in GRIB record 71 of GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> > VarInfo "WIND/Z10" in GRIB record 72 of GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
records
> > matching VarInfo "WIND/Z10" in GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading V-wind
records.
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> > VarInfo "WIND/Z10" in GRIB record 73 of GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> > VarInfo "WIND/Z10" in GRIB record 74 of GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> > VarInfo "WIND/Z10" in GRIB record 75 of GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match
for
> > VarInfo "WIND/Z10" in GRIB record 76 of GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
records
> > matching VarInfo "WIND/Z10" in GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> > DEBUG 3: Deriving wind speed from U and V wind components.
> > DEBUG 3: Deriving wind speed from U and V wind components.
> > DEBUG 3: Deriving wind speed from U and V wind components.
> > DEBUG 3: Deriving wind speed from U and V wind components.
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
records
> > matching VarInfo "WIND/Z10" in GRIB file
> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
> > DEBUG 3: Found 0 climatology fields.
> > DEBUG 3: Found 0 climatology mean and 0 climatology standard
deviation
> > field(s) for forecast PROB(WIND>15.000)/Z10/PROB.
> > DEBUG 2: Processing masking regions.
> > DEBUG 3: Processing grid mask: FULL
> >
> > So, MET is reading the climo file just fine but does not find a
match
> with
> > the FCST Probability (WIND>15).   Do I have to convert the NCEP
climo
> file
> > in a probabilistic format where winds > 15 is 1 and winds < 15 is
0?
> > I need the climo field to get a Brier Skill Score.
> >
> > Thanks
> > Bob Craig
> >
> >
> >
>
>
>
>

------------------------------------------------
Subject: Grid Stat Climo File
From: robert.craig.2 at us.af.mil
Time: Thu Dec 20 11:44:26 2018

Minna, I figured out the issue I was having - I forgot to set the
climo match_data field to false.  So that ticket can be closed.



A new issue for John or anyone else is a problem where I am degribbing
a GRIB 2 file of forecast probabilities.  I am using the notation
below in the grid stat configuration file for get 500mb HGT
probabilities.



fcst = {



   field = [

             { name = "PROB"; level = "P500";

               prob = { name = "HGT"; thresh_hi = 5400; } },

             { name = "PROB"; level = "P700";

               prob = { name = "HGT"; thresh_hi = 3000; } },

             { name = "PROB"; level = "P500";

               prob = { name = "HGT"; thresh_lo = 5800; } }



           ];

           cat_thresh=[==0.1];

}



In the grib2 file, the table of contents is listed below, there are
multiple 500MB HGT entries which include the probability I am after.
The // are data I left out for clarity.  I highlighted the record I am
after.



1:0:d=2018112600:PRES:surface:6 hour fcst:mean all members

2:552511:d=2018112600:PRES:surface:6 hour fcst:std dev (cluster mean)

3:975062:d=2018112600:HGT:200 mb:6 hour fcst:mean all members

4:1365123:d=2018112600:HGT:250 mb:6 hour fcst:mean all members

5:1755184:d=2018112600:HGT:500 mb:6 hour fcst:mean all members

6:2112755:d=2018112600:HGT:700 mb:6 hour fcst:mean all members

7:2437836:d=2018112600:HGT:850 mb:6 hour fcst:mean all members

8:2762917:d=2018112600:HGT:925 mb:6 hour fcst:mean all members

//

12:3803321:d=2018112600:HGT:500 mb:6 hour fcst:std dev (cluster mean)

//

148:29723303:d=2018112600:TMP:2 m above ground:24 hour fcst:prob
<273.15:probability forecast

149:29950925:d=2018112600:TMP:2 m above ground:24 hour fcst:prob
>305.37:probability forecast

150:30178547:d=2018112600:HGT:500 mb:24 hour fcst:prob
>5800:probability forecast

151:30406169:d=2018112600:GEOLAT:surface:24 hour fcst:

152:30666268:d=2018112600:GEOLON:surface:24 hour fcst:

153:30958857:d=2018112600:HGT:500 mb:24 hour fcst:prob
<5400:probability forecast

154:31186479:d=2018112600:HGT:700 mb:24 hour fcst:prob
<3000:probability forecast

155:31414101:d=2018112600:PRMSL:mean sea level:24 hour fcst:prob
>102500:probability forecast



Below is the MET out listing for this:



DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_Gb2".

DEBUG 3:

DEBUG 3: Multiple GRIB2 table entries match lookup criteria (parm_name
= HGT):

DEBUG 3:   parm_name: HGT, index_a = 0, grib2_mtab = 1, grib2_cntr =
0, grib2_ltab = 0, index_b = 3, index_c = 5

DEBUG 3:   parm_name: HGT, index_a = 0, grib2_mtab = 14, grib2_cntr =
57, grib2_ltab = 0, index_b = 3, index_c = 5

DEBUG 3: Using the first match found:   parm_name: HGT, index_a = 0,
grib2_mtab = 1, grib2_cntr = 0, grib2_ltab = 0, index_b = 3, index_c =
5

DEBUG 3:

DEBUG 3:

DEBUG 3: Multiple GRIB2 table entries match lookup criteria (parm_name
= HGT):

DEBUG 3:   parm_name: HGT, index_a = 0, grib2_mtab = 1, grib2_cntr =
0, grib2_ltab = 0, index_b = 3, index_c = 5

DEBUG 3:   parm_name: HGT, index_a = 0, grib2_mtab = 14, grib2_cntr =
57, grib2_ltab = 0, index_b = 3, index_c = 5

DEBUG 3: Using the first match found:   parm_name: HGT, index_a = 0,
grib2_mtab = 1, grib2_cntr = 0, grib2_ltab = 0, index_b = 3, index_c =
5

DEBUG 3:

DEBUG 3: Use the forecast grid.

DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 720 Ny: 361 lat_ll:
-90.000 lon_ll: -0.000 delta_lat: 0.500 delta_lon: 0.500

GSL_RNG_TYPE=mt19937

GSL_RNG_SEED=249568092

DEBUG 1: Forecast File:
/h/data/global/ENSEMBLE/data/post/GEPS/global/2018112600/grib2.2018112600.0024.EPOSTOUT

DEBUG 1: Observation File: /h/data/global/WXQC/data/galwem/1811270000

DEBUG 3: Reading forecast data for PROB(HGT<5400.000)/P500/PROB.

WARNING:

WARNING: MetGrib2DataFile::data_plane() - Multiple matching records
found for 'PROB(HGT<5400.000)/P500/PROB' - 5, 12, 153, 213 - using 5

WARNING:

DEBUG 3: MetGrib2DataFile::data_plane() - Found exact match for
'PROB(HGT<5400.000)/P500/PROB' in GRIB2 record 5 field 1 of GRIB2 file
'/h/data/global/ENSEMBLE/data/post/GEPS/global/2018112600/grib2.2018112600.0024.EPOSTOUT'



Notice it finds multiple records and chooses the first on it sees (5)
instead of the correct one (153).    Record 5 and 12 are not
probability forecasts so why is it considering them?  If I pull using
the record number I will get the correct record but this can’t be done
for automated runs.



Thanks for your help.



Bob









-----Original Message-----
From: Minna Win via RT <met_help at ucar.edu>
Sent: Wednesday, December 19, 2018 12:51 PM
To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
Subject: Re: [Non-DoD Source] Re: [rt.rap.ucar.edu #88200] Grid Stat
Climo File



Hi Bob,



Thanks for your patience.  I checked with a statistician and it looks
like this is a question for John.  I'm assigning this ticket to him so
he can take a look at it when he returns.



Thanks,

Minna

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

Minna Win

NCAR

Research Applications Lab

Phone: 303-497-8423

Fax:   303-497-8401







On Tue, Dec 18, 2018 at 10:37 PM  <mailto:robert.craig.2 at us.af.mil>
robert.craig.2 at us.af.mil via RT <  <mailto:met_help at ucar.edu>
met_help at ucar.edu> wrote:



>

> <URL:  <https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88200>
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88200 >

>

> Minna, That’s fine - I can wait for John if there is not an easy

> answer to this.

>

> Thanks

> Bob

>

> -----Original Message-----

> From: Minna Win via RT < <mailto:met_help at ucar.edu>
met_help at ucar.edu>

> Sent: Monday, December 17, 2018 5:39 PM

> To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN

> < <mailto:robert.craig.2 at us.af.mil> robert.craig.2 at us.af.mil>

> Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #88200] Grid Stat
Climo

> File

>

> Hi Bob,

>

> John is out of the office for the next two weeks.  I will try and
find

> another content-expert who can handle your questions.

>

> Thanks for your patience.

>

> Regards,

> Minna

> ---------------

> Minna Win

> NCAR

> Research Applications Lab

> Phone: 303-497-8423

> Fax:   303-497-8401

>

>

>

> On Mon, Dec 17, 2018 at 8:14 PM  <mailto:robert.craig.2 at us.af.mil>
robert.craig.2 at us.af.mil via RT <

>  <mailto:met_help at ucar.edu> met_help at ucar.edu> wrote:

>

> >

> > Mon Dec 17 13:13:44 2018: Request 88200 was acted upon.

> > Transaction: Ticket created by  <mailto:robert.craig.2 at us.af.mil>
robert.craig.2 at us.af.mil

> >        Queue: met_help

> >      Subject: Grid Stat Climo File

> >        Owner: Nobody

> >   Requestors:  <mailto:robert.craig.2 at us.af.mil>
robert.craig.2 at us.af.mil

> >       Status: new

> >  Ticket <URL:

> >  <https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88200>
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88200

> > >

> >

> >

> > John Gotway recently sent me a reply to Grid Stat question (
#88155).  I

> > have on additional question.   I am using NCEP climo data to
verify wind

> > probability forecasts.  When I run MET, everything works except
that

> > it found 0 climatology files to match PROB(WIND>15)/z10/PROB.  The

> > listing is

> > below:

> >

> > GSL_RNG_SEED=18446744072601605980

> > DEBUG 1: Forecast File:

> >
/h/data/global/ENSEMBLE/data/post/GEPS/global/2018112600/grib2.20181

> > 12 600.0006.EPOSTOUT DEBUG 1: Observation File:

> > /h/data/global/WXQC/data/galwem/1811260600

> > DEBUG 3: Reading forecast data for PROB(WIND>15.000)/Z10/PROB.

> > DEBUG 3: MetGrib2DataFile::data_plane() - Found exact match for

> > 'PROB(WIND>15.000)/Z10/PROB' in GRIB2 record 97 field 1 of GRIB2

> > file

> >

>
'/h/data/global/ENSEMBLE/data/post/GEPS/global/2018112600/grib2.2018112600.0006.EPOSTOUT'

> > DEBUG 3: Rescaling probabilistic field from [0,100] to [0,1].

> > DEBUG 3: Reading observation data for WIND/Z10.

> > DEBUG 3: MetGrib2DataFile::data_plane_array() - Found exact match

> > for "WIND/Z10" in GRIB2 record 99 field 1 of GRIB2 file

> > "/h/data/global/WXQC/data/galwem/1811260600"

> > DEBUG 3: MetGrib2DataFile::data_plane_array() - Found exact match

> > for "WIND/Z10" in GRIB2 record 100 field 1 of GRIB2 file

> > "/h/data/global/WXQC/data/galwem/1811260600"

> > DEBUG 3: Deriving wind speed from U and V wind components.

> > DEBUG 3: MetGrib2DataFile::data_plane() - Found derived match for

> > 'WIND/Z10' in GRIB2 file
'/h/data/global/WXQC/data/galwem/1811260600'

> > DEBUG 1: Regridding field WIND/Z10 to the verification grid.

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Attempt to derive

> > winds from U and V components.

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading U-wind
records.

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match

> > for VarInfo "WIND/Z10" in GRIB record 69 of GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match

> > for VarInfo "WIND/Z10" in GRIB record 70 of GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match

> > for VarInfo "WIND/Z10" in GRIB record 71 of GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match

> > for VarInfo "WIND/Z10" in GRIB record 72 of GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB

> > records matching VarInfo "WIND/Z10" in GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading V-wind
records.

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match

> > for VarInfo "WIND/Z10" in GRIB record 73 of GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match

> > for VarInfo "WIND/Z10" in GRIB record 74 of GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match

> > for VarInfo "WIND/Z10" in GRIB record 75 of GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match

> > for VarInfo "WIND/Z10" in GRIB record 76 of GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB

> > records matching VarInfo "WIND/Z10" in GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".

> > DEBUG 3: Deriving wind speed from U and V wind components.

> > DEBUG 3: Deriving wind speed from U and V wind components.

> > DEBUG 3: Deriving wind speed from U and V wind components.

> > DEBUG 3: Deriving wind speed from U and V wind components.

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB

> > records matching VarInfo "WIND/Z10" in GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".

> > DEBUG 3: Found 0 climatology fields.

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Attempt to derive

> > winds from U and V components.

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading U-wind
records.

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match

> > for VarInfo "WIND/Z10" in GRIB record 69 of GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match

> > for VarInfo "WIND/Z10" in GRIB record 70 of GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match

> > for VarInfo "WIND/Z10" in GRIB record 71 of GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match

> > for VarInfo "WIND/Z10" in GRIB record 72 of GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB

> > records matching VarInfo "WIND/Z10" in GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading V-wind
records.

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match

> > for VarInfo "WIND/Z10" in GRIB record 73 of GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match

> > for VarInfo "WIND/Z10" in GRIB record 74 of GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match

> > for VarInfo "WIND/Z10" in GRIB record 75 of GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact match

> > for VarInfo "WIND/Z10" in GRIB record 76 of GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB

> > records matching VarInfo "WIND/Z10" in GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".

> > DEBUG 3: Deriving wind speed from U and V wind components.

> > DEBUG 3: Deriving wind speed from U and V wind components.

> > DEBUG 3: Deriving wind speed from U and V wind components.

> > DEBUG 3: Deriving wind speed from U and V wind components.

> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB

> > records matching VarInfo "WIND/Z10" in GRIB file

> >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".

> > DEBUG 3: Found 0 climatology fields.

> > DEBUG 3: Found 0 climatology mean and 0 climatology standard

> > deviation

> > field(s) for forecast PROB(WIND>15.000)/Z10/PROB.

> > DEBUG 2: Processing masking regions.

> > DEBUG 3: Processing grid mask: FULL

> >

> > So, MET is reading the climo file just fine but does not find a

> > match

> with

> > the FCST Probability (WIND>15).   Do I have to convert the NCEP
climo

> file

> > in a probabilistic format where winds > 15 is 1 and winds < 15 is
0?

> > I need the climo field to get a Brier Skill Score.

> >

> > Thanks

> > Bob Craig

> >

> >

> >

>

>

>

>




------------------------------------------------
Subject: Grid Stat Climo File
From: Minna Win
Time: Thu Dec 20 11:46:22 2018

Hi Bob,

Thanks for the update.

I'll close this ticket.

Regards,
Minna
---------------
Minna Win
NCAR
Research Applications Lab
Phone: 303-497-8423
Fax:   303-497-8401



On Thu, Dec 20, 2018 at 6:44 PM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88200 >
>
> Minna, I figured out the issue I was having - I forgot to set the
climo
> match_data field to false.  So that ticket can be closed.
>
>
>
> A new issue for John or anyone else is a problem where I am
degribbing a
> GRIB 2 file of forecast probabilities.  I am using the notation
below in
> the grid stat configuration file for get 500mb HGT probabilities.
>
>
>
> fcst = {
>
>
>
>    field = [
>
>              { name = "PROB"; level = "P500";
>
>                prob = { name = "HGT"; thresh_hi = 5400; } },
>
>              { name = "PROB"; level = "P700";
>
>                prob = { name = "HGT"; thresh_hi = 3000; } },
>
>              { name = "PROB"; level = "P500";
>
>                prob = { name = "HGT"; thresh_lo = 5800; } }
>
>
>
>            ];
>
>            cat_thresh=[==0.1];
>
> }
>
>
>
> In the grib2 file, the table of contents is listed below, there are
> multiple 500MB HGT entries which include the probability I am after.
The
> // are data I left out for clarity.  I highlighted the record I am
after.
>
>
>
> 1:0:d=2018112600:PRES:surface:6 hour fcst:mean all members
>
> 2:552511:d=2018112600:PRES:surface:6 hour fcst:std dev (cluster
mean)
>
> 3:975062:d=2018112600:HGT:200 mb:6 hour fcst:mean all members
>
> 4:1365123:d=2018112600:HGT:250 mb:6 hour fcst:mean all members
>
> 5:1755184:d=2018112600:HGT:500 mb:6 hour fcst:mean all members
>
> 6:2112755:d=2018112600:HGT:700 mb:6 hour fcst:mean all members
>
> 7:2437836:d=2018112600:HGT:850 mb:6 hour fcst:mean all members
>
> 8:2762917:d=2018112600:HGT:925 mb:6 hour fcst:mean all members
>
> //
>
> 12:3803321:d=2018112600:HGT:500 mb:6 hour fcst:std dev (cluster
mean)
>
> //
>
> 148:29723303:d=2018112600:TMP:2 m above ground:24 hour fcst:prob
> <273.15:probability forecast
>
> 149:29950925:d=2018112600:TMP:2 m above ground:24 hour fcst:prob
> >305.37:probability forecast
>
> 150:30178547:d=2018112600:HGT:500 mb:24 hour fcst:prob
>5800:probability
> forecast
>
> 151:30406169:d=2018112600:GEOLAT:surface:24 hour fcst:
>
> 152:30666268:d=2018112600:GEOLON:surface:24 hour fcst:
>
> 153:30958857:d=2018112600:HGT:500 mb:24 hour fcst:prob
<5400:probability
> forecast
>
> 154:31186479:d=2018112600:HGT:700 mb:24 hour fcst:prob
<3000:probability
> forecast
>
> 155:31414101:d=2018112600:PRMSL:mean sea level:24 hour fcst:prob
> >102500:probability forecast
>
>
>
> Below is the MET out listing for this:
>
>
>
> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> type "FileType_Gb2".
>
> DEBUG 3:
>
> DEBUG 3: Multiple GRIB2 table entries match lookup criteria
(parm_name =
> HGT):
>
> DEBUG 3:   parm_name: HGT, index_a = 0, grib2_mtab = 1, grib2_cntr =
0,
> grib2_ltab = 0, index_b = 3, index_c = 5
>
> DEBUG 3:   parm_name: HGT, index_a = 0, grib2_mtab = 14, grib2_cntr
= 57,
> grib2_ltab = 0, index_b = 3, index_c = 5
>
> DEBUG 3: Using the first match found:   parm_name: HGT, index_a = 0,
> grib2_mtab = 1, grib2_cntr = 0, grib2_ltab = 0, index_b = 3, index_c
= 5
>
> DEBUG 3:
>
> DEBUG 3:
>
> DEBUG 3: Multiple GRIB2 table entries match lookup criteria
(parm_name =
> HGT):
>
> DEBUG 3:   parm_name: HGT, index_a = 0, grib2_mtab = 1, grib2_cntr =
0,
> grib2_ltab = 0, index_b = 3, index_c = 5
>
> DEBUG 3:   parm_name: HGT, index_a = 0, grib2_mtab = 14, grib2_cntr
= 57,
> grib2_ltab = 0, index_b = 3, index_c = 5
>
> DEBUG 3: Using the first match found:   parm_name: HGT, index_a = 0,
> grib2_mtab = 1, grib2_cntr = 0, grib2_ltab = 0, index_b = 3, index_c
= 5
>
> DEBUG 3:
>
> DEBUG 3: Use the forecast grid.
>
> DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 720 Ny: 361
lat_ll:
> -90.000 lon_ll: -0.000 delta_lat: 0.500 delta_lon: 0.500
>
> GSL_RNG_TYPE=mt19937
>
> GSL_RNG_SEED=249568092
>
> DEBUG 1: Forecast File:
>
/h/data/global/ENSEMBLE/data/post/GEPS/global/2018112600/grib2.2018112600.0024.EPOSTOUT
>
> DEBUG 1: Observation File:
/h/data/global/WXQC/data/galwem/1811270000
>
> DEBUG 3: Reading forecast data for PROB(HGT<5400.000)/P500/PROB.
>
> WARNING:
>
> WARNING: MetGrib2DataFile::data_plane() - Multiple matching records
found
> for 'PROB(HGT<5400.000)/P500/PROB' - 5, 12, 153, 213 - using 5
>
> WARNING:
>
> DEBUG 3: MetGrib2DataFile::data_plane() - Found exact match for
> 'PROB(HGT<5400.000)/P500/PROB' in GRIB2 record 5 field 1 of GRIB2
file
>
'/h/data/global/ENSEMBLE/data/post/GEPS/global/2018112600/grib2.2018112600.0024.EPOSTOUT'
>
>
>
> Notice it finds multiple records and chooses the first on it sees
(5)
> instead of the correct one (153).    Record 5 and 12 are not
probability
> forecasts so why is it considering them?  If I pull using the record
number
> I will get the correct record but this can’t be done for automated
runs.
>
>
>
> Thanks for your help.
>
>
>
> Bob
>
>
>
>
>
>
>
>
>
> -----Original Message-----
> From: Minna Win via RT <met_help at ucar.edu>
> Sent: Wednesday, December 19, 2018 12:51 PM
> To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
> Subject: Re: [Non-DoD Source] Re: [rt.rap.ucar.edu #88200] Grid Stat
> Climo File
>
>
>
> Hi Bob,
>
>
>
> Thanks for your patience.  I checked with a statistician and it
looks like
> this is a question for John.  I'm assigning this ticket to him so he
can
> take a look at it when he returns.
>
>
>
> Thanks,
>
> Minna
>
> ---------------
>
> Minna Win
>
> NCAR
>
> Research Applications Lab
>
> Phone: 303-497-8423
>
> Fax:   303-497-8401
>
>
>
>
>
>
>
> On Tue, Dec 18, 2018 at 10:37 PM  <mailto:robert.craig.2 at us.af.mil>
> robert.craig.2 at us.af.mil via RT <  <mailto:met_help at ucar.edu>
> met_help at ucar.edu> wrote:
>
>
>
> >
>
> > <URL:  <https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88200>
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88200 >
>
> >
>
> > Minna, That’s fine - I can wait for John if there is not an easy
>
> > answer to this.
>
> >
>
> > Thanks
>
> > Bob
>
> >
>
> > -----Original Message-----
>
> > From: Minna Win via RT < <mailto:met_help at ucar.edu>
met_help at ucar.edu>
>
> > Sent: Monday, December 17, 2018 5:39 PM
>
> > To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
>
> > < <mailto:robert.craig.2 at us.af.mil> robert.craig.2 at us.af.mil>
>
> > Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #88200] Grid Stat
Climo
>
> > File
>
> >
>
> > Hi Bob,
>
> >
>
> > John is out of the office for the next two weeks.  I will try and
find
>
> > another content-expert who can handle your questions.
>
> >
>
> > Thanks for your patience.
>
> >
>
> > Regards,
>
> > Minna
>
> > ---------------
>
> > Minna Win
>
> > NCAR
>
> > Research Applications Lab
>
> > Phone: 303-497-8423
>
> > Fax:   303-497-8401
>
> >
>
> >
>
> >
>
> > On Mon, Dec 17, 2018 at 8:14 PM  <mailto:robert.craig.2 at us.af.mil>
> robert.craig.2 at us.af.mil via RT <
>
> >  <mailto:met_help at ucar.edu> met_help at ucar.edu> wrote:
>
> >
>
> > >
>
> > > Mon Dec 17 13:13:44 2018: Request 88200 was acted upon.
>
> > > Transaction: Ticket created by
<mailto:robert.craig.2 at us.af.mil>
> robert.craig.2 at us.af.mil
>
> > >        Queue: met_help
>
> > >      Subject: Grid Stat Climo File
>
> > >        Owner: Nobody
>
> > >   Requestors:  <mailto:robert.craig.2 at us.af.mil>
> robert.craig.2 at us.af.mil
>
> > >       Status: new
>
> > >  Ticket <URL:
>
> > >  <https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88200>
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88200
>
> > > >
>
> > >
>
> > >
>
> > > John Gotway recently sent me a reply to Grid Stat question (
#88155).
> I
>
> > > have on additional question.   I am using NCEP climo data to
verify
> wind
>
> > > probability forecasts.  When I run MET, everything works except
that
>
> > > it found 0 climatology files to match PROB(WIND>15)/z10/PROB.
The
>
> > > listing is
>
> > > below:
>
> > >
>
> > > GSL_RNG_SEED=18446744072601605980
>
> > > DEBUG 1: Forecast File:
>
> > >
/h/data/global/ENSEMBLE/data/post/GEPS/global/2018112600/grib2.20181
>
> > > 12 600.0006.EPOSTOUT DEBUG 1: Observation File:
>
> > > /h/data/global/WXQC/data/galwem/1811260600
>
> > > DEBUG 3: Reading forecast data for PROB(WIND>15.000)/Z10/PROB.
>
> > > DEBUG 3: MetGrib2DataFile::data_plane() - Found exact match for
>
> > > 'PROB(WIND>15.000)/Z10/PROB' in GRIB2 record 97 field 1 of GRIB2
>
> > > file
>
> > >
>
> >
>
'/h/data/global/ENSEMBLE/data/post/GEPS/global/2018112600/grib2.2018112600.0006.EPOSTOUT'
>
> > > DEBUG 3: Rescaling probabilistic field from [0,100] to [0,1].
>
> > > DEBUG 3: Reading observation data for WIND/Z10.
>
> > > DEBUG 3: MetGrib2DataFile::data_plane_array() - Found exact
match
>
> > > for "WIND/Z10" in GRIB2 record 99 field 1 of GRIB2 file
>
> > > "/h/data/global/WXQC/data/galwem/1811260600"
>
> > > DEBUG 3: MetGrib2DataFile::data_plane_array() - Found exact
match
>
> > > for "WIND/Z10" in GRIB2 record 100 field 1 of GRIB2 file
>
> > > "/h/data/global/WXQC/data/galwem/1811260600"
>
> > > DEBUG 3: Deriving wind speed from U and V wind components.
>
> > > DEBUG 3: MetGrib2DataFile::data_plane() - Found derived match
for
>
> > > 'WIND/Z10' in GRIB2 file
'/h/data/global/WXQC/data/galwem/1811260600'
>
> > > DEBUG 1: Regridding field WIND/Z10 to the verification grid.
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Attempt to
derive
>
> > > winds from U and V components.
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading U-wind
> records.
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact
match
>
> > > for VarInfo "WIND/Z10" in GRIB record 69 of GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact
match
>
> > > for VarInfo "WIND/Z10" in GRIB record 70 of GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact
match
>
> > > for VarInfo "WIND/Z10" in GRIB record 71 of GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact
match
>
> > > for VarInfo "WIND/Z10" in GRIB record 72 of GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
>
> > > records matching VarInfo "WIND/Z10" in GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading V-wind
> records.
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact
match
>
> > > for VarInfo "WIND/Z10" in GRIB record 73 of GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact
match
>
> > > for VarInfo "WIND/Z10" in GRIB record 74 of GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact
match
>
> > > for VarInfo "WIND/Z10" in GRIB record 75 of GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact
match
>
> > > for VarInfo "WIND/Z10" in GRIB record 76 of GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
>
> > > records matching VarInfo "WIND/Z10" in GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
>
> > > DEBUG 3: Deriving wind speed from U and V wind components.
>
> > > DEBUG 3: Deriving wind speed from U and V wind components.
>
> > > DEBUG 3: Deriving wind speed from U and V wind components.
>
> > > DEBUG 3: Deriving wind speed from U and V wind components.
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
>
> > > records matching VarInfo "WIND/Z10" in GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_mean.19591115".
>
> > > DEBUG 3: Found 0 climatology fields.
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Attempt to
derive
>
> > > winds from U and V components.
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading U-wind
> records.
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact
match
>
> > > for VarInfo "WIND/Z10" in GRIB record 69 of GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact
match
>
> > > for VarInfo "WIND/Z10" in GRIB record 70 of GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact
match
>
> > > for VarInfo "WIND/Z10" in GRIB record 71 of GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact
match
>
> > > for VarInfo "WIND/Z10" in GRIB record 72 of GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
>
> > > records matching VarInfo "WIND/Z10" in GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Reading V-wind
> records.
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact
match
>
> > > for VarInfo "WIND/Z10" in GRIB record 73 of GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact
match
>
> > > for VarInfo "WIND/Z10" in GRIB record 74 of GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact
match
>
> > > for VarInfo "WIND/Z10" in GRIB record 75 of GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found exact
match
>
> > > for VarInfo "WIND/Z10" in GRIB record 76 of GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
>
> > > records matching VarInfo "WIND/Z10" in GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
>
> > > DEBUG 3: Deriving wind speed from U and V wind components.
>
> > > DEBUG 3: Deriving wind speed from U and V wind components.
>
> > > DEBUG 3: Deriving wind speed from U and V wind components.
>
> > > DEBUG 3: Deriving wind speed from U and V wind components.
>
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 4 GRIB
>
> > > records matching VarInfo "WIND/Z10" in GRIB file
>
> > >
"/h/data/global/WXQC/data/climo/ncep/NCEP_2.5deg/pgba_stdv.19591115".
>
> > > DEBUG 3: Found 0 climatology fields.
>
> > > DEBUG 3: Found 0 climatology mean and 0 climatology standard
>
> > > deviation
>
> > > field(s) for forecast PROB(WIND>15.000)/Z10/PROB.
>
> > > DEBUG 2: Processing masking regions.
>
> > > DEBUG 3: Processing grid mask: FULL
>
> > >
>
> > > So, MET is reading the climo file just fine but does not find a
>
> > > match
>
> > with
>
> > > the FCST Probability (WIND>15).   Do I have to convert the NCEP
climo
>
> > file
>
> > > in a probabilistic format where winds > 15 is 1 and winds < 15
is 0?
>
> > > I need the climo field to get a Brier Skill Score.
>
> > >
>
> > > Thanks
>
> > > Bob Craig
>
> > >
>
> > >
>
> > >
>
> >
>
> >
>
> >
>
> >
>
>
>
>
>

------------------------------------------------
Subject: Grid Stat Climo File
From: Minna Win
Time: Thu Dec 20 11:47:10 2018

This ticket has been closed.  The user had forgot to set the climo
match_data field to false.

------------------------------------------------
Subject: RE: [Non-DoD Source] [rt.rap.ucar.edu #88200] Grid Stat Climo File
From: robert.craig.2 at us.af.mil
Time: Thu Dec 20 12:44:23 2018

Minna, do you assign a new ticket to the new problem I sent?

thanks

-----Original Message-----
From: Minna Win via RT <met_help at ucar.edu>
Sent: Thursday, December 20, 2018 12:47 PM
To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
Subject: [Non-DoD Source] [rt.rap.ucar.edu #88200] Grid Stat Climo
File

This ticket has been closed.  The user had forgot to set the climo
match_data field to false.


------------------------------------------------
Subject: Grid Stat Climo File
From: Minna Win
Time: Thu Dec 20 12:54:43 2018

Hi Bob,

A new ticket number was automatically generated when you submitted
your
latest query. Do you need this re-opened?

Regards,
Minna
---------------
Minna Win
NCAR
Research Applications Lab
Phone: 303-497-8423
Fax:   303-497-8401



On Thu, Dec 20, 2018 at 7:44 PM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88200 >
>
> Minna, do you assign a new ticket to the new problem I sent?
>
> thanks
>
> -----Original Message-----
> From: Minna Win via RT <met_help at ucar.edu>
> Sent: Thursday, December 20, 2018 12:47 PM
> To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
> Subject: [Non-DoD Source] [rt.rap.ucar.edu #88200] Grid Stat Climo
File
>
> This ticket has been closed.  The user had forgot to set the climo
> match_data field to false.
>
>
>

------------------------------------------------
Subject: RE: [Non-DoD Source] [rt.rap.ucar.edu #88200] Grid Stat Climo File
From: robert.craig.2 at us.af.mil
Time: Thu Dec 20 13:21:17 2018

Minna

No. I didn't get a kick back from met_help on the second problem, that
is why I asked.

Bob

-----Original Message-----
From: Minna Win via RT <met_help at ucar.edu>
Sent: Thursday, December 20, 2018 1:55 PM
To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
Subject: Re: [Non-DoD Source] [rt.rap.ucar.edu #88200] Grid Stat Climo
File

Hi Bob,

A new ticket number was automatically generated when you submitted
your latest query. Do you need this re-opened?

Regards,
Minna
---------------
Minna Win
NCAR
Research Applications Lab
Phone: 303-497-8423
Fax:   303-497-8401



On Thu, Dec 20, 2018 at 7:44 PM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88200 >
>
> Minna, do you assign a new ticket to the new problem I sent?
>
> thanks
>
> -----Original Message-----
> From: Minna Win via RT <met_help at ucar.edu>
> Sent: Thursday, December 20, 2018 12:47 PM
> To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
> <robert.craig.2 at us.af.mil>
> Subject: [Non-DoD Source] [rt.rap.ucar.edu #88200] Grid Stat Climo
> File
>
> This ticket has been closed.  The user had forgot to set the climo
> match_data field to false.
>
>
>



------------------------------------------------
Subject: Grid Stat Climo File
From: Minna Win
Time: Thu Dec 20 13:26:33 2018

Hi Bob,

I forwarded the message that changed the ticket to 'closed'.
Sometimes the
ticketing system doesn't send out all the emails that we expect
(especially
from those of us who aren't as familiar with the ticketing system).

Thanks for your patience.

Regards,
Minna
---------------
Minna Win
NCAR
Research Applications Lab
Phone: 303-497-8423
Fax:   303-497-8401



On Thu, Dec 20, 2018 at 8:21 PM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88200 >
>
> Minna
>
> No. I didn't get a kick back from met_help on the second problem,
that is
> why I asked.
>
> Bob
>
> -----Original Message-----
> From: Minna Win via RT <met_help at ucar.edu>
> Sent: Thursday, December 20, 2018 1:55 PM
> To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
> Subject: Re: [Non-DoD Source] [rt.rap.ucar.edu #88200] Grid Stat
Climo
> File
>
> Hi Bob,
>
> A new ticket number was automatically generated when you submitted
your
> latest query. Do you need this re-opened?
>
> Regards,
> Minna
> ---------------
> Minna Win
> NCAR
> Research Applications Lab
> Phone: 303-497-8423
> Fax:   303-497-8401
>
>
>
> On Thu, Dec 20, 2018 at 7:44 PM robert.craig.2 at us.af.mil via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=88200 >
> >
> > Minna, do you assign a new ticket to the new problem I sent?
> >
> > thanks
> >
> > -----Original Message-----
> > From: Minna Win via RT <met_help at ucar.edu>
> > Sent: Thursday, December 20, 2018 12:47 PM
> > To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
> > <robert.craig.2 at us.af.mil>
> > Subject: [Non-DoD Source] [rt.rap.ucar.edu #88200] Grid Stat Climo
> > File
> >
> > This ticket has been closed.  The user had forgot to set the climo
> > match_data field to false.
> >
> >
> >
>
>
>
>

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

More information about the Met_help mailing list