[Met_help] [rt.rap.ucar.edu #88200] History for Grid Stat Climo File
Minna Win via RT
met_help at ucar.edu
Thu Dec 20 11:47:11 MST 2018
----------------------------------------------------------------
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.
------------------------------------------------
More information about the Met_help
mailing list