[Met_help] [rt.rap.ucar.edu #81800] History for MET6.0/Point_stat
John Halley Gotway via RT
met_help at ucar.edu
Wed Aug 30 09:32:55 MDT 2017
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
I ran into some errors running point_stat, with the command, its arguments
and output below. To describe what I was attempting:
I ran pcp_combine on 4 wpc 2.5km 6 hourly precip accumulations
to get a 24 hour accumulation. Similarly for 2.5km urma. Then I wanted to
run point_stat, with results below (config file attached):
The config file was based on the tutorial example.
One other question: can you provide an example of the simplest possible
config file for point_stat? I would like to start with that and then get
more complex as my understanding grows.
Thanks for your help, Brad Diehl
(g10a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ point_stat out/pcp_combine/
wpc_fcst_24_2017082400.nc out/pcp_combine/wpc_urma_24_2017082400.nc
config/PointStatConfig_24 -outdir out/point_stat -v 4
DEBUG 1: Default Config File:
/usrx/local/dev/MET/6.0/share/met/config/PointStatConfig_default
DEBUG 1: User Config File: config/PointStatConfig_24
DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
Met2dDataFile object of type "FileType_NcMet".
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: VarInfoFactory::new_var_info() -> created new VarInfo object of
type "FileType_NcMet".
DEBUG 4: parse_conf_sid_exc() -> Station ID exclusion list contains 0
entries.
ERROR :
ERROR : VarInfoGrib::add_grib_code() -> unrecognized GRIB1 field
abbreviation 'APCP_24' for table version 2
ERROR :
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: MET6.0/Point_stat
From: John Halley Gotway
Time: Tue Aug 29 12:39:13 2017
Hi Brad,
It sounds like you're comparing a gridded forecast of 24-hour
accumulated
precip to a gridded analysis of 24-hour accumulated precip (URMA). If
that's the case, you should be running the Grid-Stat tool and not
Point-Stat. Grid-Stat is for grid-to-grid verification while Point-
Stat is
for grid-to-point verification.
I see that these are the files you're comparing:
wpc_fcst_24_2017082400.nc
wpc_urma_24_2017082400.nc
Set up your Grid-Stat config file like this:
fcst = {
field = [
{
name = "APCP_24";
level = "(*,*)";
cat_thresh = [ >0, >5, >10 ];
}
];
}
obs = fcst;
And then run Grid-Stat to compare those two files.
Make sense?
Thanks,
John
On Tue, Aug 29, 2017 at 12:17 PM, Brad Diehl - NOAA Federal via RT <
met_help at ucar.edu> wrote:
>
> Tue Aug 29 12:17:28 2017: Request 81800 was acted upon.
> Transaction: Ticket created by brad.diehl at noaa.gov
> Queue: met_help
> Subject: MET6.0/Point_stat
> Owner: Nobody
> Requestors: brad.diehl at noaa.gov
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81800 >
>
>
> I ran into some errors running point_stat, with the command, its
arguments
> and output below. To describe what I was attempting:
>
> I ran pcp_combine on 4 wpc 2.5km 6 hourly precip accumulations
> to get a 24 hour accumulation. Similarly for 2.5km urma. Then I
wanted to
> run point_stat, with results below (config file attached):
>
> The config file was based on the tutorial example.
>
> One other question: can you provide an example of the simplest
possible
> config file for point_stat? I would like to start with that and then
get
> more complex as my understanding grows.
>
> Thanks for your help, Brad Diehl
>
>
> (g10a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ point_stat
> out/pcp_combine/
> wpc_fcst_24_2017082400.nc out/pcp_combine/wpc_urma_24_2017082400.nc
> config/PointStatConfig_24 -outdir out/point_stat -v 4
> DEBUG 1: Default Config File:
> /usrx/local/dev/MET/6.0/share/met/config/PointStatConfig_default
> DEBUG 1: User Config File: config/PointStatConfig_24
> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
> Met2dDataFile object of type "FileType_NcMet".
> 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: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> type "FileType_NcMet".
> DEBUG 4: parse_conf_sid_exc() -> Station ID exclusion list contains
0
> entries.
> ERROR :
> ERROR : VarInfoGrib::add_grib_code() -> unrecognized GRIB1 field
> abbreviation 'APCP_24' for table version 2
> ERROR :
>
>
------------------------------------------------
Subject: MET6.0/Point_stat
From: Brad Diehl - NOAA Federal
Time: Tue Aug 29 12:41:24 2017
Upon trying to run Grid_stat, its output made me realize the grid
domain of
the obs and fcst are not the same. I'll have to provide a mask to the
wpc
2.5km qpf to get it to the same grid and number of grid points as the
urma
grid. Partial output of Grid_stat follows:
ERROR : parse_vx_grid() -> The forecast and observation grids do not
match:
ERROR : Projection: Lambert Conformal Nx: 2145 Ny: 1597 Lat_LL:
20.192
Lon_LL: 121.554 Lon_orient: 95.000 Alpha: 6509.082 Cone: 0.423 Bx:
1088.0030 By: 5483.6617 !=
ERROR : Projection: Lambert Conformal Nx: 2145 Ny: 1377 Lat_LL:
20.192
Lon_LL: 121.554 Lon_orient: 95.000 Alpha: 6509.082 Cone: 0.423 Bx:
1088.0030 By: 5483.6617
ERROR : Specify regridding logic in the config file "regrid" section.
Thanks, Brad Diehl
On Tue, Aug 29, 2017 at 6:17 PM, met_help at ucar.edu via RT
<met_help at ucar.edu
> wrote:
> NOTE: THE MET TEAM HAS LIMITED AVAILABILITY BETWEEN 5/23/17 AND
5/30/17.
> RESPONSES MAY BE DELAYED.
>
> Greetings,
>
> This message has been automatically generated in response to the
> creation of a trouble ticket regarding:
> "MET6.0/Point_stat",
> a summary of which appears below.
>
> There is no need to reply to this message right now. Your ticket
has been
> assigned an ID of [rt.rap.ucar.edu #81800].
>
> Please include the string:
>
> [rt.rap.ucar.edu #81800]
>
> in the subject line of all future correspondence about this issue.
To do
> so,
> you may reply to this message.
>
> Thank you,
> met_help at ucar.edu
>
>
-------------------------------------------------------------------------
> I ran into some errors running point_stat, with the command, its
arguments
> and output below. To describe what I was attempting:
>
> I ran pcp_combine on 4 wpc 2.5km 6 hourly precip accumulations
> to get a 24 hour accumulation. Similarly for 2.5km urma. Then I
wanted to
> run point_stat, with results below (config file attached):
>
> The config file was based on the tutorial example.
>
> One other question: can you provide an example of the simplest
possible
> config file for point_stat? I would like to start with that and then
get
> more complex as my understanding grows.
>
> Thanks for your help, Brad Diehl
>
>
> (g10a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ point_stat
> out/pcp_combine/
> wpc_fcst_24_2017082400.nc out/pcp_combine/wpc_urma_24_2017082400.nc
> config/PointStatConfig_24 -outdir out/point_stat -v 4
> DEBUG 1: Default Config File:
> /usrx/local/dev/MET/6.0/share/met/config/PointStatConfig_default
> DEBUG 1: User Config File: config/PointStatConfig_24
> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
> Met2dDataFile object of type "FileType_NcMet".
> 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: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> type "FileType_NcMet".
> DEBUG 4: parse_conf_sid_exc() -> Station ID exclusion list contains
0
> entries.
> ERROR :
> ERROR : VarInfoGrib::add_grib_code() -> unrecognized GRIB1 field
> abbreviation 'APCP_24' for table version 2
> ERROR :
>
>
------------------------------------------------
Subject: MET6.0/Point_stat
From: Brad Diehl - NOAA Federal
Time: Tue Aug 29 12:46:04 2017
I see your point about point_stat, stupid me, sorry for wasting your
time.
Sent you the following about Grid_stat, and believe what I said is
true,
I'll have to provide a mask to get the obs and fcst file the same.
Thanks, Brad
Upon trying to run Grid_stat, its output made me realize the grid
domain of
the obs and fcst are not the same. I'll have to provide a mask to the
wpc
2.5km qpf to get it to the same grid and number of grid points as the
urma
grid. Partial output of Grid_stat follows:
ERROR : parse_vx_grid() -> The forecast and observation grids do not
match:
ERROR : Projection: Lambert Conformal Nx: 2145 Ny: 1597 Lat_LL:
20.192
Lon_LL: 121.554 Lon_orient: 95.000 Alpha: 6509.082 Cone: 0.423 Bx:
1088.0030 By: 5483.6617 !=
ERROR : Projection: Lambert Conformal Nx: 2145 Ny: 1377 Lat_LL:
20.192
Lon_LL: 121.554 Lon_orient: 95.000 Alpha: 6509.082 Cone: 0.423 Bx:
1088.0030 By: 5483.6617
ERROR : Specify regridding logic in the config file "regrid" section.
On Tue, Aug 29, 2017 at 6:39 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Hi Brad,
>
> It sounds like you're comparing a gridded forecast of 24-hour
accumulated
> precip to a gridded analysis of 24-hour accumulated precip (URMA).
If
> that's the case, you should be running the Grid-Stat tool and not
> Point-Stat. Grid-Stat is for grid-to-grid verification while Point-
Stat is
> for grid-to-point verification.
>
> I see that these are the files you're comparing:
> wpc_fcst_24_2017082400.nc
> wpc_urma_24_2017082400.nc
>
> Set up your Grid-Stat config file like this:
>
> fcst = {
> field = [
> {
> name = "APCP_24";
> level = "(*,*)";
> cat_thresh = [ >0, >5, >10 ];
> }
> ];
> }
> obs = fcst;
>
> And then run Grid-Stat to compare those two files.
>
> Make sense?
>
> Thanks,
> John
>
> On Tue, Aug 29, 2017 at 12:17 PM, Brad Diehl - NOAA Federal via RT <
> met_help at ucar.edu> wrote:
>
> >
> > Tue Aug 29 12:17:28 2017: Request 81800 was acted upon.
> > Transaction: Ticket created by brad.diehl at noaa.gov
> > Queue: met_help
> > Subject: MET6.0/Point_stat
> > Owner: Nobody
> > Requestors: brad.diehl at noaa.gov
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81800 >
> >
> >
> > I ran into some errors running point_stat, with the command, its
> arguments
> > and output below. To describe what I was attempting:
> >
> > I ran pcp_combine on 4 wpc 2.5km 6 hourly precip accumulations
> > to get a 24 hour accumulation. Similarly for 2.5km urma. Then I
wanted to
> > run point_stat, with results below (config file attached):
> >
> > The config file was based on the tutorial example.
> >
> > One other question: can you provide an example of the simplest
possible
> > config file for point_stat? I would like to start with that and
then get
> > more complex as my understanding grows.
> >
> > Thanks for your help, Brad Diehl
> >
> >
> > (g10a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ point_stat
> > out/pcp_combine/
> > wpc_fcst_24_2017082400.nc
out/pcp_combine/wpc_urma_24_2017082400.nc
> > config/PointStatConfig_24 -outdir out/point_stat -v 4
> > DEBUG 1: Default Config File:
> > /usrx/local/dev/MET/6.0/share/met/config/PointStatConfig_default
> > DEBUG 1: User Config File: config/PointStatConfig_24
> > DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
> > Met2dDataFile object of type "FileType_NcMet".
> > 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: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> > type "FileType_NcMet".
> > DEBUG 4: parse_conf_sid_exc() -> Station ID exclusion list
contains 0
> > entries.
> > ERROR :
> > ERROR : VarInfoGrib::add_grib_code() -> unrecognized GRIB1 field
> > abbreviation 'APCP_24' for table version 2
> > ERROR :
> >
> >
>
>
------------------------------------------------
Subject: MET6.0/Point_stat
From: John Halley Gotway
Time: Tue Aug 29 12:55:19 2017
Brad,
Good news. You can just instruct MET to do the regridding for you on
the
fly. This is configurable via the "regrid" section of the Grid-Stat
config
file.
Please read about the "regrid" option on page 49 of the MET User's
Guide:
http://www.dtcenter.org/met/users/docs/users_guide/MET_
Users_Guide_v6.0.pdf
Since you're processing precip here, it's usually recommended to use
BUDGET
interpolation, like this:
regrid = {
to_grid = FCST;
method = BUDGET;
width = 2;
vld_thresh = 0.5;
}
Above, I set "to_grid = FCST;" which tells Grid-Stat to regrid the
URMA
observations to the WPC forecast domain. If you want to go the other
way,
just set "to_grid = OBS;".
Thanks,
John
On Tue, Aug 29, 2017 at 12:46 PM, Brad Diehl - NOAA Federal via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81800 >
>
> I see your point about point_stat, stupid me, sorry for wasting your
time.
>
> Sent you the following about Grid_stat, and believe what I said is
true,
> I'll have to provide a mask to get the obs and fcst file the same.
>
> Thanks, Brad
>
> Upon trying to run Grid_stat, its output made me realize the grid
domain of
> the obs and fcst are not the same. I'll have to provide a mask to
the wpc
> 2.5km qpf to get it to the same grid and number of grid points as
the urma
> grid. Partial output of Grid_stat follows:
>
> ERROR : parse_vx_grid() -> The forecast and observation grids do
not
> match:
> ERROR : Projection: Lambert Conformal Nx: 2145 Ny: 1597 Lat_LL:
20.192
> Lon_LL: 121.554 Lon_orient: 95.000 Alpha: 6509.082 Cone: 0.423 Bx:
> 1088.0030 By: 5483.6617 !=
> ERROR : Projection: Lambert Conformal Nx: 2145 Ny: 1377 Lat_LL:
20.192
> Lon_LL: 121.554 Lon_orient: 95.000 Alpha: 6509.082 Cone: 0.423 Bx:
> 1088.0030 By: 5483.6617
> ERROR : Specify regridding logic in the config file "regrid"
section.
>
> On Tue, Aug 29, 2017 at 6:39 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
> > Hi Brad,
> >
> > It sounds like you're comparing a gridded forecast of 24-hour
accumulated
> > precip to a gridded analysis of 24-hour accumulated precip (URMA).
If
> > that's the case, you should be running the Grid-Stat tool and not
> > Point-Stat. Grid-Stat is for grid-to-grid verification while
Point-Stat
> is
> > for grid-to-point verification.
> >
> > I see that these are the files you're comparing:
> > wpc_fcst_24_2017082400.nc
> > wpc_urma_24_2017082400.nc
> >
> > Set up your Grid-Stat config file like this:
> >
> > fcst = {
> > field = [
> > {
> > name = "APCP_24";
> > level = "(*,*)";
> > cat_thresh = [ >0, >5, >10 ];
> > }
> > ];
> > }
> > obs = fcst;
> >
> > And then run Grid-Stat to compare those two files.
> >
> > Make sense?
> >
> > Thanks,
> > John
> >
> > On Tue, Aug 29, 2017 at 12:17 PM, Brad Diehl - NOAA Federal via RT
<
> > met_help at ucar.edu> wrote:
> >
> > >
> > > Tue Aug 29 12:17:28 2017: Request 81800 was acted upon.
> > > Transaction: Ticket created by brad.diehl at noaa.gov
> > > Queue: met_help
> > > Subject: MET6.0/Point_stat
> > > Owner: Nobody
> > > Requestors: brad.diehl at noaa.gov
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81800
> >
> > >
> > >
> > > I ran into some errors running point_stat, with the command, its
> > arguments
> > > and output below. To describe what I was attempting:
> > >
> > > I ran pcp_combine on 4 wpc 2.5km 6 hourly precip accumulations
> > > to get a 24 hour accumulation. Similarly for 2.5km urma. Then I
wanted
> to
> > > run point_stat, with results below (config file attached):
> > >
> > > The config file was based on the tutorial example.
> > >
> > > One other question: can you provide an example of the simplest
possible
> > > config file for point_stat? I would like to start with that and
then
> get
> > > more complex as my understanding grows.
> > >
> > > Thanks for your help, Brad Diehl
> > >
> > >
> > > (g10a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ point_stat
> > > out/pcp_combine/
> > > wpc_fcst_24_2017082400.nc
out/pcp_combine/wpc_urma_24_2017082400.nc
> > > config/PointStatConfig_24 -outdir out/point_stat -v 4
> > > DEBUG 1: Default Config File:
> > > /usrx/local/dev/MET/6.0/share/met/config/PointStatConfig_default
> > > DEBUG 1: User Config File: config/PointStatConfig_24
> > > DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
> > > Met2dDataFile object of type "FileType_NcMet".
> > > 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: VarInfoFactory::new_var_info() -> created new VarInfo
object
> of
> > > type "FileType_NcMet".
> > > DEBUG 4: parse_conf_sid_exc() -> Station ID exclusion list
contains 0
> > > entries.
> > > ERROR :
> > > ERROR : VarInfoGrib::add_grib_code() -> unrecognized GRIB1
field
> > > abbreviation 'APCP_24' for table version 2
> > > ERROR :
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: MET6.0/Point_stat
From: Brad Diehl - NOAA Federal
Time: Tue Aug 29 14:56:31 2017
Was finally able to work out the bugs, and with your recommendation
about
re-grid, am good to go. Plotted results with R script and output is
what is
expected from our previous discussions. Quite a learning curve, but
finally
getting there. Certainly more questions/concerns will come up, but
appreciate your efforts thus far. You should close the ticket.
On Tue, Aug 29, 2017 at 6:55 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Brad,
>
> Good news. You can just instruct MET to do the regridding for you
on the
> fly. This is configurable via the "regrid" section of the Grid-Stat
config
> file.
>
> Please read about the "regrid" option on page 49 of the MET User's
Guide:
> http://www.dtcenter.org/met/users/docs/users_guide/MET_
> Users_Guide_v6.0.pdf
>
> Since you're processing precip here, it's usually recommended to use
BUDGET
> interpolation, like this:
>
> regrid = {
> to_grid = FCST;
> method = BUDGET;
> width = 2;
> vld_thresh = 0.5;
> }
>
> Above, I set "to_grid = FCST;" which tells Grid-Stat to regrid the
URMA
> observations to the WPC forecast domain. If you want to go the
other way,
> just set "to_grid = OBS;".
>
> Thanks,
> John
>
>
>
>
> On Tue, Aug 29, 2017 at 12:46 PM, Brad Diehl - NOAA Federal via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=81800 >
> >
> > I see your point about point_stat, stupid me, sorry for wasting
your
> time.
> >
> > Sent you the following about Grid_stat, and believe what I said is
true,
> > I'll have to provide a mask to get the obs and fcst file the same.
> >
> > Thanks, Brad
> >
> > Upon trying to run Grid_stat, its output made me realize the grid
domain
> of
> > the obs and fcst are not the same. I'll have to provide a mask to
the wpc
> > 2.5km qpf to get it to the same grid and number of grid points as
the
> urma
> > grid. Partial output of Grid_stat follows:
> >
> > ERROR : parse_vx_grid() -> The forecast and observation grids do
not
> > match:
> > ERROR : Projection: Lambert Conformal Nx: 2145 Ny: 1597 Lat_LL:
20.192
> > Lon_LL: 121.554 Lon_orient: 95.000 Alpha: 6509.082 Cone: 0.423 Bx:
> > 1088.0030 By: 5483.6617 !=
> > ERROR : Projection: Lambert Conformal Nx: 2145 Ny: 1377 Lat_LL:
20.192
> > Lon_LL: 121.554 Lon_orient: 95.000 Alpha: 6509.082 Cone: 0.423 Bx:
> > 1088.0030 By: 5483.6617
> > ERROR : Specify regridding logic in the config file "regrid"
section.
> >
> > On Tue, Aug 29, 2017 at 6:39 PM, John Halley Gotway via RT <
> > met_help at ucar.edu> wrote:
> >
> > > Hi Brad,
> > >
> > > It sounds like you're comparing a gridded forecast of 24-hour
> accumulated
> > > precip to a gridded analysis of 24-hour accumulated precip
(URMA). If
> > > that's the case, you should be running the Grid-Stat tool and
not
> > > Point-Stat. Grid-Stat is for grid-to-grid verification while
> Point-Stat
> > is
> > > for grid-to-point verification.
> > >
> > > I see that these are the files you're comparing:
> > > wpc_fcst_24_2017082400.nc
> > > wpc_urma_24_2017082400.nc
> > >
> > > Set up your Grid-Stat config file like this:
> > >
> > > fcst = {
> > > field = [
> > > {
> > > name = "APCP_24";
> > > level = "(*,*)";
> > > cat_thresh = [ >0, >5, >10 ];
> > > }
> > > ];
> > > }
> > > obs = fcst;
> > >
> > > And then run Grid-Stat to compare those two files.
> > >
> > > Make sense?
> > >
> > > Thanks,
> > > John
> > >
> > > On Tue, Aug 29, 2017 at 12:17 PM, Brad Diehl - NOAA Federal via
RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > Tue Aug 29 12:17:28 2017: Request 81800 was acted upon.
> > > > Transaction: Ticket created by brad.diehl at noaa.gov
> > > > Queue: met_help
> > > > Subject: MET6.0/Point_stat
> > > > Owner: Nobody
> > > > Requestors: brad.diehl at noaa.gov
> > > > Status: new
> > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=81800
> > >
> > > >
> > > >
> > > > I ran into some errors running point_stat, with the command,
its
> > > arguments
> > > > and output below. To describe what I was attempting:
> > > >
> > > > I ran pcp_combine on 4 wpc 2.5km 6 hourly precip accumulations
> > > > to get a 24 hour accumulation. Similarly for 2.5km urma. Then
I
> wanted
> > to
> > > > run point_stat, with results below (config file attached):
> > > >
> > > > The config file was based on the tutorial example.
> > > >
> > > > One other question: can you provide an example of the simplest
> possible
> > > > config file for point_stat? I would like to start with that
and then
> > get
> > > > more complex as my understanding grows.
> > > >
> > > > Thanks for your help, Brad Diehl
> > > >
> > > >
> > > > (g10a2) /u/Brad.Diehl/MET_Tutorial_6.0/tutorial$ point_stat
> > > > out/pcp_combine/
> > > > wpc_fcst_24_2017082400.nc
out/pcp_combine/wpc_urma_24_2017082400.nc
> > > > config/PointStatConfig_24 -outdir out/point_stat -v 4
> > > > DEBUG 1: Default Config File:
> > > >
/usrx/local/dev/MET/6.0/share/met/config/PointStatConfig_default
> > > > DEBUG 1: User Config File: config/PointStatConfig_24
> > > > DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() ->
created new
> > > > Met2dDataFile object of type "FileType_NcMet".
> > > > 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: VarInfoFactory::new_var_info() -> created new VarInfo
object
> > of
> > > > type "FileType_NcMet".
> > > > DEBUG 4: parse_conf_sid_exc() -> Station ID exclusion list
contains 0
> > > > entries.
> > > > ERROR :
> > > > ERROR : VarInfoGrib::add_grib_code() -> unrecognized GRIB1
field
> > > > abbreviation 'APCP_24' for table version 2
> > > > ERROR :
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
More information about the Met_help
mailing list