[Met_help] [rt.rap.ucar.edu #80881] History for how set set up the level parameters in the configure file for Met gridstate
John Halley Gotway via RT
met_help at ucar.edu
Fri Jun 16 11:19:48 MDT 2017
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Dear Mr. or Mrs.,
I'm running Met grid_stat,
I have two kinds of files (fcstfile and obsfile) respectively, as shown in
the following:
1)
^^
wgrib2 -v pcp_03hr_fcst03hrpcp_hi_imp_2015103018.gb2
1:0:d=2015103018:WEASD Water Equivalent of Accumulated Snow Depth
[kg/m^2]:surface:3 hour fcst:
2:105540:d=2015103018:APCP Total Precipitation [kg/m^2]:surface:2-3 hour
acc fcst:
3:509588:d=2015103018:WEASD Water Equivalent of Accumulated Snow Depth
[kg/m^2]:surface:2-3 hour acc fcst:
VV
2)
^^
wgrib2 GaugeOnly_QPE_01H_00.00_20151208-080000.grib2
1:0:d=2015120808:var discipline=209 master_table=255 parmcat=6 parm=16:0 m
above mean sea level:anl:
VV
What should I set as for the level parameter in the configure file in the
following:
^^
fcst = {
wind_thresh = [ NA ];
field = [
{
name = "APCP";
level = [ "" ];
cat_thresh = [ >0.0, >0.254, >0.508, >1.270, >2.540, >3.810,
>6.350, >8.890, >12.700, >25.400 ];
}
];
};
V
I tried surface for the first part, however, failed. I also have no clue
what to fill in for the second file.
Your help is appreciated.
Best Wishes,
Ting
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: how set set up the level parameters in the configure file for Met gridstate
From: John Halley Gotway
Time: Thu Jun 15 14:08:40 2017
Ting,
It looks like you're trying to verify accumulated precipitation.
Here's
the important line from the output of wgrib2 that you sent:
2:105540:d=2015103018:APCP Total Precipitation [kg/m^2]:surface:2-3
hour
acc fcst:
Notices that this is an accumulation between hours 2 and 3. So that's
1
hour of accumulation.
Please try running MET's plot_data_plane utility to make sure MET is
reading your data correctly. Set the level to "A1" to indicate that
you
want a 1-hourly accumulation interval:
met-6.0/bin/plot_data_plane \
pcp_03hr_fcst03hrpcp_hi_imp_2015103018.gb2 \
apcp_01.ps \
'name="APCP"; level="A1";'
That generates a PostScript plot which you can view using a PostScript
viewer, like gv or ghostview. Alternatively, you can convert it to a
png
and view that:
convert -rotate 90 -background white apcp_01.ps apcp_01.png
Does your data look like what you'd expect?
Thanks,
John
On Thu, Jun 15, 2017 at 1:58 PM, Ting Lei - NOAA Affiliate via RT <
met_help at ucar.edu> wrote:
>
> Thu Jun 15 13:58:23 2017: Request 80881 was acted upon.
> Transaction: Ticket created by ting.lei at noaa.gov
> Queue: met_help
> Subject: how set set up the level parameters in the configure
file
> for Met gridstate
> Owner: Nobody
> Requestors: ting.lei at noaa.gov
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=80881 >
>
>
> Dear Mr. or Mrs.,
> I'm running Met grid_stat,
> I have two kinds of files (fcstfile and obsfile) respectively, as
shown in
> the following:
> 1)
> ^^
>
> wgrib2 -v pcp_03hr_fcst03hrpcp_hi_imp_2015103018.gb2
> 1:0:d=2015103018:WEASD Water Equivalent of Accumulated Snow Depth
> [kg/m^2]:surface:3 hour fcst:
> 2:105540:d=2015103018:APCP Total Precipitation [kg/m^2]:surface:2-3
hour
> acc fcst:
> 3:509588:d=2015103018:WEASD Water Equivalent of Accumulated Snow
Depth
> [kg/m^2]:surface:2-3 hour acc fcst:
>
> VV
> 2)
> ^^
> wgrib2 GaugeOnly_QPE_01H_00.00_20151208-080000.grib2
> 1:0:d=2015120808:var discipline=209 master_table=255 parmcat=6
parm=16:0 m
> above mean sea level:anl:
>
> VV
> What should I set as for the level parameter in the configure file
in the
> following:
> ^^
> fcst = {
> wind_thresh = [ NA ];
>
> field = [
> {
> name = "APCP";
> level = [ "" ];
> cat_thresh = [ >0.0, >0.254, >0.508, >1.270, >2.540, >3.810,
> >6.350, >8.890, >12.700, >25.400 ];
> }
> ];
>
> };
>
>
> V
> I tried surface for the first part, however, failed. I also have no
clue
> what to fill in for the second file.
>
> Your help is appreciated.
> Best Wishes,
> Ting
>
>
------------------------------------------------
Subject: how set set up the level parameters in the configure file for Met gridstate
From: Ting Lei - NOAA Affiliate
Time: Thu Jun 15 14:41:15 2017
John,
Thank you so much. I will follow your suggestion.
How about the second kind of data (obs),
^^
wgrib2 GaugeOnly_QPE_01H_00.00_20151208-080000.grib2
1:0:d=2015120808 <(201)%20512-0808>:var discipline=209
master_table=255
parmcat=6 parm=16:0 m above mean sea level:anl:
VV
using g2ctl.pl to generate the grads description file, I got
^^
var209616_0m 0,102,0 209,6,16 ** 0 m above mean sea level desc
[unit]
VV
For this one, how to let MET treat it correctly.
Your further help is appreciated.
Best,
Ting
On Thu, Jun 15, 2017 at 4:08 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Ting,
>
> It looks like you're trying to verify accumulated precipitation.
Here's
> the important line from the output of wgrib2 that you sent:
>
> 2:105540:d=2015103018:APCP Total Precipitation [kg/m^2]:surface:2-3
hour
> acc fcst:
>
> Notices that this is an accumulation between hours 2 and 3. So
that's 1
> hour of accumulation.
>
> Please try running MET's plot_data_plane utility to make sure MET is
> reading your data correctly. Set the level to "A1" to indicate that
you
> want a 1-hourly accumulation interval:
>
> met-6.0/bin/plot_data_plane \
> pcp_03hr_fcst03hrpcp_hi_imp_2015103018.gb2 \
> apcp_01.ps \
> 'name="APCP"; level="A1";'
>
> That generates a PostScript plot which you can view using a
PostScript
> viewer, like gv or ghostview. Alternatively, you can convert it to
a png
> and view that:
>
> convert -rotate 90 -background white apcp_01.ps apcp_01.png
>
> Does your data look like what you'd expect?
>
> Thanks,
> John
>
> On Thu, Jun 15, 2017 at 1:58 PM, Ting Lei - NOAA Affiliate via RT <
> met_help at ucar.edu> wrote:
>
> >
> > Thu Jun 15 13:58:23 2017: Request 80881 was acted upon.
> > Transaction: Ticket created by ting.lei at noaa.gov
> > Queue: met_help
> > Subject: how set set up the level parameters in the configure
file
> > for Met gridstate
> > Owner: Nobody
> > Requestors: ting.lei at noaa.gov
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=80881 >
> >
> >
> > Dear Mr. or Mrs.,
> > I'm running Met grid_stat,
> > I have two kinds of files (fcstfile and obsfile) respectively, as
shown
> in
> > the following:
> > 1)
> > ^^
> >
> > wgrib2 -v pcp_03hr_fcst03hrpcp_hi_imp_2015103018.gb2
> > 1:0:d=2015103018:WEASD Water Equivalent of Accumulated Snow Depth
> > [kg/m^2]:surface:3 hour fcst:
> > 2:105540:d=2015103018:APCP Total Precipitation [kg/m^2]:surface:2-
3 hour
> > acc fcst:
> > 3:509588:d=2015103018:WEASD Water Equivalent of Accumulated Snow
Depth
> > [kg/m^2]:surface:2-3 hour acc fcst:
> >
> > VV
> > 2)
> > ^^
> > wgrib2 GaugeOnly_QPE_01H_00.00_20151208-080000.grib2
> > 1:0:d=2015120808:var discipline=209 master_table=255 parmcat=6
parm=16:0
> m
> > above mean sea level:anl:
> >
> > VV
> > What should I set as for the level parameter in the configure file
in the
> > following:
> > ^^
> > fcst = {
> > wind_thresh = [ NA ];
> >
> > field = [
> > {
> > name = "APCP";
> > level = [ "" ];
> > cat_thresh = [ >0.0, >0.254, >0.508, >1.270, >2.540,
>3.810,
> > >6.350, >8.890, >12.700, >25.400 ];
> > }
> > ];
> >
> > };
> >
> >
> > V
> > I tried surface for the first part, however, failed. I also have
no clue
> > what to fill in for the second file.
> >
> > Your help is appreciated.
> > Best Wishes,
> > Ting
> >
> >
>
>
------------------------------------------------
Subject: how set set up the level parameters in the configure file for Met gridstate
From: John Halley Gotway
Time: Thu Jun 15 14:59:42 2017
Ting,
wgrib2 doesn't know how to handle this data. That's why it listed the
values for the discipline, master_table, parmcat, and parm instead of
listing the corresponding abbreviation:
wgrib2 GaugeOnly_QPE_01H_00.00_20151208-080000.grib2
1:0:d=2015120808:var discipline=209 master_table=255 parmcat=6
parm=16:0 m
above mean sea level:anl:
This combination of settings doesn't appear in it's GRIB tables... and
likely won't appear in the ones that MET uses either.
The quick solution is that you could just tell MET to look for first
record
number:
met-6.0/bin/plot_data_plane \
GaugeOnly_QPE_01H_00.00_20151208-080000.grib2 \
gauge.ps \
'name="APCP"; level="R1";'
Setting level to "R1" tells MET to use the first record regardless of
what
it actually contains!
Alternatively, you could define a custom user GRIB table to define
what
you'd like to call this data. Read about MET_GRIB_TABLES on page 47
of the
user's guide:
http://www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v6.0.pdf
Thanks,
John
On Thu, Jun 15, 2017 at 2:41 PM, Ting Lei - NOAA Affiliate via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=80881 >
>
> John,
> Thank you so much. I will follow your suggestion.
> How about the second kind of data (obs),
> ^^
> wgrib2 GaugeOnly_QPE_01H_00.00_20151208-080000.grib2
> 1:0:d=2015120808 <(201)%20512-0808>:var discipline=209
master_table=255
> parmcat=6 parm=16:0 m above mean sea level:anl:
> VV
> using g2ctl.pl to generate the grads description file, I got
> ^^
> var209616_0m 0,102,0 209,6,16 ** 0 m above mean sea level desc
[unit]
> VV
> For this one, how to let MET treat it correctly.
> Your further help is appreciated.
> Best,
> Ting
>
> On Thu, Jun 15, 2017 at 4:08 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
> > Ting,
> >
> > It looks like you're trying to verify accumulated precipitation.
Here's
> > the important line from the output of wgrib2 that you sent:
> >
> > 2:105540:d=2015103018:APCP Total Precipitation [kg/m^2]:surface:2-
3 hour
> > acc fcst:
> >
> > Notices that this is an accumulation between hours 2 and 3. So
that's 1
> > hour of accumulation.
> >
> > Please try running MET's plot_data_plane utility to make sure MET
is
> > reading your data correctly. Set the level to "A1" to indicate
that you
> > want a 1-hourly accumulation interval:
> >
> > met-6.0/bin/plot_data_plane \
> > pcp_03hr_fcst03hrpcp_hi_imp_2015103018.gb2 \
> > apcp_01.ps \
> > 'name="APCP"; level="A1";'
> >
> > That generates a PostScript plot which you can view using a
PostScript
> > viewer, like gv or ghostview. Alternatively, you can convert it
to a png
> > and view that:
> >
> > convert -rotate 90 -background white apcp_01.ps apcp_01.png
> >
> > Does your data look like what you'd expect?
> >
> > Thanks,
> > John
> >
> > On Thu, Jun 15, 2017 at 1:58 PM, Ting Lei - NOAA Affiliate via RT
<
> > met_help at ucar.edu> wrote:
> >
> > >
> > > Thu Jun 15 13:58:23 2017: Request 80881 was acted upon.
> > > Transaction: Ticket created by ting.lei at noaa.gov
> > > Queue: met_help
> > > Subject: how set set up the level parameters in the
configure file
> > > for Met gridstate
> > > Owner: Nobody
> > > Requestors: ting.lei at noaa.gov
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=80881
> >
> > >
> > >
> > > Dear Mr. or Mrs.,
> > > I'm running Met grid_stat,
> > > I have two kinds of files (fcstfile and obsfile) respectively,
as shown
> > in
> > > the following:
> > > 1)
> > > ^^
> > >
> > > wgrib2 -v pcp_03hr_fcst03hrpcp_hi_imp_2015103018.gb2
> > > 1:0:d=2015103018:WEASD Water Equivalent of Accumulated Snow
Depth
> > > [kg/m^2]:surface:3 hour fcst:
> > > 2:105540:d=2015103018:APCP Total Precipitation
[kg/m^2]:surface:2-3
> hour
> > > acc fcst:
> > > 3:509588:d=2015103018:WEASD Water Equivalent of Accumulated Snow
Depth
> > > [kg/m^2]:surface:2-3 hour acc fcst:
> > >
> > > VV
> > > 2)
> > > ^^
> > > wgrib2 GaugeOnly_QPE_01H_00.00_20151208-080000.grib2
> > > 1:0:d=2015120808:var discipline=209 master_table=255 parmcat=6
> parm=16:0
> > m
> > > above mean sea level:anl:
> > >
> > > VV
> > > What should I set as for the level parameter in the configure
file in
> the
> > > following:
> > > ^^
> > > fcst = {
> > > wind_thresh = [ NA ];
> > >
> > > field = [
> > > {
> > > name = "APCP";
> > > level = [ "" ];
> > > cat_thresh = [ >0.0, >0.254, >0.508, >1.270, >2.540,
>3.810,
> > > >6.350, >8.890, >12.700, >25.400 ];
> > > }
> > > ];
> > >
> > > };
> > >
> > >
> > > V
> > > I tried surface for the first part, however, failed. I also have
no
> clue
> > > what to fill in for the second file.
> > >
> > > Your help is appreciated.
> > > Best Wishes,
> > > Ting
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: how set set up the level parameters in the configure file for Met gridstate
From: Ting Lei - NOAA Affiliate
Time: Fri Jun 16 09:16:09 2017
Hi, John,
So many thanks.
I chose the quick solution you gave and it works seemingly. That is
great!.
I will dig deeper on the more formal solution you suggested later.
Just a problem on which your clarification is appreciated.
In my case, the fcst file and obs file are on different domains of
different sizes.
How can I specify which domain is used as the reference domain. Then,
for
example, the total obs number (relevant in different statistics) will
be
defined as the obs number in the reference domain.
Your further help is appreciated.
Best Wishes,
Ting
On Thu, Jun 15, 2017 at 4:59 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Ting,
>
> wgrib2 doesn't know how to handle this data. That's why it listed
the
> values for the discipline, master_table, parmcat, and parm instead
of
> listing the corresponding abbreviation:
>
> wgrib2 GaugeOnly_QPE_01H_00.00_20151208-080000.grib2
> 1:0:d=2015120808:var discipline=209 master_table=255 parmcat=6
parm=16:0 m
> above mean sea level:anl:
>
> This combination of settings doesn't appear in it's GRIB tables...
and
> likely won't appear in the ones that MET uses either.
>
> The quick solution is that you could just tell MET to look for first
record
> number:
>
> met-6.0/bin/plot_data_plane \
> GaugeOnly_QPE_01H_00.00_20151208-080000.grib2 \
> gauge.ps \
> 'name="APCP"; level="R1";'
>
> Setting level to "R1" tells MET to use the first record regardless
of what
> it actually contains!
>
> Alternatively, you could define a custom user GRIB table to define
what
> you'd like to call this data. Read about MET_GRIB_TABLES on page 47
of the
> user's guide:
>
> http://www.dtcenter.org/met/users/docs/users_guide/MET_
> Users_Guide_v6.0.pdf
>
> Thanks,
> John
>
> On Thu, Jun 15, 2017 at 2:41 PM, Ting Lei - NOAA Affiliate via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=80881 >
> >
> > John,
> > Thank you so much. I will follow your suggestion.
> > How about the second kind of data (obs),
> > ^^
> > wgrib2 GaugeOnly_QPE_01H_00.00_20151208-080000.grib2
> > 1:0:d=2015120808 <(201)%20512-0808>:var discipline=209
master_table=255
> > parmcat=6 parm=16:0 m above mean sea level:anl:
> > VV
> > using g2ctl.pl to generate the grads description file, I got
> > ^^
> > var209616_0m 0,102,0 209,6,16 ** 0 m above mean sea level desc
[unit]
> > VV
> > For this one, how to let MET treat it correctly.
> > Your further help is appreciated.
> > Best,
> > Ting
> >
> > On Thu, Jun 15, 2017 at 4:08 PM, John Halley Gotway via RT <
> > met_help at ucar.edu> wrote:
> >
> > > Ting,
> > >
> > > It looks like you're trying to verify accumulated precipitation.
> Here's
> > > the important line from the output of wgrib2 that you sent:
> > >
> > > 2:105540:d=2015103018:APCP Total Precipitation
[kg/m^2]:surface:2-3
> hour
> > > acc fcst:
> > >
> > > Notices that this is an accumulation between hours 2 and 3. So
that's
> 1
> > > hour of accumulation.
> > >
> > > Please try running MET's plot_data_plane utility to make sure
MET is
> > > reading your data correctly. Set the level to "A1" to indicate
that
> you
> > > want a 1-hourly accumulation interval:
> > >
> > > met-6.0/bin/plot_data_plane \
> > > pcp_03hr_fcst03hrpcp_hi_imp_2015103018.gb2 \
> > > apcp_01.ps \
> > > 'name="APCP"; level="A1";'
> > >
> > > That generates a PostScript plot which you can view using a
PostScript
> > > viewer, like gv or ghostview. Alternatively, you can convert it
to a
> png
> > > and view that:
> > >
> > > convert -rotate 90 -background white apcp_01.ps apcp_01.png
> > >
> > > Does your data look like what you'd expect?
> > >
> > > Thanks,
> > > John
> > >
> > > On Thu, Jun 15, 2017 at 1:58 PM, Ting Lei - NOAA Affiliate via
RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > Thu Jun 15 13:58:23 2017: Request 80881 was acted upon.
> > > > Transaction: Ticket created by ting.lei at noaa.gov
> > > > Queue: met_help
> > > > Subject: how set set up the level parameters in the
configure
> file
> > > > for Met gridstate
> > > > Owner: Nobody
> > > > Requestors: ting.lei at noaa.gov
> > > > Status: new
> > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=80881
> > >
> > > >
> > > >
> > > > Dear Mr. or Mrs.,
> > > > I'm running Met grid_stat,
> > > > I have two kinds of files (fcstfile and obsfile) respectively,
as
> shown
> > > in
> > > > the following:
> > > > 1)
> > > > ^^
> > > >
> > > > wgrib2 -v pcp_03hr_fcst03hrpcp_hi_imp_2015103018.gb2
> > > > 1:0:d=2015103018:WEASD Water Equivalent of Accumulated Snow
Depth
> > > > [kg/m^2]:surface:3 hour fcst:
> > > > 2:105540:d=2015103018:APCP Total Precipitation
[kg/m^2]:surface:2-3
> > hour
> > > > acc fcst:
> > > > 3:509588:d=2015103018:WEASD Water Equivalent of Accumulated
Snow
> Depth
> > > > [kg/m^2]:surface:2-3 hour acc fcst:
> > > >
> > > > VV
> > > > 2)
> > > > ^^
> > > > wgrib2 GaugeOnly_QPE_01H_00.00_20151208-080000.grib2
> > > > 1:0:d=2015120808:var discipline=209 master_table=255 parmcat=6
> > parm=16:0
> > > m
> > > > above mean sea level:anl:
> > > >
> > > > VV
> > > > What should I set as for the level parameter in the configure
file in
> > the
> > > > following:
> > > > ^^
> > > > fcst = {
> > > > wind_thresh = [ NA ];
> > > >
> > > > field = [
> > > > {
> > > > name = "APCP";
> > > > level = [ "" ];
> > > > cat_thresh = [ >0.0, >0.254, >0.508, >1.270, >2.540,
>3.810,
> > > > >6.350, >8.890, >12.700, >25.400 ];
> > > > }
> > > > ];
> > > >
> > > > };
> > > >
> > > >
> > > > V
> > > > I tried surface for the first part, however, failed. I also
have no
> > clue
> > > > what to fill in for the second file.
> > > >
> > > > Your help is appreciated.
> > > > Best Wishes,
> > > > Ting
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: how set set up the level parameters in the configure file for Met gridstate
From: John Halley Gotway
Time: Fri Jun 16 10:59:39 2017
Ting,
Please take a look in the Grid-Stat configuration file at the section
named
"regrid". The is where you specify the grid over which verification
should
be performed. Setting "to_grid = FCST;" tells Grid-Stat to regrid the
observation data to the forecast grid. Setting "to_grid = OBS;" tells
Grid-Stat to regrid the forecast data to the observation grid. Or you
could specify a 3rd grid that differs from the forecast and
observation
domains... and both will be regridded to that.
The method, width, vld_thresh, and shape options control how the
regridding
is done. I'll assume you want to regrid the forecast to the
observation
grid, and since you're processing precip, "budget" interpolation is
the
recommended method.
These options are shown below:
//
// Verification grid
//
regrid = {
to_grid = OBS;
method = BUDGET;
width = 2;
vld_thresh = 0.5;
shape = SQUARE;
}
Hope that helps.
Thanks,
John
On Fri, Jun 16, 2017 at 9:16 AM, Ting Lei - NOAA Affiliate via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=80881 >
>
> Hi, John,
> So many thanks.
> I chose the quick solution you gave and it works seemingly. That is
great!.
> I will dig deeper on the more formal solution you suggested later.
> Just a problem on which your clarification is appreciated.
> In my case, the fcst file and obs file are on different domains of
> different sizes.
> How can I specify which domain is used as the reference domain.
Then, for
> example, the total obs number (relevant in different statistics)
will be
> defined as the obs number in the reference domain.
> Your further help is appreciated.
> Best Wishes,
> Ting
>
>
> On Thu, Jun 15, 2017 at 4:59 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
> > Ting,
> >
> > wgrib2 doesn't know how to handle this data. That's why it listed
the
> > values for the discipline, master_table, parmcat, and parm instead
of
> > listing the corresponding abbreviation:
> >
> > wgrib2 GaugeOnly_QPE_01H_00.00_20151208-080000.grib2
> > 1:0:d=2015120808:var discipline=209 master_table=255 parmcat=6
parm=16:0
> m
> > above mean sea level:anl:
> >
> > This combination of settings doesn't appear in it's GRIB tables...
and
> > likely won't appear in the ones that MET uses either.
> >
> > The quick solution is that you could just tell MET to look for
first
> record
> > number:
> >
> > met-6.0/bin/plot_data_plane \
> > GaugeOnly_QPE_01H_00.00_20151208-080000.grib2 \
> > gauge.ps \
> > 'name="APCP"; level="R1";'
> >
> > Setting level to "R1" tells MET to use the first record regardless
of
> what
> > it actually contains!
> >
> > Alternatively, you could define a custom user GRIB table to define
what
> > you'd like to call this data. Read about MET_GRIB_TABLES on page
47 of
> the
> > user's guide:
> >
> > http://www.dtcenter.org/met/users/docs/users_guide/MET_
> > Users_Guide_v6.0.pdf
> >
> > Thanks,
> > John
> >
> > On Thu, Jun 15, 2017 at 2:41 PM, Ting Lei - NOAA Affiliate via RT
<
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=80881 >
> > >
> > > John,
> > > Thank you so much. I will follow your suggestion.
> > > How about the second kind of data (obs),
> > > ^^
> > > wgrib2 GaugeOnly_QPE_01H_00.00_20151208-080000.grib2
> > > 1:0:d=2015120808 <(201)%20512-0808>:var discipline=209
> master_table=255
> > > parmcat=6 parm=16:0 m above mean sea level:anl:
> > > VV
> > > using g2ctl.pl to generate the grads description file, I got
> > > ^^
> > > var209616_0m 0,102,0 209,6,16 ** 0 m above mean sea level
desc
> [unit]
> > > VV
> > > For this one, how to let MET treat it correctly.
> > > Your further help is appreciated.
> > > Best,
> > > Ting
> > >
> > > On Thu, Jun 15, 2017 at 4:08 PM, John Halley Gotway via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > > Ting,
> > > >
> > > > It looks like you're trying to verify accumulated
precipitation.
> > Here's
> > > > the important line from the output of wgrib2 that you sent:
> > > >
> > > > 2:105540:d=2015103018:APCP Total Precipitation
[kg/m^2]:surface:2-3
> > hour
> > > > acc fcst:
> > > >
> > > > Notices that this is an accumulation between hours 2 and 3.
So
> that's
> > 1
> > > > hour of accumulation.
> > > >
> > > > Please try running MET's plot_data_plane utility to make sure
MET is
> > > > reading your data correctly. Set the level to "A1" to
indicate that
> > you
> > > > want a 1-hourly accumulation interval:
> > > >
> > > > met-6.0/bin/plot_data_plane \
> > > > pcp_03hr_fcst03hrpcp_hi_imp_2015103018.gb2 \
> > > > apcp_01.ps \
> > > > 'name="APCP"; level="A1";'
> > > >
> > > > That generates a PostScript plot which you can view using a
> PostScript
> > > > viewer, like gv or ghostview. Alternatively, you can convert
it to a
> > png
> > > > and view that:
> > > >
> > > > convert -rotate 90 -background white apcp_01.ps apcp_01.png
> > > >
> > > > Does your data look like what you'd expect?
> > > >
> > > > Thanks,
> > > > John
> > > >
> > > > On Thu, Jun 15, 2017 at 1:58 PM, Ting Lei - NOAA Affiliate via
RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > Thu Jun 15 13:58:23 2017: Request 80881 was acted upon.
> > > > > Transaction: Ticket created by ting.lei at noaa.gov
> > > > > Queue: met_help
> > > > > Subject: how set set up the level parameters in the
configure
> > file
> > > > > for Met gridstate
> > > > > Owner: Nobody
> > > > > Requestors: ting.lei at noaa.gov
> > > > > Status: new
> > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > Ticket/Display.html?id=80881
> > > >
> > > > >
> > > > >
> > > > > Dear Mr. or Mrs.,
> > > > > I'm running Met grid_stat,
> > > > > I have two kinds of files (fcstfile and obsfile)
respectively, as
> > shown
> > > > in
> > > > > the following:
> > > > > 1)
> > > > > ^^
> > > > >
> > > > > wgrib2 -v pcp_03hr_fcst03hrpcp_hi_imp_2015103018.gb2
> > > > > 1:0:d=2015103018:WEASD Water Equivalent of Accumulated Snow
Depth
> > > > > [kg/m^2]:surface:3 hour fcst:
> > > > > 2:105540:d=2015103018:APCP Total Precipitation
[kg/m^2]:surface:2-3
> > > hour
> > > > > acc fcst:
> > > > > 3:509588:d=2015103018:WEASD Water Equivalent of Accumulated
Snow
> > Depth
> > > > > [kg/m^2]:surface:2-3 hour acc fcst:
> > > > >
> > > > > VV
> > > > > 2)
> > > > > ^^
> > > > > wgrib2 GaugeOnly_QPE_01H_00.00_20151208-080000.grib2
> > > > > 1:0:d=2015120808:var discipline=209 master_table=255
parmcat=6
> > > parm=16:0
> > > > m
> > > > > above mean sea level:anl:
> > > > >
> > > > > VV
> > > > > What should I set as for the level parameter in the
configure file
> in
> > > the
> > > > > following:
> > > > > ^^
> > > > > fcst = {
> > > > > wind_thresh = [ NA ];
> > > > >
> > > > > field = [
> > > > > {
> > > > > name = "APCP";
> > > > > level = [ "" ];
> > > > > cat_thresh = [ >0.0, >0.254, >0.508, >1.270, >2.540,
> >3.810,
> > > > > >6.350, >8.890, >12.700, >25.400 ];
> > > > > }
> > > > > ];
> > > > >
> > > > > };
> > > > >
> > > > >
> > > > > V
> > > > > I tried surface for the first part, however, failed. I also
have no
> > > clue
> > > > > what to fill in for the second file.
> > > > >
> > > > > Your help is appreciated.
> > > > > Best Wishes,
> > > > > Ting
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: how set set up the level parameters in the configure file for Met gridstate
From: Ting Lei - NOAA Affiliate
Time: Fri Jun 16 11:10:27 2017
John,
Thank you so much for your perfect explanation as well as the
solutions.
I will follow your suggestions to check my configure file.
Your help is really appreciated.
Have a good weekend,
Ting
On Fri, Jun 16, 2017 at 12:59 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Ting,
>
> Please take a look in the Grid-Stat configuration file at the
section named
> "regrid". The is where you specify the grid over which verification
should
> be performed. Setting "to_grid = FCST;" tells Grid-Stat to regrid
the
> observation data to the forecast grid. Setting "to_grid = OBS;"
tells
> Grid-Stat to regrid the forecast data to the observation grid. Or
you
> could specify a 3rd grid that differs from the forecast and
observation
> domains... and both will be regridded to that.
>
> The method, width, vld_thresh, and shape options control how the
regridding
> is done. I'll assume you want to regrid the forecast to the
observation
> grid, and since you're processing precip, "budget" interpolation is
the
> recommended method.
>
> These options are shown below:
>
> //
> // Verification grid
> //
> regrid = {
> to_grid = OBS;
> method = BUDGET;
> width = 2;
> vld_thresh = 0.5;
> shape = SQUARE;
> }
>
> Hope that helps.
>
> Thanks,
> John
>
> On Fri, Jun 16, 2017 at 9:16 AM, Ting Lei - NOAA Affiliate via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=80881 >
> >
> > Hi, John,
> > So many thanks.
> > I chose the quick solution you gave and it works seemingly. That
is
> great!.
> > I will dig deeper on the more formal solution you suggested later.
> > Just a problem on which your clarification is appreciated.
> > In my case, the fcst file and obs file are on different domains of
> > different sizes.
> > How can I specify which domain is used as the reference domain.
Then, for
> > example, the total obs number (relevant in different statistics)
will be
> > defined as the obs number in the reference domain.
> > Your further help is appreciated.
> > Best Wishes,
> > Ting
> >
> >
> > On Thu, Jun 15, 2017 at 4:59 PM, John Halley Gotway via RT <
> > met_help at ucar.edu> wrote:
> >
> > > Ting,
> > >
> > > wgrib2 doesn't know how to handle this data. That's why it
listed the
> > > values for the discipline, master_table, parmcat, and parm
instead of
> > > listing the corresponding abbreviation:
> > >
> > > wgrib2 GaugeOnly_QPE_01H_00.00_20151208-080000.grib2
> > > 1:0:d=2015120808:var discipline=209 master_table=255 parmcat=6
> parm=16:0
> > m
> > > above mean sea level:anl:
> > >
> > > This combination of settings doesn't appear in it's GRIB
tables... and
> > > likely won't appear in the ones that MET uses either.
> > >
> > > The quick solution is that you could just tell MET to look for
first
> > record
> > > number:
> > >
> > > met-6.0/bin/plot_data_plane \
> > > GaugeOnly_QPE_01H_00.00_20151208-080000.grib2 \
> > > gauge.ps \
> > > 'name="APCP"; level="R1";'
> > >
> > > Setting level to "R1" tells MET to use the first record
regardless of
> > what
> > > it actually contains!
> > >
> > > Alternatively, you could define a custom user GRIB table to
define what
> > > you'd like to call this data. Read about MET_GRIB_TABLES on
page 47 of
> > the
> > > user's guide:
> > >
> > > http://www.dtcenter.org/met/users/docs/users_guide/MET_
> > > Users_Guide_v6.0.pdf
> > >
> > > Thanks,
> > > John
> > >
> > > On Thu, Jun 15, 2017 at 2:41 PM, Ting Lei - NOAA Affiliate via
RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=80881
>
> > > >
> > > > John,
> > > > Thank you so much. I will follow your suggestion.
> > > > How about the second kind of data (obs),
> > > > ^^
> > > > wgrib2 GaugeOnly_QPE_01H_00.00_20151208-080000.grib2
> > > > 1:0:d=2015120808 <(201)%20512-0808>:var discipline=209
> > master_table=255
> > > > parmcat=6 parm=16:0 m above mean sea level:anl:
> > > > VV
> > > > using g2ctl.pl to generate the grads description file, I got
> > > > ^^
> > > > var209616_0m 0,102,0 209,6,16 ** 0 m above mean sea level
desc
> > [unit]
> > > > VV
> > > > For this one, how to let MET treat it correctly.
> > > > Your further help is appreciated.
> > > > Best,
> > > > Ting
> > > >
> > > > On Thu, Jun 15, 2017 at 4:08 PM, John Halley Gotway via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > > Ting,
> > > > >
> > > > > It looks like you're trying to verify accumulated
precipitation.
> > > Here's
> > > > > the important line from the output of wgrib2 that you sent:
> > > > >
> > > > > 2:105540:d=2015103018:APCP Total Precipitation
> [kg/m^2]:surface:2-3
> > > hour
> > > > > acc fcst:
> > > > >
> > > > > Notices that this is an accumulation between hours 2 and 3.
So
> > that's
> > > 1
> > > > > hour of accumulation.
> > > > >
> > > > > Please try running MET's plot_data_plane utility to make
sure MET
> is
> > > > > reading your data correctly. Set the level to "A1" to
indicate
> that
> > > you
> > > > > want a 1-hourly accumulation interval:
> > > > >
> > > > > met-6.0/bin/plot_data_plane \
> > > > > pcp_03hr_fcst03hrpcp_hi_imp_2015103018.gb2 \
> > > > > apcp_01.ps \
> > > > > 'name="APCP"; level="A1";'
> > > > >
> > > > > That generates a PostScript plot which you can view using a
> > PostScript
> > > > > viewer, like gv or ghostview. Alternatively, you can
convert it
> to a
> > > png
> > > > > and view that:
> > > > >
> > > > > convert -rotate 90 -background white apcp_01.ps apcp_01.png
> > > > >
> > > > > Does your data look like what you'd expect?
> > > > >
> > > > > Thanks,
> > > > > John
> > > > >
> > > > > On Thu, Jun 15, 2017 at 1:58 PM, Ting Lei - NOAA Affiliate
via RT <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > > >
> > > > > > Thu Jun 15 13:58:23 2017: Request 80881 was acted upon.
> > > > > > Transaction: Ticket created by ting.lei at noaa.gov
> > > > > > Queue: met_help
> > > > > > Subject: how set set up the level parameters in the
> configure
> > > file
> > > > > > for Met gridstate
> > > > > > Owner: Nobody
> > > > > > Requestors: ting.lei at noaa.gov
> > > > > > Status: new
> > > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > > Ticket/Display.html?id=80881
> > > > >
> > > > > >
> > > > > >
> > > > > > Dear Mr. or Mrs.,
> > > > > > I'm running Met grid_stat,
> > > > > > I have two kinds of files (fcstfile and obsfile)
respectively, as
> > > shown
> > > > > in
> > > > > > the following:
> > > > > > 1)
> > > > > > ^^
> > > > > >
> > > > > > wgrib2 -v pcp_03hr_fcst03hrpcp_hi_imp_2015103018.gb2
> > > > > > 1:0:d=2015103018:WEASD Water Equivalent of Accumulated
Snow
> Depth
> > > > > > [kg/m^2]:surface:3 hour fcst:
> > > > > > 2:105540:d=2015103018:APCP Total Precipitation
> [kg/m^2]:surface:2-3
> > > > hour
> > > > > > acc fcst:
> > > > > > 3:509588:d=2015103018:WEASD Water Equivalent of
Accumulated Snow
> > > Depth
> > > > > > [kg/m^2]:surface:2-3 hour acc fcst:
> > > > > >
> > > > > > VV
> > > > > > 2)
> > > > > > ^^
> > > > > > wgrib2 GaugeOnly_QPE_01H_00.00_20151208-080000.grib2
> > > > > > 1:0:d=2015120808:var discipline=209 master_table=255
parmcat=6
> > > > parm=16:0
> > > > > m
> > > > > > above mean sea level:anl:
> > > > > >
> > > > > > VV
> > > > > > What should I set as for the level parameter in the
configure
> file
> > in
> > > > the
> > > > > > following:
> > > > > > ^^
> > > > > > fcst = {
> > > > > > wind_thresh = [ NA ];
> > > > > >
> > > > > > field = [
> > > > > > {
> > > > > > name = "APCP";
> > > > > > level = [ "" ];
> > > > > > cat_thresh = [ >0.0, >0.254, >0.508, >1.270,
>2.540,
> > >3.810,
> > > > > > >6.350, >8.890, >12.700, >25.400 ];
> > > > > > }
> > > > > > ];
> > > > > >
> > > > > > };
> > > > > >
> > > > > >
> > > > > > V
> > > > > > I tried surface for the first part, however, failed. I
also have
> no
> > > > clue
> > > > > > what to fill in for the second file.
> > > > > >
> > > > > > Your help is appreciated.
> > > > > > Best Wishes,
> > > > > > Ting
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
More information about the Met_help
mailing list