[Met_help] [rt.rap.ucar.edu #81829] History for grib2 files as input into MET grid_stat
John Halley Gotway via RT
met_help at ucar.edu
Fri Sep 1 13:16:47 MDT 2017
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
I tried to use grid_stat with two grib2 files as input, and received an
error regarding the GRIB2 field abbreviation 'APCP06' (see output at bottom
of message). Yet when I run degrib on the fcst and obs files that's what I
get:
(g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
p06m_2017082312f006.grb2 -I
MsgNum, Byte, GRIB-Version, elem, level, reference(UTC), valid(UTC),
Proj(hr)
1.0, 0, 2, APCP06="06 hr Total precipitation [kg/(m^2)]", 0-SFC, 08/23/2017
12:00, 08/23/2017 18:00, 6.00
(g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
urma2p5.2017082318.pcp_06h.184.grb2 -I
MsgNum, Byte, GRIB-Version, elem, level, reference(UTC), valid(UTC),
Proj(hr)
1.0, 0, 2, APCP06="06 hr Total precipitation [kg/(m^2)]", 0-SFC, 08/23/2017
12:00, 08/23/2017 18:00, 6.00
(g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
If I put both the fcst and obs files through pcp_combine and use the
resultant .nc files, grid_stat runs fine, with the caveat that the field
name becomes 'APCP_06'.
Thanks for any comments/suggestions, Brad Diehl
(g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ grid_stat
p06m_2017082312f006.grb2 urma2p5.2017082318.pcp_06h.184.grb2
config/GridStatConfig_APCP_06_grib -outdir out/grid_stat -v 4DEBUG 1:
Default Config File:
/usrx/local/dev/MET/6.0/share/met/config/GridStatConfig_default
DEBUG 1: User Config File: config/GridStatConfig_APCP_06_grib
DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
Met2dDataFile object of type "FileType_Gb2".
WARNING:
WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not currently
support Lambert Conformal grids where dx (2.5397) != dy (2.5397) and may
produce unexpected results!
WARNING:
DEBUG 4:
DEBUG 4: Lambert Conformal Grid Data:
DEBUG 4: scale_lat_1: 25
DEBUG 4: scale_lat_2: 25
DEBUG 4: lat_pin: 20.192
DEBUG 4: lon_pin: 121.554
DEBUG 4: x_pin: 0
DEBUG 4: y_pin: 0
DEBUG 4: lon_orient: 95
DEBUG 4: d_km: 2.5397
DEBUG 4: r_km: 6371.2
DEBUG 4: nx: 2145
DEBUG 4: ny: 1597
DEBUG 4:
DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
Met2dDataFile object of type "FileType_Gb2".
DEBUG 4:
DEBUG 4: Lambert Conformal Grid Data:
DEBUG 4: scale_lat_1: 25
DEBUG 4: scale_lat_2: 25
DEBUG 4: lat_pin: 20.192
DEBUG 4: lon_pin: 121.554
DEBUG 4: x_pin: 0
DEBUG 4: y_pin: 0
DEBUG 4: lon_orient: 95
DEBUG 4: d_km: 2.5397
DEBUG 4: r_km: 6371.2
DEBUG 4: nx: 2145
DEBUG 4: ny: 1377
DEBUG 4:
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object of
type "FileType_Gb2".
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object of
type "FileType_Gb2".
ERROR :
ERROR : VarInfoGrib2::set_dict() -> unrecognized GRIB2 field abbreviation
'APCP06'
ERROR :
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: grib2 files as input into MET grid_stat
From: Minna Win
Time: Thu Aug 31 08:51:22 2017
HI Brad,
wgrib is an NDFD tool and I don't know how it is interpreting the APCP
variable in your grib file. Try using wgrib2 on your grib file and
grep for
APCP and use that variable abbreviation in you MET Grid_stat config
file. I
hope that helps. Please let us know if you still have problems using
Grid_stat.
Regards,
Minna
---------------
Minna Win
NCAR
Research Applications Lab
Phone: 303-497-8423
Fax: 302-497-8401
On Thu, Aug 31, 2017 at 2:33 PM, Brad Diehl - NOAA Federal via RT <
met_help at ucar.edu> wrote:
>
> Thu Aug 31 08:33:00 2017: Request 81829 was acted upon.
> Transaction: Ticket created by brad.diehl at noaa.gov
> Queue: met_help
> Subject: grib2 files as input into MET grid_stat
> Owner: Nobody
> Requestors: brad.diehl at noaa.gov
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
>
>
> I tried to use grid_stat with two grib2 files as input, and received
an
> error regarding the GRIB2 field abbreviation 'APCP06' (see output at
bottom
> of message). Yet when I run degrib on the fcst and obs files that's
what I
> get:
>
> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
> p06m_2017082312f006.grb2 -I
> MsgNum, Byte, GRIB-Version, elem, level, reference(UTC), valid(UTC),
> Proj(hr)
> 1.0, 0, 2, APCP06="06 hr Total precipitation [kg/(m^2)]", 0-SFC,
08/23/2017
> 12:00, 08/23/2017 18:00, 6.00
>
> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
> urma2p5.2017082318.pcp_06h.184.grb2 -I
> MsgNum, Byte, GRIB-Version, elem, level, reference(UTC), valid(UTC),
> Proj(hr)
> 1.0, 0, 2, APCP06="06 hr Total precipitation [kg/(m^2)]", 0-SFC,
08/23/2017
> 12:00, 08/23/2017 18:00, 6.00
> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
>
> If I put both the fcst and obs files through pcp_combine and use the
> resultant .nc files, grid_stat runs fine, with the caveat that the
field
> name becomes 'APCP_06'.
>
> Thanks for any comments/suggestions, Brad Diehl
>
> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ grid_stat
> p06m_2017082312f006.grb2 urma2p5.2017082318.pcp_06h.184.grb2
> config/GridStatConfig_APCP_06_grib -outdir out/grid_stat -v 4DEBUG
1:
> Default Config File:
> /usrx/local/dev/MET/6.0/share/met/config/GridStatConfig_default
> DEBUG 1: User Config File: config/GridStatConfig_APCP_06_grib
> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
> Met2dDataFile object of type "FileType_Gb2".
> WARNING:
> WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
currently
> support Lambert Conformal grids where dx (2.5397) != dy (2.5397) and
may
> produce unexpected results!
> WARNING:
> DEBUG 4:
> DEBUG 4: Lambert Conformal Grid Data:
> DEBUG 4: scale_lat_1: 25
> DEBUG 4: scale_lat_2: 25
> DEBUG 4: lat_pin: 20.192
> DEBUG 4: lon_pin: 121.554
> DEBUG 4: x_pin: 0
> DEBUG 4: y_pin: 0
> DEBUG 4: lon_orient: 95
> DEBUG 4: d_km: 2.5397
> DEBUG 4: r_km: 6371.2
> DEBUG 4: nx: 2145
> DEBUG 4: ny: 1597
> DEBUG 4:
> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
> Met2dDataFile object of type "FileType_Gb2".
> DEBUG 4:
> DEBUG 4: Lambert Conformal Grid Data:
> DEBUG 4: scale_lat_1: 25
> DEBUG 4: scale_lat_2: 25
> DEBUG 4: lat_pin: 20.192
> DEBUG 4: lon_pin: 121.554
> DEBUG 4: x_pin: 0
> DEBUG 4: y_pin: 0
> DEBUG 4: lon_orient: 95
> DEBUG 4: d_km: 2.5397
> DEBUG 4: r_km: 6371.2
> DEBUG 4: nx: 2145
> DEBUG 4: ny: 1377
> DEBUG 4:
> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> type "FileType_Gb2".
> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> type "FileType_Gb2".
> ERROR :
> ERROR : VarInfoGrib2::set_dict() -> unrecognized GRIB2 field
abbreviation
> 'APCP06'
> ERROR :
>
>
------------------------------------------------
Subject: grib2 files as input into MET grid_stat
From: Minna Win
Time: Thu Aug 31 08:58:32 2017
Hi Brad,
Please disregard my previous response. I didn't see that you had run
pcp_combine-that *should* have worked. I'll look a little further or
see
if someone else has seen this before.
Regards,
Minna
---------------
Minna Win
NCAR
Research Applications Lab
Phone: 303-497-8423
Fax: 302-497-8401
On Thu, Aug 31, 2017 at 2:50 PM, Minna Win-Gildenmeister
<minnawin at ucar.edu>
wrote:
> HI Brad,
>
> wgrib is an NDFD tool and I don't know how it is interpreting the
APCP
> variable in your grib file. Try using wgrib2 on your grib file and
grep for
> APCP and use that variable abbreviation in you MET Grid_stat config
file. I
> hope that helps. Please let us know if you still have problems
using
> Grid_stat.
>
> Regards,
> Minna
>
> ---------------
> Minna Win
> NCAR
> Research Applications Lab
> Phone: 303-497-8423 <(303)%20497-8423>
> Fax: 302-497-8401 <(302)%20497-8401>
>
>
> On Thu, Aug 31, 2017 at 2:33 PM, Brad Diehl - NOAA Federal via RT <
> met_help at ucar.edu> wrote:
>
>>
>> Thu Aug 31 08:33:00 2017: Request 81829 was acted upon.
>> Transaction: Ticket created by brad.diehl at noaa.gov
>> Queue: met_help
>> Subject: grib2 files as input into MET grid_stat
>> Owner: Nobody
>> Requestors: brad.diehl at noaa.gov
>> Status: new
>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
>>
>>
>> I tried to use grid_stat with two grib2 files as input, and
received an
>> error regarding the GRIB2 field abbreviation 'APCP06' (see output
at
>> bottom
>> of message). Yet when I run degrib on the fcst and obs files that's
what I
>> get:
>>
>> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
>> p06m_2017082312f006.grb2 -I
>> MsgNum, Byte, GRIB-Version, elem, level, reference(UTC),
valid(UTC),
>> Proj(hr)
>> 1.0, 0, 2, APCP06="06 hr Total precipitation [kg/(m^2)]", 0-SFC,
>> 08/23/2017
>> 12:00, 08/23/2017 18:00, 6.00
>>
>> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
>> urma2p5.2017082318 <(201)%20708-2318>.pcp_06h.184.grb2 -I
>> MsgNum, Byte, GRIB-Version, elem, level, reference(UTC),
valid(UTC),
>> Proj(hr)
>> 1.0, 0, 2, APCP06="06 hr Total precipitation [kg/(m^2)]", 0-SFC,
>> 08/23/2017
>> 12:00, 08/23/2017 18:00, 6.00
>> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
>>
>> If I put both the fcst and obs files through pcp_combine and use
the
>> resultant .nc files, grid_stat runs fine, with the caveat that the
field
>> name becomes 'APCP_06'.
>>
>> Thanks for any comments/suggestions, Brad Diehl
>>
>> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ grid_stat
>> p06m_2017082312f006.grb2 urma2p5.2017082318 <(201)%20708-2318>
>> .pcp_06h.184.grb2
>> config/GridStatConfig_APCP_06_grib -outdir out/grid_stat -v 4DEBUG
1:
>> Default Config File:
>> /usrx/local/dev/MET/6.0/share/met/config/GridStatConfig_default
>> DEBUG 1: User Config File: config/GridStatConfig_APCP_06_grib
>> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
>> Met2dDataFile object of type "FileType_Gb2".
>> WARNING:
>> WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
currently
>> support Lambert Conformal grids where dx (2.5397) != dy (2.5397)
and may
>> produce unexpected results!
>> WARNING:
>> DEBUG 4:
>> DEBUG 4: Lambert Conformal Grid Data:
>> DEBUG 4: scale_lat_1: 25
>> DEBUG 4: scale_lat_2: 25
>> DEBUG 4: lat_pin: 20.192
>> DEBUG 4: lon_pin: 121.554
>> DEBUG 4: x_pin: 0
>> DEBUG 4: y_pin: 0
>> DEBUG 4: lon_orient: 95
>> DEBUG 4: d_km: 2.5397
>> DEBUG 4: r_km: 6371.2
>> DEBUG 4: nx: 2145
>> DEBUG 4: ny: 1597
>> DEBUG 4:
>> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
>> Met2dDataFile object of type "FileType_Gb2".
>> DEBUG 4:
>> DEBUG 4: Lambert Conformal Grid Data:
>> DEBUG 4: scale_lat_1: 25
>> DEBUG 4: scale_lat_2: 25
>> DEBUG 4: lat_pin: 20.192
>> DEBUG 4: lon_pin: 121.554
>> DEBUG 4: x_pin: 0
>> DEBUG 4: y_pin: 0
>> DEBUG 4: lon_orient: 95
>> DEBUG 4: d_km: 2.5397
>> DEBUG 4: r_km: 6371.2
>> DEBUG 4: nx: 2145
>> DEBUG 4: ny: 1377
>> DEBUG 4:
>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
>> type "FileType_Gb2".
>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
>> type "FileType_Gb2".
>> ERROR :
>> ERROR : VarInfoGrib2::set_dict() -> unrecognized GRIB2 field
abbreviation
>> 'APCP06'
>> ERROR :
>>
>>
>
------------------------------------------------
Subject: grib2 files as input into MET grid_stat
From: Minna Win
Time: Thu Aug 31 09:05:26 2017
Hi Brad,
My apologies, I thought your pcp_combine/grid_stat combination failed
but
it succeeded. Please try using wgrib2 on your grib2 file and grep for
APCP
to see what wgrib2 returns for the variable abbreviation and use that
in
grid_stat, as originally posted.
Thanks,
Minna
---------------
Minna Win
NCAR
Research Applications Lab
Phone: 303-497-8423
Fax: 302-497-8401
On Thu, Aug 31, 2017 at 2:57 PM, Minna Win-Gildenmeister
<minnawin at ucar.edu>
wrote:
> Hi Brad,
>
> Please disregard my previous response. I didn't see that you had
run
> pcp_combine-that *should* have worked. I'll look a little further
or see
> if someone else has seen this before.
>
>
>
> Regards,
> Minna
>
> ---------------
> Minna Win
> NCAR
> Research Applications Lab
> Phone: 303-497-8423 <(303)%20497-8423>
> Fax: 302-497-8401 <(302)%20497-8401>
>
>
> On Thu, Aug 31, 2017 at 2:50 PM, Minna Win-Gildenmeister <
> minnawin at ucar.edu> wrote:
>
>> HI Brad,
>>
>> wgrib is an NDFD tool and I don't know how it is interpreting the
APCP
>> variable in your grib file. Try using wgrib2 on your grib file and
grep for
>> APCP and use that variable abbreviation in you MET Grid_stat config
file. I
>> hope that helps. Please let us know if you still have problems
using
>> Grid_stat.
>>
>> Regards,
>> Minna
>>
>> ---------------
>> Minna Win
>> NCAR
>> Research Applications Lab
>> Phone: 303-497-8423 <(303)%20497-8423>
>> Fax: 302-497-8401 <(302)%20497-8401>
>>
>>
>> On Thu, Aug 31, 2017 at 2:33 PM, Brad Diehl - NOAA Federal via RT <
>> met_help at ucar.edu> wrote:
>>
>>>
>>> Thu Aug 31 08:33:00 2017: Request 81829 was acted upon.
>>> Transaction: Ticket created by brad.diehl at noaa.gov
>>> Queue: met_help
>>> Subject: grib2 files as input into MET grid_stat
>>> Owner: Nobody
>>> Requestors: brad.diehl at noaa.gov
>>> Status: new
>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
>>>
>>>
>>> I tried to use grid_stat with two grib2 files as input, and
received an
>>> error regarding the GRIB2 field abbreviation 'APCP06' (see output
at
>>> bottom
>>> of message). Yet when I run degrib on the fcst and obs files
that's what
>>> I
>>> get:
>>>
>>> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
>>> p06m_2017082312f006.grb2 -I
>>> MsgNum, Byte, GRIB-Version, elem, level, reference(UTC),
valid(UTC),
>>> Proj(hr)
>>> 1.0, 0, 2, APCP06="06 hr Total precipitation [kg/(m^2)]", 0-SFC,
>>> 08/23/2017
>>> 12:00, 08/23/2017 18:00, 6.00
>>>
>>> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
>>> urma2p5.2017082318 <(201)%20708-2318>.pcp_06h.184.grb2 -I
>>> MsgNum, Byte, GRIB-Version, elem, level, reference(UTC),
valid(UTC),
>>> Proj(hr)
>>> 1.0, 0, 2, APCP06="06 hr Total precipitation [kg/(m^2)]", 0-SFC,
>>> 08/23/2017
>>> 12:00, 08/23/2017 18:00, 6.00
>>> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
>>>
>>> If I put both the fcst and obs files through pcp_combine and use
the
>>> resultant .nc files, grid_stat runs fine, with the caveat that the
field
>>> name becomes 'APCP_06'.
>>>
>>> Thanks for any comments/suggestions, Brad Diehl
>>>
>>> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ grid_stat
>>> p06m_2017082312f006.grb2 urma2p5.2017082318 <(201)%20708-2318>
>>> .pcp_06h.184.grb2
>>> config/GridStatConfig_APCP_06_grib -outdir out/grid_stat -v 4DEBUG
1:
>>> Default Config File:
>>> /usrx/local/dev/MET/6.0/share/met/config/GridStatConfig_default
>>> DEBUG 1: User Config File: config/GridStatConfig_APCP_06_grib
>>> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
>>> Met2dDataFile object of type "FileType_Gb2".
>>> WARNING:
>>> WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
currently
>>> support Lambert Conformal grids where dx (2.5397) != dy (2.5397)
and may
>>> produce unexpected results!
>>> WARNING:
>>> DEBUG 4:
>>> DEBUG 4: Lambert Conformal Grid Data:
>>> DEBUG 4: scale_lat_1: 25
>>> DEBUG 4: scale_lat_2: 25
>>> DEBUG 4: lat_pin: 20.192
>>> DEBUG 4: lon_pin: 121.554
>>> DEBUG 4: x_pin: 0
>>> DEBUG 4: y_pin: 0
>>> DEBUG 4: lon_orient: 95
>>> DEBUG 4: d_km: 2.5397
>>> DEBUG 4: r_km: 6371.2
>>> DEBUG 4: nx: 2145
>>> DEBUG 4: ny: 1597
>>> DEBUG 4:
>>> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
>>> Met2dDataFile object of type "FileType_Gb2".
>>> DEBUG 4:
>>> DEBUG 4: Lambert Conformal Grid Data:
>>> DEBUG 4: scale_lat_1: 25
>>> DEBUG 4: scale_lat_2: 25
>>> DEBUG 4: lat_pin: 20.192
>>> DEBUG 4: lon_pin: 121.554
>>> DEBUG 4: x_pin: 0
>>> DEBUG 4: y_pin: 0
>>> DEBUG 4: lon_orient: 95
>>> DEBUG 4: d_km: 2.5397
>>> DEBUG 4: r_km: 6371.2
>>> DEBUG 4: nx: 2145
>>> DEBUG 4: ny: 1377
>>> DEBUG 4:
>>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
>>> type "FileType_Gb2".
>>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
>>> type "FileType_Gb2".
>>> ERROR :
>>> ERROR : VarInfoGrib2::set_dict() -> unrecognized GRIB2 field
>>> abbreviation
>>> 'APCP06'
>>> ERROR :
>>>
>>>
>>
>
------------------------------------------------
Subject: grib2 files as input into MET grid_stat
From: Brad Diehl - NOAA Federal
Time: Thu Aug 31 09:09:18 2017
After running pcp_combine, everything did work correctly, it is the
grib2
fcts/obs part I'm having trouble with. I think you put me on the right
track by suggesting wgrib2. Here's what I get using wgrib2 on the fcst
file:
(g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ wgrib2
p06m_2017082312f006.grb2
1:21:d=2017082312:APCP:surface - surface:CodeTable 4.11=missing:
in the config file I'm using:
fcst = {
field = [
{
name = "APCP";
level = "Z0" ;
cat_thresh = [ >0.0, >=5.0 , >=10.0];
}
];
for level, should I be using surface or Z0?, neither of which seemed
to
work.
On Thu, Aug 31, 2017 at 2:58 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
> Hi Brad,
>
> Please disregard my previous response. I didn't see that you had
run
> pcp_combine-that *should* have worked. I'll look a little further
or see
> if someone else has seen this before.
>
>
>
> Regards,
> Minna
>
> ---------------
> Minna Win
> NCAR
> Research Applications Lab
> Phone: 303-497-8423
> Fax: 302-497-8401
>
>
> On Thu, Aug 31, 2017 at 2:50 PM, Minna Win-Gildenmeister <
> minnawin at ucar.edu>
> wrote:
>
> > HI Brad,
> >
> > wgrib is an NDFD tool and I don't know how it is interpreting the
APCP
> > variable in your grib file. Try using wgrib2 on your grib file and
grep
> for
> > APCP and use that variable abbreviation in you MET Grid_stat
config
> file. I
> > hope that helps. Please let us know if you still have problems
using
> > Grid_stat.
> >
> > Regards,
> > Minna
> >
> > ---------------
> > Minna Win
> > NCAR
> > Research Applications Lab
> > Phone: 303-497-8423 <(303)%20497-8423>
> > Fax: 302-497-8401 <(302)%20497-8401>
> >
> >
> > On Thu, Aug 31, 2017 at 2:33 PM, Brad Diehl - NOAA Federal via RT
<
> > met_help at ucar.edu> wrote:
> >
> >>
> >> Thu Aug 31 08:33:00 2017: Request 81829 was acted upon.
> >> Transaction: Ticket created by brad.diehl at noaa.gov
> >> Queue: met_help
> >> Subject: grib2 files as input into MET grid_stat
> >> Owner: Nobody
> >> Requestors: brad.diehl at noaa.gov
> >> Status: new
> >> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> >>
> >>
> >> I tried to use grid_stat with two grib2 files as input, and
received an
> >> error regarding the GRIB2 field abbreviation 'APCP06' (see output
at
> >> bottom
> >> of message). Yet when I run degrib on the fcst and obs files
that's
> what I
> >> get:
> >>
> >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
> >> p06m_2017082312f006.grb2 -I
> >> MsgNum, Byte, GRIB-Version, elem, level, reference(UTC),
valid(UTC),
> >> Proj(hr)
> >> 1.0, 0, 2, APCP06="06 hr Total precipitation [kg/(m^2)]", 0-SFC,
> >> 08/23/2017
> >> 12:00, 08/23/2017 18:00, 6.00
> >>
> >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
> >> urma2p5.2017082318 <(201)%20708-2318>.pcp_06h.184.grb2 -I
> >> MsgNum, Byte, GRIB-Version, elem, level, reference(UTC),
valid(UTC),
> >> Proj(hr)
> >> 1.0, 0, 2, APCP06="06 hr Total precipitation [kg/(m^2)]", 0-SFC,
> >> 08/23/2017
> >> 12:00, 08/23/2017 18:00, 6.00
> >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> >>
> >> If I put both the fcst and obs files through pcp_combine and use
the
> >> resultant .nc files, grid_stat runs fine, with the caveat that
the field
> >> name becomes 'APCP_06'.
> >>
> >> Thanks for any comments/suggestions, Brad Diehl
> >>
> >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ grid_stat
> >> p06m_2017082312f006.grb2 urma2p5.2017082318 <(201)%20708-2318>
> >> .pcp_06h.184.grb2
> >> config/GridStatConfig_APCP_06_grib -outdir out/grid_stat -v
4DEBUG 1:
> >> Default Config File:
> >> /usrx/local/dev/MET/6.0/share/met/config/GridStatConfig_default
> >> DEBUG 1: User Config File: config/GridStatConfig_APCP_06_grib
> >> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
> >> Met2dDataFile object of type "FileType_Gb2".
> >> WARNING:
> >> WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
currently
> >> support Lambert Conformal grids where dx (2.5397) != dy (2.5397)
and may
> >> produce unexpected results!
> >> WARNING:
> >> DEBUG 4:
> >> DEBUG 4: Lambert Conformal Grid Data:
> >> DEBUG 4: scale_lat_1: 25
> >> DEBUG 4: scale_lat_2: 25
> >> DEBUG 4: lat_pin: 20.192
> >> DEBUG 4: lon_pin: 121.554
> >> DEBUG 4: x_pin: 0
> >> DEBUG 4: y_pin: 0
> >> DEBUG 4: lon_orient: 95
> >> DEBUG 4: d_km: 2.5397
> >> DEBUG 4: r_km: 6371.2
> >> DEBUG 4: nx: 2145
> >> DEBUG 4: ny: 1597
> >> DEBUG 4:
> >> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
> >> Met2dDataFile object of type "FileType_Gb2".
> >> DEBUG 4:
> >> DEBUG 4: Lambert Conformal Grid Data:
> >> DEBUG 4: scale_lat_1: 25
> >> DEBUG 4: scale_lat_2: 25
> >> DEBUG 4: lat_pin: 20.192
> >> DEBUG 4: lon_pin: 121.554
> >> DEBUG 4: x_pin: 0
> >> DEBUG 4: y_pin: 0
> >> DEBUG 4: lon_orient: 95
> >> DEBUG 4: d_km: 2.5397
> >> DEBUG 4: r_km: 6371.2
> >> DEBUG 4: nx: 2145
> >> DEBUG 4: ny: 1377
> >> DEBUG 4:
> >> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> >> type "FileType_Gb2".
> >> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> >> type "FileType_Gb2".
> >> ERROR :
> >> ERROR : VarInfoGrib2::set_dict() -> unrecognized GRIB2 field
> abbreviation
> >> 'APCP06'
> >> ERROR :
> >>
> >>
> >
>
>
------------------------------------------------
Subject: grib2 files as input into MET grid_stat
From: Minna Win
Time: Thu Aug 31 10:01:16 2017
Hi Brad,
Can you try "L0", and "P0"? If that doesn't work, could you please
send me
your grib file? Instructions are as follows:
Non-Windows FTP client:
* ftp -p ftp.rap.ucar.edu <http://ftp.rap.ucar.edu> username=anonymous
password="your email address" cd incoming/irap/met_help mkdir
"surname_data" # where surname is your last name cd "surname_data"
put
"your files" bye *
Windows FTP client:
* ftp ftp.rap.ucar.edu <http://ftp.rap.ucar.edu> username=anonymous
password="your email address" quote pasv cd incoming/irap/met_help
mkdir
"surname_data" # where surname is your last name cd "surname_data" put
"your files" bye*
I can take a look at the level type and value, you might need to use
special syntax in your config file to get grid_stat to work. Let me
know
if you ftp'd a grib file so I can grab it before the file gets
scrubbed.
Thanks,
Minna
---------------
Minna Win
NCAR
Research Applications Lab
Phone: 303-497-8423
Fax: 302-497-8401
On Thu, Aug 31, 2017 at 3:09 PM, Brad Diehl - NOAA Federal via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
>
> After running pcp_combine, everything did work correctly, it is the
grib2
> fcts/obs part I'm having trouble with. I think you put me on the
right
> track by suggesting wgrib2. Here's what I get using wgrib2 on the
fcst
> file:
>
> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ wgrib2
> p06m_2017082312f006.grb2
> 1:21:d=2017082312:APCP:surface - surface:CodeTable 4.11=missing:
>
> in the config file I'm using:
>
>
> fcst = {
>
> field = [
> {
> name = "APCP";
> level = "Z0" ;
> cat_thresh = [ >0.0, >=5.0 , >=10.0];
> }
> ];
>
> for level, should I be using surface or Z0?, neither of which seemed
to
> work.
>
> On Thu, Aug 31, 2017 at 2:58 PM, Minna Win via RT
<met_help at ucar.edu>
> wrote:
>
> > Hi Brad,
> >
> > Please disregard my previous response. I didn't see that you had
run
> > pcp_combine-that *should* have worked. I'll look a little further
or see
> > if someone else has seen this before.
> >
> >
> >
> > Regards,
> > Minna
> >
> > ---------------
> > Minna Win
> > NCAR
> > Research Applications Lab
> > Phone: 303-497-8423
> > Fax: 302-497-8401
> >
> >
> > On Thu, Aug 31, 2017 at 2:50 PM, Minna Win-Gildenmeister <
> > minnawin at ucar.edu>
> > wrote:
> >
> > > HI Brad,
> > >
> > > wgrib is an NDFD tool and I don't know how it is interpreting
the APCP
> > > variable in your grib file. Try using wgrib2 on your grib file
and grep
> > for
> > > APCP and use that variable abbreviation in you MET Grid_stat
config
> > file. I
> > > hope that helps. Please let us know if you still have problems
using
> > > Grid_stat.
> > >
> > > Regards,
> > > Minna
> > >
> > > ---------------
> > > Minna Win
> > > NCAR
> > > Research Applications Lab
> > > Phone: 303-497-8423 <(303)%20497-8423>
> > > Fax: 302-497-8401 <(302)%20497-8401>
> > >
> > >
> > > On Thu, Aug 31, 2017 at 2:33 PM, Brad Diehl - NOAA Federal via
RT <
> > > met_help at ucar.edu> wrote:
> > >
> > >>
> > >> Thu Aug 31 08:33:00 2017: Request 81829 was acted upon.
> > >> Transaction: Ticket created by brad.diehl at noaa.gov
> > >> Queue: met_help
> > >> Subject: grib2 files as input into MET grid_stat
> > >> Owner: Nobody
> > >> Requestors: brad.diehl at noaa.gov
> > >> Status: new
> > >> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829
> >
> > >>
> > >>
> > >> I tried to use grid_stat with two grib2 files as input, and
received
> an
> > >> error regarding the GRIB2 field abbreviation 'APCP06' (see
output at
> > >> bottom
> > >> of message). Yet when I run degrib on the fcst and obs files
that's
> > what I
> > >> get:
> > >>
> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
> > >> p06m_2017082312f006.grb2 -I
> > >> MsgNum, Byte, GRIB-Version, elem, level, reference(UTC),
valid(UTC),
> > >> Proj(hr)
> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation [kg/(m^2)]", 0-
SFC,
> > >> 08/23/2017
> > >> 12:00, 08/23/2017 18:00, 6.00
> > >>
> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
> > >> urma2p5.2017082318 <(201)%20708-2318>.pcp_06h.184.grb2 -I
> > >> MsgNum, Byte, GRIB-Version, elem, level, reference(UTC),
valid(UTC),
> > >> Proj(hr)
> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation [kg/(m^2)]", 0-
SFC,
> > >> 08/23/2017
> > >> 12:00, 08/23/2017 18:00, 6.00
> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > >>
> > >> If I put both the fcst and obs files through pcp_combine and
use the
> > >> resultant .nc files, grid_stat runs fine, with the caveat that
the
> field
> > >> name becomes 'APCP_06'.
> > >>
> > >> Thanks for any comments/suggestions, Brad Diehl
> > >>
> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ grid_stat
> > >> p06m_2017082312f006.grb2 urma2p5.2017082318 <(201)%20708-2318>
> > >> .pcp_06h.184.grb2
> > >> config/GridStatConfig_APCP_06_grib -outdir out/grid_stat -v
4DEBUG 1:
> > >> Default Config File:
> > >> /usrx/local/dev/MET/6.0/share/met/config/GridStatConfig_default
> > >> DEBUG 1: User Config File: config/GridStatConfig_APCP_06_grib
> > >> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() ->
created new
> > >> Met2dDataFile object of type "FileType_Gb2".
> > >> WARNING:
> > >> WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
> currently
> > >> support Lambert Conformal grids where dx (2.5397) != dy
(2.5397) and
> may
> > >> produce unexpected results!
> > >> WARNING:
> > >> DEBUG 4:
> > >> DEBUG 4: Lambert Conformal Grid Data:
> > >> DEBUG 4: scale_lat_1: 25
> > >> DEBUG 4: scale_lat_2: 25
> > >> DEBUG 4: lat_pin: 20.192
> > >> DEBUG 4: lon_pin: 121.554
> > >> DEBUG 4: x_pin: 0
> > >> DEBUG 4: y_pin: 0
> > >> DEBUG 4: lon_orient: 95
> > >> DEBUG 4: d_km: 2.5397
> > >> DEBUG 4: r_km: 6371.2
> > >> DEBUG 4: nx: 2145
> > >> DEBUG 4: ny: 1597
> > >> DEBUG 4:
> > >> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() ->
created new
> > >> Met2dDataFile object of type "FileType_Gb2".
> > >> DEBUG 4:
> > >> DEBUG 4: Lambert Conformal Grid Data:
> > >> DEBUG 4: scale_lat_1: 25
> > >> DEBUG 4: scale_lat_2: 25
> > >> DEBUG 4: lat_pin: 20.192
> > >> DEBUG 4: lon_pin: 121.554
> > >> DEBUG 4: x_pin: 0
> > >> DEBUG 4: y_pin: 0
> > >> DEBUG 4: lon_orient: 95
> > >> DEBUG 4: d_km: 2.5397
> > >> DEBUG 4: r_km: 6371.2
> > >> DEBUG 4: nx: 2145
> > >> DEBUG 4: ny: 1377
> > >> DEBUG 4:
> > >> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
> of
> > >> type "FileType_Gb2".
> > >> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
> of
> > >> type "FileType_Gb2".
> > >> ERROR :
> > >> ERROR : VarInfoGrib2::set_dict() -> unrecognized GRIB2 field
> > abbreviation
> > >> 'APCP06'
> > >> ERROR :
> > >>
> > >>
> > >
> >
> >
>
>
------------------------------------------------
Subject: grib2 files as input into MET grid_stat
From: Minna Win
Time: Thu Aug 31 10:11:18 2017
Hi Brad,
Tara Jensen brought this to my attention, for verify APCP, you need to
set
the level to the value pre-pended with the letter 'A'. For instance,
for 6
hour accumulation level, set level="A06".
Thanks,
Minna
---------------
Minna Win
NCAR
Research Applications Lab
Phone: 303-497-8423
Fax: 302-497-8401
On Thu, Aug 31, 2017 at 3:52 PM, Minna Win-Gildenmeister
<minnawin at ucar.edu>
wrote:
> Hi Brad,
>
> Can you try "L0", and "P0"? If that doesn't work, could you please
send
> me your grib file? Instructions are as follows:
>
> Non-Windows FTP client:
>
>
>
>
>
>
>
>
>
> * ftp -p ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
username=anonymous
> password="your email address" cd incoming/irap/met_help mkdir
> "surname_data" # where surname is your last name cd
"surname_data" put
> "your files" bye *
>
> Windows FTP client:
>
>
>
>
>
>
>
>
> * ftp ftp.rap.ucar.edu <http://ftp.rap.ucar.edu> username=anonymous
> password="your email address" quote pasv cd incoming/irap/met_help
mkdir
> "surname_data" # where surname is your last name cd "surname_data"
put
> "your files" bye*
>
> I can take a look at the level type and value, you might need to use
> special syntax in your config file to get grid_stat to work. Let me
know
> if you ftp'd a grib file so I can grab it before the file gets
scrubbed.
>
> Thanks,
> Minna
>
>
> ---------------
> Minna Win
> NCAR
> Research Applications Lab
> Phone: 303-497-8423 <(303)%20497-8423>
> Fax: 302-497-8401 <(302)%20497-8401>
>
>
> On Thu, Aug 31, 2017 at 3:09 PM, Brad Diehl - NOAA Federal via RT <
> met_help at ucar.edu> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
>>
>> After running pcp_combine, everything did work correctly, it is the
grib2
>> fcts/obs part I'm having trouble with. I think you put me on the
right
>> track by suggesting wgrib2. Here's what I get using wgrib2 on the
fcst
>> file:
>>
>> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ wgrib2
>> p06m_2017082312f006.grb2
>> 1:21:d=2017082312 <(201)%20708-2312>:APCP:surface -
surface:CodeTable
>> 4.11=missing:
>>
>> in the config file I'm using:
>>
>>
>> fcst = {
>>
>> field = [
>> {
>> name = "APCP";
>> level = "Z0" ;
>> cat_thresh = [ >0.0, >=5.0 , >=10.0];
>> }
>> ];
>>
>> for level, should I be using surface or Z0?, neither of which
seemed to
>> work.
>>
>> On Thu, Aug 31, 2017 at 2:58 PM, Minna Win via RT
<met_help at ucar.edu>
>> wrote:
>>
>> > Hi Brad,
>> >
>> > Please disregard my previous response. I didn't see that you had
run
>> > pcp_combine-that *should* have worked. I'll look a little
further or
>> see
>> > if someone else has seen this before.
>> >
>> >
>> >
>> > Regards,
>> > Minna
>> >
>> > ---------------
>> > Minna Win
>> > NCAR
>> > Research Applications Lab
>> > Phone: 303-497-8423
>> > Fax: 302-497-8401
>> >
>> >
>> > On Thu, Aug 31, 2017 at 2:50 PM, Minna Win-Gildenmeister <
>> > minnawin at ucar.edu>
>> > wrote:
>> >
>> > > HI Brad,
>> > >
>> > > wgrib is an NDFD tool and I don't know how it is interpreting
the APCP
>> > > variable in your grib file. Try using wgrib2 on your grib file
and
>> grep
>> > for
>> > > APCP and use that variable abbreviation in you MET Grid_stat
config
>> > file. I
>> > > hope that helps. Please let us know if you still have problems
using
>> > > Grid_stat.
>> > >
>> > > Regards,
>> > > Minna
>> > >
>> > > ---------------
>> > > Minna Win
>> > > NCAR
>> > > Research Applications Lab
>> > > Phone: 303-497-8423 <(303)%20497-8423>
>> > > Fax: 302-497-8401 <(302)%20497-8401>
>> > >
>> > >
>> > > On Thu, Aug 31, 2017 at 2:33 PM, Brad Diehl - NOAA Federal via
RT <
>> > > met_help at ucar.edu> wrote:
>> > >
>> > >>
>> > >> Thu Aug 31 08:33:00 2017: Request 81829 was acted upon.
>> > >> Transaction: Ticket created by brad.diehl at noaa.gov
>> > >> Queue: met_help
>> > >> Subject: grib2 files as input into MET grid_stat
>> > >> Owner: Nobody
>> > >> Requestors: brad.diehl at noaa.gov
>> > >> Status: new
>> > >> Ticket <URL: https://rt.rap.ucar.edu/rt/Tic
>> ket/Display.html?id=81829 >
>> > >>
>> > >>
>> > >> I tried to use grid_stat with two grib2 files as input, and
received
>> an
>> > >> error regarding the GRIB2 field abbreviation 'APCP06' (see
output at
>> > >> bottom
>> > >> of message). Yet when I run degrib on the fcst and obs files
that's
>> > what I
>> > >> get:
>> > >>
>> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
>> > >> p06m_2017082312f006.grb2 -I
>> > >> MsgNum, Byte, GRIB-Version, elem, level, reference(UTC),
valid(UTC),
>> > >> Proj(hr)
>> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation [kg/(m^2)]", 0-
SFC,
>> > >> 08/23/2017
>> > >> 12:00, 08/23/2017 18:00, 6.00
>> > >>
>> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
>> > >> urma2p5.2017082318 <(201)%20708-2318>.pcp_06h.184.grb2 -I
>> > >> MsgNum, Byte, GRIB-Version, elem, level, reference(UTC),
valid(UTC),
>> > >> Proj(hr)
>> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation [kg/(m^2)]", 0-
SFC,
>> > >> 08/23/2017
>> > >> 12:00, 08/23/2017 18:00, 6.00
>> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
>> > >>
>> > >> If I put both the fcst and obs files through pcp_combine and
use the
>> > >> resultant .nc files, grid_stat runs fine, with the caveat that
the
>> field
>> > >> name becomes 'APCP_06'.
>> > >>
>> > >> Thanks for any comments/suggestions, Brad Diehl
>> > >>
>> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ grid_stat
>> > >> p06m_2017082312f006.grb2 urma2p5.2017082318 <(201)%20708-2318>
>> > >> .pcp_06h.184.grb2
>> > >> config/GridStatConfig_APCP_06_grib -outdir out/grid_stat -v
4DEBUG
>> 1:
>> > >> Default Config File:
>> > >>
/usrx/local/dev/MET/6.0/share/met/config/GridStatConfig_default
>> > >> DEBUG 1: User Config File: config/GridStatConfig_APCP_06_grib
>> > >> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() ->
created new
>> > >> Met2dDataFile object of type "FileType_Gb2".
>> > >> WARNING:
>> > >> WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
>> currently
>> > >> support Lambert Conformal grids where dx (2.5397) != dy
(2.5397) and
>> may
>> > >> produce unexpected results!
>> > >> WARNING:
>> > >> DEBUG 4:
>> > >> DEBUG 4: Lambert Conformal Grid Data:
>> > >> DEBUG 4: scale_lat_1: 25
>> > >> DEBUG 4: scale_lat_2: 25
>> > >> DEBUG 4: lat_pin: 20.192
>> > >> DEBUG 4: lon_pin: 121.554
>> > >> DEBUG 4: x_pin: 0
>> > >> DEBUG 4: y_pin: 0
>> > >> DEBUG 4: lon_orient: 95
>> > >> DEBUG 4: d_km: 2.5397
>> > >> DEBUG 4: r_km: 6371.2
>> > >> DEBUG 4: nx: 2145
>> > >> DEBUG 4: ny: 1597
>> > >> DEBUG 4:
>> > >> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() ->
created new
>> > >> Met2dDataFile object of type "FileType_Gb2".
>> > >> DEBUG 4:
>> > >> DEBUG 4: Lambert Conformal Grid Data:
>> > >> DEBUG 4: scale_lat_1: 25
>> > >> DEBUG 4: scale_lat_2: 25
>> > >> DEBUG 4: lat_pin: 20.192
>> > >> DEBUG 4: lon_pin: 121.554
>> > >> DEBUG 4: x_pin: 0
>> > >> DEBUG 4: y_pin: 0
>> > >> DEBUG 4: lon_orient: 95
>> > >> DEBUG 4: d_km: 2.5397
>> > >> DEBUG 4: r_km: 6371.2
>> > >> DEBUG 4: nx: 2145
>> > >> DEBUG 4: ny: 1377
>> > >> DEBUG 4:
>> > >> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
>> object of
>> > >> type "FileType_Gb2".
>> > >> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
>> object of
>> > >> type "FileType_Gb2".
>> > >> ERROR :
>> > >> ERROR : VarInfoGrib2::set_dict() -> unrecognized GRIB2 field
>> > abbreviation
>> > >> 'APCP06'
>> > >> ERROR :
>> > >>
>> > >>
>> > >
>> >
>> >
>>
>>
>
------------------------------------------------
Subject: grib2 files as input into MET grid_stat
From: Brad Diehl - NOAA Federal
Time: Thu Aug 31 10:21:06 2017
the A06 worked. Prior to trying that I sent you via mailx the fcst/obs
grib2 files, you can ignore that now.
But it does beg the question, how does one find the level names, is it
in
the grib file itself? Very confusing, at least to me. And what does
*,*
mean when I see that in some examples?
Thanks for all your help........Brad
On Thu, Aug 31, 2017 at 4:11 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
> Hi Brad,
>
> Tara Jensen brought this to my attention, for verify APCP, you need
to set
> the level to the value pre-pended with the letter 'A'. For
instance, for 6
> hour accumulation level, set level="A06".
>
> Thanks,
> Minna
>
> ---------------
> Minna Win
> NCAR
> Research Applications Lab
> Phone: 303-497-8423
> Fax: 302-497-8401
>
>
> On Thu, Aug 31, 2017 at 3:52 PM, Minna Win-Gildenmeister <
> minnawin at ucar.edu>
> wrote:
>
> > Hi Brad,
> >
> > Can you try "L0", and "P0"? If that doesn't work, could you
please send
> > me your grib file? Instructions are as follows:
> >
> > Non-Windows FTP client:
> >
> >
> >
> >
> >
> >
> >
> >
> >
> > * ftp -p ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
username=anonymous
> > password="your email address" cd incoming/irap/met_help mkdir
> > "surname_data" # where surname is your last name cd
"surname_data" put
> > "your files" bye *
> >
> > Windows FTP client:
> >
> >
> >
> >
> >
> >
> >
> >
> > * ftp ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
username=anonymous
> > password="your email address" quote pasv cd incoming/irap/met_help
mkdir
> > "surname_data" # where surname is your last name cd "surname_data"
put
> > "your files" bye*
> >
> > I can take a look at the level type and value, you might need to
use
> > special syntax in your config file to get grid_stat to work. Let
me know
> > if you ftp'd a grib file so I can grab it before the file gets
scrubbed.
> >
> > Thanks,
> > Minna
> >
> >
> > ---------------
> > Minna Win
> > NCAR
> > Research Applications Lab
> > Phone: 303-497-8423 <(303)%20497-8423>
> > Fax: 302-497-8401 <(302)%20497-8401>
> >
> >
> > On Thu, Aug 31, 2017 at 3:09 PM, Brad Diehl - NOAA Federal via RT
<
> > met_help at ucar.edu> wrote:
> >
> >>
> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> >>
> >> After running pcp_combine, everything did work correctly, it is
the
> grib2
> >> fcts/obs part I'm having trouble with. I think you put me on the
right
> >> track by suggesting wgrib2. Here's what I get using wgrib2 on the
fcst
> >> file:
> >>
> >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ wgrib2
> >> p06m_2017082312f006.grb2
> >> 1:21:d=2017082312 <(201)%20708-2312>:APCP:surface -
surface:CodeTable
> >> 4.11=missing:
> >>
> >> in the config file I'm using:
> >>
> >>
> >> fcst = {
> >>
> >> field = [
> >> {
> >> name = "APCP";
> >> level = "Z0" ;
> >> cat_thresh = [ >0.0, >=5.0 , >=10.0];
> >> }
> >> ];
> >>
> >> for level, should I be using surface or Z0?, neither of which
seemed to
> >> work.
> >>
> >> On Thu, Aug 31, 2017 at 2:58 PM, Minna Win via RT
<met_help at ucar.edu>
> >> wrote:
> >>
> >> > Hi Brad,
> >> >
> >> > Please disregard my previous response. I didn't see that you
had run
> >> > pcp_combine-that *should* have worked. I'll look a little
further or
> >> see
> >> > if someone else has seen this before.
> >> >
> >> >
> >> >
> >> > Regards,
> >> > Minna
> >> >
> >> > ---------------
> >> > Minna Win
> >> > NCAR
> >> > Research Applications Lab
> >> > Phone: 303-497-8423
> >> > Fax: 302-497-8401
> >> >
> >> >
> >> > On Thu, Aug 31, 2017 at 2:50 PM, Minna Win-Gildenmeister <
> >> > minnawin at ucar.edu>
> >> > wrote:
> >> >
> >> > > HI Brad,
> >> > >
> >> > > wgrib is an NDFD tool and I don't know how it is interpreting
the
> APCP
> >> > > variable in your grib file. Try using wgrib2 on your grib
file and
> >> grep
> >> > for
> >> > > APCP and use that variable abbreviation in you MET Grid_stat
config
> >> > file. I
> >> > > hope that helps. Please let us know if you still have
problems
> using
> >> > > Grid_stat.
> >> > >
> >> > > Regards,
> >> > > Minna
> >> > >
> >> > > ---------------
> >> > > Minna Win
> >> > > NCAR
> >> > > Research Applications Lab
> >> > > Phone: 303-497-8423 <(303)%20497-8423>
> >> > > Fax: 302-497-8401 <(302)%20497-8401>
> >> > >
> >> > >
> >> > > On Thu, Aug 31, 2017 at 2:33 PM, Brad Diehl - NOAA Federal
via RT <
> >> > > met_help at ucar.edu> wrote:
> >> > >
> >> > >>
> >> > >> Thu Aug 31 08:33:00 2017: Request 81829 was acted upon.
> >> > >> Transaction: Ticket created by brad.diehl at noaa.gov
> >> > >> Queue: met_help
> >> > >> Subject: grib2 files as input into MET grid_stat
> >> > >> Owner: Nobody
> >> > >> Requestors: brad.diehl at noaa.gov
> >> > >> Status: new
> >> > >> Ticket <URL: https://rt.rap.ucar.edu/rt/Tic
> >> ket/Display.html?id=81829 >
> >> > >>
> >> > >>
> >> > >> I tried to use grid_stat with two grib2 files as input, and
> received
> >> an
> >> > >> error regarding the GRIB2 field abbreviation 'APCP06' (see
output
> at
> >> > >> bottom
> >> > >> of message). Yet when I run degrib on the fcst and obs files
that's
> >> > what I
> >> > >> get:
> >> > >>
> >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
> >> > >> p06m_2017082312f006.grb2 -I
> >> > >> MsgNum, Byte, GRIB-Version, elem, level, reference(UTC),
> valid(UTC),
> >> > >> Proj(hr)
> >> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation [kg/(m^2)]", 0-
SFC,
> >> > >> 08/23/2017
> >> > >> 12:00, 08/23/2017 18:00, 6.00
> >> > >>
> >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
> >> > >> urma2p5.2017082318 <(201)%20708-2318>.pcp_06h.184.grb2 -I
> >> > >> MsgNum, Byte, GRIB-Version, elem, level, reference(UTC),
> valid(UTC),
> >> > >> Proj(hr)
> >> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation [kg/(m^2)]", 0-
SFC,
> >> > >> 08/23/2017
> >> > >> 12:00, 08/23/2017 18:00, 6.00
> >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> >> > >>
> >> > >> If I put both the fcst and obs files through pcp_combine and
use
> the
> >> > >> resultant .nc files, grid_stat runs fine, with the caveat
that the
> >> field
> >> > >> name becomes 'APCP_06'.
> >> > >>
> >> > >> Thanks for any comments/suggestions, Brad Diehl
> >> > >>
> >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ grid_stat
> >> > >> p06m_2017082312f006.grb2 urma2p5.2017082318 <(201)%20708-
2318>
> >> > >> .pcp_06h.184.grb2
> >> > >> config/GridStatConfig_APCP_06_grib -outdir out/grid_stat -v
4DEBUG
> >> 1:
> >> > >> Default Config File:
> >> > >>
/usrx/local/dev/MET/6.0/share/met/config/GridStatConfig_default
> >> > >> DEBUG 1: User Config File:
config/GridStatConfig_APCP_06_grib
> >> > >> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() ->
created
> new
> >> > >> Met2dDataFile object of type "FileType_Gb2".
> >> > >> WARNING:
> >> > >> WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does not
> >> currently
> >> > >> support Lambert Conformal grids where dx (2.5397) != dy
(2.5397)
> and
> >> may
> >> > >> produce unexpected results!
> >> > >> WARNING:
> >> > >> DEBUG 4:
> >> > >> DEBUG 4: Lambert Conformal Grid Data:
> >> > >> DEBUG 4: scale_lat_1: 25
> >> > >> DEBUG 4: scale_lat_2: 25
> >> > >> DEBUG 4: lat_pin: 20.192
> >> > >> DEBUG 4: lon_pin: 121.554
> >> > >> DEBUG 4: x_pin: 0
> >> > >> DEBUG 4: y_pin: 0
> >> > >> DEBUG 4: lon_orient: 95
> >> > >> DEBUG 4: d_km: 2.5397
> >> > >> DEBUG 4: r_km: 6371.2
> >> > >> DEBUG 4: nx: 2145
> >> > >> DEBUG 4: ny: 1597
> >> > >> DEBUG 4:
> >> > >> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() ->
created
> new
> >> > >> Met2dDataFile object of type "FileType_Gb2".
> >> > >> DEBUG 4:
> >> > >> DEBUG 4: Lambert Conformal Grid Data:
> >> > >> DEBUG 4: scale_lat_1: 25
> >> > >> DEBUG 4: scale_lat_2: 25
> >> > >> DEBUG 4: lat_pin: 20.192
> >> > >> DEBUG 4: lon_pin: 121.554
> >> > >> DEBUG 4: x_pin: 0
> >> > >> DEBUG 4: y_pin: 0
> >> > >> DEBUG 4: lon_orient: 95
> >> > >> DEBUG 4: d_km: 2.5397
> >> > >> DEBUG 4: r_km: 6371.2
> >> > >> DEBUG 4: nx: 2145
> >> > >> DEBUG 4: ny: 1377
> >> > >> DEBUG 4:
> >> > >> DEBUG 4: VarInfoFactory::new_var_info() -> created new
VarInfo
> >> object of
> >> > >> type "FileType_Gb2".
> >> > >> DEBUG 4: VarInfoFactory::new_var_info() -> created new
VarInfo
> >> object of
> >> > >> type "FileType_Gb2".
> >> > >> ERROR :
> >> > >> ERROR : VarInfoGrib2::set_dict() -> unrecognized GRIB2
field
> >> > abbreviation
> >> > >> 'APCP06'
> >> > >> ERROR :
> >> > >>
> >> > >>
> >> > >
> >> >
> >> >
> >>
> >>
> >
>
>
------------------------------------------------
Subject: grib2 files as input into MET grid_stat
From: Minna Win
Time: Thu Aug 31 12:25:44 2017
Hi Brad,
For GRIB2, the level (if not an accumulated precip type) is typically
"Z0"
for surface (sometimes "L0" works, too), but for other levels such as
entire atmosphere, you will need to refer to the GRIB2 table "Table 2.
Parameters and Units":
http://www.wmo.int/pages/prog/www/WDM/Guides/Guide-binary-2.html#GRIB
.
>From that table, find the value that corresponds to your level type
and set
it to GRIB_lvl_typ and then set GRIB_lvl_val1 to the appropriate
level.
For 'entire atmosphere', there is no numerical value for the level, so
it
would be set to 0 or 255 (you might have to try both to see what
works, it
depends on the grib file). For something like a pressure difference,
you
would need to set GRIB_lvl_typ and the GRIB_lvl_va1 and
GRIB_lvl_val2,
i.e. GRIB_lvl_val2=top level GRIB_lvl_val1=bottom level, which you
would
obtain from performing a wgrib2/grep combination on your grib file and
variable of interest. You can use the MET tool plot_data_plane to
check
that your GRIB_lvl_typ and GRIB_lvl_val1 (and GRIB_lvl_val2 if used)
produced a match to the record of interest in the grib file.
As for your question about the (*,*) you see in some examples, there
is an
explanation in Chapter 3 of the MET User's Manual (excerpted):
> // - The "level" entry specifies the dimensions to be used:
> // - (i,...,j,*,*) for a single field, where i,...,j specifies
fixed
> // dimension values and *,* specifies the two dimensions for
the
> // gridded field.
>
field = [
{
name = "QVAPOR";
level = "( = "(0,5,*,*)"
},
{
name = "TMP_P850_ENS_MEAN;
level = [ "(*,*)" ];
}
];
I hope that clarifies things. Please don't hesitate to contact us if
you
have further questions.
Regards,
Minna
---------------
Minna Win
NCAR
Research Applications Lab
Phone: 303-497-8423
Fax: 302-497-8401
On Thu, Aug 31, 2017 at 4:21 PM, Brad Diehl - NOAA Federal via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
>
> the A06 worked. Prior to trying that I sent you via mailx the
fcst/obs
> grib2 files, you can ignore that now.
>
> But it does beg the question, how does one find the level names, is
it in
> the grib file itself? Very confusing, at least to me. And what does
*,*
> mean when I see that in some examples?
>
> Thanks for all your help........Brad
>
> On Thu, Aug 31, 2017 at 4:11 PM, Minna Win via RT
<met_help at ucar.edu>
> wrote:
>
> > Hi Brad,
> >
> > Tara Jensen brought this to my attention, for verify APCP, you
need to
> set
> > the level to the value pre-pended with the letter 'A'. For
instance,
> for 6
> > hour accumulation level, set level="A06".
> >
> > Thanks,
> > Minna
> >
> > ---------------
> > Minna Win
> > NCAR
> > Research Applications Lab
> > Phone: 303-497-8423
> > Fax: 302-497-8401
> >
> >
> > On Thu, Aug 31, 2017 at 3:52 PM, Minna Win-Gildenmeister <
> > minnawin at ucar.edu>
> > wrote:
> >
> > > Hi Brad,
> > >
> > > Can you try "L0", and "P0"? If that doesn't work, could you
please
> send
> > > me your grib file? Instructions are as follows:
> > >
> > > Non-Windows FTP client:
> > >
> > >
> > >
> > >
> > >
> > >
> > >
> > >
> > >
> > > * ftp -p ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
username=anonymous
> > > password="your email address" cd incoming/irap/met_help mkdir
> > > "surname_data" # where surname is your last name cd
"surname_data"
> put
> > > "your files" bye *
> > >
> > > Windows FTP client:
> > >
> > >
> > >
> > >
> > >
> > >
> > >
> > >
> > > * ftp ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
username=anonymous
> > > password="your email address" quote pasv cd
incoming/irap/met_help
> mkdir
> > > "surname_data" # where surname is your last name cd
"surname_data" put
> > > "your files" bye*
> > >
> > > I can take a look at the level type and value, you might need to
use
> > > special syntax in your config file to get grid_stat to work.
Let me
> know
> > > if you ftp'd a grib file so I can grab it before the file gets
> scrubbed.
> > >
> > > Thanks,
> > > Minna
> > >
> > >
> > > ---------------
> > > Minna Win
> > > NCAR
> > > Research Applications Lab
> > > Phone: 303-497-8423 <(303)%20497-8423>
> > > Fax: 302-497-8401 <(302)%20497-8401>
> > >
> > >
> > > On Thu, Aug 31, 2017 at 3:09 PM, Brad Diehl - NOAA Federal via
RT <
> > > met_help at ucar.edu> wrote:
> > >
> > >>
> > >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> > >>
> > >> After running pcp_combine, everything did work correctly, it is
the
> > grib2
> > >> fcts/obs part I'm having trouble with. I think you put me on
the right
> > >> track by suggesting wgrib2. Here's what I get using wgrib2 on
the fcst
> > >> file:
> > >>
> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ wgrib2
> > >> p06m_2017082312f006.grb2
> > >> 1:21:d=2017082312 <(201)%20708-2312>:APCP:surface -
surface:CodeTable
> > >> 4.11=missing:
> > >>
> > >> in the config file I'm using:
> > >>
> > >>
> > >> fcst = {
> > >>
> > >> field = [
> > >> {
> > >> name = "APCP";
> > >> level = "Z0" ;
> > >> cat_thresh = [ >0.0, >=5.0 , >=10.0];
> > >> }
> > >> ];
> > >>
> > >> for level, should I be using surface or Z0?, neither of which
seemed
> to
> > >> work.
> > >>
> > >> On Thu, Aug 31, 2017 at 2:58 PM, Minna Win via RT
<met_help at ucar.edu>
> > >> wrote:
> > >>
> > >> > Hi Brad,
> > >> >
> > >> > Please disregard my previous response. I didn't see that you
had
> run
> > >> > pcp_combine-that *should* have worked. I'll look a little
further
> or
> > >> see
> > >> > if someone else has seen this before.
> > >> >
> > >> >
> > >> >
> > >> > Regards,
> > >> > Minna
> > >> >
> > >> > ---------------
> > >> > Minna Win
> > >> > NCAR
> > >> > Research Applications Lab
> > >> > Phone: 303-497-8423
> > >> > Fax: 302-497-8401
> > >> >
> > >> >
> > >> > On Thu, Aug 31, 2017 at 2:50 PM, Minna Win-Gildenmeister <
> > >> > minnawin at ucar.edu>
> > >> > wrote:
> > >> >
> > >> > > HI Brad,
> > >> > >
> > >> > > wgrib is an NDFD tool and I don't know how it is
interpreting the
> > APCP
> > >> > > variable in your grib file. Try using wgrib2 on your grib
file and
> > >> grep
> > >> > for
> > >> > > APCP and use that variable abbreviation in you MET
Grid_stat
> config
> > >> > file. I
> > >> > > hope that helps. Please let us know if you still have
problems
> > using
> > >> > > Grid_stat.
> > >> > >
> > >> > > Regards,
> > >> > > Minna
> > >> > >
> > >> > > ---------------
> > >> > > Minna Win
> > >> > > NCAR
> > >> > > Research Applications Lab
> > >> > > Phone: 303-497-8423 <(303)%20497-8423>
> > >> > > Fax: 302-497-8401 <(302)%20497-8401>
> > >> > >
> > >> > >
> > >> > > On Thu, Aug 31, 2017 at 2:33 PM, Brad Diehl - NOAA Federal
via RT
> <
> > >> > > met_help at ucar.edu> wrote:
> > >> > >
> > >> > >>
> > >> > >> Thu Aug 31 08:33:00 2017: Request 81829 was acted upon.
> > >> > >> Transaction: Ticket created by brad.diehl at noaa.gov
> > >> > >> Queue: met_help
> > >> > >> Subject: grib2 files as input into MET grid_stat
> > >> > >> Owner: Nobody
> > >> > >> Requestors: brad.diehl at noaa.gov
> > >> > >> Status: new
> > >> > >> Ticket <URL: https://rt.rap.ucar.edu/rt/Tic
> > >> ket/Display.html?id=81829 >
> > >> > >>
> > >> > >>
> > >> > >> I tried to use grid_stat with two grib2 files as input,
and
> > received
> > >> an
> > >> > >> error regarding the GRIB2 field abbreviation 'APCP06' (see
output
> > at
> > >> > >> bottom
> > >> > >> of message). Yet when I run degrib on the fcst and obs
files
> that's
> > >> > what I
> > >> > >> get:
> > >> > >>
> > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
> > >> > >> p06m_2017082312f006.grb2 -I
> > >> > >> MsgNum, Byte, GRIB-Version, elem, level, reference(UTC),
> > valid(UTC),
> > >> > >> Proj(hr)
> > >> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation [kg/(m^2)]",
0-SFC,
> > >> > >> 08/23/2017
> > >> > >> 12:00, 08/23/2017 18:00, 6.00
> > >> > >>
> > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
> > >> > >> urma2p5.2017082318 <(201)%20708-2318>.pcp_06h.184.grb2 -I
> > >> > >> MsgNum, Byte, GRIB-Version, elem, level, reference(UTC),
> > valid(UTC),
> > >> > >> Proj(hr)
> > >> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation [kg/(m^2)]",
0-SFC,
> > >> > >> 08/23/2017
> > >> > >> 12:00, 08/23/2017 18:00, 6.00
> > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > >> > >>
> > >> > >> If I put both the fcst and obs files through pcp_combine
and use
> > the
> > >> > >> resultant .nc files, grid_stat runs fine, with the caveat
that
> the
> > >> field
> > >> > >> name becomes 'APCP_06'.
> > >> > >>
> > >> > >> Thanks for any comments/suggestions, Brad Diehl
> > >> > >>
> > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ grid_stat
> > >> > >> p06m_2017082312f006.grb2 urma2p5.2017082318 <(201)%20708-
2318>
> > >> > >> .pcp_06h.184.grb2
> > >> > >> config/GridStatConfig_APCP_06_grib -outdir out/grid_stat
-v
> 4DEBUG
> > >> 1:
> > >> > >> Default Config File:
> > >> > >>
/usrx/local/dev/MET/6.0/share/met/config/GridStatConfig_default
> > >> > >> DEBUG 1: User Config File:
config/GridStatConfig_APCP_06_grib
> > >> > >> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() ->
created
> > new
> > >> > >> Met2dDataFile object of type "FileType_Gb2".
> > >> > >> WARNING:
> > >> > >> WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does
not
> > >> currently
> > >> > >> support Lambert Conformal grids where dx (2.5397) != dy
(2.5397)
> > and
> > >> may
> > >> > >> produce unexpected results!
> > >> > >> WARNING:
> > >> > >> DEBUG 4:
> > >> > >> DEBUG 4: Lambert Conformal Grid Data:
> > >> > >> DEBUG 4: scale_lat_1: 25
> > >> > >> DEBUG 4: scale_lat_2: 25
> > >> > >> DEBUG 4: lat_pin: 20.192
> > >> > >> DEBUG 4: lon_pin: 121.554
> > >> > >> DEBUG 4: x_pin: 0
> > >> > >> DEBUG 4: y_pin: 0
> > >> > >> DEBUG 4: lon_orient: 95
> > >> > >> DEBUG 4: d_km: 2.5397
> > >> > >> DEBUG 4: r_km: 6371.2
> > >> > >> DEBUG 4: nx: 2145
> > >> > >> DEBUG 4: ny: 1597
> > >> > >> DEBUG 4:
> > >> > >> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() ->
created
> > new
> > >> > >> Met2dDataFile object of type "FileType_Gb2".
> > >> > >> DEBUG 4:
> > >> > >> DEBUG 4: Lambert Conformal Grid Data:
> > >> > >> DEBUG 4: scale_lat_1: 25
> > >> > >> DEBUG 4: scale_lat_2: 25
> > >> > >> DEBUG 4: lat_pin: 20.192
> > >> > >> DEBUG 4: lon_pin: 121.554
> > >> > >> DEBUG 4: x_pin: 0
> > >> > >> DEBUG 4: y_pin: 0
> > >> > >> DEBUG 4: lon_orient: 95
> > >> > >> DEBUG 4: d_km: 2.5397
> > >> > >> DEBUG 4: r_km: 6371.2
> > >> > >> DEBUG 4: nx: 2145
> > >> > >> DEBUG 4: ny: 1377
> > >> > >> DEBUG 4:
> > >> > >> DEBUG 4: VarInfoFactory::new_var_info() -> created new
VarInfo
> > >> object of
> > >> > >> type "FileType_Gb2".
> > >> > >> DEBUG 4: VarInfoFactory::new_var_info() -> created new
VarInfo
> > >> object of
> > >> > >> type "FileType_Gb2".
> > >> > >> ERROR :
> > >> > >> ERROR : VarInfoGrib2::set_dict() -> unrecognized GRIB2
field
> > >> > abbreviation
> > >> > >> 'APCP06'
> > >> > >> ERROR :
> > >> > >>
> > >> > >>
> > >> > >
> > >> >
> > >> >
> > >>
> > >>
> > >
> >
> >
>
>
------------------------------------------------
Subject: grib2 files as input into MET grid_stat
From: Julie Prestopnik
Time: Thu Aug 31 12:35:19 2017
Hi Brad,
I just thought I'd add that you can find more information for file-
format
specific settings for the "field" entry (including level) for MET on
pg. 53
of the MET User's Guide for met-6.0:
http://www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v6.0.pdf
or in the README file in the met-6.0/data/config area of your install.
For example:
// File-format specific settings for the "field" entry:
//
// - GRIB1 and GRIB2:
// - For custom GRIB tables, see note about MET_GRIB_TABLES.
// - The "name" entry specifies a GRIB code number or
abbreviation.
// - GRIB1 Product Definition Section:
// http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
// - GRIB2 Product Definition Section:
//
http://www.nco.ncep.noaa.gov/pmb/docs/grib2/grib2_doc.shtml
// - The "level" entry specifies a level type and value:
// - ANNN for accumulation interval NNN
// - ZNNN for vertical level NNN
// - ZNNN-NNN for a range of vertical levels
// - PNNN for pressure level NNN in hPa
// - PNNN-NNN for a range of pressure levels in hPa
// - LNNN for a generic level type
// - RNNN for a specific GRIB record number
I hope that helps!
Julie
On Thu, Aug 31, 2017 at 12:25 PM, Minna Win via RT <met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
>
> Hi Brad,
>
> For GRIB2, the level (if not an accumulated precip type) is
typically "Z0"
> for surface (sometimes "L0" works, too), but for other levels such
as
> entire atmosphere, you will need to refer to the GRIB2 table "Table
2.
> Parameters and Units":
> http://www.wmo.int/pages/prog/www/WDM/Guides/Guide-binary-
2.html#GRIB .
> From that table, find the value that corresponds to your level type
and set
> it to GRIB_lvl_typ and then set GRIB_lvl_val1 to the appropriate
level.
> For 'entire atmosphere', there is no numerical value for the level,
so it
> would be set to 0 or 255 (you might have to try both to see what
works, it
> depends on the grib file). For something like a pressure
difference, you
> would need to set GRIB_lvl_typ and the GRIB_lvl_va1 and
GRIB_lvl_val2,
> i.e. GRIB_lvl_val2=top level GRIB_lvl_val1=bottom level, which you
would
> obtain from performing a wgrib2/grep combination on your grib file
and
> variable of interest. You can use the MET tool plot_data_plane to
check
> that your GRIB_lvl_typ and GRIB_lvl_val1 (and GRIB_lvl_val2 if used)
> produced a match to the record of interest in the grib file.
>
> As for your question about the (*,*) you see in some examples, there
is an
> explanation in Chapter 3 of the MET User's Manual (excerpted):
>
>
> > // - The "level" entry specifies the dimensions to be used:
> > // - (i,...,j,*,*) for a single field, where i,...,j specifies
fixed
> > // dimension values and *,* specifies the two dimensions for
the
> > // gridded field.
> >
>
>
> field = [
> {
> name = "QVAPOR";
> level = "( = "(0,5,*,*)"
>
> },
> {
>
> name = "TMP_P850_ENS_MEAN;
> level = [ "(*,*)" ];
> }
> ];
>
>
> I hope that clarifies things. Please don't hesitate to contact us
if you
> have further questions.
>
> Regards,
> Minna
>
> ---------------
> Minna Win
> NCAR
> Research Applications Lab
> Phone: 303-497-8423
> Fax: 302-497-8401
>
>
> On Thu, Aug 31, 2017 at 4:21 PM, Brad Diehl - NOAA Federal via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> >
> > the A06 worked. Prior to trying that I sent you via mailx the
fcst/obs
> > grib2 files, you can ignore that now.
> >
> > But it does beg the question, how does one find the level names,
is it in
> > the grib file itself? Very confusing, at least to me. And what
does *,*
> > mean when I see that in some examples?
> >
> > Thanks for all your help........Brad
> >
> > On Thu, Aug 31, 2017 at 4:11 PM, Minna Win via RT
<met_help at ucar.edu>
> > wrote:
> >
> > > Hi Brad,
> > >
> > > Tara Jensen brought this to my attention, for verify APCP, you
need to
> > set
> > > the level to the value pre-pended with the letter 'A'. For
instance,
> > for 6
> > > hour accumulation level, set level="A06".
> > >
> > > Thanks,
> > > Minna
> > >
> > > ---------------
> > > Minna Win
> > > NCAR
> > > Research Applications Lab
> > > Phone: 303-497-8423
> > > Fax: 302-497-8401
> > >
> > >
> > > On Thu, Aug 31, 2017 at 3:52 PM, Minna Win-Gildenmeister <
> > > minnawin at ucar.edu>
> > > wrote:
> > >
> > > > Hi Brad,
> > > >
> > > > Can you try "L0", and "P0"? If that doesn't work, could you
please
> > send
> > > > me your grib file? Instructions are as follows:
> > > >
> > > > Non-Windows FTP client:
> > > >
> > > >
> > > >
> > > >
> > > >
> > > >
> > > >
> > > >
> > > >
> > > > * ftp -p ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
> username=anonymous
> > > > password="your email address" cd incoming/irap/met_help
mkdir
> > > > "surname_data" # where surname is your last name cd
"surname_data"
> > put
> > > > "your files" bye *
> > > >
> > > > Windows FTP client:
> > > >
> > > >
> > > >
> > > >
> > > >
> > > >
> > > >
> > > >
> > > > * ftp ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
username=anonymous
> > > > password="your email address" quote pasv cd
incoming/irap/met_help
> > mkdir
> > > > "surname_data" # where surname is your last name cd
"surname_data"
> put
> > > > "your files" bye*
> > > >
> > > > I can take a look at the level type and value, you might need
to use
> > > > special syntax in your config file to get grid_stat to work.
Let me
> > know
> > > > if you ftp'd a grib file so I can grab it before the file gets
> > scrubbed.
> > > >
> > > > Thanks,
> > > > Minna
> > > >
> > > >
> > > > ---------------
> > > > Minna Win
> > > > NCAR
> > > > Research Applications Lab
> > > > Phone: 303-497-8423 <(303)%20497-8423>
> > > > Fax: 302-497-8401 <(302)%20497-8401>
> > > >
> > > >
> > > > On Thu, Aug 31, 2017 at 3:09 PM, Brad Diehl - NOAA Federal via
RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > >>
> > > >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829
>
> > > >>
> > > >> After running pcp_combine, everything did work correctly, it
is the
> > > grib2
> > > >> fcts/obs part I'm having trouble with. I think you put me on
the
> right
> > > >> track by suggesting wgrib2. Here's what I get using wgrib2 on
the
> fcst
> > > >> file:
> > > >>
> > > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ wgrib2
> > > >> p06m_2017082312f006.grb2
> > > >> 1:21:d=2017082312 <(201)%20708-2312>:APCP:surface -
> surface:CodeTable
> > > >> 4.11=missing:
> > > >>
> > > >> in the config file I'm using:
> > > >>
> > > >>
> > > >> fcst = {
> > > >>
> > > >> field = [
> > > >> {
> > > >> name = "APCP";
> > > >> level = "Z0" ;
> > > >> cat_thresh = [ >0.0, >=5.0 , >=10.0];
> > > >> }
> > > >> ];
> > > >>
> > > >> for level, should I be using surface or Z0?, neither of which
seemed
> > to
> > > >> work.
> > > >>
> > > >> On Thu, Aug 31, 2017 at 2:58 PM, Minna Win via RT <
> met_help at ucar.edu>
> > > >> wrote:
> > > >>
> > > >> > Hi Brad,
> > > >> >
> > > >> > Please disregard my previous response. I didn't see that
you had
> > run
> > > >> > pcp_combine-that *should* have worked. I'll look a little
further
> > or
> > > >> see
> > > >> > if someone else has seen this before.
> > > >> >
> > > >> >
> > > >> >
> > > >> > Regards,
> > > >> > Minna
> > > >> >
> > > >> > ---------------
> > > >> > Minna Win
> > > >> > NCAR
> > > >> > Research Applications Lab
> > > >> > Phone: 303-497-8423
> > > >> > Fax: 302-497-8401
> > > >> >
> > > >> >
> > > >> > On Thu, Aug 31, 2017 at 2:50 PM, Minna Win-Gildenmeister <
> > > >> > minnawin at ucar.edu>
> > > >> > wrote:
> > > >> >
> > > >> > > HI Brad,
> > > >> > >
> > > >> > > wgrib is an NDFD tool and I don't know how it is
interpreting
> the
> > > APCP
> > > >> > > variable in your grib file. Try using wgrib2 on your grib
file
> and
> > > >> grep
> > > >> > for
> > > >> > > APCP and use that variable abbreviation in you MET
Grid_stat
> > config
> > > >> > file. I
> > > >> > > hope that helps. Please let us know if you still have
problems
> > > using
> > > >> > > Grid_stat.
> > > >> > >
> > > >> > > Regards,
> > > >> > > Minna
> > > >> > >
> > > >> > > ---------------
> > > >> > > Minna Win
> > > >> > > NCAR
> > > >> > > Research Applications Lab
> > > >> > > Phone: 303-497-8423 <(303)%20497-8423>
> > > >> > > Fax: 302-497-8401 <(302)%20497-8401>
> > > >> > >
> > > >> > >
> > > >> > > On Thu, Aug 31, 2017 at 2:33 PM, Brad Diehl - NOAA
Federal via
> RT
> > <
> > > >> > > met_help at ucar.edu> wrote:
> > > >> > >
> > > >> > >>
> > > >> > >> Thu Aug 31 08:33:00 2017: Request 81829 was acted upon.
> > > >> > >> Transaction: Ticket created by brad.diehl at noaa.gov
> > > >> > >> Queue: met_help
> > > >> > >> Subject: grib2 files as input into MET grid_stat
> > > >> > >> Owner: Nobody
> > > >> > >> Requestors: brad.diehl at noaa.gov
> > > >> > >> Status: new
> > > >> > >> Ticket <URL: https://rt.rap.ucar.edu/rt/Tic
> > > >> ket/Display.html?id=81829 >
> > > >> > >>
> > > >> > >>
> > > >> > >> I tried to use grid_stat with two grib2 files as input,
and
> > > received
> > > >> an
> > > >> > >> error regarding the GRIB2 field abbreviation 'APCP06'
(see
> output
> > > at
> > > >> > >> bottom
> > > >> > >> of message). Yet when I run degrib on the fcst and obs
files
> > that's
> > > >> > what I
> > > >> > >> get:
> > > >> > >>
> > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
> > > >> > >> p06m_2017082312f006.grb2 -I
> > > >> > >> MsgNum, Byte, GRIB-Version, elem, level, reference(UTC),
> > > valid(UTC),
> > > >> > >> Proj(hr)
> > > >> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation
[kg/(m^2)]",
> 0-SFC,
> > > >> > >> 08/23/2017
> > > >> > >> 12:00, 08/23/2017 18:00, 6.00
> > > >> > >>
> > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ degrib
> > > >> > >> urma2p5.2017082318 <(201)%20708-2318>.pcp_06h.184.grb2
-I
> > > >> > >> MsgNum, Byte, GRIB-Version, elem, level, reference(UTC),
> > > valid(UTC),
> > > >> > >> Proj(hr)
> > > >> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation
[kg/(m^2)]",
> 0-SFC,
> > > >> > >> 08/23/2017
> > > >> > >> 12:00, 08/23/2017 18:00, 6.00
> > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > > >> > >>
> > > >> > >> If I put both the fcst and obs files through pcp_combine
and
> use
> > > the
> > > >> > >> resultant .nc files, grid_stat runs fine, with the
caveat that
> > the
> > > >> field
> > > >> > >> name becomes 'APCP_06'.
> > > >> > >>
> > > >> > >> Thanks for any comments/suggestions, Brad Diehl
> > > >> > >>
> > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
grid_stat
> > > >> > >> p06m_2017082312f006.grb2 urma2p5.2017082318
<(201)%20708-2318>
> > > >> > >> .pcp_06h.184.grb2
> > > >> > >> config/GridStatConfig_APCP_06_grib -outdir out/grid_stat
-v
> > 4DEBUG
> > > >> 1:
> > > >> > >> Default Config File:
> > > >> > >> /usrx/local/dev/MET/6.0/share/met/config/GridStatConfig_
> default
> > > >> > >> DEBUG 1: User Config File:
config/GridStatConfig_APCP_06_grib
> > > >> > >> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() ->
> created
> > > new
> > > >> > >> Met2dDataFile object of type "FileType_Gb2".
> > > >> > >> WARNING:
> > > >> > >> WARNING: MetGrib2DataFile::read_grib2_grid() -> MET does
not
> > > >> currently
> > > >> > >> support Lambert Conformal grids where dx (2.5397) != dy
> (2.5397)
> > > and
> > > >> may
> > > >> > >> produce unexpected results!
> > > >> > >> WARNING:
> > > >> > >> DEBUG 4:
> > > >> > >> DEBUG 4: Lambert Conformal Grid Data:
> > > >> > >> DEBUG 4: scale_lat_1: 25
> > > >> > >> DEBUG 4: scale_lat_2: 25
> > > >> > >> DEBUG 4: lat_pin: 20.192
> > > >> > >> DEBUG 4: lon_pin: 121.554
> > > >> > >> DEBUG 4: x_pin: 0
> > > >> > >> DEBUG 4: y_pin: 0
> > > >> > >> DEBUG 4: lon_orient: 95
> > > >> > >> DEBUG 4: d_km: 2.5397
> > > >> > >> DEBUG 4: r_km: 6371.2
> > > >> > >> DEBUG 4: nx: 2145
> > > >> > >> DEBUG 4: ny: 1597
> > > >> > >> DEBUG 4:
> > > >> > >> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() ->
> created
> > > new
> > > >> > >> Met2dDataFile object of type "FileType_Gb2".
> > > >> > >> DEBUG 4:
> > > >> > >> DEBUG 4: Lambert Conformal Grid Data:
> > > >> > >> DEBUG 4: scale_lat_1: 25
> > > >> > >> DEBUG 4: scale_lat_2: 25
> > > >> > >> DEBUG 4: lat_pin: 20.192
> > > >> > >> DEBUG 4: lon_pin: 121.554
> > > >> > >> DEBUG 4: x_pin: 0
> > > >> > >> DEBUG 4: y_pin: 0
> > > >> > >> DEBUG 4: lon_orient: 95
> > > >> > >> DEBUG 4: d_km: 2.5397
> > > >> > >> DEBUG 4: r_km: 6371.2
> > > >> > >> DEBUG 4: nx: 2145
> > > >> > >> DEBUG 4: ny: 1377
> > > >> > >> DEBUG 4:
> > > >> > >> DEBUG 4: VarInfoFactory::new_var_info() -> created new
VarInfo
> > > >> object of
> > > >> > >> type "FileType_Gb2".
> > > >> > >> DEBUG 4: VarInfoFactory::new_var_info() -> created new
VarInfo
> > > >> object of
> > > >> > >> type "FileType_Gb2".
> > > >> > >> ERROR :
> > > >> > >> ERROR : VarInfoGrib2::set_dict() -> unrecognized GRIB2
field
> > > >> > abbreviation
> > > >> > >> 'APCP06'
> > > >> > >> ERROR :
> > > >> > >>
> > > >> > >>
> > > >> > >
> > > >> >
> > > >> >
> > > >>
> > > >>
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: grib2 files as input into MET grid_stat
From: Brad Diehl - NOAA Federal
Time: Thu Aug 31 13:15:05 2017
Thanks for your help and suggestions today. A question came up from
one of
our developers regarding netcdf and CF compliance. Is the output of
MET/6.0
netcdf CF compliant?
Thanks in advance.......Brad
On Thu, Aug 31, 2017 at 6:35 PM, Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> Hi Brad,
>
> I just thought I'd add that you can find more information for file-
format
> specific settings for the "field" entry (including level) for MET on
pg. 53
> of the MET User's Guide for met-6.0:
>
> http://www.dtcenter.org/met/users/docs/users_guide/MET_
> Users_Guide_v6.0.pdf
>
> or in the README file in the met-6.0/data/config area of your
install.
>
> For example:
> // File-format specific settings for the "field" entry:
> //
> // - GRIB1 and GRIB2:
> // - For custom GRIB tables, see note about MET_GRIB_TABLES.
> // - The "name" entry specifies a GRIB code number or
abbreviation.
> // - GRIB1 Product Definition Section:
> // http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
> // - GRIB2 Product Definition Section:
> //
http://www.nco.ncep.noaa.gov/pmb/docs/grib2/grib2_doc.shtml
> // - The "level" entry specifies a level type and value:
> // - ANNN for accumulation interval NNN
> // - ZNNN for vertical level NNN
> // - ZNNN-NNN for a range of vertical levels
> // - PNNN for pressure level NNN in hPa
> // - PNNN-NNN for a range of pressure levels in hPa
> // - LNNN for a generic level type
> // - RNNN for a specific GRIB record number
>
> I hope that helps!
>
> Julie
>
> On Thu, Aug 31, 2017 at 12:25 PM, Minna Win via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> >
> > Hi Brad,
> >
> > For GRIB2, the level (if not an accumulated precip type) is
typically
> "Z0"
> > for surface (sometimes "L0" works, too), but for other levels such
as
> > entire atmosphere, you will need to refer to the GRIB2 table
"Table 2.
> > Parameters and Units":
> > http://www.wmo.int/pages/prog/www/WDM/Guides/Guide-binary-
2.html#GRIB .
> > From that table, find the value that corresponds to your level
type and
> set
> > it to GRIB_lvl_typ and then set GRIB_lvl_val1 to the appropriate
level.
> > For 'entire atmosphere', there is no numerical value for the
level, so it
> > would be set to 0 or 255 (you might have to try both to see what
works,
> it
> > depends on the grib file). For something like a pressure
difference, you
> > would need to set GRIB_lvl_typ and the GRIB_lvl_va1 and
GRIB_lvl_val2,
> > i.e. GRIB_lvl_val2=top level GRIB_lvl_val1=bottom level, which you
would
> > obtain from performing a wgrib2/grep combination on your grib file
and
> > variable of interest. You can use the MET tool plot_data_plane to
check
> > that your GRIB_lvl_typ and GRIB_lvl_val1 (and GRIB_lvl_val2 if
used)
> > produced a match to the record of interest in the grib file.
> >
> > As for your question about the (*,*) you see in some examples,
there is
> an
> > explanation in Chapter 3 of the MET User's Manual (excerpted):
> >
> >
> > > // - The "level" entry specifies the dimensions to be used:
> > > // - (i,...,j,*,*) for a single field, where i,...,j
specifies fixed
> > > // dimension values and *,* specifies the two dimensions
for the
> > > // gridded field.
> > >
> >
> >
> > field = [
> > {
> > name = "QVAPOR";
> > level = "( = "(0,5,*,*)"
> >
> > },
> > {
> >
> > name = "TMP_P850_ENS_MEAN;
> > level = [ "(*,*)" ];
> > }
> > ];
> >
> >
> > I hope that clarifies things. Please don't hesitate to contact us
if you
> > have further questions.
> >
> > Regards,
> > Minna
> >
> > ---------------
> > Minna Win
> > NCAR
> > Research Applications Lab
> > Phone: 303-497-8423
> > Fax: 302-497-8401
> >
> >
> > On Thu, Aug 31, 2017 at 4:21 PM, Brad Diehl - NOAA Federal via RT
<
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> > >
> > > the A06 worked. Prior to trying that I sent you via mailx the
fcst/obs
> > > grib2 files, you can ignore that now.
> > >
> > > But it does beg the question, how does one find the level names,
is it
> in
> > > the grib file itself? Very confusing, at least to me. And what
does *,*
> > > mean when I see that in some examples?
> > >
> > > Thanks for all your help........Brad
> > >
> > > On Thu, Aug 31, 2017 at 4:11 PM, Minna Win via RT
<met_help at ucar.edu>
> > > wrote:
> > >
> > > > Hi Brad,
> > > >
> > > > Tara Jensen brought this to my attention, for verify APCP, you
need
> to
> > > set
> > > > the level to the value pre-pended with the letter 'A'. For
instance,
> > > for 6
> > > > hour accumulation level, set level="A06".
> > > >
> > > > Thanks,
> > > > Minna
> > > >
> > > > ---------------
> > > > Minna Win
> > > > NCAR
> > > > Research Applications Lab
> > > > Phone: 303-497-8423
> > > > Fax: 302-497-8401
> > > >
> > > >
> > > > On Thu, Aug 31, 2017 at 3:52 PM, Minna Win-Gildenmeister <
> > > > minnawin at ucar.edu>
> > > > wrote:
> > > >
> > > > > Hi Brad,
> > > > >
> > > > > Can you try "L0", and "P0"? If that doesn't work, could you
please
> > > send
> > > > > me your grib file? Instructions are as follows:
> > > > >
> > > > > Non-Windows FTP client:
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > > * ftp -p ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
> > username=anonymous
> > > > > password="your email address" cd incoming/irap/met_help
mkdir
> > > > > "surname_data" # where surname is your last name cd
> "surname_data"
> > > put
> > > > > "your files" bye *
> > > > >
> > > > > Windows FTP client:
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > >
> > > > > * ftp ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
> username=anonymous
> > > > > password="your email address" quote pasv cd
incoming/irap/met_help
> > > mkdir
> > > > > "surname_data" # where surname is your last name cd
"surname_data"
> > put
> > > > > "your files" bye*
> > > > >
> > > > > I can take a look at the level type and value, you might
need to
> use
> > > > > special syntax in your config file to get grid_stat to work.
Let
> me
> > > know
> > > > > if you ftp'd a grib file so I can grab it before the file
gets
> > > scrubbed.
> > > > >
> > > > > Thanks,
> > > > > Minna
> > > > >
> > > > >
> > > > > ---------------
> > > > > Minna Win
> > > > > NCAR
> > > > > Research Applications Lab
> > > > > Phone: 303-497-8423 <(303)%20497-8423>
> > > > > Fax: 302-497-8401 <(302)%20497-8401>
> > > > >
> > > > >
> > > > > On Thu, Aug 31, 2017 at 3:09 PM, Brad Diehl - NOAA Federal
via RT <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > >>
> > > > >> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> > > > >>
> > > > >> After running pcp_combine, everything did work correctly,
it is
> the
> > > > grib2
> > > > >> fcts/obs part I'm having trouble with. I think you put me
on the
> > right
> > > > >> track by suggesting wgrib2. Here's what I get using wgrib2
on the
> > fcst
> > > > >> file:
> > > > >>
> > > > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ wgrib2
> > > > >> p06m_2017082312f006.grb2
> > > > >> 1:21:d=2017082312 <(201)%20708-2312>:APCP:surface -
> > surface:CodeTable
> > > > >> 4.11=missing:
> > > > >>
> > > > >> in the config file I'm using:
> > > > >>
> > > > >>
> > > > >> fcst = {
> > > > >>
> > > > >> field = [
> > > > >> {
> > > > >> name = "APCP";
> > > > >> level = "Z0" ;
> > > > >> cat_thresh = [ >0.0, >=5.0 , >=10.0];
> > > > >> }
> > > > >> ];
> > > > >>
> > > > >> for level, should I be using surface or Z0?, neither of
which
> seemed
> > > to
> > > > >> work.
> > > > >>
> > > > >> On Thu, Aug 31, 2017 at 2:58 PM, Minna Win via RT <
> > met_help at ucar.edu>
> > > > >> wrote:
> > > > >>
> > > > >> > Hi Brad,
> > > > >> >
> > > > >> > Please disregard my previous response. I didn't see that
you
> had
> > > run
> > > > >> > pcp_combine-that *should* have worked. I'll look a
little
> further
> > > or
> > > > >> see
> > > > >> > if someone else has seen this before.
> > > > >> >
> > > > >> >
> > > > >> >
> > > > >> > Regards,
> > > > >> > Minna
> > > > >> >
> > > > >> > ---------------
> > > > >> > Minna Win
> > > > >> > NCAR
> > > > >> > Research Applications Lab
> > > > >> > Phone: 303-497-8423
> > > > >> > Fax: 302-497-8401
> > > > >> >
> > > > >> >
> > > > >> > On Thu, Aug 31, 2017 at 2:50 PM, Minna Win-Gildenmeister
<
> > > > >> > minnawin at ucar.edu>
> > > > >> > wrote:
> > > > >> >
> > > > >> > > HI Brad,
> > > > >> > >
> > > > >> > > wgrib is an NDFD tool and I don't know how it is
interpreting
> > the
> > > > APCP
> > > > >> > > variable in your grib file. Try using wgrib2 on your
grib file
> > and
> > > > >> grep
> > > > >> > for
> > > > >> > > APCP and use that variable abbreviation in you MET
Grid_stat
> > > config
> > > > >> > file. I
> > > > >> > > hope that helps. Please let us know if you still have
> problems
> > > > using
> > > > >> > > Grid_stat.
> > > > >> > >
> > > > >> > > Regards,
> > > > >> > > Minna
> > > > >> > >
> > > > >> > > ---------------
> > > > >> > > Minna Win
> > > > >> > > NCAR
> > > > >> > > Research Applications Lab
> > > > >> > > Phone: 303-497-8423 <(303)%20497-8423>
> > > > >> > > Fax: 302-497-8401 <(302)%20497-8401>
> > > > >> > >
> > > > >> > >
> > > > >> > > On Thu, Aug 31, 2017 at 2:33 PM, Brad Diehl - NOAA
Federal via
> > RT
> > > <
> > > > >> > > met_help at ucar.edu> wrote:
> > > > >> > >
> > > > >> > >>
> > > > >> > >> Thu Aug 31 08:33:00 2017: Request 81829 was acted
upon.
> > > > >> > >> Transaction: Ticket created by brad.diehl at noaa.gov
> > > > >> > >> Queue: met_help
> > > > >> > >> Subject: grib2 files as input into MET grid_stat
> > > > >> > >> Owner: Nobody
> > > > >> > >> Requestors: brad.diehl at noaa.gov
> > > > >> > >> Status: new
> > > > >> > >> Ticket <URL: https://rt.rap.ucar.edu/rt/Tic
> > > > >> ket/Display.html?id=81829 >
> > > > >> > >>
> > > > >> > >>
> > > > >> > >> I tried to use grid_stat with two grib2 files as
input, and
> > > > received
> > > > >> an
> > > > >> > >> error regarding the GRIB2 field abbreviation 'APCP06'
(see
> > output
> > > > at
> > > > >> > >> bottom
> > > > >> > >> of message). Yet when I run degrib on the fcst and obs
files
> > > that's
> > > > >> > what I
> > > > >> > >> get:
> > > > >> > >>
> > > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
degrib
> > > > >> > >> p06m_2017082312f006.grb2 -I
> > > > >> > >> MsgNum, Byte, GRIB-Version, elem, level,
reference(UTC),
> > > > valid(UTC),
> > > > >> > >> Proj(hr)
> > > > >> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation
[kg/(m^2)]",
> > 0-SFC,
> > > > >> > >> 08/23/2017
> > > > >> > >> 12:00, 08/23/2017 18:00, 6.00
> > > > >> > >>
> > > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
degrib
> > > > >> > >> urma2p5.2017082318 <(201)%20708-2318>.pcp_06h.184.grb2
-I
> > > > >> > >> MsgNum, Byte, GRIB-Version, elem, level,
reference(UTC),
> > > > valid(UTC),
> > > > >> > >> Proj(hr)
> > > > >> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation
[kg/(m^2)]",
> > 0-SFC,
> > > > >> > >> 08/23/2017
> > > > >> > >> 12:00, 08/23/2017 18:00, 6.00
> > > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > > > >> > >>
> > > > >> > >> If I put both the fcst and obs files through
pcp_combine and
> > use
> > > > the
> > > > >> > >> resultant .nc files, grid_stat runs fine, with the
caveat
> that
> > > the
> > > > >> field
> > > > >> > >> name becomes 'APCP_06'.
> > > > >> > >>
> > > > >> > >> Thanks for any comments/suggestions, Brad Diehl
> > > > >> > >>
> > > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
grid_stat
> > > > >> > >> p06m_2017082312f006.grb2 urma2p5.2017082318
> <(201)%20708-2318>
> > > > >> > >> .pcp_06h.184.grb2
> > > > >> > >> config/GridStatConfig_APCP_06_grib -outdir
out/grid_stat -v
> > > 4DEBUG
> > > > >> 1:
> > > > >> > >> Default Config File:
> > > > >> > >>
/usrx/local/dev/MET/6.0/share/met/config/GridStatConfig_
> > default
> > > > >> > >> DEBUG 1: User Config File:
config/GridStatConfig_APCP_06_
> grib
> > > > >> > >> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file()
->
> > created
> > > > new
> > > > >> > >> Met2dDataFile object of type "FileType_Gb2".
> > > > >> > >> WARNING:
> > > > >> > >> WARNING: MetGrib2DataFile::read_grib2_grid() -> MET
does not
> > > > >> currently
> > > > >> > >> support Lambert Conformal grids where dx (2.5397) !=
dy
> > (2.5397)
> > > > and
> > > > >> may
> > > > >> > >> produce unexpected results!
> > > > >> > >> WARNING:
> > > > >> > >> DEBUG 4:
> > > > >> > >> DEBUG 4: Lambert Conformal Grid Data:
> > > > >> > >> DEBUG 4: scale_lat_1: 25
> > > > >> > >> DEBUG 4: scale_lat_2: 25
> > > > >> > >> DEBUG 4: lat_pin: 20.192
> > > > >> > >> DEBUG 4: lon_pin: 121.554
> > > > >> > >> DEBUG 4: x_pin: 0
> > > > >> > >> DEBUG 4: y_pin: 0
> > > > >> > >> DEBUG 4: lon_orient: 95
> > > > >> > >> DEBUG 4: d_km: 2.5397
> > > > >> > >> DEBUG 4: r_km: 6371.2
> > > > >> > >> DEBUG 4: nx: 2145
> > > > >> > >> DEBUG 4: ny: 1597
> > > > >> > >> DEBUG 4:
> > > > >> > >> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file()
->
> > created
> > > > new
> > > > >> > >> Met2dDataFile object of type "FileType_Gb2".
> > > > >> > >> DEBUG 4:
> > > > >> > >> DEBUG 4: Lambert Conformal Grid Data:
> > > > >> > >> DEBUG 4: scale_lat_1: 25
> > > > >> > >> DEBUG 4: scale_lat_2: 25
> > > > >> > >> DEBUG 4: lat_pin: 20.192
> > > > >> > >> DEBUG 4: lon_pin: 121.554
> > > > >> > >> DEBUG 4: x_pin: 0
> > > > >> > >> DEBUG 4: y_pin: 0
> > > > >> > >> DEBUG 4: lon_orient: 95
> > > > >> > >> DEBUG 4: d_km: 2.5397
> > > > >> > >> DEBUG 4: r_km: 6371.2
> > > > >> > >> DEBUG 4: nx: 2145
> > > > >> > >> DEBUG 4: ny: 1377
> > > > >> > >> DEBUG 4:
> > > > >> > >> DEBUG 4: VarInfoFactory::new_var_info() -> created new
> VarInfo
> > > > >> object of
> > > > >> > >> type "FileType_Gb2".
> > > > >> > >> DEBUG 4: VarInfoFactory::new_var_info() -> created new
> VarInfo
> > > > >> object of
> > > > >> > >> type "FileType_Gb2".
> > > > >> > >> ERROR :
> > > > >> > >> ERROR : VarInfoGrib2::set_dict() -> unrecognized
GRIB2 field
> > > > >> > abbreviation
> > > > >> > >> 'APCP06'
> > > > >> > >> ERROR :
> > > > >> > >>
> > > > >> > >>
> > > > >> > >
> > > > >> >
> > > > >> >
> > > > >>
> > > > >>
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: grib2 files as input into MET grid_stat
From: John Halley Gotway
Time: Fri Sep 01 09:35:19 2017
Hi Brad,
Unfortunately, no, the NetCDF output generated by the MET tools
themselves
are not CF-compliant.
Since we always suggest that others use that convention, we should
too!
And we do plan to make that switch over eventually, but haven't done
so yet.
Thanks,
John
On Thu, Aug 31, 2017 at 1:15 PM, Brad Diehl - NOAA Federal via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
>
> Thanks for your help and suggestions today. A question came up from
one of
> our developers regarding netcdf and CF compliance. Is the output of
MET/6.0
> netcdf CF compliant?
>
> Thanks in advance.......Brad
>
> On Thu, Aug 31, 2017 at 6:35 PM, Julie Prestopnik via RT <
> met_help at ucar.edu>
> wrote:
>
> > Hi Brad,
> >
> > I just thought I'd add that you can find more information for
file-format
> > specific settings for the "field" entry (including level) for MET
on pg.
> 53
> > of the MET User's Guide for met-6.0:
> >
> > http://www.dtcenter.org/met/users/docs/users_guide/MET_
> > Users_Guide_v6.0.pdf
> >
> > or in the README file in the met-6.0/data/config area of your
install.
> >
> > For example:
> > // File-format specific settings for the "field" entry:
> > //
> > // - GRIB1 and GRIB2:
> > // - For custom GRIB tables, see note about MET_GRIB_TABLES.
> > // - The "name" entry specifies a GRIB code number or
abbreviation.
> > // - GRIB1 Product Definition Section:
> > //
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
> > // - GRIB2 Product Definition Section:
> > //
http://www.nco.ncep.noaa.gov/pmb/docs/grib2/grib2_doc.shtml
> > // - The "level" entry specifies a level type and value:
> > // - ANNN for accumulation interval NNN
> > // - ZNNN for vertical level NNN
> > // - ZNNN-NNN for a range of vertical levels
> > // - PNNN for pressure level NNN in hPa
> > // - PNNN-NNN for a range of pressure levels in hPa
> > // - LNNN for a generic level type
> > // - RNNN for a specific GRIB record number
> >
> > I hope that helps!
> >
> > Julie
> >
> > On Thu, Aug 31, 2017 at 12:25 PM, Minna Win via RT
<met_help at ucar.edu>
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> > >
> > > Hi Brad,
> > >
> > > For GRIB2, the level (if not an accumulated precip type) is
typically
> > "Z0"
> > > for surface (sometimes "L0" works, too), but for other levels
such as
> > > entire atmosphere, you will need to refer to the GRIB2 table
"Table 2.
> > > Parameters and Units":
> > > http://www.wmo.int/pages/prog/www/WDM/Guides/Guide-binary-
2.html#GRIB
> .
> > > From that table, find the value that corresponds to your level
type and
> > set
> > > it to GRIB_lvl_typ and then set GRIB_lvl_val1 to the appropriate
level.
> > > For 'entire atmosphere', there is no numerical value for the
level, so
> it
> > > would be set to 0 or 255 (you might have to try both to see what
works,
> > it
> > > depends on the grib file). For something like a pressure
difference,
> you
> > > would need to set GRIB_lvl_typ and the GRIB_lvl_va1 and
GRIB_lvl_val2,
> > > i.e. GRIB_lvl_val2=top level GRIB_lvl_val1=bottom level, which
you
> would
> > > obtain from performing a wgrib2/grep combination on your grib
file and
> > > variable of interest. You can use the MET tool plot_data_plane
to check
> > > that your GRIB_lvl_typ and GRIB_lvl_val1 (and GRIB_lvl_val2 if
used)
> > > produced a match to the record of interest in the grib file.
> > >
> > > As for your question about the (*,*) you see in some examples,
there is
> > an
> > > explanation in Chapter 3 of the MET User's Manual (excerpted):
> > >
> > >
> > > > // - The "level" entry specifies the dimensions to be used:
> > > > // - (i,...,j,*,*) for a single field, where i,...,j
specifies
> fixed
> > > > // dimension values and *,* specifies the two dimensions
for
> the
> > > > // gridded field.
> > > >
> > >
> > >
> > > field = [
> > > {
> > > name = "QVAPOR";
> > > level = "( = "(0,5,*,*)"
> > >
> > > },
> > > {
> > >
> > > name = "TMP_P850_ENS_MEAN;
> > > level = [ "(*,*)" ];
> > > }
> > > ];
> > >
> > >
> > > I hope that clarifies things. Please don't hesitate to contact
us if
> you
> > > have further questions.
> > >
> > > Regards,
> > > Minna
> > >
> > > ---------------
> > > Minna Win
> > > NCAR
> > > Research Applications Lab
> > > Phone: 303-497-8423
> > > Fax: 302-497-8401
> > >
> > >
> > > On Thu, Aug 31, 2017 at 4:21 PM, Brad Diehl - NOAA Federal via
RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829
>
> > > >
> > > > the A06 worked. Prior to trying that I sent you via mailx the
> fcst/obs
> > > > grib2 files, you can ignore that now.
> > > >
> > > > But it does beg the question, how does one find the level
names, is
> it
> > in
> > > > the grib file itself? Very confusing, at least to me. And what
does
> *,*
> > > > mean when I see that in some examples?
> > > >
> > > > Thanks for all your help........Brad
> > > >
> > > > On Thu, Aug 31, 2017 at 4:11 PM, Minna Win via RT
<met_help at ucar.edu
> >
> > > > wrote:
> > > >
> > > > > Hi Brad,
> > > > >
> > > > > Tara Jensen brought this to my attention, for verify APCP,
you need
> > to
> > > > set
> > > > > the level to the value pre-pended with the letter 'A'. For
> instance,
> > > > for 6
> > > > > hour accumulation level, set level="A06".
> > > > >
> > > > > Thanks,
> > > > > Minna
> > > > >
> > > > > ---------------
> > > > > Minna Win
> > > > > NCAR
> > > > > Research Applications Lab
> > > > > Phone: 303-497-8423
> > > > > Fax: 302-497-8401
> > > > >
> > > > >
> > > > > On Thu, Aug 31, 2017 at 3:52 PM, Minna Win-Gildenmeister <
> > > > > minnawin at ucar.edu>
> > > > > wrote:
> > > > >
> > > > > > Hi Brad,
> > > > > >
> > > > > > Can you try "L0", and "P0"? If that doesn't work, could
you
> please
> > > > send
> > > > > > me your grib file? Instructions are as follows:
> > > > > >
> > > > > > Non-Windows FTP client:
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > > * ftp -p ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
> > > username=anonymous
> > > > > > password="your email address" cd incoming/irap/met_help
mkdir
> > > > > > "surname_data" # where surname is your last name cd
> > "surname_data"
> > > > put
> > > > > > "your files" bye *
> > > > > >
> > > > > > Windows FTP client:
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > >
> > > > > > * ftp ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
> > username=anonymous
> > > > > > password="your email address" quote pasv cd
> incoming/irap/met_help
> > > > mkdir
> > > > > > "surname_data" # where surname is your last name cd
> "surname_data"
> > > put
> > > > > > "your files" bye*
> > > > > >
> > > > > > I can take a look at the level type and value, you might
need to
> > use
> > > > > > special syntax in your config file to get grid_stat to
work. Let
> > me
> > > > know
> > > > > > if you ftp'd a grib file so I can grab it before the file
gets
> > > > scrubbed.
> > > > > >
> > > > > > Thanks,
> > > > > > Minna
> > > > > >
> > > > > >
> > > > > > ---------------
> > > > > > Minna Win
> > > > > > NCAR
> > > > > > Research Applications Lab
> > > > > > Phone: 303-497-8423 <(303)%20497-8423>
> > > > > > Fax: 302-497-8401 <(302)%20497-8401>
> > > > > >
> > > > > >
> > > > > > On Thu, Aug 31, 2017 at 3:09 PM, Brad Diehl - NOAA Federal
via
> RT <
> > > > > > met_help at ucar.edu> wrote:
> > > > > >
> > > > > >>
> > > > > >> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> > > > > >>
> > > > > >> After running pcp_combine, everything did work correctly,
it is
> > the
> > > > > grib2
> > > > > >> fcts/obs part I'm having trouble with. I think you put me
on the
> > > right
> > > > > >> track by suggesting wgrib2. Here's what I get using
wgrib2 on
> the
> > > fcst
> > > > > >> file:
> > > > > >>
> > > > > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ wgrib2
> > > > > >> p06m_2017082312f006.grb2
> > > > > >> 1:21:d=2017082312 <(201)%20708-2312>:APCP:surface -
> > > surface:CodeTable
> > > > > >> 4.11=missing:
> > > > > >>
> > > > > >> in the config file I'm using:
> > > > > >>
> > > > > >>
> > > > > >> fcst = {
> > > > > >>
> > > > > >> field = [
> > > > > >> {
> > > > > >> name = "APCP";
> > > > > >> level = "Z0" ;
> > > > > >> cat_thresh = [ >0.0, >=5.0 , >=10.0];
> > > > > >> }
> > > > > >> ];
> > > > > >>
> > > > > >> for level, should I be using surface or Z0?, neither of
which
> > seemed
> > > > to
> > > > > >> work.
> > > > > >>
> > > > > >> On Thu, Aug 31, 2017 at 2:58 PM, Minna Win via RT <
> > > met_help at ucar.edu>
> > > > > >> wrote:
> > > > > >>
> > > > > >> > Hi Brad,
> > > > > >> >
> > > > > >> > Please disregard my previous response. I didn't see
that you
> > had
> > > > run
> > > > > >> > pcp_combine-that *should* have worked. I'll look a
little
> > further
> > > > or
> > > > > >> see
> > > > > >> > if someone else has seen this before.
> > > > > >> >
> > > > > >> >
> > > > > >> >
> > > > > >> > Regards,
> > > > > >> > Minna
> > > > > >> >
> > > > > >> > ---------------
> > > > > >> > Minna Win
> > > > > >> > NCAR
> > > > > >> > Research Applications Lab
> > > > > >> > Phone: 303-497-8423
> > > > > >> > Fax: 302-497-8401
> > > > > >> >
> > > > > >> >
> > > > > >> > On Thu, Aug 31, 2017 at 2:50 PM, Minna Win-
Gildenmeister <
> > > > > >> > minnawin at ucar.edu>
> > > > > >> > wrote:
> > > > > >> >
> > > > > >> > > HI Brad,
> > > > > >> > >
> > > > > >> > > wgrib is an NDFD tool and I don't know how it is
> interpreting
> > > the
> > > > > APCP
> > > > > >> > > variable in your grib file. Try using wgrib2 on your
grib
> file
> > > and
> > > > > >> grep
> > > > > >> > for
> > > > > >> > > APCP and use that variable abbreviation in you MET
Grid_stat
> > > > config
> > > > > >> > file. I
> > > > > >> > > hope that helps. Please let us know if you still
have
> > problems
> > > > > using
> > > > > >> > > Grid_stat.
> > > > > >> > >
> > > > > >> > > Regards,
> > > > > >> > > Minna
> > > > > >> > >
> > > > > >> > > ---------------
> > > > > >> > > Minna Win
> > > > > >> > > NCAR
> > > > > >> > > Research Applications Lab
> > > > > >> > > Phone: 303-497-8423 <(303)%20497-8423>
> > > > > >> > > Fax: 302-497-8401 <(302)%20497-8401>
> > > > > >> > >
> > > > > >> > >
> > > > > >> > > On Thu, Aug 31, 2017 at 2:33 PM, Brad Diehl - NOAA
Federal
> via
> > > RT
> > > > <
> > > > > >> > > met_help at ucar.edu> wrote:
> > > > > >> > >
> > > > > >> > >>
> > > > > >> > >> Thu Aug 31 08:33:00 2017: Request 81829 was acted
upon.
> > > > > >> > >> Transaction: Ticket created by brad.diehl at noaa.gov
> > > > > >> > >> Queue: met_help
> > > > > >> > >> Subject: grib2 files as input into MET
grid_stat
> > > > > >> > >> Owner: Nobody
> > > > > >> > >> Requestors: brad.diehl at noaa.gov
> > > > > >> > >> Status: new
> > > > > >> > >> Ticket <URL: https://rt.rap.ucar.edu/rt/Tic
> > > > > >> ket/Display.html?id=81829 >
> > > > > >> > >>
> > > > > >> > >>
> > > > > >> > >> I tried to use grid_stat with two grib2 files as
input, and
> > > > > received
> > > > > >> an
> > > > > >> > >> error regarding the GRIB2 field abbreviation
'APCP06' (see
> > > output
> > > > > at
> > > > > >> > >> bottom
> > > > > >> > >> of message). Yet when I run degrib on the fcst and
obs
> files
> > > > that's
> > > > > >> > what I
> > > > > >> > >> get:
> > > > > >> > >>
> > > > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
degrib
> > > > > >> > >> p06m_2017082312f006.grb2 -I
> > > > > >> > >> MsgNum, Byte, GRIB-Version, elem, level,
reference(UTC),
> > > > > valid(UTC),
> > > > > >> > >> Proj(hr)
> > > > > >> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation
[kg/(m^2)]",
> > > 0-SFC,
> > > > > >> > >> 08/23/2017
> > > > > >> > >> 12:00, 08/23/2017 18:00, 6.00
> > > > > >> > >>
> > > > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
degrib
> > > > > >> > >> urma2p5.2017082318 <(201)%20708-
2318>.pcp_06h.184.grb2 -I
> > > > > >> > >> MsgNum, Byte, GRIB-Version, elem, level,
reference(UTC),
> > > > > valid(UTC),
> > > > > >> > >> Proj(hr)
> > > > > >> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation
[kg/(m^2)]",
> > > 0-SFC,
> > > > > >> > >> 08/23/2017
> > > > > >> > >> 12:00, 08/23/2017 18:00, 6.00
> > > > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > > > > >> > >>
> > > > > >> > >> If I put both the fcst and obs files through
pcp_combine
> and
> > > use
> > > > > the
> > > > > >> > >> resultant .nc files, grid_stat runs fine, with the
caveat
> > that
> > > > the
> > > > > >> field
> > > > > >> > >> name becomes 'APCP_06'.
> > > > > >> > >>
> > > > > >> > >> Thanks for any comments/suggestions, Brad Diehl
> > > > > >> > >>
> > > > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
grid_stat
> > > > > >> > >> p06m_2017082312f006.grb2 urma2p5.2017082318
> > <(201)%20708-2318>
> > > > > >> > >> .pcp_06h.184.grb2
> > > > > >> > >> config/GridStatConfig_APCP_06_grib -outdir
out/grid_stat
> -v
> > > > 4DEBUG
> > > > > >> 1:
> > > > > >> > >> Default Config File:
> > > > > >> > >>
/usrx/local/dev/MET/6.0/share/met/config/GridStatConfig_
> > > default
> > > > > >> > >> DEBUG 1: User Config File:
config/GridStatConfig_APCP_06_
> > grib
> > > > > >> > >> DEBUG 4:
Met2dDataFileFactory::new_met_2d_data_file() ->
> > > created
> > > > > new
> > > > > >> > >> Met2dDataFile object of type "FileType_Gb2".
> > > > > >> > >> WARNING:
> > > > > >> > >> WARNING: MetGrib2DataFile::read_grib2_grid() -> MET
does
> not
> > > > > >> currently
> > > > > >> > >> support Lambert Conformal grids where dx (2.5397) !=
dy
> > > (2.5397)
> > > > > and
> > > > > >> may
> > > > > >> > >> produce unexpected results!
> > > > > >> > >> WARNING:
> > > > > >> > >> DEBUG 4:
> > > > > >> > >> DEBUG 4: Lambert Conformal Grid Data:
> > > > > >> > >> DEBUG 4: scale_lat_1: 25
> > > > > >> > >> DEBUG 4: scale_lat_2: 25
> > > > > >> > >> DEBUG 4: lat_pin: 20.192
> > > > > >> > >> DEBUG 4: lon_pin: 121.554
> > > > > >> > >> DEBUG 4: x_pin: 0
> > > > > >> > >> DEBUG 4: y_pin: 0
> > > > > >> > >> DEBUG 4: lon_orient: 95
> > > > > >> > >> DEBUG 4: d_km: 2.5397
> > > > > >> > >> DEBUG 4: r_km: 6371.2
> > > > > >> > >> DEBUG 4: nx: 2145
> > > > > >> > >> DEBUG 4: ny: 1597
> > > > > >> > >> DEBUG 4:
> > > > > >> > >> DEBUG 4:
Met2dDataFileFactory::new_met_2d_data_file() ->
> > > created
> > > > > new
> > > > > >> > >> Met2dDataFile object of type "FileType_Gb2".
> > > > > >> > >> DEBUG 4:
> > > > > >> > >> DEBUG 4: Lambert Conformal Grid Data:
> > > > > >> > >> DEBUG 4: scale_lat_1: 25
> > > > > >> > >> DEBUG 4: scale_lat_2: 25
> > > > > >> > >> DEBUG 4: lat_pin: 20.192
> > > > > >> > >> DEBUG 4: lon_pin: 121.554
> > > > > >> > >> DEBUG 4: x_pin: 0
> > > > > >> > >> DEBUG 4: y_pin: 0
> > > > > >> > >> DEBUG 4: lon_orient: 95
> > > > > >> > >> DEBUG 4: d_km: 2.5397
> > > > > >> > >> DEBUG 4: r_km: 6371.2
> > > > > >> > >> DEBUG 4: nx: 2145
> > > > > >> > >> DEBUG 4: ny: 1377
> > > > > >> > >> DEBUG 4:
> > > > > >> > >> DEBUG 4: VarInfoFactory::new_var_info() -> created
new
> > VarInfo
> > > > > >> object of
> > > > > >> > >> type "FileType_Gb2".
> > > > > >> > >> DEBUG 4: VarInfoFactory::new_var_info() -> created
new
> > VarInfo
> > > > > >> object of
> > > > > >> > >> type "FileType_Gb2".
> > > > > >> > >> ERROR :
> > > > > >> > >> ERROR : VarInfoGrib2::set_dict() -> unrecognized
GRIB2
> field
> > > > > >> > abbreviation
> > > > > >> > >> 'APCP06'
> > > > > >> > >> ERROR :
> > > > > >> > >>
> > > > > >> > >>
> > > > > >> > >
> > > > > >> >
> > > > > >> >
> > > > > >>
> > > > > >>
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: grib2 files as input into MET grid_stat
From: Brad Diehl - NOAA Federal
Time: Fri Sep 01 10:29:21 2017
Is there any kind of time frame to make that switch over to
compliance?
Brad Diehl
On Fri, Sep 1, 2017 at 3:35 PM, John Halley Gotway via RT
<met_help at ucar.edu
> wrote:
> Hi Brad,
>
> Unfortunately, no, the NetCDF output generated by the MET tools
themselves
> are not CF-compliant.
>
> Since we always suggest that others use that convention, we should
too!
> And we do plan to make that switch over eventually, but haven't done
so
> yet.
>
> Thanks,
> John
>
> On Thu, Aug 31, 2017 at 1:15 PM, Brad Diehl - NOAA Federal via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> >
> > Thanks for your help and suggestions today. A question came up
from one
> of
> > our developers regarding netcdf and CF compliance. Is the output
of
> MET/6.0
> > netcdf CF compliant?
> >
> > Thanks in advance.......Brad
> >
> > On Thu, Aug 31, 2017 at 6:35 PM, Julie Prestopnik via RT <
> > met_help at ucar.edu>
> > wrote:
> >
> > > Hi Brad,
> > >
> > > I just thought I'd add that you can find more information for
> file-format
> > > specific settings for the "field" entry (including level) for
MET on
> pg.
> > 53
> > > of the MET User's Guide for met-6.0:
> > >
> > > http://www.dtcenter.org/met/users/docs/users_guide/MET_
> > > Users_Guide_v6.0.pdf
> > >
> > > or in the README file in the met-6.0/data/config area of your
install.
> > >
> > > For example:
> > > // File-format specific settings for the "field" entry:
> > > //
> > > // - GRIB1 and GRIB2:
> > > // - For custom GRIB tables, see note about
MET_GRIB_TABLES.
> > > // - The "name" entry specifies a GRIB code number or
> abbreviation.
> > > // - GRIB1 Product Definition Section:
> > > //
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
> > > // - GRIB2 Product Definition Section:
> > > // http://www.nco.ncep.noaa.gov/
> pmb/docs/grib2/grib2_doc.shtml
> > > // - The "level" entry specifies a level type and value:
> > > // - ANNN for accumulation interval NNN
> > > // - ZNNN for vertical level NNN
> > > // - ZNNN-NNN for a range of vertical levels
> > > // - PNNN for pressure level NNN in hPa
> > > // - PNNN-NNN for a range of pressure levels in hPa
> > > // - LNNN for a generic level type
> > > // - RNNN for a specific GRIB record number
> > >
> > > I hope that helps!
> > >
> > > Julie
> > >
> > > On Thu, Aug 31, 2017 at 12:25 PM, Minna Win via RT
<met_help at ucar.edu>
> > > wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829
>
> > > >
> > > > Hi Brad,
> > > >
> > > > For GRIB2, the level (if not an accumulated precip type) is
typically
> > > "Z0"
> > > > for surface (sometimes "L0" works, too), but for other levels
such as
> > > > entire atmosphere, you will need to refer to the GRIB2 table
"Table
> 2.
> > > > Parameters and Units":
> > > > http://www.wmo.int/pages/prog/www/WDM/Guides/Guide-binary-2.
> html#GRIB
> > .
> > > > From that table, find the value that corresponds to your level
type
> and
> > > set
> > > > it to GRIB_lvl_typ and then set GRIB_lvl_val1 to the
appropriate
> level.
> > > > For 'entire atmosphere', there is no numerical value for the
level,
> so
> > it
> > > > would be set to 0 or 255 (you might have to try both to see
what
> works,
> > > it
> > > > depends on the grib file). For something like a pressure
difference,
> > you
> > > > would need to set GRIB_lvl_typ and the GRIB_lvl_va1 and
> GRIB_lvl_val2,
> > > > i.e. GRIB_lvl_val2=top level GRIB_lvl_val1=bottom level, which
you
> > would
> > > > obtain from performing a wgrib2/grep combination on your grib
file
> and
> > > > variable of interest. You can use the MET tool plot_data_plane
to
> check
> > > > that your GRIB_lvl_typ and GRIB_lvl_val1 (and GRIB_lvl_val2 if
used)
> > > > produced a match to the record of interest in the grib file.
> > > >
> > > > As for your question about the (*,*) you see in some examples,
there
> is
> > > an
> > > > explanation in Chapter 3 of the MET User's Manual (excerpted):
> > > >
> > > >
> > > > > // - The "level" entry specifies the dimensions to be used:
> > > > > // - (i,...,j,*,*) for a single field, where i,...,j
specifies
> > fixed
> > > > > // dimension values and *,* specifies the two
dimensions for
> > the
> > > > > // gridded field.
> > > > >
> > > >
> > > >
> > > > field = [
> > > > {
> > > > name = "QVAPOR";
> > > > level = "( = "(0,5,*,*)"
> > > >
> > > > },
> > > > {
> > > >
> > > > name = "TMP_P850_ENS_MEAN;
> > > > level = [ "(*,*)" ];
> > > > }
> > > > ];
> > > >
> > > >
> > > > I hope that clarifies things. Please don't hesitate to
contact us if
> > you
> > > > have further questions.
> > > >
> > > > Regards,
> > > > Minna
> > > >
> > > > ---------------
> > > > Minna Win
> > > > NCAR
> > > > Research Applications Lab
> > > > Phone: 303-497-8423
> > > > Fax: 302-497-8401
> > > >
> > > >
> > > > On Thu, Aug 31, 2017 at 4:21 PM, Brad Diehl - NOAA Federal via
RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> > > > >
> > > > > the A06 worked. Prior to trying that I sent you via mailx
the
> > fcst/obs
> > > > > grib2 files, you can ignore that now.
> > > > >
> > > > > But it does beg the question, how does one find the level
names, is
> > it
> > > in
> > > > > the grib file itself? Very confusing, at least to me. And
what does
> > *,*
> > > > > mean when I see that in some examples?
> > > > >
> > > > > Thanks for all your help........Brad
> > > > >
> > > > > On Thu, Aug 31, 2017 at 4:11 PM, Minna Win via RT <
> met_help at ucar.edu
> > >
> > > > > wrote:
> > > > >
> > > > > > Hi Brad,
> > > > > >
> > > > > > Tara Jensen brought this to my attention, for verify APCP,
you
> need
> > > to
> > > > > set
> > > > > > the level to the value pre-pended with the letter 'A'.
For
> > instance,
> > > > > for 6
> > > > > > hour accumulation level, set level="A06".
> > > > > >
> > > > > > Thanks,
> > > > > > Minna
> > > > > >
> > > > > > ---------------
> > > > > > Minna Win
> > > > > > NCAR
> > > > > > Research Applications Lab
> > > > > > Phone: 303-497-8423
> > > > > > Fax: 302-497-8401
> > > > > >
> > > > > >
> > > > > > On Thu, Aug 31, 2017 at 3:52 PM, Minna Win-Gildenmeister <
> > > > > > minnawin at ucar.edu>
> > > > > > wrote:
> > > > > >
> > > > > > > Hi Brad,
> > > > > > >
> > > > > > > Can you try "L0", and "P0"? If that doesn't work, could
you
> > please
> > > > > send
> > > > > > > me your grib file? Instructions are as follows:
> > > > > > >
> > > > > > > Non-Windows FTP client:
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > * ftp -p ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
> > > > username=anonymous
> > > > > > > password="your email address" cd
incoming/irap/met_help
> mkdir
> > > > > > > "surname_data" # where surname is your last name cd
> > > "surname_data"
> > > > > put
> > > > > > > "your files" bye *
> > > > > > >
> > > > > > > Windows FTP client:
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > * ftp ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
> > > username=anonymous
> > > > > > > password="your email address" quote pasv cd
> > incoming/irap/met_help
> > > > > mkdir
> > > > > > > "surname_data" # where surname is your last name cd
> > "surname_data"
> > > > put
> > > > > > > "your files" bye*
> > > > > > >
> > > > > > > I can take a look at the level type and value, you might
need
> to
> > > use
> > > > > > > special syntax in your config file to get grid_stat to
work.
> Let
> > > me
> > > > > know
> > > > > > > if you ftp'd a grib file so I can grab it before the
file gets
> > > > > scrubbed.
> > > > > > >
> > > > > > > Thanks,
> > > > > > > Minna
> > > > > > >
> > > > > > >
> > > > > > > ---------------
> > > > > > > Minna Win
> > > > > > > NCAR
> > > > > > > Research Applications Lab
> > > > > > > Phone: 303-497-8423 <(303)%20497-8423>
> > > > > > > Fax: 302-497-8401 <(302)%20497-8401>
> > > > > > >
> > > > > > >
> > > > > > > On Thu, Aug 31, 2017 at 3:09 PM, Brad Diehl - NOAA
Federal via
> > RT <
> > > > > > > met_help at ucar.edu> wrote:
> > > > > > >
> > > > > > >>
> > > > > > >> <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829
> >
> > > > > > >>
> > > > > > >> After running pcp_combine, everything did work
correctly, it
> is
> > > the
> > > > > > grib2
> > > > > > >> fcts/obs part I'm having trouble with. I think you put
me on
> the
> > > > right
> > > > > > >> track by suggesting wgrib2. Here's what I get using
wgrib2 on
> > the
> > > > fcst
> > > > > > >> file:
> > > > > > >>
> > > > > > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
wgrib2
> > > > > > >> p06m_2017082312f006.grb2
> > > > > > >> 1:21:d=2017082312 <(201)%20708-2312>:APCP:surface -
> > > > surface:CodeTable
> > > > > > >> 4.11=missing:
> > > > > > >>
> > > > > > >> in the config file I'm using:
> > > > > > >>
> > > > > > >>
> > > > > > >> fcst = {
> > > > > > >>
> > > > > > >> field = [
> > > > > > >> {
> > > > > > >> name = "APCP";
> > > > > > >> level = "Z0" ;
> > > > > > >> cat_thresh = [ >0.0, >=5.0 , >=10.0];
> > > > > > >> }
> > > > > > >> ];
> > > > > > >>
> > > > > > >> for level, should I be using surface or Z0?, neither of
which
> > > seemed
> > > > > to
> > > > > > >> work.
> > > > > > >>
> > > > > > >> On Thu, Aug 31, 2017 at 2:58 PM, Minna Win via RT <
> > > > met_help at ucar.edu>
> > > > > > >> wrote:
> > > > > > >>
> > > > > > >> > Hi Brad,
> > > > > > >> >
> > > > > > >> > Please disregard my previous response. I didn't see
that
> you
> > > had
> > > > > run
> > > > > > >> > pcp_combine-that *should* have worked. I'll look a
little
> > > further
> > > > > or
> > > > > > >> see
> > > > > > >> > if someone else has seen this before.
> > > > > > >> >
> > > > > > >> >
> > > > > > >> >
> > > > > > >> > Regards,
> > > > > > >> > Minna
> > > > > > >> >
> > > > > > >> > ---------------
> > > > > > >> > Minna Win
> > > > > > >> > NCAR
> > > > > > >> > Research Applications Lab
> > > > > > >> > Phone: 303-497-8423
> > > > > > >> > Fax: 302-497-8401
> > > > > > >> >
> > > > > > >> >
> > > > > > >> > On Thu, Aug 31, 2017 at 2:50 PM, Minna Win-
Gildenmeister <
> > > > > > >> > minnawin at ucar.edu>
> > > > > > >> > wrote:
> > > > > > >> >
> > > > > > >> > > HI Brad,
> > > > > > >> > >
> > > > > > >> > > wgrib is an NDFD tool and I don't know how it is
> > interpreting
> > > > the
> > > > > > APCP
> > > > > > >> > > variable in your grib file. Try using wgrib2 on
your grib
> > file
> > > > and
> > > > > > >> grep
> > > > > > >> > for
> > > > > > >> > > APCP and use that variable abbreviation in you MET
> Grid_stat
> > > > > config
> > > > > > >> > file. I
> > > > > > >> > > hope that helps. Please let us know if you still
have
> > > problems
> > > > > > using
> > > > > > >> > > Grid_stat.
> > > > > > >> > >
> > > > > > >> > > Regards,
> > > > > > >> > > Minna
> > > > > > >> > >
> > > > > > >> > > ---------------
> > > > > > >> > > Minna Win
> > > > > > >> > > NCAR
> > > > > > >> > > Research Applications Lab
> > > > > > >> > > Phone: 303-497-8423 <(303)%20497-8423>
> > > > > > >> > > Fax: 302-497-8401 <(302)%20497-8401>
> > > > > > >> > >
> > > > > > >> > >
> > > > > > >> > > On Thu, Aug 31, 2017 at 2:33 PM, Brad Diehl - NOAA
Federal
> > via
> > > > RT
> > > > > <
> > > > > > >> > > met_help at ucar.edu> wrote:
> > > > > > >> > >
> > > > > > >> > >>
> > > > > > >> > >> Thu Aug 31 08:33:00 2017: Request 81829 was acted
upon.
> > > > > > >> > >> Transaction: Ticket created by brad.diehl at noaa.gov
> > > > > > >> > >> Queue: met_help
> > > > > > >> > >> Subject: grib2 files as input into MET
grid_stat
> > > > > > >> > >> Owner: Nobody
> > > > > > >> > >> Requestors: brad.diehl at noaa.gov
> > > > > > >> > >> Status: new
> > > > > > >> > >> Ticket <URL: https://rt.rap.ucar.edu/rt/Tic
> > > > > > >> ket/Display.html?id=81829 >
> > > > > > >> > >>
> > > > > > >> > >>
> > > > > > >> > >> I tried to use grid_stat with two grib2 files as
input,
> and
> > > > > > received
> > > > > > >> an
> > > > > > >> > >> error regarding the GRIB2 field abbreviation
'APCP06'
> (see
> > > > output
> > > > > > at
> > > > > > >> > >> bottom
> > > > > > >> > >> of message). Yet when I run degrib on the fcst and
obs
> > files
> > > > > that's
> > > > > > >> > what I
> > > > > > >> > >> get:
> > > > > > >> > >>
> > > > > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
degrib
> > > > > > >> > >> p06m_2017082312f006.grb2 -I
> > > > > > >> > >> MsgNum, Byte, GRIB-Version, elem, level,
reference(UTC),
> > > > > > valid(UTC),
> > > > > > >> > >> Proj(hr)
> > > > > > >> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation
[kg/(m^2)]",
> > > > 0-SFC,
> > > > > > >> > >> 08/23/2017
> > > > > > >> > >> 12:00, 08/23/2017 18:00, 6.00
> > > > > > >> > >>
> > > > > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
degrib
> > > > > > >> > >> urma2p5.2017082318 <(201)%20708-
2318>.pcp_06h.184.grb2
> -I
> > > > > > >> > >> MsgNum, Byte, GRIB-Version, elem, level,
reference(UTC),
> > > > > > valid(UTC),
> > > > > > >> > >> Proj(hr)
> > > > > > >> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation
[kg/(m^2)]",
> > > > 0-SFC,
> > > > > > >> > >> 08/23/2017
> > > > > > >> > >> 12:00, 08/23/2017 18:00, 6.00
> > > > > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > > > > > >> > >>
> > > > > > >> > >> If I put both the fcst and obs files through
pcp_combine
> > and
> > > > use
> > > > > > the
> > > > > > >> > >> resultant .nc files, grid_stat runs fine, with the
caveat
> > > that
> > > > > the
> > > > > > >> field
> > > > > > >> > >> name becomes 'APCP_06'.
> > > > > > >> > >>
> > > > > > >> > >> Thanks for any comments/suggestions, Brad Diehl
> > > > > > >> > >>
> > > > > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> grid_stat
> > > > > > >> > >> p06m_2017082312f006.grb2 urma2p5.2017082318
> > > <(201)%20708-2318>
> > > > > > >> > >> .pcp_06h.184.grb2
> > > > > > >> > >> config/GridStatConfig_APCP_06_grib -outdir
out/grid_stat
> > -v
> > > > > 4DEBUG
> > > > > > >> 1:
> > > > > > >> > >> Default Config File:
> > > > > > >> > >>
/usrx/local/dev/MET/6.0/share/met/config/GridStatConfig_
> > > > default
> > > > > > >> > >> DEBUG 1: User Config File:
config/GridStatConfig_APCP_06_
> > > grib
> > > > > > >> > >> DEBUG 4:
Met2dDataFileFactory::new_met_2d_data_file() ->
> > > > created
> > > > > > new
> > > > > > >> > >> Met2dDataFile object of type "FileType_Gb2".
> > > > > > >> > >> WARNING:
> > > > > > >> > >> WARNING: MetGrib2DataFile::read_grib2_grid() ->
MET does
> > not
> > > > > > >> currently
> > > > > > >> > >> support Lambert Conformal grids where dx (2.5397)
!= dy
> > > > (2.5397)
> > > > > > and
> > > > > > >> may
> > > > > > >> > >> produce unexpected results!
> > > > > > >> > >> WARNING:
> > > > > > >> > >> DEBUG 4:
> > > > > > >> > >> DEBUG 4: Lambert Conformal Grid Data:
> > > > > > >> > >> DEBUG 4: scale_lat_1: 25
> > > > > > >> > >> DEBUG 4: scale_lat_2: 25
> > > > > > >> > >> DEBUG 4: lat_pin: 20.192
> > > > > > >> > >> DEBUG 4: lon_pin: 121.554
> > > > > > >> > >> DEBUG 4: x_pin: 0
> > > > > > >> > >> DEBUG 4: y_pin: 0
> > > > > > >> > >> DEBUG 4: lon_orient: 95
> > > > > > >> > >> DEBUG 4: d_km: 2.5397
> > > > > > >> > >> DEBUG 4: r_km: 6371.2
> > > > > > >> > >> DEBUG 4: nx: 2145
> > > > > > >> > >> DEBUG 4: ny: 1597
> > > > > > >> > >> DEBUG 4:
> > > > > > >> > >> DEBUG 4:
Met2dDataFileFactory::new_met_2d_data_file() ->
> > > > created
> > > > > > new
> > > > > > >> > >> Met2dDataFile object of type "FileType_Gb2".
> > > > > > >> > >> DEBUG 4:
> > > > > > >> > >> DEBUG 4: Lambert Conformal Grid Data:
> > > > > > >> > >> DEBUG 4: scale_lat_1: 25
> > > > > > >> > >> DEBUG 4: scale_lat_2: 25
> > > > > > >> > >> DEBUG 4: lat_pin: 20.192
> > > > > > >> > >> DEBUG 4: lon_pin: 121.554
> > > > > > >> > >> DEBUG 4: x_pin: 0
> > > > > > >> > >> DEBUG 4: y_pin: 0
> > > > > > >> > >> DEBUG 4: lon_orient: 95
> > > > > > >> > >> DEBUG 4: d_km: 2.5397
> > > > > > >> > >> DEBUG 4: r_km: 6371.2
> > > > > > >> > >> DEBUG 4: nx: 2145
> > > > > > >> > >> DEBUG 4: ny: 1377
> > > > > > >> > >> DEBUG 4:
> > > > > > >> > >> DEBUG 4: VarInfoFactory::new_var_info() -> created
new
> > > VarInfo
> > > > > > >> object of
> > > > > > >> > >> type "FileType_Gb2".
> > > > > > >> > >> DEBUG 4: VarInfoFactory::new_var_info() -> created
new
> > > VarInfo
> > > > > > >> object of
> > > > > > >> > >> type "FileType_Gb2".
> > > > > > >> > >> ERROR :
> > > > > > >> > >> ERROR : VarInfoGrib2::set_dict() -> unrecognized
GRIB2
> > field
> > > > > > >> > abbreviation
> > > > > > >> > >> 'APCP06'
> > > > > > >> > >> ERROR :
> > > > > > >> > >>
> > > > > > >> > >>
> > > > > > >> > >
> > > > > > >> >
> > > > > > >> >
> > > > > > >>
> > > > > > >>
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: grib2 files as input into MET grid_stat
From: John Halley Gotway
Time: Fri Sep 01 10:47:38 2017
Brad,
It won't be in the next version, 6.1, and we don't currently have it
assigned to version 6.2 either. So no specific timeline.
How does this affect the decisions you're making? In our application
of
MET to projects in the Development Testbed Center (DTC), the primary
output
of MET are the ASCII statistics files, to which NetCDF CF conventions
do
not apply. We load them into a METViewer database, aggregate through
time,
and make plots to illustrate model performance.
Which NetCDF files produced by the MET tools are of the most concern?
Thanks,
John
On Fri, Sep 1, 2017 at 10:29 AM, Brad Diehl - NOAA Federal via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
>
> Is there any kind of time frame to make that switch over to
compliance?
>
> Brad Diehl
>
> On Fri, Sep 1, 2017 at 3:35 PM, John Halley Gotway via RT <
> met_help at ucar.edu
> > wrote:
>
> > Hi Brad,
> >
> > Unfortunately, no, the NetCDF output generated by the MET tools
> themselves
> > are not CF-compliant.
> >
> > Since we always suggest that others use that convention, we should
too!
> > And we do plan to make that switch over eventually, but haven't
done so
> > yet.
> >
> > Thanks,
> > John
> >
> > On Thu, Aug 31, 2017 at 1:15 PM, Brad Diehl - NOAA Federal via RT
<
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> > >
> > > Thanks for your help and suggestions today. A question came up
from one
> > of
> > > our developers regarding netcdf and CF compliance. Is the output
of
> > MET/6.0
> > > netcdf CF compliant?
> > >
> > > Thanks in advance.......Brad
> > >
> > > On Thu, Aug 31, 2017 at 6:35 PM, Julie Prestopnik via RT <
> > > met_help at ucar.edu>
> > > wrote:
> > >
> > > > Hi Brad,
> > > >
> > > > I just thought I'd add that you can find more information for
> > file-format
> > > > specific settings for the "field" entry (including level) for
MET on
> > pg.
> > > 53
> > > > of the MET User's Guide for met-6.0:
> > > >
> > > > http://www.dtcenter.org/met/users/docs/users_guide/MET_
> > > > Users_Guide_v6.0.pdf
> > > >
> > > > or in the README file in the met-6.0/data/config area of your
> install.
> > > >
> > > > For example:
> > > > // File-format specific settings for the "field" entry:
> > > > //
> > > > // - GRIB1 and GRIB2:
> > > > // - For custom GRIB tables, see note about
MET_GRIB_TABLES.
> > > > // - The "name" entry specifies a GRIB code number or
> > abbreviation.
> > > > // - GRIB1 Product Definition Section:
> > > > //
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table2.html
> > > > // - GRIB2 Product Definition Section:
> > > > // http://www.nco.ncep.noaa.gov/
> > pmb/docs/grib2/grib2_doc.shtml
> > > > // - The "level" entry specifies a level type and value:
> > > > // - ANNN for accumulation interval NNN
> > > > // - ZNNN for vertical level NNN
> > > > // - ZNNN-NNN for a range of vertical levels
> > > > // - PNNN for pressure level NNN in hPa
> > > > // - PNNN-NNN for a range of pressure levels in hPa
> > > > // - LNNN for a generic level type
> > > > // - RNNN for a specific GRIB record number
> > > >
> > > > I hope that helps!
> > > >
> > > > Julie
> > > >
> > > > On Thu, Aug 31, 2017 at 12:25 PM, Minna Win via RT <
> met_help at ucar.edu>
> > > > wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> > > > >
> > > > > Hi Brad,
> > > > >
> > > > > For GRIB2, the level (if not an accumulated precip type) is
> typically
> > > > "Z0"
> > > > > for surface (sometimes "L0" works, too), but for other
levels such
> as
> > > > > entire atmosphere, you will need to refer to the GRIB2 table
"Table
> > 2.
> > > > > Parameters and Units":
> > > > > http://www.wmo.int/pages/prog/www/WDM/Guides/Guide-binary-2.
> > html#GRIB
> > > .
> > > > > From that table, find the value that corresponds to your
level type
> > and
> > > > set
> > > > > it to GRIB_lvl_typ and then set GRIB_lvl_val1 to the
appropriate
> > level.
> > > > > For 'entire atmosphere', there is no numerical value for the
level,
> > so
> > > it
> > > > > would be set to 0 or 255 (you might have to try both to see
what
> > works,
> > > > it
> > > > > depends on the grib file). For something like a pressure
> difference,
> > > you
> > > > > would need to set GRIB_lvl_typ and the GRIB_lvl_va1 and
> > GRIB_lvl_val2,
> > > > > i.e. GRIB_lvl_val2=top level GRIB_lvl_val1=bottom level,
which you
> > > would
> > > > > obtain from performing a wgrib2/grep combination on your
grib file
> > and
> > > > > variable of interest. You can use the MET tool
plot_data_plane to
> > check
> > > > > that your GRIB_lvl_typ and GRIB_lvl_val1 (and GRIB_lvl_val2
if
> used)
> > > > > produced a match to the record of interest in the grib file.
> > > > >
> > > > > As for your question about the (*,*) you see in some
examples,
> there
> > is
> > > > an
> > > > > explanation in Chapter 3 of the MET User's Manual
(excerpted):
> > > > >
> > > > >
> > > > > > // - The "level" entry specifies the dimensions to be
used:
> > > > > > // - (i,...,j,*,*) for a single field, where i,...,j
specifies
> > > fixed
> > > > > > // dimension values and *,* specifies the two
dimensions
> for
> > > the
> > > > > > // gridded field.
> > > > > >
> > > > >
> > > > >
> > > > > field = [
> > > > > {
> > > > > name = "QVAPOR";
> > > > > level = "( = "(0,5,*,*)"
> > > > >
> > > > > },
> > > > > {
> > > > >
> > > > > name = "TMP_P850_ENS_MEAN;
> > > > > level = [ "(*,*)" ];
> > > > > }
> > > > > ];
> > > > >
> > > > >
> > > > > I hope that clarifies things. Please don't hesitate to
contact us
> if
> > > you
> > > > > have further questions.
> > > > >
> > > > > Regards,
> > > > > Minna
> > > > >
> > > > > ---------------
> > > > > Minna Win
> > > > > NCAR
> > > > > Research Applications Lab
> > > > > Phone: 303-497-8423
> > > > > Fax: 302-497-8401
> > > > >
> > > > >
> > > > > On Thu, Aug 31, 2017 at 4:21 PM, Brad Diehl - NOAA Federal
via RT <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > > >
> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> > > > > >
> > > > > > the A06 worked. Prior to trying that I sent you via mailx
the
> > > fcst/obs
> > > > > > grib2 files, you can ignore that now.
> > > > > >
> > > > > > But it does beg the question, how does one find the level
names,
> is
> > > it
> > > > in
> > > > > > the grib file itself? Very confusing, at least to me. And
what
> does
> > > *,*
> > > > > > mean when I see that in some examples?
> > > > > >
> > > > > > Thanks for all your help........Brad
> > > > > >
> > > > > > On Thu, Aug 31, 2017 at 4:11 PM, Minna Win via RT <
> > met_help at ucar.edu
> > > >
> > > > > > wrote:
> > > > > >
> > > > > > > Hi Brad,
> > > > > > >
> > > > > > > Tara Jensen brought this to my attention, for verify
APCP, you
> > need
> > > > to
> > > > > > set
> > > > > > > the level to the value pre-pended with the letter 'A'.
For
> > > instance,
> > > > > > for 6
> > > > > > > hour accumulation level, set level="A06".
> > > > > > >
> > > > > > > Thanks,
> > > > > > > Minna
> > > > > > >
> > > > > > > ---------------
> > > > > > > Minna Win
> > > > > > > NCAR
> > > > > > > Research Applications Lab
> > > > > > > Phone: 303-497-8423
> > > > > > > Fax: 302-497-8401
> > > > > > >
> > > > > > >
> > > > > > > On Thu, Aug 31, 2017 at 3:52 PM, Minna Win-Gildenmeister
<
> > > > > > > minnawin at ucar.edu>
> > > > > > > wrote:
> > > > > > >
> > > > > > > > Hi Brad,
> > > > > > > >
> > > > > > > > Can you try "L0", and "P0"? If that doesn't work,
could you
> > > please
> > > > > > send
> > > > > > > > me your grib file? Instructions are as follows:
> > > > > > > >
> > > > > > > > Non-Windows FTP client:
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > * ftp -p ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
> > > > > username=anonymous
> > > > > > > > password="your email address" cd
incoming/irap/met_help
> > mkdir
> > > > > > > > "surname_data" # where surname is your last name cd
> > > > "surname_data"
> > > > > > put
> > > > > > > > "your files" bye *
> > > > > > > >
> > > > > > > > Windows FTP client:
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > * ftp ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
> > > > username=anonymous
> > > > > > > > password="your email address" quote pasv cd
> > > incoming/irap/met_help
> > > > > > mkdir
> > > > > > > > "surname_data" # where surname is your last name cd
> > > "surname_data"
> > > > > put
> > > > > > > > "your files" bye*
> > > > > > > >
> > > > > > > > I can take a look at the level type and value, you
might need
> > to
> > > > use
> > > > > > > > special syntax in your config file to get grid_stat to
work.
> > Let
> > > > me
> > > > > > know
> > > > > > > > if you ftp'd a grib file so I can grab it before the
file
> gets
> > > > > > scrubbed.
> > > > > > > >
> > > > > > > > Thanks,
> > > > > > > > Minna
> > > > > > > >
> > > > > > > >
> > > > > > > > ---------------
> > > > > > > > Minna Win
> > > > > > > > NCAR
> > > > > > > > Research Applications Lab
> > > > > > > > Phone: 303-497-8423 <(303)%20497-8423>
> > > > > > > > Fax: 302-497-8401 <(302)%20497-8401>
> > > > > > > >
> > > > > > > >
> > > > > > > > On Thu, Aug 31, 2017 at 3:09 PM, Brad Diehl - NOAA
Federal
> via
> > > RT <
> > > > > > > > met_help at ucar.edu> wrote:
> > > > > > > >
> > > > > > > >>
> > > > > > > >> <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=81829
> > >
> > > > > > > >>
> > > > > > > >> After running pcp_combine, everything did work
correctly, it
> > is
> > > > the
> > > > > > > grib2
> > > > > > > >> fcts/obs part I'm having trouble with. I think you
put me on
> > the
> > > > > right
> > > > > > > >> track by suggesting wgrib2. Here's what I get using
wgrib2
> on
> > > the
> > > > > fcst
> > > > > > > >> file:
> > > > > > > >>
> > > > > > > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
wgrib2
> > > > > > > >> p06m_2017082312f006.grb2
> > > > > > > >> 1:21:d=2017082312 <(201)%20708-2312>:APCP:surface -
> > > > > surface:CodeTable
> > > > > > > >> 4.11=missing:
> > > > > > > >>
> > > > > > > >> in the config file I'm using:
> > > > > > > >>
> > > > > > > >>
> > > > > > > >> fcst = {
> > > > > > > >>
> > > > > > > >> field = [
> > > > > > > >> {
> > > > > > > >> name = "APCP";
> > > > > > > >> level = "Z0" ;
> > > > > > > >> cat_thresh = [ >0.0, >=5.0 , >=10.0];
> > > > > > > >> }
> > > > > > > >> ];
> > > > > > > >>
> > > > > > > >> for level, should I be using surface or Z0?, neither
of
> which
> > > > seemed
> > > > > > to
> > > > > > > >> work.
> > > > > > > >>
> > > > > > > >> On Thu, Aug 31, 2017 at 2:58 PM, Minna Win via RT <
> > > > > met_help at ucar.edu>
> > > > > > > >> wrote:
> > > > > > > >>
> > > > > > > >> > Hi Brad,
> > > > > > > >> >
> > > > > > > >> > Please disregard my previous response. I didn't
see that
> > you
> > > > had
> > > > > > run
> > > > > > > >> > pcp_combine-that *should* have worked. I'll look a
little
> > > > further
> > > > > > or
> > > > > > > >> see
> > > > > > > >> > if someone else has seen this before.
> > > > > > > >> >
> > > > > > > >> >
> > > > > > > >> >
> > > > > > > >> > Regards,
> > > > > > > >> > Minna
> > > > > > > >> >
> > > > > > > >> > ---------------
> > > > > > > >> > Minna Win
> > > > > > > >> > NCAR
> > > > > > > >> > Research Applications Lab
> > > > > > > >> > Phone: 303-497-8423
> > > > > > > >> > Fax: 302-497-8401
> > > > > > > >> >
> > > > > > > >> >
> > > > > > > >> > On Thu, Aug 31, 2017 at 2:50 PM, Minna Win-
Gildenmeister <
> > > > > > > >> > minnawin at ucar.edu>
> > > > > > > >> > wrote:
> > > > > > > >> >
> > > > > > > >> > > HI Brad,
> > > > > > > >> > >
> > > > > > > >> > > wgrib is an NDFD tool and I don't know how it is
> > > interpreting
> > > > > the
> > > > > > > APCP
> > > > > > > >> > > variable in your grib file. Try using wgrib2 on
your
> grib
> > > file
> > > > > and
> > > > > > > >> grep
> > > > > > > >> > for
> > > > > > > >> > > APCP and use that variable abbreviation in you
MET
> > Grid_stat
> > > > > > config
> > > > > > > >> > file. I
> > > > > > > >> > > hope that helps. Please let us know if you still
have
> > > > problems
> > > > > > > using
> > > > > > > >> > > Grid_stat.
> > > > > > > >> > >
> > > > > > > >> > > Regards,
> > > > > > > >> > > Minna
> > > > > > > >> > >
> > > > > > > >> > > ---------------
> > > > > > > >> > > Minna Win
> > > > > > > >> > > NCAR
> > > > > > > >> > > Research Applications Lab
> > > > > > > >> > > Phone: 303-497-8423 <(303)%20497-8423>
> > > > > > > >> > > Fax: 302-497-8401 <(302)%20497-8401>
> > > > > > > >> > >
> > > > > > > >> > >
> > > > > > > >> > > On Thu, Aug 31, 2017 at 2:33 PM, Brad Diehl -
NOAA
> Federal
> > > via
> > > > > RT
> > > > > > <
> > > > > > > >> > > met_help at ucar.edu> wrote:
> > > > > > > >> > >
> > > > > > > >> > >>
> > > > > > > >> > >> Thu Aug 31 08:33:00 2017: Request 81829 was
acted upon.
> > > > > > > >> > >> Transaction: Ticket created by
brad.diehl at noaa.gov
> > > > > > > >> > >> Queue: met_help
> > > > > > > >> > >> Subject: grib2 files as input into MET
grid_stat
> > > > > > > >> > >> Owner: Nobody
> > > > > > > >> > >> Requestors: brad.diehl at noaa.gov
> > > > > > > >> > >> Status: new
> > > > > > > >> > >> Ticket <URL: https://rt.rap.ucar.edu/rt/Tic
> > > > > > > >> ket/Display.html?id=81829 >
> > > > > > > >> > >>
> > > > > > > >> > >>
> > > > > > > >> > >> I tried to use grid_stat with two grib2 files as
input,
> > and
> > > > > > > received
> > > > > > > >> an
> > > > > > > >> > >> error regarding the GRIB2 field abbreviation
'APCP06'
> > (see
> > > > > output
> > > > > > > at
> > > > > > > >> > >> bottom
> > > > > > > >> > >> of message). Yet when I run degrib on the fcst
and obs
> > > files
> > > > > > that's
> > > > > > > >> > what I
> > > > > > > >> > >> get:
> > > > > > > >> > >>
> > > > > > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> degrib
> > > > > > > >> > >> p06m_2017082312f006.grb2 -I
> > > > > > > >> > >> MsgNum, Byte, GRIB-Version, elem, level,
> reference(UTC),
> > > > > > > valid(UTC),
> > > > > > > >> > >> Proj(hr)
> > > > > > > >> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation
> [kg/(m^2)]",
> > > > > 0-SFC,
> > > > > > > >> > >> 08/23/2017
> > > > > > > >> > >> 12:00, 08/23/2017 18:00, 6.00
> > > > > > > >> > >>
> > > > > > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> degrib
> > > > > > > >> > >> urma2p5.2017082318 <(201)%20708-
2318>.pcp_06h.184.grb2
> > -I
> > > > > > > >> > >> MsgNum, Byte, GRIB-Version, elem, level,
> reference(UTC),
> > > > > > > valid(UTC),
> > > > > > > >> > >> Proj(hr)
> > > > > > > >> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation
> [kg/(m^2)]",
> > > > > 0-SFC,
> > > > > > > >> > >> 08/23/2017
> > > > > > > >> > >> 12:00, 08/23/2017 18:00, 6.00
> > > > > > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > > > > > > >> > >>
> > > > > > > >> > >> If I put both the fcst and obs files through
> pcp_combine
> > > and
> > > > > use
> > > > > > > the
> > > > > > > >> > >> resultant .nc files, grid_stat runs fine, with
the
> caveat
> > > > that
> > > > > > the
> > > > > > > >> field
> > > > > > > >> > >> name becomes 'APCP_06'.
> > > > > > > >> > >>
> > > > > > > >> > >> Thanks for any comments/suggestions, Brad Diehl
> > > > > > > >> > >>
> > > > > > > >> > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > grid_stat
> > > > > > > >> > >> p06m_2017082312f006.grb2 urma2p5.2017082318
> > > > <(201)%20708-2318>
> > > > > > > >> > >> .pcp_06h.184.grb2
> > > > > > > >> > >> config/GridStatConfig_APCP_06_grib -outdir
> out/grid_stat
> > > -v
> > > > > > 4DEBUG
> > > > > > > >> 1:
> > > > > > > >> > >> Default Config File:
> > > > > > > >> > >> /usrx/local/dev/MET/6.0/share/
> met/config/GridStatConfig_
> > > > > default
> > > > > > > >> > >> DEBUG 1: User Config File:
> config/GridStatConfig_APCP_06_
> > > > grib
> > > > > > > >> > >> DEBUG 4:
Met2dDataFileFactory::new_met_2d_data_file()
> ->
> > > > > created
> > > > > > > new
> > > > > > > >> > >> Met2dDataFile object of type "FileType_Gb2".
> > > > > > > >> > >> WARNING:
> > > > > > > >> > >> WARNING: MetGrib2DataFile::read_grib2_grid() ->
MET
> does
> > > not
> > > > > > > >> currently
> > > > > > > >> > >> support Lambert Conformal grids where dx
(2.5397) != dy
> > > > > (2.5397)
> > > > > > > and
> > > > > > > >> may
> > > > > > > >> > >> produce unexpected results!
> > > > > > > >> > >> WARNING:
> > > > > > > >> > >> DEBUG 4:
> > > > > > > >> > >> DEBUG 4: Lambert Conformal Grid Data:
> > > > > > > >> > >> DEBUG 4: scale_lat_1: 25
> > > > > > > >> > >> DEBUG 4: scale_lat_2: 25
> > > > > > > >> > >> DEBUG 4: lat_pin: 20.192
> > > > > > > >> > >> DEBUG 4: lon_pin: 121.554
> > > > > > > >> > >> DEBUG 4: x_pin: 0
> > > > > > > >> > >> DEBUG 4: y_pin: 0
> > > > > > > >> > >> DEBUG 4: lon_orient: 95
> > > > > > > >> > >> DEBUG 4: d_km: 2.5397
> > > > > > > >> > >> DEBUG 4: r_km: 6371.2
> > > > > > > >> > >> DEBUG 4: nx: 2145
> > > > > > > >> > >> DEBUG 4: ny: 1597
> > > > > > > >> > >> DEBUG 4:
> > > > > > > >> > >> DEBUG 4:
Met2dDataFileFactory::new_met_2d_data_file()
> ->
> > > > > created
> > > > > > > new
> > > > > > > >> > >> Met2dDataFile object of type "FileType_Gb2".
> > > > > > > >> > >> DEBUG 4:
> > > > > > > >> > >> DEBUG 4: Lambert Conformal Grid Data:
> > > > > > > >> > >> DEBUG 4: scale_lat_1: 25
> > > > > > > >> > >> DEBUG 4: scale_lat_2: 25
> > > > > > > >> > >> DEBUG 4: lat_pin: 20.192
> > > > > > > >> > >> DEBUG 4: lon_pin: 121.554
> > > > > > > >> > >> DEBUG 4: x_pin: 0
> > > > > > > >> > >> DEBUG 4: y_pin: 0
> > > > > > > >> > >> DEBUG 4: lon_orient: 95
> > > > > > > >> > >> DEBUG 4: d_km: 2.5397
> > > > > > > >> > >> DEBUG 4: r_km: 6371.2
> > > > > > > >> > >> DEBUG 4: nx: 2145
> > > > > > > >> > >> DEBUG 4: ny: 1377
> > > > > > > >> > >> DEBUG 4:
> > > > > > > >> > >> DEBUG 4: VarInfoFactory::new_var_info() ->
created new
> > > > VarInfo
> > > > > > > >> object of
> > > > > > > >> > >> type "FileType_Gb2".
> > > > > > > >> > >> DEBUG 4: VarInfoFactory::new_var_info() ->
created new
> > > > VarInfo
> > > > > > > >> object of
> > > > > > > >> > >> type "FileType_Gb2".
> > > > > > > >> > >> ERROR :
> > > > > > > >> > >> ERROR : VarInfoGrib2::set_dict() ->
unrecognized GRIB2
> > > field
> > > > > > > >> > abbreviation
> > > > > > > >> > >> 'APCP06'
> > > > > > > >> > >> ERROR :
> > > > > > > >> > >>
> > > > > > > >> > >>
> > > > > > > >> > >
> > > > > > > >> >
> > > > > > > >> >
> > > > > > > >>
> > > > > > > >>
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: grib2 files as input into MET grid_stat
From: Brad Diehl - NOAA Federal
Time: Fri Sep 01 12:02:18 2017
Thanks for your reply again. As you stated, the primary output is text
files. But we are working on the very early development of the
Quantitative
Precipitation Forecast Verification System for the NWS, and will want
to
give the users the ability to look at the difference fields, for
example
the ability to see the difference fields geo-referenced, possibly
using the
netcdf output. So we're in the early stages of assessing how we are
going
to do all this.
Anyway, thanks for your help earlier, and have a nice long
weekend....Brad
On Fri, Sep 1, 2017 at 4:47 PM, John Halley Gotway via RT
<met_help at ucar.edu
> wrote:
> Brad,
>
> It won't be in the next version, 6.1, and we don't currently have it
> assigned to version 6.2 either. So no specific timeline.
>
> How does this affect the decisions you're making? In our
application of
> MET to projects in the Development Testbed Center (DTC), the primary
output
> of MET are the ASCII statistics files, to which NetCDF CF
conventions do
> not apply. We load them into a METViewer database, aggregate
through time,
> and make plots to illustrate model performance.
>
> Which NetCDF files produced by the MET tools are of the most
concern?
>
> Thanks,
> John
>
>
>
> On Fri, Sep 1, 2017 at 10:29 AM, Brad Diehl - NOAA Federal via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> >
> > Is there any kind of time frame to make that switch over to
compliance?
> >
> > Brad Diehl
> >
> > On Fri, Sep 1, 2017 at 3:35 PM, John Halley Gotway via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > > Hi Brad,
> > >
> > > Unfortunately, no, the NetCDF output generated by the MET tools
> > themselves
> > > are not CF-compliant.
> > >
> > > Since we always suggest that others use that convention, we
should too!
> > > And we do plan to make that switch over eventually, but haven't
done so
> > > yet.
> > >
> > > Thanks,
> > > John
> > >
> > > On Thu, Aug 31, 2017 at 1:15 PM, Brad Diehl - NOAA Federal via
RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829
>
> > > >
> > > > Thanks for your help and suggestions today. A question came up
from
> one
> > > of
> > > > our developers regarding netcdf and CF compliance. Is the
output of
> > > MET/6.0
> > > > netcdf CF compliant?
> > > >
> > > > Thanks in advance.......Brad
> > > >
> > > > On Thu, Aug 31, 2017 at 6:35 PM, Julie Prestopnik via RT <
> > > > met_help at ucar.edu>
> > > > wrote:
> > > >
> > > > > Hi Brad,
> > > > >
> > > > > I just thought I'd add that you can find more information
for
> > > file-format
> > > > > specific settings for the "field" entry (including level)
for MET
> on
> > > pg.
> > > > 53
> > > > > of the MET User's Guide for met-6.0:
> > > > >
> > > > > http://www.dtcenter.org/met/users/docs/users_guide/MET_
> > > > > Users_Guide_v6.0.pdf
> > > > >
> > > > > or in the README file in the met-6.0/data/config area of
your
> > install.
> > > > >
> > > > > For example:
> > > > > // File-format specific settings for the "field" entry:
> > > > > //
> > > > > // - GRIB1 and GRIB2:
> > > > > // - For custom GRIB tables, see note about
MET_GRIB_TABLES.
> > > > > // - The "name" entry specifies a GRIB code number or
> > > abbreviation.
> > > > > // - GRIB1 Product Definition Section:
> > > > > // http://www.nco.ncep.noaa.gov/
> pmb/docs/on388/table2.html
> > > > > // - GRIB2 Product Definition Section:
> > > > > // http://www.nco.ncep.noaa.gov/
> > > pmb/docs/grib2/grib2_doc.shtml
> > > > > // - The "level" entry specifies a level type and
value:
> > > > > // - ANNN for accumulation interval NNN
> > > > > // - ZNNN for vertical level NNN
> > > > > // - ZNNN-NNN for a range of vertical levels
> > > > > // - PNNN for pressure level NNN in hPa
> > > > > // - PNNN-NNN for a range of pressure levels in hPa
> > > > > // - LNNN for a generic level type
> > > > > // - RNNN for a specific GRIB record number
> > > > >
> > > > > I hope that helps!
> > > > >
> > > > > Julie
> > > > >
> > > > > On Thu, Aug 31, 2017 at 12:25 PM, Minna Win via RT <
> > met_help at ucar.edu>
> > > > > wrote:
> > > > >
> > > > > >
> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> > > > > >
> > > > > > Hi Brad,
> > > > > >
> > > > > > For GRIB2, the level (if not an accumulated precip type)
is
> > typically
> > > > > "Z0"
> > > > > > for surface (sometimes "L0" works, too), but for other
levels
> such
> > as
> > > > > > entire atmosphere, you will need to refer to the GRIB2
table
> "Table
> > > 2.
> > > > > > Parameters and Units":
> > > > > > http://www.wmo.int/pages/prog/www/WDM/Guides/Guide-binary-
2.
> > > html#GRIB
> > > > .
> > > > > > From that table, find the value that corresponds to your
level
> type
> > > and
> > > > > set
> > > > > > it to GRIB_lvl_typ and then set GRIB_lvl_val1 to the
appropriate
> > > level.
> > > > > > For 'entire atmosphere', there is no numerical value for
the
> level,
> > > so
> > > > it
> > > > > > would be set to 0 or 255 (you might have to try both to
see what
> > > works,
> > > > > it
> > > > > > depends on the grib file). For something like a pressure
> > difference,
> > > > you
> > > > > > would need to set GRIB_lvl_typ and the GRIB_lvl_va1 and
> > > GRIB_lvl_val2,
> > > > > > i.e. GRIB_lvl_val2=top level GRIB_lvl_val1=bottom level,
which
> you
> > > > would
> > > > > > obtain from performing a wgrib2/grep combination on your
grib
> file
> > > and
> > > > > > variable of interest. You can use the MET tool
plot_data_plane to
> > > check
> > > > > > that your GRIB_lvl_typ and GRIB_lvl_val1 (and
GRIB_lvl_val2 if
> > used)
> > > > > > produced a match to the record of interest in the grib
file.
> > > > > >
> > > > > > As for your question about the (*,*) you see in some
examples,
> > there
> > > is
> > > > > an
> > > > > > explanation in Chapter 3 of the MET User's Manual
(excerpted):
> > > > > >
> > > > > >
> > > > > > > // - The "level" entry specifies the dimensions to be
used:
> > > > > > > // - (i,...,j,*,*) for a single field, where i,...,j
> specifies
> > > > fixed
> > > > > > > // dimension values and *,* specifies the two
dimensions
> > for
> > > > the
> > > > > > > // gridded field.
> > > > > > >
> > > > > >
> > > > > >
> > > > > > field = [
> > > > > > {
> > > > > > name = "QVAPOR";
> > > > > > level = "( = "(0,5,*,*)"
> > > > > >
> > > > > > },
> > > > > > {
> > > > > >
> > > > > > name = "TMP_P850_ENS_MEAN;
> > > > > > level = [ "(*,*)" ];
> > > > > > }
> > > > > > ];
> > > > > >
> > > > > >
> > > > > > I hope that clarifies things. Please don't hesitate to
contact
> us
> > if
> > > > you
> > > > > > have further questions.
> > > > > >
> > > > > > Regards,
> > > > > > Minna
> > > > > >
> > > > > > ---------------
> > > > > > Minna Win
> > > > > > NCAR
> > > > > > Research Applications Lab
> > > > > > Phone: 303-497-8423
> > > > > > Fax: 302-497-8401
> > > > > >
> > > > > >
> > > > > > On Thu, Aug 31, 2017 at 4:21 PM, Brad Diehl - NOAA Federal
via
> RT <
> > > > > > met_help at ucar.edu> wrote:
> > > > > >
> > > > > > >
> > > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829
> >
> > > > > > >
> > > > > > > the A06 worked. Prior to trying that I sent you via
mailx the
> > > > fcst/obs
> > > > > > > grib2 files, you can ignore that now.
> > > > > > >
> > > > > > > But it does beg the question, how does one find the
level
> names,
> > is
> > > > it
> > > > > in
> > > > > > > the grib file itself? Very confusing, at least to me.
And what
> > does
> > > > *,*
> > > > > > > mean when I see that in some examples?
> > > > > > >
> > > > > > > Thanks for all your help........Brad
> > > > > > >
> > > > > > > On Thu, Aug 31, 2017 at 4:11 PM, Minna Win via RT <
> > > met_help at ucar.edu
> > > > >
> > > > > > > wrote:
> > > > > > >
> > > > > > > > Hi Brad,
> > > > > > > >
> > > > > > > > Tara Jensen brought this to my attention, for verify
APCP,
> you
> > > need
> > > > > to
> > > > > > > set
> > > > > > > > the level to the value pre-pended with the letter 'A'.
For
> > > > instance,
> > > > > > > for 6
> > > > > > > > hour accumulation level, set level="A06".
> > > > > > > >
> > > > > > > > Thanks,
> > > > > > > > Minna
> > > > > > > >
> > > > > > > > ---------------
> > > > > > > > Minna Win
> > > > > > > > NCAR
> > > > > > > > Research Applications Lab
> > > > > > > > Phone: 303-497-8423
> > > > > > > > Fax: 302-497-8401
> > > > > > > >
> > > > > > > >
> > > > > > > > On Thu, Aug 31, 2017 at 3:52 PM, Minna Win-
Gildenmeister <
> > > > > > > > minnawin at ucar.edu>
> > > > > > > > wrote:
> > > > > > > >
> > > > > > > > > Hi Brad,
> > > > > > > > >
> > > > > > > > > Can you try "L0", and "P0"? If that doesn't work,
could
> you
> > > > please
> > > > > > > send
> > > > > > > > > me your grib file? Instructions are as follows:
> > > > > > > > >
> > > > > > > > > Non-Windows FTP client:
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > > * ftp -p ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
> > > > > > username=anonymous
> > > > > > > > > password="your email address" cd
incoming/irap/met_help
> > > mkdir
> > > > > > > > > "surname_data" # where surname is your last name
cd
> > > > > "surname_data"
> > > > > > > put
> > > > > > > > > "your files" bye *
> > > > > > > > >
> > > > > > > > > Windows FTP client:
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > > > * ftp ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
> > > > > username=anonymous
> > > > > > > > > password="your email address" quote pasv cd
> > > > incoming/irap/met_help
> > > > > > > mkdir
> > > > > > > > > "surname_data" # where surname is your last name cd
> > > > "surname_data"
> > > > > > put
> > > > > > > > > "your files" bye*
> > > > > > > > >
> > > > > > > > > I can take a look at the level type and value, you
might
> need
> > > to
> > > > > use
> > > > > > > > > special syntax in your config file to get grid_stat
to
> work.
> > > Let
> > > > > me
> > > > > > > know
> > > > > > > > > if you ftp'd a grib file so I can grab it before the
file
> > gets
> > > > > > > scrubbed.
> > > > > > > > >
> > > > > > > > > Thanks,
> > > > > > > > > Minna
> > > > > > > > >
> > > > > > > > >
> > > > > > > > > ---------------
> > > > > > > > > Minna Win
> > > > > > > > > NCAR
> > > > > > > > > Research Applications Lab
> > > > > > > > > Phone: 303-497-8423 <(303)%20497-8423>
> > > > > > > > > Fax: 302-497-8401 <(302)%20497-8401>
> > > > > > > > >
> > > > > > > > >
> > > > > > > > > On Thu, Aug 31, 2017 at 3:09 PM, Brad Diehl - NOAA
Federal
> > via
> > > > RT <
> > > > > > > > > met_help at ucar.edu> wrote:
> > > > > > > > >
> > > > > > > > >>
> > > > > > > > >> <URL: https://rt.rap.ucar.edu/rt/
> > Ticket/Display.html?id=81829
> > > >
> > > > > > > > >>
> > > > > > > > >> After running pcp_combine, everything did work
correctly,
> it
> > > is
> > > > > the
> > > > > > > > grib2
> > > > > > > > >> fcts/obs part I'm having trouble with. I think you
put me
> on
> > > the
> > > > > > right
> > > > > > > > >> track by suggesting wgrib2. Here's what I get using
wgrib2
> > on
> > > > the
> > > > > > fcst
> > > > > > > > >> file:
> > > > > > > > >>
> > > > > > > > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
wgrib2
> > > > > > > > >> p06m_2017082312f006.grb2
> > > > > > > > >> 1:21:d=2017082312 <(201)%20708-2312>:APCP:surface -
> > > > > > surface:CodeTable
> > > > > > > > >> 4.11=missing:
> > > > > > > > >>
> > > > > > > > >> in the config file I'm using:
> > > > > > > > >>
> > > > > > > > >>
> > > > > > > > >> fcst = {
> > > > > > > > >>
> > > > > > > > >> field = [
> > > > > > > > >> {
> > > > > > > > >> name = "APCP";
> > > > > > > > >> level = "Z0" ;
> > > > > > > > >> cat_thresh = [ >0.0, >=5.0 , >=10.0];
> > > > > > > > >> }
> > > > > > > > >> ];
> > > > > > > > >>
> > > > > > > > >> for level, should I be using surface or Z0?,
neither of
> > which
> > > > > seemed
> > > > > > > to
> > > > > > > > >> work.
> > > > > > > > >>
> > > > > > > > >> On Thu, Aug 31, 2017 at 2:58 PM, Minna Win via RT <
> > > > > > met_help at ucar.edu>
> > > > > > > > >> wrote:
> > > > > > > > >>
> > > > > > > > >> > Hi Brad,
> > > > > > > > >> >
> > > > > > > > >> > Please disregard my previous response. I didn't
see
> that
> > > you
> > > > > had
> > > > > > > run
> > > > > > > > >> > pcp_combine-that *should* have worked. I'll look
a
> little
> > > > > further
> > > > > > > or
> > > > > > > > >> see
> > > > > > > > >> > if someone else has seen this before.
> > > > > > > > >> >
> > > > > > > > >> >
> > > > > > > > >> >
> > > > > > > > >> > Regards,
> > > > > > > > >> > Minna
> > > > > > > > >> >
> > > > > > > > >> > ---------------
> > > > > > > > >> > Minna Win
> > > > > > > > >> > NCAR
> > > > > > > > >> > Research Applications Lab
> > > > > > > > >> > Phone: 303-497-8423
> > > > > > > > >> > Fax: 302-497-8401
> > > > > > > > >> >
> > > > > > > > >> >
> > > > > > > > >> > On Thu, Aug 31, 2017 at 2:50 PM, Minna
> Win-Gildenmeister <
> > > > > > > > >> > minnawin at ucar.edu>
> > > > > > > > >> > wrote:
> > > > > > > > >> >
> > > > > > > > >> > > HI Brad,
> > > > > > > > >> > >
> > > > > > > > >> > > wgrib is an NDFD tool and I don't know how it
is
> > > > interpreting
> > > > > > the
> > > > > > > > APCP
> > > > > > > > >> > > variable in your grib file. Try using wgrib2 on
your
> > grib
> > > > file
> > > > > > and
> > > > > > > > >> grep
> > > > > > > > >> > for
> > > > > > > > >> > > APCP and use that variable abbreviation in you
MET
> > > Grid_stat
> > > > > > > config
> > > > > > > > >> > file. I
> > > > > > > > >> > > hope that helps. Please let us know if you
still have
> > > > > problems
> > > > > > > > using
> > > > > > > > >> > > Grid_stat.
> > > > > > > > >> > >
> > > > > > > > >> > > Regards,
> > > > > > > > >> > > Minna
> > > > > > > > >> > >
> > > > > > > > >> > > ---------------
> > > > > > > > >> > > Minna Win
> > > > > > > > >> > > NCAR
> > > > > > > > >> > > Research Applications Lab
> > > > > > > > >> > > Phone: 303-497-8423 <(303)%20497-8423>
> > > > > > > > >> > > Fax: 302-497-8401 <(302)%20497-8401>
> > > > > > > > >> > >
> > > > > > > > >> > >
> > > > > > > > >> > > On Thu, Aug 31, 2017 at 2:33 PM, Brad Diehl -
NOAA
> > Federal
> > > > via
> > > > > > RT
> > > > > > > <
> > > > > > > > >> > > met_help at ucar.edu> wrote:
> > > > > > > > >> > >
> > > > > > > > >> > >>
> > > > > > > > >> > >> Thu Aug 31 08:33:00 2017: Request 81829 was
acted
> upon.
> > > > > > > > >> > >> Transaction: Ticket created by
brad.diehl at noaa.gov
> > > > > > > > >> > >> Queue: met_help
> > > > > > > > >> > >> Subject: grib2 files as input into MET
grid_stat
> > > > > > > > >> > >> Owner: Nobody
> > > > > > > > >> > >> Requestors: brad.diehl at noaa.gov
> > > > > > > > >> > >> Status: new
> > > > > > > > >> > >> Ticket <URL: https://rt.rap.ucar.edu/rt/Tic
> > > > > > > > >> ket/Display.html?id=81829 >
> > > > > > > > >> > >>
> > > > > > > > >> > >>
> > > > > > > > >> > >> I tried to use grid_stat with two grib2 files
as
> input,
> > > and
> > > > > > > > received
> > > > > > > > >> an
> > > > > > > > >> > >> error regarding the GRIB2 field abbreviation
'APCP06'
> > > (see
> > > > > > output
> > > > > > > > at
> > > > > > > > >> > >> bottom
> > > > > > > > >> > >> of message). Yet when I run degrib on the fcst
and
> obs
> > > > files
> > > > > > > that's
> > > > > > > > >> > what I
> > > > > > > > >> > >> get:
> > > > > > > > >> > >>
> > > > > > > > >> > >> (g14a2)
/u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > degrib
> > > > > > > > >> > >> p06m_2017082312f006.grb2 -I
> > > > > > > > >> > >> MsgNum, Byte, GRIB-Version, elem, level,
> > reference(UTC),
> > > > > > > > valid(UTC),
> > > > > > > > >> > >> Proj(hr)
> > > > > > > > >> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation
> > [kg/(m^2)]",
> > > > > > 0-SFC,
> > > > > > > > >> > >> 08/23/2017
> > > > > > > > >> > >> 12:00, 08/23/2017 18:00, 6.00
> > > > > > > > >> > >>
> > > > > > > > >> > >> (g14a2)
/u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > degrib
> > > > > > > > >> > >> urma2p5.2017082318 <(201)%20708-2318>.pcp_06h.
> 184.grb2
> > > -I
> > > > > > > > >> > >> MsgNum, Byte, GRIB-Version, elem, level,
> > reference(UTC),
> > > > > > > > valid(UTC),
> > > > > > > > >> > >> Proj(hr)
> > > > > > > > >> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation
> > [kg/(m^2)]",
> > > > > > 0-SFC,
> > > > > > > > >> > >> 08/23/2017
> > > > > > > > >> > >> 12:00, 08/23/2017 18:00, 6.00
> > > > > > > > >> > >> (g14a2)
/u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > > > > > > > >> > >>
> > > > > > > > >> > >> If I put both the fcst and obs files through
> > pcp_combine
> > > > and
> > > > > > use
> > > > > > > > the
> > > > > > > > >> > >> resultant .nc files, grid_stat runs fine, with
the
> > caveat
> > > > > that
> > > > > > > the
> > > > > > > > >> field
> > > > > > > > >> > >> name becomes 'APCP_06'.
> > > > > > > > >> > >>
> > > > > > > > >> > >> Thanks for any comments/suggestions, Brad
Diehl
> > > > > > > > >> > >>
> > > > > > > > >> > >> (g14a2)
/u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > > grid_stat
> > > > > > > > >> > >> p06m_2017082312f006.grb2 urma2p5.2017082318
> > > > > <(201)%20708-2318>
> > > > > > > > >> > >> .pcp_06h.184.grb2
> > > > > > > > >> > >> config/GridStatConfig_APCP_06_grib -outdir
> > out/grid_stat
> > > > -v
> > > > > > > 4DEBUG
> > > > > > > > >> 1:
> > > > > > > > >> > >> Default Config File:
> > > > > > > > >> > >> /usrx/local/dev/MET/6.0/share/
> > met/config/GridStatConfig_
> > > > > > default
> > > > > > > > >> > >> DEBUG 1: User Config File:
> > config/GridStatConfig_APCP_06_
> > > > > grib
> > > > > > > > >> > >> DEBUG 4: Met2dDataFileFactory::new_met_
> 2d_data_file()
> > ->
> > > > > > created
> > > > > > > > new
> > > > > > > > >> > >> Met2dDataFile object of type "FileType_Gb2".
> > > > > > > > >> > >> WARNING:
> > > > > > > > >> > >> WARNING: MetGrib2DataFile::read_grib2_grid()
-> MET
> > does
> > > > not
> > > > > > > > >> currently
> > > > > > > > >> > >> support Lambert Conformal grids where dx
(2.5397) !=
> dy
> > > > > > (2.5397)
> > > > > > > > and
> > > > > > > > >> may
> > > > > > > > >> > >> produce unexpected results!
> > > > > > > > >> > >> WARNING:
> > > > > > > > >> > >> DEBUG 4:
> > > > > > > > >> > >> DEBUG 4: Lambert Conformal Grid Data:
> > > > > > > > >> > >> DEBUG 4: scale_lat_1: 25
> > > > > > > > >> > >> DEBUG 4: scale_lat_2: 25
> > > > > > > > >> > >> DEBUG 4: lat_pin: 20.192
> > > > > > > > >> > >> DEBUG 4: lon_pin: 121.554
> > > > > > > > >> > >> DEBUG 4: x_pin: 0
> > > > > > > > >> > >> DEBUG 4: y_pin: 0
> > > > > > > > >> > >> DEBUG 4: lon_orient: 95
> > > > > > > > >> > >> DEBUG 4: d_km: 2.5397
> > > > > > > > >> > >> DEBUG 4: r_km: 6371.2
> > > > > > > > >> > >> DEBUG 4: nx: 2145
> > > > > > > > >> > >> DEBUG 4: ny: 1597
> > > > > > > > >> > >> DEBUG 4:
> > > > > > > > >> > >> DEBUG 4: Met2dDataFileFactory::new_met_
> 2d_data_file()
> > ->
> > > > > > created
> > > > > > > > new
> > > > > > > > >> > >> Met2dDataFile object of type "FileType_Gb2".
> > > > > > > > >> > >> DEBUG 4:
> > > > > > > > >> > >> DEBUG 4: Lambert Conformal Grid Data:
> > > > > > > > >> > >> DEBUG 4: scale_lat_1: 25
> > > > > > > > >> > >> DEBUG 4: scale_lat_2: 25
> > > > > > > > >> > >> DEBUG 4: lat_pin: 20.192
> > > > > > > > >> > >> DEBUG 4: lon_pin: 121.554
> > > > > > > > >> > >> DEBUG 4: x_pin: 0
> > > > > > > > >> > >> DEBUG 4: y_pin: 0
> > > > > > > > >> > >> DEBUG 4: lon_orient: 95
> > > > > > > > >> > >> DEBUG 4: d_km: 2.5397
> > > > > > > > >> > >> DEBUG 4: r_km: 6371.2
> > > > > > > > >> > >> DEBUG 4: nx: 2145
> > > > > > > > >> > >> DEBUG 4: ny: 1377
> > > > > > > > >> > >> DEBUG 4:
> > > > > > > > >> > >> DEBUG 4: VarInfoFactory::new_var_info() ->
created
> new
> > > > > VarInfo
> > > > > > > > >> object of
> > > > > > > > >> > >> type "FileType_Gb2".
> > > > > > > > >> > >> DEBUG 4: VarInfoFactory::new_var_info() ->
created
> new
> > > > > VarInfo
> > > > > > > > >> object of
> > > > > > > > >> > >> type "FileType_Gb2".
> > > > > > > > >> > >> ERROR :
> > > > > > > > >> > >> ERROR : VarInfoGrib2::set_dict() ->
unrecognized
> GRIB2
> > > > field
> > > > > > > > >> > abbreviation
> > > > > > > > >> > >> 'APCP06'
> > > > > > > > >> > >> ERROR :
> > > > > > > > >> > >>
> > > > > > > > >> > >>
> > > > > > > > >> > >
> > > > > > > > >> >
> > > > > > > > >> >
> > > > > > > > >>
> > > > > > > > >>
> > > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: grib2 files as input into MET grid_stat
From: John Halley Gotway
Time: Fri Sep 01 12:41:20 2017
Brad,
OK, thanks for the info. We would certainly be interested in working
closely with you in your development of QPF verification for NWS. If
MET
proves useful, I wouldn't let the file formats stand in the way.
Using a
script to reformat the existing MET NetCDF output to follow CF
conventions
would be pretty straight-forward with the longer term goal of MET
producing
them directly.
I'll be in DC visiting NCEP 9/19 - 9/22 and would be happy to discuss
any
issues or questions and/or you could talk directly with the MET
project
lead, Tara Jensen.
Thanks,
John
On Fri, Sep 1, 2017 at 12:02 PM, Brad Diehl - NOAA Federal via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
>
> Thanks for your reply again. As you stated, the primary output is
text
> files. But we are working on the very early development of the
Quantitative
> Precipitation Forecast Verification System for the NWS, and will
want to
> give the users the ability to look at the difference fields, for
example
> the ability to see the difference fields geo-referenced, possibly
using the
> netcdf output. So we're in the early stages of assessing how we are
going
> to do all this.
>
> Anyway, thanks for your help earlier, and have a nice long
weekend....Brad
>
> On Fri, Sep 1, 2017 at 4:47 PM, John Halley Gotway via RT <
> met_help at ucar.edu
> > wrote:
>
> > Brad,
> >
> > It won't be in the next version, 6.1, and we don't currently have
it
> > assigned to version 6.2 either. So no specific timeline.
> >
> > How does this affect the decisions you're making? In our
application of
> > MET to projects in the Development Testbed Center (DTC), the
primary
> output
> > of MET are the ASCII statistics files, to which NetCDF CF
conventions do
> > not apply. We load them into a METViewer database, aggregate
through
> time,
> > and make plots to illustrate model performance.
> >
> > Which NetCDF files produced by the MET tools are of the most
concern?
> >
> > Thanks,
> > John
> >
> >
> >
> > On Fri, Sep 1, 2017 at 10:29 AM, Brad Diehl - NOAA Federal via RT
<
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> > >
> > > Is there any kind of time frame to make that switch over to
compliance?
> > >
> > > Brad Diehl
> > >
> > > On Fri, Sep 1, 2017 at 3:35 PM, John Halley Gotway via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > > Hi Brad,
> > > >
> > > > Unfortunately, no, the NetCDF output generated by the MET
tools
> > > themselves
> > > > are not CF-compliant.
> > > >
> > > > Since we always suggest that others use that convention, we
should
> too!
> > > > And we do plan to make that switch over eventually, but
haven't done
> so
> > > > yet.
> > > >
> > > > Thanks,
> > > > John
> > > >
> > > > On Thu, Aug 31, 2017 at 1:15 PM, Brad Diehl - NOAA Federal via
RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> > > > >
> > > > > Thanks for your help and suggestions today. A question came
up from
> > one
> > > > of
> > > > > our developers regarding netcdf and CF compliance. Is the
output of
> > > > MET/6.0
> > > > > netcdf CF compliant?
> > > > >
> > > > > Thanks in advance.......Brad
> > > > >
> > > > > On Thu, Aug 31, 2017 at 6:35 PM, Julie Prestopnik via RT <
> > > > > met_help at ucar.edu>
> > > > > wrote:
> > > > >
> > > > > > Hi Brad,
> > > > > >
> > > > > > I just thought I'd add that you can find more information
for
> > > > file-format
> > > > > > specific settings for the "field" entry (including level)
for MET
> > on
> > > > pg.
> > > > > 53
> > > > > > of the MET User's Guide for met-6.0:
> > > > > >
> > > > > > http://www.dtcenter.org/met/users/docs/users_guide/MET_
> > > > > > Users_Guide_v6.0.pdf
> > > > > >
> > > > > > or in the README file in the met-6.0/data/config area of
your
> > > install.
> > > > > >
> > > > > > For example:
> > > > > > // File-format specific settings for the "field" entry:
> > > > > > //
> > > > > > // - GRIB1 and GRIB2:
> > > > > > // - For custom GRIB tables, see note about
> MET_GRIB_TABLES.
> > > > > > // - The "name" entry specifies a GRIB code number
or
> > > > abbreviation.
> > > > > > // - GRIB1 Product Definition Section:
> > > > > > // http://www.nco.ncep.noaa.gov/
> > pmb/docs/on388/table2.html
> > > > > > // - GRIB2 Product Definition Section:
> > > > > > // http://www.nco.ncep.noaa.gov/
> > > > pmb/docs/grib2/grib2_doc.shtml
> > > > > > // - The "level" entry specifies a level type and
value:
> > > > > > // - ANNN for accumulation interval NNN
> > > > > > // - ZNNN for vertical level NNN
> > > > > > // - ZNNN-NNN for a range of vertical levels
> > > > > > // - PNNN for pressure level NNN in hPa
> > > > > > // - PNNN-NNN for a range of pressure levels in
hPa
> > > > > > // - LNNN for a generic level type
> > > > > > // - RNNN for a specific GRIB record number
> > > > > >
> > > > > > I hope that helps!
> > > > > >
> > > > > > Julie
> > > > > >
> > > > > > On Thu, Aug 31, 2017 at 12:25 PM, Minna Win via RT <
> > > met_help at ucar.edu>
> > > > > > wrote:
> > > > > >
> > > > > > >
> > > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829
> >
> > > > > > >
> > > > > > > Hi Brad,
> > > > > > >
> > > > > > > For GRIB2, the level (if not an accumulated precip type)
is
> > > typically
> > > > > > "Z0"
> > > > > > > for surface (sometimes "L0" works, too), but for other
levels
> > such
> > > as
> > > > > > > entire atmosphere, you will need to refer to the GRIB2
table
> > "Table
> > > > 2.
> > > > > > > Parameters and Units":
> > > > > > > http://www.wmo.int/pages/prog/www/WDM/Guides/Guide-
binary-2.
> > > > html#GRIB
> > > > > .
> > > > > > > From that table, find the value that corresponds to your
level
> > type
> > > > and
> > > > > > set
> > > > > > > it to GRIB_lvl_typ and then set GRIB_lvl_val1 to the
> appropriate
> > > > level.
> > > > > > > For 'entire atmosphere', there is no numerical value for
the
> > level,
> > > > so
> > > > > it
> > > > > > > would be set to 0 or 255 (you might have to try both to
see
> what
> > > > works,
> > > > > > it
> > > > > > > depends on the grib file). For something like a
pressure
> > > difference,
> > > > > you
> > > > > > > would need to set GRIB_lvl_typ and the GRIB_lvl_va1 and
> > > > GRIB_lvl_val2,
> > > > > > > i.e. GRIB_lvl_val2=top level GRIB_lvl_val1=bottom level,
which
> > you
> > > > > would
> > > > > > > obtain from performing a wgrib2/grep combination on your
grib
> > file
> > > > and
> > > > > > > variable of interest. You can use the MET tool
plot_data_plane
> to
> > > > check
> > > > > > > that your GRIB_lvl_typ and GRIB_lvl_val1 (and
GRIB_lvl_val2 if
> > > used)
> > > > > > > produced a match to the record of interest in the grib
file.
> > > > > > >
> > > > > > > As for your question about the (*,*) you see in some
examples,
> > > there
> > > > is
> > > > > > an
> > > > > > > explanation in Chapter 3 of the MET User's Manual
(excerpted):
> > > > > > >
> > > > > > >
> > > > > > > > // - The "level" entry specifies the dimensions to be
used:
> > > > > > > > // - (i,...,j,*,*) for a single field, where
i,...,j
> > specifies
> > > > > fixed
> > > > > > > > // dimension values and *,* specifies the two
> dimensions
> > > for
> > > > > the
> > > > > > > > // gridded field.
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > field = [
> > > > > > > {
> > > > > > > name = "QVAPOR";
> > > > > > > level = "( = "(0,5,*,*)"
> > > > > > >
> > > > > > > },
> > > > > > > {
> > > > > > >
> > > > > > > name = "TMP_P850_ENS_MEAN;
> > > > > > > level = [ "(*,*)" ];
> > > > > > > }
> > > > > > > ];
> > > > > > >
> > > > > > >
> > > > > > > I hope that clarifies things. Please don't hesitate to
contact
> > us
> > > if
> > > > > you
> > > > > > > have further questions.
> > > > > > >
> > > > > > > Regards,
> > > > > > > Minna
> > > > > > >
> > > > > > > ---------------
> > > > > > > Minna Win
> > > > > > > NCAR
> > > > > > > Research Applications Lab
> > > > > > > Phone: 303-497-8423
> > > > > > > Fax: 302-497-8401
> > > > > > >
> > > > > > >
> > > > > > > On Thu, Aug 31, 2017 at 4:21 PM, Brad Diehl - NOAA
Federal via
> > RT <
> > > > > > > met_help at ucar.edu> wrote:
> > > > > > >
> > > > > > > >
> > > > > > > > <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=81829
> > >
> > > > > > > >
> > > > > > > > the A06 worked. Prior to trying that I sent you via
mailx the
> > > > > fcst/obs
> > > > > > > > grib2 files, you can ignore that now.
> > > > > > > >
> > > > > > > > But it does beg the question, how does one find the
level
> > names,
> > > is
> > > > > it
> > > > > > in
> > > > > > > > the grib file itself? Very confusing, at least to me.
And
> what
> > > does
> > > > > *,*
> > > > > > > > mean when I see that in some examples?
> > > > > > > >
> > > > > > > > Thanks for all your help........Brad
> > > > > > > >
> > > > > > > > On Thu, Aug 31, 2017 at 4:11 PM, Minna Win via RT <
> > > > met_help at ucar.edu
> > > > > >
> > > > > > > > wrote:
> > > > > > > >
> > > > > > > > > Hi Brad,
> > > > > > > > >
> > > > > > > > > Tara Jensen brought this to my attention, for verify
APCP,
> > you
> > > > need
> > > > > > to
> > > > > > > > set
> > > > > > > > > the level to the value pre-pended with the letter
'A'. For
> > > > > instance,
> > > > > > > > for 6
> > > > > > > > > hour accumulation level, set level="A06".
> > > > > > > > >
> > > > > > > > > Thanks,
> > > > > > > > > Minna
> > > > > > > > >
> > > > > > > > > ---------------
> > > > > > > > > Minna Win
> > > > > > > > > NCAR
> > > > > > > > > Research Applications Lab
> > > > > > > > > Phone: 303-497-8423
> > > > > > > > > Fax: 302-497-8401
> > > > > > > > >
> > > > > > > > >
> > > > > > > > > On Thu, Aug 31, 2017 at 3:52 PM, Minna Win-
Gildenmeister <
> > > > > > > > > minnawin at ucar.edu>
> > > > > > > > > wrote:
> > > > > > > > >
> > > > > > > > > > Hi Brad,
> > > > > > > > > >
> > > > > > > > > > Can you try "L0", and "P0"? If that doesn't work,
could
> > you
> > > > > please
> > > > > > > > send
> > > > > > > > > > me your grib file? Instructions are as follows:
> > > > > > > > > >
> > > > > > > > > > Non-Windows FTP client:
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > > * ftp -p ftp.rap.ucar.edu
<http://ftp.rap.ucar.edu>
> > > > > > > username=anonymous
> > > > > > > > > > password="your email address" cd
> incoming/irap/met_help
> > > > mkdir
> > > > > > > > > > "surname_data" # where surname is your last
name cd
> > > > > > "surname_data"
> > > > > > > > put
> > > > > > > > > > "your files" bye *
> > > > > > > > > >
> > > > > > > > > > Windows FTP client:
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > > * ftp ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
> > > > > > username=anonymous
> > > > > > > > > > password="your email address" quote pasv cd
> > > > > incoming/irap/met_help
> > > > > > > > mkdir
> > > > > > > > > > "surname_data" # where surname is your last name
cd
> > > > > "surname_data"
> > > > > > > put
> > > > > > > > > > "your files" bye*
> > > > > > > > > >
> > > > > > > > > > I can take a look at the level type and value, you
might
> > need
> > > > to
> > > > > > use
> > > > > > > > > > special syntax in your config file to get
grid_stat to
> > work.
> > > > Let
> > > > > > me
> > > > > > > > know
> > > > > > > > > > if you ftp'd a grib file so I can grab it before
the file
> > > gets
> > > > > > > > scrubbed.
> > > > > > > > > >
> > > > > > > > > > Thanks,
> > > > > > > > > > Minna
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > > ---------------
> > > > > > > > > > Minna Win
> > > > > > > > > > NCAR
> > > > > > > > > > Research Applications Lab
> > > > > > > > > > Phone: 303-497-8423 <(303)%20497-8423>
> > > > > > > > > > Fax: 302-497-8401 <(302)%20497-8401>
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > > On Thu, Aug 31, 2017 at 3:09 PM, Brad Diehl - NOAA
> Federal
> > > via
> > > > > RT <
> > > > > > > > > > met_help at ucar.edu> wrote:
> > > > > > > > > >
> > > > > > > > > >>
> > > > > > > > > >> <URL: https://rt.rap.ucar.edu/rt/
> > > Ticket/Display.html?id=81829
> > > > >
> > > > > > > > > >>
> > > > > > > > > >> After running pcp_combine, everything did work
> correctly,
> > it
> > > > is
> > > > > > the
> > > > > > > > > grib2
> > > > > > > > > >> fcts/obs part I'm having trouble with. I think
you put
> me
> > on
> > > > the
> > > > > > > right
> > > > > > > > > >> track by suggesting wgrib2. Here's what I get
using
> wgrib2
> > > on
> > > > > the
> > > > > > > fcst
> > > > > > > > > >> file:
> > > > > > > > > >>
> > > > > > > > > >> (g14a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> wgrib2
> > > > > > > > > >> p06m_2017082312f006.grb2
> > > > > > > > > >> 1:21:d=2017082312 <(201)%20708-2312>:APCP:surface
-
> > > > > > > surface:CodeTable
> > > > > > > > > >> 4.11=missing:
> > > > > > > > > >>
> > > > > > > > > >> in the config file I'm using:
> > > > > > > > > >>
> > > > > > > > > >>
> > > > > > > > > >> fcst = {
> > > > > > > > > >>
> > > > > > > > > >> field = [
> > > > > > > > > >> {
> > > > > > > > > >> name = "APCP";
> > > > > > > > > >> level = "Z0" ;
> > > > > > > > > >> cat_thresh = [ >0.0, >=5.0 , >=10.0];
> > > > > > > > > >> }
> > > > > > > > > >> ];
> > > > > > > > > >>
> > > > > > > > > >> for level, should I be using surface or Z0?,
neither of
> > > which
> > > > > > seemed
> > > > > > > > to
> > > > > > > > > >> work.
> > > > > > > > > >>
> > > > > > > > > >> On Thu, Aug 31, 2017 at 2:58 PM, Minna Win via RT
<
> > > > > > > met_help at ucar.edu>
> > > > > > > > > >> wrote:
> > > > > > > > > >>
> > > > > > > > > >> > Hi Brad,
> > > > > > > > > >> >
> > > > > > > > > >> > Please disregard my previous response. I
didn't see
> > that
> > > > you
> > > > > > had
> > > > > > > > run
> > > > > > > > > >> > pcp_combine-that *should* have worked. I'll
look a
> > little
> > > > > > further
> > > > > > > > or
> > > > > > > > > >> see
> > > > > > > > > >> > if someone else has seen this before.
> > > > > > > > > >> >
> > > > > > > > > >> >
> > > > > > > > > >> >
> > > > > > > > > >> > Regards,
> > > > > > > > > >> > Minna
> > > > > > > > > >> >
> > > > > > > > > >> > ---------------
> > > > > > > > > >> > Minna Win
> > > > > > > > > >> > NCAR
> > > > > > > > > >> > Research Applications Lab
> > > > > > > > > >> > Phone: 303-497-8423
> > > > > > > > > >> > Fax: 302-497-8401
> > > > > > > > > >> >
> > > > > > > > > >> >
> > > > > > > > > >> > On Thu, Aug 31, 2017 at 2:50 PM, Minna
> > Win-Gildenmeister <
> > > > > > > > > >> > minnawin at ucar.edu>
> > > > > > > > > >> > wrote:
> > > > > > > > > >> >
> > > > > > > > > >> > > HI Brad,
> > > > > > > > > >> > >
> > > > > > > > > >> > > wgrib is an NDFD tool and I don't know how it
is
> > > > > interpreting
> > > > > > > the
> > > > > > > > > APCP
> > > > > > > > > >> > > variable in your grib file. Try using wgrib2
on your
> > > grib
> > > > > file
> > > > > > > and
> > > > > > > > > >> grep
> > > > > > > > > >> > for
> > > > > > > > > >> > > APCP and use that variable abbreviation in
you MET
> > > > Grid_stat
> > > > > > > > config
> > > > > > > > > >> > file. I
> > > > > > > > > >> > > hope that helps. Please let us know if you
still
> have
> > > > > > problems
> > > > > > > > > using
> > > > > > > > > >> > > Grid_stat.
> > > > > > > > > >> > >
> > > > > > > > > >> > > Regards,
> > > > > > > > > >> > > Minna
> > > > > > > > > >> > >
> > > > > > > > > >> > > ---------------
> > > > > > > > > >> > > Minna Win
> > > > > > > > > >> > > NCAR
> > > > > > > > > >> > > Research Applications Lab
> > > > > > > > > >> > > Phone: 303-497-8423 <(303)%20497-8423>
> > > > > > > > > >> > > Fax: 302-497-8401 <(302)%20497-8401>
> > > > > > > > > >> > >
> > > > > > > > > >> > >
> > > > > > > > > >> > > On Thu, Aug 31, 2017 at 2:33 PM, Brad Diehl -
NOAA
> > > Federal
> > > > > via
> > > > > > > RT
> > > > > > > > <
> > > > > > > > > >> > > met_help at ucar.edu> wrote:
> > > > > > > > > >> > >
> > > > > > > > > >> > >>
> > > > > > > > > >> > >> Thu Aug 31 08:33:00 2017: Request 81829 was
acted
> > upon.
> > > > > > > > > >> > >> Transaction: Ticket created by
brad.diehl at noaa.gov
> > > > > > > > > >> > >> Queue: met_help
> > > > > > > > > >> > >> Subject: grib2 files as input into MET
> grid_stat
> > > > > > > > > >> > >> Owner: Nobody
> > > > > > > > > >> > >> Requestors: brad.diehl at noaa.gov
> > > > > > > > > >> > >> Status: new
> > > > > > > > > >> > >> Ticket <URL: https://rt.rap.ucar.edu/rt/Tic
> > > > > > > > > >> ket/Display.html?id=81829 >
> > > > > > > > > >> > >>
> > > > > > > > > >> > >>
> > > > > > > > > >> > >> I tried to use grid_stat with two grib2
files as
> > input,
> > > > and
> > > > > > > > > received
> > > > > > > > > >> an
> > > > > > > > > >> > >> error regarding the GRIB2 field abbreviation
> 'APCP06'
> > > > (see
> > > > > > > output
> > > > > > > > > at
> > > > > > > > > >> > >> bottom
> > > > > > > > > >> > >> of message). Yet when I run degrib on the
fcst and
> > obs
> > > > > files
> > > > > > > > that's
> > > > > > > > > >> > what I
> > > > > > > > > >> > >> get:
> > > > > > > > > >> > >>
> > > > > > > > > >> > >> (g14a2)
/u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > > degrib
> > > > > > > > > >> > >> p06m_2017082312f006.grb2 -I
> > > > > > > > > >> > >> MsgNum, Byte, GRIB-Version, elem, level,
> > > reference(UTC),
> > > > > > > > > valid(UTC),
> > > > > > > > > >> > >> Proj(hr)
> > > > > > > > > >> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation
> > > [kg/(m^2)]",
> > > > > > > 0-SFC,
> > > > > > > > > >> > >> 08/23/2017
> > > > > > > > > >> > >> 12:00, 08/23/2017 18:00, 6.00
> > > > > > > > > >> > >>
> > > > > > > > > >> > >> (g14a2)
/u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > > degrib
> > > > > > > > > >> > >> urma2p5.2017082318 <(201)%20708-
2318>.pcp_06h.
> > 184.grb2
> > > > -I
> > > > > > > > > >> > >> MsgNum, Byte, GRIB-Version, elem, level,
> > > reference(UTC),
> > > > > > > > > valid(UTC),
> > > > > > > > > >> > >> Proj(hr)
> > > > > > > > > >> > >> 1.0, 0, 2, APCP06="06 hr Total precipitation
> > > [kg/(m^2)]",
> > > > > > > 0-SFC,
> > > > > > > > > >> > >> 08/23/2017
> > > > > > > > > >> > >> 12:00, 08/23/2017 18:00, 6.00
> > > > > > > > > >> > >> (g14a2)
/u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > > > > > > > > >> > >>
> > > > > > > > > >> > >> If I put both the fcst and obs files through
> > > pcp_combine
> > > > > and
> > > > > > > use
> > > > > > > > > the
> > > > > > > > > >> > >> resultant .nc files, grid_stat runs fine,
with the
> > > caveat
> > > > > > that
> > > > > > > > the
> > > > > > > > > >> field
> > > > > > > > > >> > >> name becomes 'APCP_06'.
> > > > > > > > > >> > >>
> > > > > > > > > >> > >> Thanks for any comments/suggestions, Brad
Diehl
> > > > > > > > > >> > >>
> > > > > > > > > >> > >> (g14a2)
/u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > > > grid_stat
> > > > > > > > > >> > >> p06m_2017082312f006.grb2 urma2p5.2017082318
> > > > > > <(201)%20708-2318>
> > > > > > > > > >> > >> .pcp_06h.184.grb2
> > > > > > > > > >> > >> config/GridStatConfig_APCP_06_grib -outdir
> > > out/grid_stat
> > > > > -v
> > > > > > > > 4DEBUG
> > > > > > > > > >> 1:
> > > > > > > > > >> > >> Default Config File:
> > > > > > > > > >> > >> /usrx/local/dev/MET/6.0/share/
> > > met/config/GridStatConfig_
> > > > > > > default
> > > > > > > > > >> > >> DEBUG 1: User Config File:
> > > config/GridStatConfig_APCP_06_
> > > > > > grib
> > > > > > > > > >> > >> DEBUG 4: Met2dDataFileFactory::new_met_
> > 2d_data_file()
> > > ->
> > > > > > > created
> > > > > > > > > new
> > > > > > > > > >> > >> Met2dDataFile object of type "FileType_Gb2".
> > > > > > > > > >> > >> WARNING:
> > > > > > > > > >> > >> WARNING: MetGrib2DataFile::read_grib2_grid()
->
> MET
> > > does
> > > > > not
> > > > > > > > > >> currently
> > > > > > > > > >> > >> support Lambert Conformal grids where dx
(2.5397)
> !=
> > dy
> > > > > > > (2.5397)
> > > > > > > > > and
> > > > > > > > > >> may
> > > > > > > > > >> > >> produce unexpected results!
> > > > > > > > > >> > >> WARNING:
> > > > > > > > > >> > >> DEBUG 4:
> > > > > > > > > >> > >> DEBUG 4: Lambert Conformal Grid Data:
> > > > > > > > > >> > >> DEBUG 4: scale_lat_1: 25
> > > > > > > > > >> > >> DEBUG 4: scale_lat_2: 25
> > > > > > > > > >> > >> DEBUG 4: lat_pin: 20.192
> > > > > > > > > >> > >> DEBUG 4: lon_pin: 121.554
> > > > > > > > > >> > >> DEBUG 4: x_pin: 0
> > > > > > > > > >> > >> DEBUG 4: y_pin: 0
> > > > > > > > > >> > >> DEBUG 4: lon_orient: 95
> > > > > > > > > >> > >> DEBUG 4: d_km: 2.5397
> > > > > > > > > >> > >> DEBUG 4: r_km: 6371.2
> > > > > > > > > >> > >> DEBUG 4: nx: 2145
> > > > > > > > > >> > >> DEBUG 4: ny: 1597
> > > > > > > > > >> > >> DEBUG 4:
> > > > > > > > > >> > >> DEBUG 4: Met2dDataFileFactory::new_met_
> > 2d_data_file()
> > > ->
> > > > > > > created
> > > > > > > > > new
> > > > > > > > > >> > >> Met2dDataFile object of type "FileType_Gb2".
> > > > > > > > > >> > >> DEBUG 4:
> > > > > > > > > >> > >> DEBUG 4: Lambert Conformal Grid Data:
> > > > > > > > > >> > >> DEBUG 4: scale_lat_1: 25
> > > > > > > > > >> > >> DEBUG 4: scale_lat_2: 25
> > > > > > > > > >> > >> DEBUG 4: lat_pin: 20.192
> > > > > > > > > >> > >> DEBUG 4: lon_pin: 121.554
> > > > > > > > > >> > >> DEBUG 4: x_pin: 0
> > > > > > > > > >> > >> DEBUG 4: y_pin: 0
> > > > > > > > > >> > >> DEBUG 4: lon_orient: 95
> > > > > > > > > >> > >> DEBUG 4: d_km: 2.5397
> > > > > > > > > >> > >> DEBUG 4: r_km: 6371.2
> > > > > > > > > >> > >> DEBUG 4: nx: 2145
> > > > > > > > > >> > >> DEBUG 4: ny: 1377
> > > > > > > > > >> > >> DEBUG 4:
> > > > > > > > > >> > >> DEBUG 4: VarInfoFactory::new_var_info() ->
created
> > new
> > > > > > VarInfo
> > > > > > > > > >> object of
> > > > > > > > > >> > >> type "FileType_Gb2".
> > > > > > > > > >> > >> DEBUG 4: VarInfoFactory::new_var_info() ->
created
> > new
> > > > > > VarInfo
> > > > > > > > > >> object of
> > > > > > > > > >> > >> type "FileType_Gb2".
> > > > > > > > > >> > >> ERROR :
> > > > > > > > > >> > >> ERROR : VarInfoGrib2::set_dict() ->
unrecognized
> > GRIB2
> > > > > field
> > > > > > > > > >> > abbreviation
> > > > > > > > > >> > >> 'APCP06'
> > > > > > > > > >> > >> ERROR :
> > > > > > > > > >> > >>
> > > > > > > > > >> > >>
> > > > > > > > > >> > >
> > > > > > > > > >> >
> > > > > > > > > >> >
> > > > > > > > > >>
> > > > > > > > > >>
> > > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: grib2 files as input into MET grid_stat
From: Brad Diehl - NOAA Federal
Time: Fri Sep 01 12:46:37 2017
Very good, let me try and set something up with our Task Lead and
we'll be
in touch....Brad
On Fri, Sep 1, 2017 at 6:41 PM, John Halley Gotway via RT
<met_help at ucar.edu
> wrote:
> Brad,
>
> OK, thanks for the info. We would certainly be interested in
working
> closely with you in your development of QPF verification for NWS.
If MET
> proves useful, I wouldn't let the file formats stand in the way.
Using a
> script to reformat the existing MET NetCDF output to follow CF
conventions
> would be pretty straight-forward with the longer term goal of MET
producing
> them directly.
>
> I'll be in DC visiting NCEP 9/19 - 9/22 and would be happy to
discuss any
> issues or questions and/or you could talk directly with the MET
project
> lead, Tara Jensen.
>
> Thanks,
> John
>
> On Fri, Sep 1, 2017 at 12:02 PM, Brad Diehl - NOAA Federal via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> >
> > Thanks for your reply again. As you stated, the primary output is
text
> > files. But we are working on the very early development of the
> Quantitative
> > Precipitation Forecast Verification System for the NWS, and will
want to
> > give the users the ability to look at the difference fields, for
example
> > the ability to see the difference fields geo-referenced, possibly
using
> the
> > netcdf output. So we're in the early stages of assessing how we
are going
> > to do all this.
> >
> > Anyway, thanks for your help earlier, and have a nice long
> weekend....Brad
> >
> > On Fri, Sep 1, 2017 at 4:47 PM, John Halley Gotway via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > > Brad,
> > >
> > > It won't be in the next version, 6.1, and we don't currently
have it
> > > assigned to version 6.2 either. So no specific timeline.
> > >
> > > How does this affect the decisions you're making? In our
application
> of
> > > MET to projects in the Development Testbed Center (DTC), the
primary
> > output
> > > of MET are the ASCII statistics files, to which NetCDF CF
conventions
> do
> > > not apply. We load them into a METViewer database, aggregate
through
> > time,
> > > and make plots to illustrate model performance.
> > >
> > > Which NetCDF files produced by the MET tools are of the most
concern?
> > >
> > > Thanks,
> > > John
> > >
> > >
> > >
> > > On Fri, Sep 1, 2017 at 10:29 AM, Brad Diehl - NOAA Federal via
RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829
>
> > > >
> > > > Is there any kind of time frame to make that switch over to
> compliance?
> > > >
> > > > Brad Diehl
> > > >
> > > > On Fri, Sep 1, 2017 at 3:35 PM, John Halley Gotway via RT <
> > > > met_help at ucar.edu
> > > > > wrote:
> > > >
> > > > > Hi Brad,
> > > > >
> > > > > Unfortunately, no, the NetCDF output generated by the MET
tools
> > > > themselves
> > > > > are not CF-compliant.
> > > > >
> > > > > Since we always suggest that others use that convention, we
should
> > too!
> > > > > And we do plan to make that switch over eventually, but
haven't
> done
> > so
> > > > > yet.
> > > > >
> > > > > Thanks,
> > > > > John
> > > > >
> > > > > On Thu, Aug 31, 2017 at 1:15 PM, Brad Diehl - NOAA Federal
via RT <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > > >
> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81829 >
> > > > > >
> > > > > > Thanks for your help and suggestions today. A question
came up
> from
> > > one
> > > > > of
> > > > > > our developers regarding netcdf and CF compliance. Is the
output
> of
> > > > > MET/6.0
> > > > > > netcdf CF compliant?
> > > > > >
> > > > > > Thanks in advance.......Brad
> > > > > >
> > > > > > On Thu, Aug 31, 2017 at 6:35 PM, Julie Prestopnik via RT <
> > > > > > met_help at ucar.edu>
> > > > > > wrote:
> > > > > >
> > > > > > > Hi Brad,
> > > > > > >
> > > > > > > I just thought I'd add that you can find more
information for
> > > > > file-format
> > > > > > > specific settings for the "field" entry (including
level) for
> MET
> > > on
> > > > > pg.
> > > > > > 53
> > > > > > > of the MET User's Guide for met-6.0:
> > > > > > >
> > > > > > > http://www.dtcenter.org/met/users/docs/users_guide/MET_
> > > > > > > Users_Guide_v6.0.pdf
> > > > > > >
> > > > > > > or in the README file in the met-6.0/data/config area of
your
> > > > install.
> > > > > > >
> > > > > > > For example:
> > > > > > > // File-format specific settings for the "field" entry:
> > > > > > > //
> > > > > > > // - GRIB1 and GRIB2:
> > > > > > > // - For custom GRIB tables, see note about
> > MET_GRIB_TABLES.
> > > > > > > // - The "name" entry specifies a GRIB code number
or
> > > > > abbreviation.
> > > > > > > // - GRIB1 Product Definition Section:
> > > > > > > // http://www.nco.ncep.noaa.gov/
> > > pmb/docs/on388/table2.html
> > > > > > > // - GRIB2 Product Definition Section:
> > > > > > > // http://www.nco.ncep.noaa.gov/
> > > > > pmb/docs/grib2/grib2_doc.shtml
> > > > > > > // - The "level" entry specifies a level type and
value:
> > > > > > > // - ANNN for accumulation interval NNN
> > > > > > > // - ZNNN for vertical level NNN
> > > > > > > // - ZNNN-NNN for a range of vertical levels
> > > > > > > // - PNNN for pressure level NNN in hPa
> > > > > > > // - PNNN-NNN for a range of pressure levels in
hPa
> > > > > > > // - LNNN for a generic level type
> > > > > > > // - RNNN for a specific GRIB record number
> > > > > > >
> > > > > > > I hope that helps!
> > > > > > >
> > > > > > > Julie
> > > > > > >
> > > > > > > On Thu, Aug 31, 2017 at 12:25 PM, Minna Win via RT <
> > > > met_help at ucar.edu>
> > > > > > > wrote:
> > > > > > >
> > > > > > > >
> > > > > > > > <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=81829
> > >
> > > > > > > >
> > > > > > > > Hi Brad,
> > > > > > > >
> > > > > > > > For GRIB2, the level (if not an accumulated precip
type) is
> > > > typically
> > > > > > > "Z0"
> > > > > > > > for surface (sometimes "L0" works, too), but for other
levels
> > > such
> > > > as
> > > > > > > > entire atmosphere, you will need to refer to the GRIB2
table
> > > "Table
> > > > > 2.
> > > > > > > > Parameters and Units":
> > > > > > > > http://www.wmo.int/pages/prog/www/WDM/Guides/Guide-
binary-2.
> > > > > html#GRIB
> > > > > > .
> > > > > > > > From that table, find the value that corresponds to
your
> level
> > > type
> > > > > and
> > > > > > > set
> > > > > > > > it to GRIB_lvl_typ and then set GRIB_lvl_val1 to the
> > appropriate
> > > > > level.
> > > > > > > > For 'entire atmosphere', there is no numerical value
for the
> > > level,
> > > > > so
> > > > > > it
> > > > > > > > would be set to 0 or 255 (you might have to try both
to see
> > what
> > > > > works,
> > > > > > > it
> > > > > > > > depends on the grib file). For something like a
pressure
> > > > difference,
> > > > > > you
> > > > > > > > would need to set GRIB_lvl_typ and the GRIB_lvl_va1
and
> > > > > GRIB_lvl_val2,
> > > > > > > > i.e. GRIB_lvl_val2=top level GRIB_lvl_val1=bottom
level,
> which
> > > you
> > > > > > would
> > > > > > > > obtain from performing a wgrib2/grep combination on
your grib
> > > file
> > > > > and
> > > > > > > > variable of interest. You can use the MET tool
> plot_data_plane
> > to
> > > > > check
> > > > > > > > that your GRIB_lvl_typ and GRIB_lvl_val1 (and
GRIB_lvl_val2
> if
> > > > used)
> > > > > > > > produced a match to the record of interest in the grib
file.
> > > > > > > >
> > > > > > > > As for your question about the (*,*) you see in some
> examples,
> > > > there
> > > > > is
> > > > > > > an
> > > > > > > > explanation in Chapter 3 of the MET User's Manual
> (excerpted):
> > > > > > > >
> > > > > > > >
> > > > > > > > > // - The "level" entry specifies the dimensions to
be used:
> > > > > > > > > // - (i,...,j,*,*) for a single field, where
i,...,j
> > > specifies
> > > > > > fixed
> > > > > > > > > // dimension values and *,* specifies the two
> > dimensions
> > > > for
> > > > > > the
> > > > > > > > > // gridded field.
> > > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > field = [
> > > > > > > > {
> > > > > > > > name = "QVAPOR";
> > > > > > > > level = "( = "(0,5,*,*)"
> > > > > > > >
> > > > > > > > },
> > > > > > > > {
> > > > > > > >
> > > > > > > > name = "TMP_P850_ENS_MEAN;
> > > > > > > > level = [ "(*,*)" ];
> > > > > > > > }
> > > > > > > > ];
> > > > > > > >
> > > > > > > >
> > > > > > > > I hope that clarifies things. Please don't hesitate
to
> contact
> > > us
> > > > if
> > > > > > you
> > > > > > > > have further questions.
> > > > > > > >
> > > > > > > > Regards,
> > > > > > > > Minna
> > > > > > > >
> > > > > > > > ---------------
> > > > > > > > Minna Win
> > > > > > > > NCAR
> > > > > > > > Research Applications Lab
> > > > > > > > Phone: 303-497-8423
> > > > > > > > Fax: 302-497-8401
> > > > > > > >
> > > > > > > >
> > > > > > > > On Thu, Aug 31, 2017 at 4:21 PM, Brad Diehl - NOAA
Federal
> via
> > > RT <
> > > > > > > > met_help at ucar.edu> wrote:
> > > > > > > >
> > > > > > > > >
> > > > > > > > > <URL: https://rt.rap.ucar.edu/rt/
> > Ticket/Display.html?id=81829
> > > >
> > > > > > > > >
> > > > > > > > > the A06 worked. Prior to trying that I sent you via
mailx
> the
> > > > > > fcst/obs
> > > > > > > > > grib2 files, you can ignore that now.
> > > > > > > > >
> > > > > > > > > But it does beg the question, how does one find the
level
> > > names,
> > > > is
> > > > > > it
> > > > > > > in
> > > > > > > > > the grib file itself? Very confusing, at least to
me. And
> > what
> > > > does
> > > > > > *,*
> > > > > > > > > mean when I see that in some examples?
> > > > > > > > >
> > > > > > > > > Thanks for all your help........Brad
> > > > > > > > >
> > > > > > > > > On Thu, Aug 31, 2017 at 4:11 PM, Minna Win via RT <
> > > > > met_help at ucar.edu
> > > > > > >
> > > > > > > > > wrote:
> > > > > > > > >
> > > > > > > > > > Hi Brad,
> > > > > > > > > >
> > > > > > > > > > Tara Jensen brought this to my attention, for
verify
> APCP,
> > > you
> > > > > need
> > > > > > > to
> > > > > > > > > set
> > > > > > > > > > the level to the value pre-pended with the letter
'A'.
> For
> > > > > > instance,
> > > > > > > > > for 6
> > > > > > > > > > hour accumulation level, set level="A06".
> > > > > > > > > >
> > > > > > > > > > Thanks,
> > > > > > > > > > Minna
> > > > > > > > > >
> > > > > > > > > > ---------------
> > > > > > > > > > Minna Win
> > > > > > > > > > NCAR
> > > > > > > > > > Research Applications Lab
> > > > > > > > > > Phone: 303-497-8423
> > > > > > > > > > Fax: 302-497-8401
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > > On Thu, Aug 31, 2017 at 3:52 PM, Minna Win-
Gildenmeister
> <
> > > > > > > > > > minnawin at ucar.edu>
> > > > > > > > > > wrote:
> > > > > > > > > >
> > > > > > > > > > > Hi Brad,
> > > > > > > > > > >
> > > > > > > > > > > Can you try "L0", and "P0"? If that doesn't
work,
> could
> > > you
> > > > > > please
> > > > > > > > > send
> > > > > > > > > > > me your grib file? Instructions are as follows:
> > > > > > > > > > >
> > > > > > > > > > > Non-Windows FTP client:
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > > * ftp -p ftp.rap.ucar.edu
<http://ftp.rap.ucar.edu>
> > > > > > > > username=anonymous
> > > > > > > > > > > password="your email address" cd
> > incoming/irap/met_help
> > > > > mkdir
> > > > > > > > > > > "surname_data" # where surname is your last
name cd
> > > > > > > "surname_data"
> > > > > > > > > put
> > > > > > > > > > > "your files" bye *
> > > > > > > > > > >
> > > > > > > > > > > Windows FTP client:
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > > * ftp ftp.rap.ucar.edu <http://ftp.rap.ucar.edu>
> > > > > > > username=anonymous
> > > > > > > > > > > password="your email address" quote pasv cd
> > > > > > incoming/irap/met_help
> > > > > > > > > mkdir
> > > > > > > > > > > "surname_data" # where surname is your last name
cd
> > > > > > "surname_data"
> > > > > > > > put
> > > > > > > > > > > "your files" bye*
> > > > > > > > > > >
> > > > > > > > > > > I can take a look at the level type and value,
you
> might
> > > need
> > > > > to
> > > > > > > use
> > > > > > > > > > > special syntax in your config file to get
grid_stat to
> > > work.
> > > > > Let
> > > > > > > me
> > > > > > > > > know
> > > > > > > > > > > if you ftp'd a grib file so I can grab it before
the
> file
> > > > gets
> > > > > > > > > scrubbed.
> > > > > > > > > > >
> > > > > > > > > > > Thanks,
> > > > > > > > > > > Minna
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > > ---------------
> > > > > > > > > > > Minna Win
> > > > > > > > > > > NCAR
> > > > > > > > > > > Research Applications Lab
> > > > > > > > > > > Phone: 303-497-8423 <(303)%20497-8423>
> > > > > > > > > > > Fax: 302-497-8401 <(302)%20497-8401>
> > > > > > > > > > >
> > > > > > > > > > >
> > > > > > > > > > > On Thu, Aug 31, 2017 at 3:09 PM, Brad Diehl -
NOAA
> > Federal
> > > > via
> > > > > > RT <
> > > > > > > > > > > met_help at ucar.edu> wrote:
> > > > > > > > > > >
> > > > > > > > > > >>
> > > > > > > > > > >> <URL: https://rt.rap.ucar.edu/rt/
> > > > Ticket/Display.html?id=81829
> > > > > >
> > > > > > > > > > >>
> > > > > > > > > > >> After running pcp_combine, everything did work
> > correctly,
> > > it
> > > > > is
> > > > > > > the
> > > > > > > > > > grib2
> > > > > > > > > > >> fcts/obs part I'm having trouble with. I think
you put
> > me
> > > on
> > > > > the
> > > > > > > > right
> > > > > > > > > > >> track by suggesting wgrib2. Here's what I get
using
> > wgrib2
> > > > on
> > > > > > the
> > > > > > > > fcst
> > > > > > > > > > >> file:
> > > > > > > > > > >>
> > > > > > > > > > >> (g14a2)
/u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > wgrib2
> > > > > > > > > > >> p06m_2017082312f006.grb2
> > > > > > > > > > >> 1:21:d=2017082312 <(201)%20708-
2312>:APCP:surface -
> > > > > > > > surface:CodeTable
> > > > > > > > > > >> 4.11=missing:
> > > > > > > > > > >>
> > > > > > > > > > >> in the config file I'm using:
> > > > > > > > > > >>
> > > > > > > > > > >>
> > > > > > > > > > >> fcst = {
> > > > > > > > > > >>
> > > > > > > > > > >> field = [
> > > > > > > > > > >> {
> > > > > > > > > > >> name = "APCP";
> > > > > > > > > > >> level = "Z0" ;
> > > > > > > > > > >> cat_thresh = [ >0.0, >=5.0 , >=10.0];
> > > > > > > > > > >> }
> > > > > > > > > > >> ];
> > > > > > > > > > >>
> > > > > > > > > > >> for level, should I be using surface or Z0?,
neither
> of
> > > > which
> > > > > > > seemed
> > > > > > > > > to
> > > > > > > > > > >> work.
> > > > > > > > > > >>
> > > > > > > > > > >> On Thu, Aug 31, 2017 at 2:58 PM, Minna Win via
RT <
> > > > > > > > met_help at ucar.edu>
> > > > > > > > > > >> wrote:
> > > > > > > > > > >>
> > > > > > > > > > >> > Hi Brad,
> > > > > > > > > > >> >
> > > > > > > > > > >> > Please disregard my previous response. I
didn't see
> > > that
> > > > > you
> > > > > > > had
> > > > > > > > > run
> > > > > > > > > > >> > pcp_combine-that *should* have worked. I'll
look a
> > > little
> > > > > > > further
> > > > > > > > > or
> > > > > > > > > > >> see
> > > > > > > > > > >> > if someone else has seen this before.
> > > > > > > > > > >> >
> > > > > > > > > > >> >
> > > > > > > > > > >> >
> > > > > > > > > > >> > Regards,
> > > > > > > > > > >> > Minna
> > > > > > > > > > >> >
> > > > > > > > > > >> > ---------------
> > > > > > > > > > >> > Minna Win
> > > > > > > > > > >> > NCAR
> > > > > > > > > > >> > Research Applications Lab
> > > > > > > > > > >> > Phone: 303-497-8423
> > > > > > > > > > >> > Fax: 302-497-8401
> > > > > > > > > > >> >
> > > > > > > > > > >> >
> > > > > > > > > > >> > On Thu, Aug 31, 2017 at 2:50 PM, Minna
> > > Win-Gildenmeister <
> > > > > > > > > > >> > minnawin at ucar.edu>
> > > > > > > > > > >> > wrote:
> > > > > > > > > > >> >
> > > > > > > > > > >> > > HI Brad,
> > > > > > > > > > >> > >
> > > > > > > > > > >> > > wgrib is an NDFD tool and I don't know how
it is
> > > > > > interpreting
> > > > > > > > the
> > > > > > > > > > APCP
> > > > > > > > > > >> > > variable in your grib file. Try using
wgrib2 on
> your
> > > > grib
> > > > > > file
> > > > > > > > and
> > > > > > > > > > >> grep
> > > > > > > > > > >> > for
> > > > > > > > > > >> > > APCP and use that variable abbreviation in
you MET
> > > > > Grid_stat
> > > > > > > > > config
> > > > > > > > > > >> > file. I
> > > > > > > > > > >> > > hope that helps. Please let us know if you
still
> > have
> > > > > > > problems
> > > > > > > > > > using
> > > > > > > > > > >> > > Grid_stat.
> > > > > > > > > > >> > >
> > > > > > > > > > >> > > Regards,
> > > > > > > > > > >> > > Minna
> > > > > > > > > > >> > >
> > > > > > > > > > >> > > ---------------
> > > > > > > > > > >> > > Minna Win
> > > > > > > > > > >> > > NCAR
> > > > > > > > > > >> > > Research Applications Lab
> > > > > > > > > > >> > > Phone: 303-497-8423 <(303)%20497-8423>
> > > > > > > > > > >> > > Fax: 302-497-8401 <(302)%20497-8401>
> > > > > > > > > > >> > >
> > > > > > > > > > >> > >
> > > > > > > > > > >> > > On Thu, Aug 31, 2017 at 2:33 PM, Brad Diehl
- NOAA
> > > > Federal
> > > > > > via
> > > > > > > > RT
> > > > > > > > > <
> > > > > > > > > > >> > > met_help at ucar.edu> wrote:
> > > > > > > > > > >> > >
> > > > > > > > > > >> > >>
> > > > > > > > > > >> > >> Thu Aug 31 08:33:00 2017: Request 81829
was acted
> > > upon.
> > > > > > > > > > >> > >> Transaction: Ticket created by
> brad.diehl at noaa.gov
> > > > > > > > > > >> > >> Queue: met_help
> > > > > > > > > > >> > >> Subject: grib2 files as input into
MET
> > grid_stat
> > > > > > > > > > >> > >> Owner: Nobody
> > > > > > > > > > >> > >> Requestors: brad.diehl at noaa.gov
> > > > > > > > > > >> > >> Status: new
> > > > > > > > > > >> > >> Ticket <URL:
https://rt.rap.ucar.edu/rt/Tic
> > > > > > > > > > >> ket/Display.html?id=81829 >
> > > > > > > > > > >> > >>
> > > > > > > > > > >> > >>
> > > > > > > > > > >> > >> I tried to use grid_stat with two grib2
files as
> > > input,
> > > > > and
> > > > > > > > > > received
> > > > > > > > > > >> an
> > > > > > > > > > >> > >> error regarding the GRIB2 field
abbreviation
> > 'APCP06'
> > > > > (see
> > > > > > > > output
> > > > > > > > > > at
> > > > > > > > > > >> > >> bottom
> > > > > > > > > > >> > >> of message). Yet when I run degrib on the
fcst
> and
> > > obs
> > > > > > files
> > > > > > > > > that's
> > > > > > > > > > >> > what I
> > > > > > > > > > >> > >> get:
> > > > > > > > > > >> > >>
> > > > > > > > > > >> > >> (g14a2)
/u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > > > degrib
> > > > > > > > > > >> > >> p06m_2017082312f006.grb2 -I
> > > > > > > > > > >> > >> MsgNum, Byte, GRIB-Version, elem, level,
> > > > reference(UTC),
> > > > > > > > > > valid(UTC),
> > > > > > > > > > >> > >> Proj(hr)
> > > > > > > > > > >> > >> 1.0, 0, 2, APCP06="06 hr Total
precipitation
> > > > [kg/(m^2)]",
> > > > > > > > 0-SFC,
> > > > > > > > > > >> > >> 08/23/2017
> > > > > > > > > > >> > >> 12:00, 08/23/2017 18:00, 6.00
> > > > > > > > > > >> > >>
> > > > > > > > > > >> > >> (g14a2)
/u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > > > degrib
> > > > > > > > > > >> > >> urma2p5.2017082318 <(201)%20708-
2318>.pcp_06h.
> > > 184.grb2
> > > > > -I
> > > > > > > > > > >> > >> MsgNum, Byte, GRIB-Version, elem, level,
> > > > reference(UTC),
> > > > > > > > > > valid(UTC),
> > > > > > > > > > >> > >> Proj(hr)
> > > > > > > > > > >> > >> 1.0, 0, 2, APCP06="06 hr Total
precipitation
> > > > [kg/(m^2)]",
> > > > > > > > 0-SFC,
> > > > > > > > > > >> > >> 08/23/2017
> > > > > > > > > > >> > >> 12:00, 08/23/2017 18:00, 6.00
> > > > > > > > > > >> > >> (g14a2)
/u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > > > > > > > > > >> > >>
> > > > > > > > > > >> > >> If I put both the fcst and obs files
through
> > > > pcp_combine
> > > > > > and
> > > > > > > > use
> > > > > > > > > > the
> > > > > > > > > > >> > >> resultant .nc files, grid_stat runs fine,
with
> the
> > > > caveat
> > > > > > > that
> > > > > > > > > the
> > > > > > > > > > >> field
> > > > > > > > > > >> > >> name becomes 'APCP_06'.
> > > > > > > > > > >> > >>
> > > > > > > > > > >> > >> Thanks for any comments/suggestions, Brad
Diehl
> > > > > > > > > > >> > >>
> > > > > > > > > > >> > >> (g14a2)
/u/Brad.Diehl/MET_Tutorial_6.0/tutorial$
> > > > > grid_stat
> > > > > > > > > > >> > >> p06m_2017082312f006.grb2
urma2p5.2017082318
> > > > > > > <(201)%20708-2318>
> > > > > > > > > > >> > >> .pcp_06h.184.grb2
> > > > > > > > > > >> > >> config/GridStatConfig_APCP_06_grib -outdir
> > > > out/grid_stat
> > > > > > -v
> > > > > > > > > 4DEBUG
> > > > > > > > > > >> 1:
> > > > > > > > > > >> > >> Default Config File:
> > > > > > > > > > >> > >> /usrx/local/dev/MET/6.0/share/
> > > > met/config/GridStatConfig_
> > > > > > > > default
> > > > > > > > > > >> > >> DEBUG 1: User Config File:
> > > > config/GridStatConfig_APCP_06_
> > > > > > > grib
> > > > > > > > > > >> > >> DEBUG 4: Met2dDataFileFactory::new_met_
> > > 2d_data_file()
> > > > ->
> > > > > > > > created
> > > > > > > > > > new
> > > > > > > > > > >> > >> Met2dDataFile object of type
"FileType_Gb2".
> > > > > > > > > > >> > >> WARNING:
> > > > > > > > > > >> > >> WARNING:
MetGrib2DataFile::read_grib2_grid() ->
> > MET
> > > > does
> > > > > > not
> > > > > > > > > > >> currently
> > > > > > > > > > >> > >> support Lambert Conformal grids where dx
(2.5397)
> > !=
> > > dy
> > > > > > > > (2.5397)
> > > > > > > > > > and
> > > > > > > > > > >> may
> > > > > > > > > > >> > >> produce unexpected results!
> > > > > > > > > > >> > >> WARNING:
> > > > > > > > > > >> > >> DEBUG 4:
> > > > > > > > > > >> > >> DEBUG 4: Lambert Conformal Grid Data:
> > > > > > > > > > >> > >> DEBUG 4: scale_lat_1: 25
> > > > > > > > > > >> > >> DEBUG 4: scale_lat_2: 25
> > > > > > > > > > >> > >> DEBUG 4: lat_pin: 20.192
> > > > > > > > > > >> > >> DEBUG 4: lon_pin: 121.554
> > > > > > > > > > >> > >> DEBUG 4: x_pin: 0
> > > > > > > > > > >> > >> DEBUG 4: y_pin: 0
> > > > > > > > > > >> > >> DEBUG 4: lon_orient: 95
> > > > > > > > > > >> > >> DEBUG 4: d_km: 2.5397
> > > > > > > > > > >> > >> DEBUG 4: r_km: 6371.2
> > > > > > > > > > >> > >> DEBUG 4: nx: 2145
> > > > > > > > > > >> > >> DEBUG 4: ny: 1597
> > > > > > > > > > >> > >> DEBUG 4:
> > > > > > > > > > >> > >> DEBUG 4: Met2dDataFileFactory::new_met_
> > > 2d_data_file()
> > > > ->
> > > > > > > > created
> > > > > > > > > > new
> > > > > > > > > > >> > >> Met2dDataFile object of type
"FileType_Gb2".
> > > > > > > > > > >> > >> DEBUG 4:
> > > > > > > > > > >> > >> DEBUG 4: Lambert Conformal Grid Data:
> > > > > > > > > > >> > >> DEBUG 4: scale_lat_1: 25
> > > > > > > > > > >> > >> DEBUG 4: scale_lat_2: 25
> > > > > > > > > > >> > >> DEBUG 4: lat_pin: 20.192
> > > > > > > > > > >> > >> DEBUG 4: lon_pin: 121.554
> > > > > > > > > > >> > >> DEBUG 4: x_pin: 0
> > > > > > > > > > >> > >> DEBUG 4: y_pin: 0
> > > > > > > > > > >> > >> DEBUG 4: lon_orient: 95
> > > > > > > > > > >> > >> DEBUG 4: d_km: 2.5397
> > > > > > > > > > >> > >> DEBUG 4: r_km: 6371.2
> > > > > > > > > > >> > >> DEBUG 4: nx: 2145
> > > > > > > > > > >> > >> DEBUG 4: ny: 1377
> > > > > > > > > > >> > >> DEBUG 4:
> > > > > > > > > > >> > >> DEBUG 4: VarInfoFactory::new_var_info() ->
> created
> > > new
> > > > > > > VarInfo
> > > > > > > > > > >> object of
> > > > > > > > > > >> > >> type "FileType_Gb2".
> > > > > > > > > > >> > >> DEBUG 4: VarInfoFactory::new_var_info() ->
> created
> > > new
> > > > > > > VarInfo
> > > > > > > > > > >> object of
> > > > > > > > > > >> > >> type "FileType_Gb2".
> > > > > > > > > > >> > >> ERROR :
> > > > > > > > > > >> > >> ERROR : VarInfoGrib2::set_dict() ->
unrecognized
> > > GRIB2
> > > > > > field
> > > > > > > > > > >> > abbreviation
> > > > > > > > > > >> > >> 'APCP06'
> > > > > > > > > > >> > >> ERROR :
> > > > > > > > > > >> > >>
> > > > > > > > > > >> > >>
> > > > > > > > > > >> > >
> > > > > > > > > > >> >
> > > > > > > > > > >> >
> > > > > > > > > > >>
> > > > > > > > > > >>
> > > > > > > > > > >
> > > > > > > > > >
> > > > > > > > > >
> > > > > > > > >
> > > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
More information about the Met_help
mailing list