[Met_help] [rt.rap.ucar.edu #74966] History for Troubles with point_stat
John Halley Gotway via RT
met_help at ucar.edu
Tue Feb 9 09:29:16 MST 2016
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Dear MET users,
I have a problem with, I think, the point_stat tool.
Here a brief history of what I did, just to check if I made some former
mistake.
I made a 24 h simulation with ARW initialised with GFS and starting
2016020100.
I setted the io_form_history=5 to get a grib1 wrfout file (called
wrfout.grb)
I created an ascii file containing the accumulated precipitation in three
hours for
the same dates filling the eleven fields described in table 4.4 of
Met_user_guide_v5.1
just for two stations.
Here is the file (ascii_obs.txt)
NA AA000A001 20160201_030000 9.458 39.961 813 061 3 813 NA 2.6
NA AA000A001 20160201_060000 9.458 39.961 813 061 3 813 NA 2.6
NA AA000A001 20160201_090000 9.458 39.961 813 061 3 813 NA 2.9
NA AA000A001 20160201_120000 9.458 39.961 813 061 3 813 NA 3.6
NA AA000A001 20160201_150000 9.458 39.961 813 061 3 813 NA 7.1
NA AA000A001 20160201_180000 9.458 39.961 813 061 3 813 NA 12.0
NA AA000A001 20160201_210000 9.458 39.961 813 061 3 813 NA 12.5
NA AA000A001 20160202_000000 9.458 39.961 813 061 3 813 NA 14.2
NA AA000A002 20160201_030000 8.832 40.941 6 061 3 6 NA 4.2
NA AA000A002 20160201_060000 8.832 40.941 6 061 3 6 NA 5.6
NA AA000A002 20160201_090000 8.832 40.941 6 061 3 6 NA 12.9
NA AA000A002 20160201_120000 8.832 40.941 6 061 3 6 NA 13.1
NA AA000A002 20160201_150000 8.832 40.941 6 061 3 6 NA 13.1
NA AA000A002 20160201_180000 8.832 40.941 6 061 3 6 NA 14.0
NA AA000A002 20160201_210000 8.832 40.941 6 061 3 6 NA 15.5
NA AA000A002 20160202_000000 8.832 40.941 6 061 3 6 NA 16.1
Afterward I created an nc file using the tool
./ascii2nc ascii_obs.txt ascii_obs.nc -v 2
The output of "ncdump -v obs_arr ascii_obs.nc" looks fine.
I modified PointStatConfig_default like this
fcst = {
sid_exc = [];
field = [
{
name = "APCP";
level = [ "A03" ];
cat_thresh = [ >1.0 ];
file_type = [ 1 ] ;
init_time = [ "20160201_000000" ];
valid_time = [ "20160202_000000" ];
}
];
}
obs = fcst;
and using "./point_stat wrfout.grb ascii_obs.nc PointStatConfig_modified
-outdir ./ -v 2"
the program just give me the headers of the statistic files and the
following message error
DEBUG 1: User Config File: PointStatConfig_modified
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=1087746390
DEBUG 1: Forecast File: wrfout_gfs.grib1
DEBUG 1: Observation File: prova_ascii_grib1.nc
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for APCP/A03.
WARNING:
WARNING: process_fcst_climo_files() -> no fields matching APCP/A03 found in
file: wrfout_gfs.grib1
WARNING:
.....
I also tried to modify the input file like this
fcst = {
sid_exc = [];
field = [
{
name = "APCP_03";
level = [ "(*,*)" ];
cat_thresh = [ >1.0 ];
file_type = [ 1 ] ;
init_time = [ "20160201_000000" ];
valid_time = [ "20160202_000000" ];
}
];
}
obs = {
message_type = [ "ADPSFC" ];
sid_exc = [];
field = [
{
name = "APCP";
level = [ "A03" ];
cat_thresh = [ >1.0 ];
file_type = [ 1 ] ;
init_time = [ "20160201_000000" ];
valid_time = [ "20160202_000000" ];
}
];
}
but in this case the message error is even worse and I get
DEBUG 1: User Config File: PointStatConfig_modified
ERROR :
ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field abbreviation
'APCP_03' for table version 2
ERROR :
without any statistic output file.
Any advise is really welcome.
Thanks,
Samuele
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Troubles with point_stat
From: John Halley Gotway
Time: Wed Feb 03 10:48:28 2016
Samuele,
I see you are having trouble running the Point-Stat tool. The
relevant
error here is this:
ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation
'APCP_03' for table version 2
You are passing Point-Stat a forecast GRIB file (wrfout.grb) but are
specifying in the config file as if it were a gridded NetCDF file.
Let's
back up and use wgrib to inventory your GRIB file and see what
accumulations intervals are present in there:
wgrib wrfout.grb | grep APCP
What do you see? Do you actually have 3-hourly accumulated precip?
If
not, what accumulation interval do you have?
Also, please edit your point observations by specifying a message type
in
the first column. I'd suggest replacing the NA string with ADPSFC,
which
is PREPBUFR message type convention for surface observations.
Assuming you actually have 3-hourly APCP in your GRIB file, use this
in
your Point-Stat config file:
fcst = {
sid_exc = [];
message_type = [ "ADPSFC" ];
field = [
{
name = "APCP";
level = [ "A03" ];
cat_thresh = [ >1.0 ];
}
];
}
obs = fcst;
Thanks,
John Halley Gotway
On Wed, Feb 3, 2016 at 3:43 AM, Samuele Salis via RT
<met_help at ucar.edu>
wrote:
>
> Wed Feb 03 03:43:55 2016: Request 74966 was acted upon.
> Transaction: Ticket created by samuelesalis at gmail.com
> Queue: met_help
> Subject: Troubles with point_stat
> Owner: Nobody
> Requestors: samuelesalis at gmail.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966 >
>
>
> Dear MET users,
>
> I have a problem with, I think, the point_stat tool.
>
> Here a brief history of what I did, just to check if I made some
former
> mistake.
>
> I made a 24 h simulation with ARW initialised with GFS and starting
> 2016020100.
> I setted the io_form_history=5 to get a grib1 wrfout file (called
> wrfout.grb)
>
> I created an ascii file containing the accumulated precipitation in
three
> hours for
> the same dates filling the eleven fields described in table 4.4 of
> Met_user_guide_v5.1
> just for two stations.
> Here is the file (ascii_obs.txt)
> NA AA000A001 20160201_030000 9.458 39.961 813 061 3 813 NA 2.6
> NA AA000A001 20160201_060000 9.458 39.961 813 061 3 813 NA 2.6
> NA AA000A001 20160201_090000 9.458 39.961 813 061 3 813 NA 2.9
> NA AA000A001 20160201_120000 9.458 39.961 813 061 3 813 NA 3.6
> NA AA000A001 20160201_150000 9.458 39.961 813 061 3 813 NA 7.1
> NA AA000A001 20160201_180000 9.458 39.961 813 061 3 813 NA 12.0
> NA AA000A001 20160201_210000 9.458 39.961 813 061 3 813 NA 12.5
> NA AA000A001 20160202_000000 9.458 39.961 813 061 3 813 NA 14.2
> NA AA000A002 20160201_030000 8.832 40.941 6 061 3 6 NA 4.2
> NA AA000A002 20160201_060000 8.832 40.941 6 061 3 6 NA 5.6
> NA AA000A002 20160201_090000 8.832 40.941 6 061 3 6 NA 12.9
> NA AA000A002 20160201_120000 8.832 40.941 6 061 3 6 NA 13.1
> NA AA000A002 20160201_150000 8.832 40.941 6 061 3 6 NA 13.1
> NA AA000A002 20160201_180000 8.832 40.941 6 061 3 6 NA 14.0
> NA AA000A002 20160201_210000 8.832 40.941 6 061 3 6 NA 15.5
> NA AA000A002 20160202_000000 8.832 40.941 6 061 3 6 NA 16.1
>
> Afterward I created an nc file using the tool
> ./ascii2nc ascii_obs.txt ascii_obs.nc -v 2
>
> The output of "ncdump -v obs_arr ascii_obs.nc" looks fine.
>
> I modified PointStatConfig_default like this
>
> fcst = {
> sid_exc = [];
>
> field = [
> {
> name = "APCP";
> level = [ "A03" ];
> cat_thresh = [ >1.0 ];
> file_type = [ 1 ] ;
> init_time = [ "20160201_000000" ];
> valid_time = [ "20160202_000000" ];
> }
> ];
>
> }
> obs = fcst;
>
> and using "./point_stat wrfout.grb ascii_obs.nc
PointStatConfig_modified
> -outdir ./ -v 2"
> the program just give me the headers of the statistic files and the
> following message error
>
> DEBUG 1: User Config File: PointStatConfig_modified
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=1087746390
> DEBUG 1: Forecast File: wrfout_gfs.grib1
> DEBUG 1: Observation File: prova_ascii_grib1.nc
> DEBUG 2:
> DEBUG 2:
>
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for APCP/A03.
> WARNING:
> WARNING: process_fcst_climo_files() -> no fields matching APCP/A03
found in
> file: wrfout_gfs.grib1
> WARNING:
> .....
>
>
> I also tried to modify the input file like this
>
> fcst = {
> sid_exc = [];
>
> field = [
> {
> name = "APCP_03";
> level = [ "(*,*)" ];
> cat_thresh = [ >1.0 ];
> file_type = [ 1 ] ;
> init_time = [ "20160201_000000" ];
> valid_time = [ "20160202_000000" ];
> }
> ];
>
> }
>
> obs = {
> message_type = [ "ADPSFC" ];
> sid_exc = [];
>
> field = [
> {
> name = "APCP";
> level = [ "A03" ];
> cat_thresh = [ >1.0 ];
> file_type = [ 1 ] ;
> init_time = [ "20160201_000000" ];
> valid_time = [ "20160202_000000" ];
> }
> ];
>
> }
>
> but in this case the message error is even worse and I get
>
> DEBUG 1: User Config File: PointStatConfig_modified
> ERROR :
> ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation
> 'APCP_03' for table version 2
> ERROR :
>
> without any statistic output file.
>
>
> Any advise is really welcome.
>
> Thanks,
> Samuele
>
>
------------------------------------------------
Subject: Troubles with point_stat
From: Samuele Salis
Time: Thu Feb 04 02:58:45 2016
Dear John,
first of all thank you for your kind reply.
The main error I found is that the APCP field is not included in my
wrfout.grb file, but I saw the ACPCP parameter, so I changed my
ascii_obs.txt file like this
ADPSFC AA000A001 20160201_030000 9.458 39.961 813 063 3 813 NA 2.6
ADPSFC AA000A001 20160201_060000 9.458 39.961 813 063 6 813 NA 2.6
ADPSFC AA000A001 20160201_090000 9.458 39.961 813 063 9 813 NA 2.9
ADPSFC AA000A001 20160201_120000 9.458 39.961 813 063 12 813 NA 3.6
ADPSFC AA000A001 20160201_150000 9.458 39.961 813 063 15 813 NA 7.1
ADPSFC AA000A001 20160201_180000 9.458 39.961 813 063 18 813 NA 12.0
ADPSFC AA000A001 20160201_210000 9.458 39.961 813 063 21 813 NA 12.5
ADPSFC AA000A001 20160202_000000 9.458 39.961 813 063 24 813 NA 14.2
ADPSFC AA000A002 20160201_030000 8.832 40.941 6 063 3 6 NA 4.2
ADPSFC AA000A002 20160201_060000 8.832 40.941 6 063 6 6 NA 5.6
ADPSFC AA000A002 20160201_090000 8.832 40.941 6 063 9 6 NA 12.9
ADPSFC AA000A002 20160201_120000 8.832 40.941 6 063 12 6 NA 13.1
ADPSFC AA000A002 20160201_150000 8.832 40.941 6 063 15 6 NA 13.1
ADPSFC AA000A002 20160201_180000 8.832 40.941 6 063 18 6 NA 14.0
ADPSFC AA000A002 20160201_210000 8.832 40.941 6 063 21 6 NA 15.5
ADPSFC AA000A002 20160202_000000 8.832 40.941 6 063 24 6 NA 16.1
and I generated again the ascii_obs.nc using ascii2nc.
I modified the PointStatConfig_default file in this way
fcst = {
wind_thresh = [ NA ];
message_type = [ "ADPSFC" ];
sid_exc = [];
field = [
{
name = "ACPCP";
level = [ "L000" ];
cat_thresh = [ >1.0 ];
}
];
}
obs = fcst;
and with the command "./point_stat wrfout.grb ascii_obs.nc
PointStatConfig_modified -outdir ./ -v 2"
I obtained apparently no error messages
DEBUG 1: User Config File: PointStatConfig_modified
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=18446744072318923606
DEBUG 1: Forecast File: wrfout.grb
DEBUG 1: Observation File: ascii_obs.nc
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for ACPCP/L0.
DEBUG 2: For ACPCP/L0 found 25 forecast levels and 0 climatology
levels.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Searching 16 observations from 16 messages.
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Processing ACPCP/L0 versus ACPCP/L0, for observation type
ADPSFC,
over region FULL, for interpolation method NEAREST(1), using 0 pairs.
DEBUG 2:
The problem is that the produced statistical point_stat files are all
empty
except the headers.
I am worry the problem is in my wrfout.grb file, so I will try with
another
one.
Thanks,
Samuele
On 3 February 2016 at 18:48, John Halley Gotway via RT
<met_help at ucar.edu>
wrote:
> Samuele,
>
> I see you are having trouble running the Point-Stat tool. The
relevant
> error here is this:
>
> ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation
> 'APCP_03' for table version 2
>
> You are passing Point-Stat a forecast GRIB file (wrfout.grb) but are
> specifying in the config file as if it were a gridded NetCDF file.
Let's
> back up and use wgrib to inventory your GRIB file and see what
> accumulations intervals are present in there:
>
> wgrib wrfout.grb | grep APCP
>
> What do you see? Do you actually have 3-hourly accumulated precip?
If
> not, what accumulation interval do you have?
>
> Also, please edit your point observations by specifying a message
type in
> the first column. I'd suggest replacing the NA string with ADPSFC,
which
> is PREPBUFR message type convention for surface observations.
>
> Assuming you actually have 3-hourly APCP in your GRIB file, use this
in
> your Point-Stat config file:
>
> fcst = {
> sid_exc = [];
> message_type = [ "ADPSFC" ];
>
> field = [
> {
> name = "APCP";
> level = [ "A03" ];
> cat_thresh = [ >1.0 ];
> }
> ];
> }
> obs = fcst;
>
>
> Thanks,
> John Halley Gotway
>
> On Wed, Feb 3, 2016 at 3:43 AM, Samuele Salis via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > Wed Feb 03 03:43:55 2016: Request 74966 was acted upon.
> > Transaction: Ticket created by samuelesalis at gmail.com
> > Queue: met_help
> > Subject: Troubles with point_stat
> > Owner: Nobody
> > Requestors: samuelesalis at gmail.com
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966 >
> >
> >
> > Dear MET users,
> >
> > I have a problem with, I think, the point_stat tool.
> >
> > Here a brief history of what I did, just to check if I made some
former
> > mistake.
> >
> > I made a 24 h simulation with ARW initialised with GFS and
starting
> > 2016020100.
> > I setted the io_form_history=5 to get a grib1 wrfout file (called
> > wrfout.grb)
> >
> > I created an ascii file containing the accumulated precipitation
in three
> > hours for
> > the same dates filling the eleven fields described in table 4.4 of
> > Met_user_guide_v5.1
> > just for two stations.
> > Here is the file (ascii_obs.txt)
> > NA AA000A001 20160201_030000 9.458 39.961 813 061 3 813 NA 2.6
> > NA AA000A001 20160201_060000 9.458 39.961 813 061 3 813 NA 2.6
> > NA AA000A001 20160201_090000 9.458 39.961 813 061 3 813 NA 2.9
> > NA AA000A001 20160201_120000 9.458 39.961 813 061 3 813 NA 3.6
> > NA AA000A001 20160201_150000 9.458 39.961 813 061 3 813 NA 7.1
> > NA AA000A001 20160201_180000 9.458 39.961 813 061 3 813 NA 12.0
> > NA AA000A001 20160201_210000 9.458 39.961 813 061 3 813 NA 12.5
> > NA AA000A001 20160202_000000 9.458 39.961 813 061 3 813 NA 14.2
> > NA AA000A002 20160201_030000 8.832 40.941 6 061 3 6 NA 4.2
> > NA AA000A002 20160201_060000 8.832 40.941 6 061 3 6 NA 5.6
> > NA AA000A002 20160201_090000 8.832 40.941 6 061 3 6 NA 12.9
> > NA AA000A002 20160201_120000 8.832 40.941 6 061 3 6 NA 13.1
> > NA AA000A002 20160201_150000 8.832 40.941 6 061 3 6 NA 13.1
> > NA AA000A002 20160201_180000 8.832 40.941 6 061 3 6 NA 14.0
> > NA AA000A002 20160201_210000 8.832 40.941 6 061 3 6 NA 15.5
> > NA AA000A002 20160202_000000 8.832 40.941 6 061 3 6 NA 16.1
> >
> > Afterward I created an nc file using the tool
> > ./ascii2nc ascii_obs.txt ascii_obs.nc -v 2
> >
> > The output of "ncdump -v obs_arr ascii_obs.nc" looks fine.
> >
> > I modified PointStatConfig_default like this
> >
> > fcst = {
> > sid_exc = [];
> >
> > field = [
> > {
> > name = "APCP";
> > level = [ "A03" ];
> > cat_thresh = [ >1.0 ];
> > file_type = [ 1 ] ;
> > init_time = [ "20160201_000000" ];
> > valid_time = [ "20160202_000000" ];
> > }
> > ];
> >
> > }
> > obs = fcst;
> >
> > and using "./point_stat wrfout.grb ascii_obs.nc
> PointStatConfig_modified
> > -outdir ./ -v 2"
> > the program just give me the headers of the statistic files and
the
> > following message error
> >
> > DEBUG 1: User Config File: PointStatConfig_modified
> > GSL_RNG_TYPE=mt19937
> > GSL_RNG_SEED=1087746390
> > DEBUG 1: Forecast File: wrfout_gfs.grib1
> > DEBUG 1: Observation File: prova_ascii_grib1.nc
> > DEBUG 2:
> > DEBUG 2:
> >
> >
>
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Reading data for APCP/A03.
> > WARNING:
> > WARNING: process_fcst_climo_files() -> no fields matching APCP/A03
found
> in
> > file: wrfout_gfs.grib1
> > WARNING:
> > .....
> >
> >
> > I also tried to modify the input file like this
> >
> > fcst = {
> > sid_exc = [];
> >
> > field = [
> > {
> > name = "APCP_03";
> > level = [ "(*,*)" ];
> > cat_thresh = [ >1.0 ];
> > file_type = [ 1 ] ;
> > init_time = [ "20160201_000000" ];
> > valid_time = [ "20160202_000000" ];
> > }
> > ];
> >
> > }
> >
> > obs = {
> > message_type = [ "ADPSFC" ];
> > sid_exc = [];
> >
> > field = [
> > {
> > name = "APCP";
> > level = [ "A03" ];
> > cat_thresh = [ >1.0 ];
> > file_type = [ 1 ] ;
> > init_time = [ "20160201_000000" ];
> > valid_time = [ "20160202_000000" ];
> > }
> > ];
> >
> > }
> >
> > but in this case the message error is even worse and I get
> >
> > DEBUG 1: User Config File: PointStatConfig_modified
> > ERROR :
> > ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation
> > 'APCP_03' for table version 2
> > ERROR :
> >
> > without any statistic output file.
> >
> >
> > Any advise is really welcome.
> >
> > Thanks,
> > Samuele
> >
> >
>
>
------------------------------------------------
Subject: Troubles with point_stat
From: John Halley Gotway
Time: Thu Feb 04 09:36:16 2016
Samuele,
When you run Point-Stat, try increasing the verbosity level to 3 or 4
(-v
4) to get more diagnostic information about why counts were or were
not
used.
At this point, I think it'd be easiest to just have you send me some
sample
data. I'll run it here and provide you some suggestions.
Please post ascii_obs.txt, wrfout.grb, and your Point-Stat
configuration
files to our anonymous ftp site following these directions:
http://www.dtcenter.org/met/users/support/met_help.php#ftp
Please let me know when you've posted these files, and I'll go grab
them.
Thanks,
John
On Thu, Feb 4, 2016 at 2:58 AM, Samuele Salis via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966 >
>
> Dear John,
>
> first of all thank you for your kind reply.
>
> The main error I found is that the APCP field is not included in my
> wrfout.grb file, but I saw the ACPCP parameter, so I changed my
> ascii_obs.txt file like this
>
> ADPSFC AA000A001 20160201_030000 9.458 39.961 813 063 3 813 NA 2.6
> ADPSFC AA000A001 20160201_060000 9.458 39.961 813 063 6 813 NA 2.6
> ADPSFC AA000A001 20160201_090000 9.458 39.961 813 063 9 813 NA 2.9
> ADPSFC AA000A001 20160201_120000 9.458 39.961 813 063 12 813 NA 3.6
> ADPSFC AA000A001 20160201_150000 9.458 39.961 813 063 15 813 NA 7.1
> ADPSFC AA000A001 20160201_180000 9.458 39.961 813 063 18 813 NA 12.0
> ADPSFC AA000A001 20160201_210000 9.458 39.961 813 063 21 813 NA 12.5
> ADPSFC AA000A001 20160202_000000 9.458 39.961 813 063 24 813 NA 14.2
> ADPSFC AA000A002 20160201_030000 8.832 40.941 6 063 3 6 NA 4.2
> ADPSFC AA000A002 20160201_060000 8.832 40.941 6 063 6 6 NA 5.6
> ADPSFC AA000A002 20160201_090000 8.832 40.941 6 063 9 6 NA 12.9
> ADPSFC AA000A002 20160201_120000 8.832 40.941 6 063 12 6 NA 13.1
> ADPSFC AA000A002 20160201_150000 8.832 40.941 6 063 15 6 NA 13.1
> ADPSFC AA000A002 20160201_180000 8.832 40.941 6 063 18 6 NA 14.0
> ADPSFC AA000A002 20160201_210000 8.832 40.941 6 063 21 6 NA 15.5
> ADPSFC AA000A002 20160202_000000 8.832 40.941 6 063 24 6 NA 16.1
>
> and I generated again the ascii_obs.nc using ascii2nc.
>
> I modified the PointStatConfig_default file in this way
>
> fcst = {
> wind_thresh = [ NA ];
> message_type = [ "ADPSFC" ];
> sid_exc = [];
>
> field = [
> {
> name = "ACPCP";
> level = [ "L000" ];
> cat_thresh = [ >1.0 ];
> }
> ];
>
> }
> obs = fcst;
>
> and with the command "./point_stat wrfout.grb ascii_obs.nc
> PointStatConfig_modified -outdir ./ -v 2"
> I obtained apparently no error messages
>
> DEBUG 1: User Config File: PointStatConfig_modified
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=18446744072318923606
> DEBUG 1: Forecast File: wrfout.grb
> DEBUG 1: Observation File: ascii_obs.nc
> DEBUG 2:
> DEBUG 2:
>
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for ACPCP/L0.
> DEBUG 2: For ACPCP/L0 found 25 forecast levels and 0 climatology
levels.
> DEBUG 2:
> DEBUG 2:
>
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Searching 16 observations from 16 messages.
> DEBUG 2:
> DEBUG 2:
>
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Processing ACPCP/L0 versus ACPCP/L0, for observation type
ADPSFC,
> over region FULL, for interpolation method NEAREST(1), using 0
pairs.
> DEBUG 2:
>
> The problem is that the produced statistical point_stat files are
all empty
> except the headers.
>
> I am worry the problem is in my wrfout.grb file, so I will try with
another
> one.
>
>
> Thanks,
> Samuele
>
>
>
> On 3 February 2016 at 18:48, John Halley Gotway via RT
<met_help at ucar.edu>
> wrote:
>
> > Samuele,
> >
> > I see you are having trouble running the Point-Stat tool. The
relevant
> > error here is this:
> >
> > ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
abbreviation
> > 'APCP_03' for table version 2
> >
> > You are passing Point-Stat a forecast GRIB file (wrfout.grb) but
are
> > specifying in the config file as if it were a gridded NetCDF file.
Let's
> > back up and use wgrib to inventory your GRIB file and see what
> > accumulations intervals are present in there:
> >
> > wgrib wrfout.grb | grep APCP
> >
> > What do you see? Do you actually have 3-hourly accumulated
precip? If
> > not, what accumulation interval do you have?
> >
> > Also, please edit your point observations by specifying a message
type in
> > the first column. I'd suggest replacing the NA string with
ADPSFC, which
> > is PREPBUFR message type convention for surface observations.
> >
> > Assuming you actually have 3-hourly APCP in your GRIB file, use
this in
> > your Point-Stat config file:
> >
> > fcst = {
> > sid_exc = [];
> > message_type = [ "ADPSFC" ];
> >
> > field = [
> > {
> > name = "APCP";
> > level = [ "A03" ];
> > cat_thresh = [ >1.0 ];
> > }
> > ];
> > }
> > obs = fcst;
> >
> >
> > Thanks,
> > John Halley Gotway
> >
> > On Wed, Feb 3, 2016 at 3:43 AM, Samuele Salis via RT
<met_help at ucar.edu>
> > wrote:
> >
> > >
> > > Wed Feb 03 03:43:55 2016: Request 74966 was acted upon.
> > > Transaction: Ticket created by samuelesalis at gmail.com
> > > Queue: met_help
> > > Subject: Troubles with point_stat
> > > Owner: Nobody
> > > Requestors: samuelesalis at gmail.com
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966
> >
> > >
> > >
> > > Dear MET users,
> > >
> > > I have a problem with, I think, the point_stat tool.
> > >
> > > Here a brief history of what I did, just to check if I made some
former
> > > mistake.
> > >
> > > I made a 24 h simulation with ARW initialised with GFS and
starting
> > > 2016020100.
> > > I setted the io_form_history=5 to get a grib1 wrfout file
(called
> > > wrfout.grb)
> > >
> > > I created an ascii file containing the accumulated precipitation
in
> three
> > > hours for
> > > the same dates filling the eleven fields described in table 4.4
of
> > > Met_user_guide_v5.1
> > > just for two stations.
> > > Here is the file (ascii_obs.txt)
> > > NA AA000A001 20160201_030000 9.458 39.961 813 061 3 813 NA 2.6
> > > NA AA000A001 20160201_060000 9.458 39.961 813 061 3 813 NA 2.6
> > > NA AA000A001 20160201_090000 9.458 39.961 813 061 3 813 NA 2.9
> > > NA AA000A001 20160201_120000 9.458 39.961 813 061 3 813 NA 3.6
> > > NA AA000A001 20160201_150000 9.458 39.961 813 061 3 813 NA 7.1
> > > NA AA000A001 20160201_180000 9.458 39.961 813 061 3 813 NA 12.0
> > > NA AA000A001 20160201_210000 9.458 39.961 813 061 3 813 NA 12.5
> > > NA AA000A001 20160202_000000 9.458 39.961 813 061 3 813 NA 14.2
> > > NA AA000A002 20160201_030000 8.832 40.941 6 061 3 6 NA 4.2
> > > NA AA000A002 20160201_060000 8.832 40.941 6 061 3 6 NA 5.6
> > > NA AA000A002 20160201_090000 8.832 40.941 6 061 3 6 NA 12.9
> > > NA AA000A002 20160201_120000 8.832 40.941 6 061 3 6 NA 13.1
> > > NA AA000A002 20160201_150000 8.832 40.941 6 061 3 6 NA 13.1
> > > NA AA000A002 20160201_180000 8.832 40.941 6 061 3 6 NA 14.0
> > > NA AA000A002 20160201_210000 8.832 40.941 6 061 3 6 NA 15.5
> > > NA AA000A002 20160202_000000 8.832 40.941 6 061 3 6 NA 16.1
> > >
> > > Afterward I created an nc file using the tool
> > > ./ascii2nc ascii_obs.txt ascii_obs.nc -v 2
> > >
> > > The output of "ncdump -v obs_arr ascii_obs.nc" looks fine.
> > >
> > > I modified PointStatConfig_default like this
> > >
> > > fcst = {
> > > sid_exc = [];
> > >
> > > field = [
> > > {
> > > name = "APCP";
> > > level = [ "A03" ];
> > > cat_thresh = [ >1.0 ];
> > > file_type = [ 1 ] ;
> > > init_time = [ "20160201_000000" ];
> > > valid_time = [ "20160202_000000" ];
> > > }
> > > ];
> > >
> > > }
> > > obs = fcst;
> > >
> > > and using "./point_stat wrfout.grb ascii_obs.nc
> > PointStatConfig_modified
> > > -outdir ./ -v 2"
> > > the program just give me the headers of the statistic files and
the
> > > following message error
> > >
> > > DEBUG 1: User Config File: PointStatConfig_modified
> > > GSL_RNG_TYPE=mt19937
> > > GSL_RNG_SEED=1087746390
> > > DEBUG 1: Forecast File: wrfout_gfs.grib1
> > > DEBUG 1: Observation File: prova_ascii_grib1.nc
> > > DEBUG 2:
> > > DEBUG 2:
> > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > DEBUG 2:
> > > DEBUG 2: Reading data for APCP/A03.
> > > WARNING:
> > > WARNING: process_fcst_climo_files() -> no fields matching
APCP/A03
> found
> > in
> > > file: wrfout_gfs.grib1
> > > WARNING:
> > > .....
> > >
> > >
> > > I also tried to modify the input file like this
> > >
> > > fcst = {
> > > sid_exc = [];
> > >
> > > field = [
> > > {
> > > name = "APCP_03";
> > > level = [ "(*,*)" ];
> > > cat_thresh = [ >1.0 ];
> > > file_type = [ 1 ] ;
> > > init_time = [ "20160201_000000" ];
> > > valid_time = [ "20160202_000000" ];
> > > }
> > > ];
> > >
> > > }
> > >
> > > obs = {
> > > message_type = [ "ADPSFC" ];
> > > sid_exc = [];
> > >
> > > field = [
> > > {
> > > name = "APCP";
> > > level = [ "A03" ];
> > > cat_thresh = [ >1.0 ];
> > > file_type = [ 1 ] ;
> > > init_time = [ "20160201_000000" ];
> > > valid_time = [ "20160202_000000" ];
> > > }
> > > ];
> > >
> > > }
> > >
> > > but in this case the message error is even worse and I get
> > >
> > > DEBUG 1: User Config File: PointStatConfig_modified
> > > ERROR :
> > > ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
> abbreviation
> > > 'APCP_03' for table version 2
> > > ERROR :
> > >
> > > without any statistic output file.
> > >
> > >
> > > Any advise is really welcome.
> > >
> > > Thanks,
> > > Samuele
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: Troubles with point_stat
From: Samuele Salis
Time: Thu Feb 04 10:10:49 2016
Ok John,
I uploaded my files.
Running a level 4 for verbosity I get
DEBUG 2: Processing ACPCP/L0 versus ACPCP/L0, for observation type
ADPSFC,
over region FULL, for interpolation method NEAREST(1), using 0 pairs.
DEBUG 3: Number of matched pairs = 0
DEBUG 3: Observations processed = 16
DEBUG 3: Rejected: SID exclusion = 0
DEBUG 3: Rejected: GRIB code = 0
DEBUG 3: Rejected: valid time = 16
DEBUG 3: Rejected: bad obs value = 0
DEBUG 3: Rejected: off the grid = 0
DEBUG 3: Rejected: level mismatch = 0
DEBUG 3: Rejected: quality marker = 0
DEBUG 3: Rejected: message type = 0
DEBUG 3: Rejected: masking region = 0
DEBUG 3: Rejected: bad fcst value = 0
DEBUG 3: Rejected: duplicates = 0
I see all valid times are rejected, I guess something is wrong there,
but
at the moment I cannot solve the problem.
Thanks,
Samuele
On 4 February 2016 at 17:36, John Halley Gotway via RT
<met_help at ucar.edu>
wrote:
> Samuele,
>
> When you run Point-Stat, try increasing the verbosity level to 3 or
4 (-v
> 4) to get more diagnostic information about why counts were or were
not
> used.
>
> At this point, I think it'd be easiest to just have you send me some
sample
> data. I'll run it here and provide you some suggestions.
>
> Please post ascii_obs.txt, wrfout.grb, and your Point-Stat
configuration
> files to our anonymous ftp site following these directions:
> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> Please let me know when you've posted these files, and I'll go grab
them.
>
> Thanks,
> John
>
>
> On Thu, Feb 4, 2016 at 2:58 AM, Samuele Salis via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966 >
> >
> > Dear John,
> >
> > first of all thank you for your kind reply.
> >
> > The main error I found is that the APCP field is not included in
my
> > wrfout.grb file, but I saw the ACPCP parameter, so I changed my
> > ascii_obs.txt file like this
> >
> > ADPSFC AA000A001 20160201_030000 9.458 39.961 813 063 3 813 NA
2.6
> > ADPSFC AA000A001 20160201_060000 9.458 39.961 813 063 6 813 NA
2.6
> > ADPSFC AA000A001 20160201_090000 9.458 39.961 813 063 9 813 NA
2.9
> > ADPSFC AA000A001 20160201_120000 9.458 39.961 813 063 12 813 NA
3.6
> > ADPSFC AA000A001 20160201_150000 9.458 39.961 813 063 15 813 NA
7.1
> > ADPSFC AA000A001 20160201_180000 9.458 39.961 813 063 18 813 NA
12.0
> > ADPSFC AA000A001 20160201_210000 9.458 39.961 813 063 21 813 NA
12.5
> > ADPSFC AA000A001 20160202_000000 9.458 39.961 813 063 24 813 NA
14.2
> > ADPSFC AA000A002 20160201_030000 8.832 40.941 6 063 3 6 NA
4.2
> > ADPSFC AA000A002 20160201_060000 8.832 40.941 6 063 6 6 NA
5.6
> > ADPSFC AA000A002 20160201_090000 8.832 40.941 6 063 9 6 NA
12.9
> > ADPSFC AA000A002 20160201_120000 8.832 40.941 6 063 12 6 NA
13.1
> > ADPSFC AA000A002 20160201_150000 8.832 40.941 6 063 15 6 NA
13.1
> > ADPSFC AA000A002 20160201_180000 8.832 40.941 6 063 18 6 NA
14.0
> > ADPSFC AA000A002 20160201_210000 8.832 40.941 6 063 21 6 NA
15.5
> > ADPSFC AA000A002 20160202_000000 8.832 40.941 6 063 24 6 NA
16.1
> >
> > and I generated again the ascii_obs.nc using ascii2nc.
> >
> > I modified the PointStatConfig_default file in this way
> >
> > fcst = {
> > wind_thresh = [ NA ];
> > message_type = [ "ADPSFC" ];
> > sid_exc = [];
> >
> > field = [
> > {
> > name = "ACPCP";
> > level = [ "L000" ];
> > cat_thresh = [ >1.0 ];
> > }
> > ];
> >
> > }
> > obs = fcst;
> >
> > and with the command "./point_stat wrfout.grb ascii_obs.nc
> > PointStatConfig_modified -outdir ./ -v 2"
> > I obtained apparently no error messages
> >
> > DEBUG 1: User Config File: PointStatConfig_modified
> > GSL_RNG_TYPE=mt19937
> > GSL_RNG_SEED=18446744072318923606
> > DEBUG 1: Forecast File: wrfout.grb
> > DEBUG 1: Observation File: ascii_obs.nc
> > DEBUG 2:
> > DEBUG 2:
> >
> >
>
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Reading data for ACPCP/L0.
> > DEBUG 2: For ACPCP/L0 found 25 forecast levels and 0 climatology
levels.
> > DEBUG 2:
> > DEBUG 2:
> >
> >
>
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Searching 16 observations from 16 messages.
> > DEBUG 2:
> > DEBUG 2:
> >
> >
>
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Processing ACPCP/L0 versus ACPCP/L0, for observation type
> ADPSFC,
> > over region FULL, for interpolation method NEAREST(1), using 0
pairs.
> > DEBUG 2:
> >
> > The problem is that the produced statistical point_stat files are
all
> empty
> > except the headers.
> >
> > I am worry the problem is in my wrfout.grb file, so I will try
with
> another
> > one.
> >
> >
> > Thanks,
> > Samuele
> >
> >
> >
> > On 3 February 2016 at 18:48, John Halley Gotway via RT <
> met_help at ucar.edu>
> > wrote:
> >
> > > Samuele,
> > >
> > > I see you are having trouble running the Point-Stat tool. The
relevant
> > > error here is this:
> > >
> > > ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
> abbreviation
> > > 'APCP_03' for table version 2
> > >
> > > You are passing Point-Stat a forecast GRIB file (wrfout.grb) but
are
> > > specifying in the config file as if it were a gridded NetCDF
file.
> Let's
> > > back up and use wgrib to inventory your GRIB file and see what
> > > accumulations intervals are present in there:
> > >
> > > wgrib wrfout.grb | grep APCP
> > >
> > > What do you see? Do you actually have 3-hourly accumulated
precip? If
> > > not, what accumulation interval do you have?
> > >
> > > Also, please edit your point observations by specifying a
message type
> in
> > > the first column. I'd suggest replacing the NA string with
ADPSFC,
> which
> > > is PREPBUFR message type convention for surface observations.
> > >
> > > Assuming you actually have 3-hourly APCP in your GRIB file, use
this in
> > > your Point-Stat config file:
> > >
> > > fcst = {
> > > sid_exc = [];
> > > message_type = [ "ADPSFC" ];
> > >
> > > field = [
> > > {
> > > name = "APCP";
> > > level = [ "A03" ];
> > > cat_thresh = [ >1.0 ];
> > > }
> > > ];
> > > }
> > > obs = fcst;
> > >
> > >
> > > Thanks,
> > > John Halley Gotway
> > >
> > > On Wed, Feb 3, 2016 at 3:43 AM, Samuele Salis via RT <
> met_help at ucar.edu>
> > > wrote:
> > >
> > > >
> > > > Wed Feb 03 03:43:55 2016: Request 74966 was acted upon.
> > > > Transaction: Ticket created by samuelesalis at gmail.com
> > > > Queue: met_help
> > > > Subject: Troubles with point_stat
> > > > Owner: Nobody
> > > > Requestors: samuelesalis at gmail.com
> > > > Status: new
> > > > Ticket <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966
> > >
> > > >
> > > >
> > > > Dear MET users,
> > > >
> > > > I have a problem with, I think, the point_stat tool.
> > > >
> > > > Here a brief history of what I did, just to check if I made
some
> former
> > > > mistake.
> > > >
> > > > I made a 24 h simulation with ARW initialised with GFS and
starting
> > > > 2016020100.
> > > > I setted the io_form_history=5 to get a grib1 wrfout file
(called
> > > > wrfout.grb)
> > > >
> > > > I created an ascii file containing the accumulated
precipitation in
> > three
> > > > hours for
> > > > the same dates filling the eleven fields described in table
4.4 of
> > > > Met_user_guide_v5.1
> > > > just for two stations.
> > > > Here is the file (ascii_obs.txt)
> > > > NA AA000A001 20160201_030000 9.458 39.961 813 061 3 813 NA
2.6
> > > > NA AA000A001 20160201_060000 9.458 39.961 813 061 3 813 NA
2.6
> > > > NA AA000A001 20160201_090000 9.458 39.961 813 061 3 813 NA
2.9
> > > > NA AA000A001 20160201_120000 9.458 39.961 813 061 3 813 NA
3.6
> > > > NA AA000A001 20160201_150000 9.458 39.961 813 061 3 813 NA
7.1
> > > > NA AA000A001 20160201_180000 9.458 39.961 813 061 3 813 NA
12.0
> > > > NA AA000A001 20160201_210000 9.458 39.961 813 061 3 813 NA
12.5
> > > > NA AA000A001 20160202_000000 9.458 39.961 813 061 3 813 NA
14.2
> > > > NA AA000A002 20160201_030000 8.832 40.941 6 061 3 6 NA
4.2
> > > > NA AA000A002 20160201_060000 8.832 40.941 6 061 3 6 NA
5.6
> > > > NA AA000A002 20160201_090000 8.832 40.941 6 061 3 6 NA
12.9
> > > > NA AA000A002 20160201_120000 8.832 40.941 6 061 3 6 NA
13.1
> > > > NA AA000A002 20160201_150000 8.832 40.941 6 061 3 6 NA
13.1
> > > > NA AA000A002 20160201_180000 8.832 40.941 6 061 3 6 NA
14.0
> > > > NA AA000A002 20160201_210000 8.832 40.941 6 061 3 6 NA
15.5
> > > > NA AA000A002 20160202_000000 8.832 40.941 6 061 3 6 NA
16.1
> > > >
> > > > Afterward I created an nc file using the tool
> > > > ./ascii2nc ascii_obs.txt ascii_obs.nc -v 2
> > > >
> > > > The output of "ncdump -v obs_arr ascii_obs.nc" looks fine.
> > > >
> > > > I modified PointStatConfig_default like this
> > > >
> > > > fcst = {
> > > > sid_exc = [];
> > > >
> > > > field = [
> > > > {
> > > > name = "APCP";
> > > > level = [ "A03" ];
> > > > cat_thresh = [ >1.0 ];
> > > > file_type = [ 1 ] ;
> > > > init_time = [ "20160201_000000" ];
> > > > valid_time = [ "20160202_000000" ];
> > > > }
> > > > ];
> > > >
> > > > }
> > > > obs = fcst;
> > > >
> > > > and using "./point_stat wrfout.grb ascii_obs.nc
> > > PointStatConfig_modified
> > > > -outdir ./ -v 2"
> > > > the program just give me the headers of the statistic files
and the
> > > > following message error
> > > >
> > > > DEBUG 1: User Config File: PointStatConfig_modified
> > > > GSL_RNG_TYPE=mt19937
> > > > GSL_RNG_SEED=1087746390
> > > > DEBUG 1: Forecast File: wrfout_gfs.grib1
> > > > DEBUG 1: Observation File: prova_ascii_grib1.nc
> > > > DEBUG 2:
> > > > DEBUG 2:
> > > >
> > > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > > DEBUG 2:
> > > > DEBUG 2: Reading data for APCP/A03.
> > > > WARNING:
> > > > WARNING: process_fcst_climo_files() -> no fields matching
APCP/A03
> > found
> > > in
> > > > file: wrfout_gfs.grib1
> > > > WARNING:
> > > > .....
> > > >
> > > >
> > > > I also tried to modify the input file like this
> > > >
> > > > fcst = {
> > > > sid_exc = [];
> > > >
> > > > field = [
> > > > {
> > > > name = "APCP_03";
> > > > level = [ "(*,*)" ];
> > > > cat_thresh = [ >1.0 ];
> > > > file_type = [ 1 ] ;
> > > > init_time = [ "20160201_000000" ];
> > > > valid_time = [ "20160202_000000" ];
> > > > }
> > > > ];
> > > >
> > > > }
> > > >
> > > > obs = {
> > > > message_type = [ "ADPSFC" ];
> > > > sid_exc = [];
> > > >
> > > > field = [
> > > > {
> > > > name = "APCP";
> > > > level = [ "A03" ];
> > > > cat_thresh = [ >1.0 ];
> > > > file_type = [ 1 ] ;
> > > > init_time = [ "20160201_000000" ];
> > > > valid_time = [ "20160202_000000" ];
> > > > }
> > > > ];
> > > >
> > > > }
> > > >
> > > > but in this case the message error is even worse and I get
> > > >
> > > > DEBUG 1: User Config File: PointStatConfig_modified
> > > > ERROR :
> > > > ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
> > abbreviation
> > > > 'APCP_03' for table version 2
> > > > ERROR :
> > > >
> > > > without any statistic output file.
> > > >
> > > >
> > > > Any advise is really welcome.
> > > >
> > > > Thanks,
> > > > Samuele
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: Troubles with point_stat
From: John Halley Gotway
Time: Thu Feb 04 10:54:31 2016
Samuele,
Great! I'm glad you were able to diagnose that behavior. I see that
the
timestamps of your observations range from 20160201 to 20160202.
I don't know the actual timestamp of your forecast GRIB data, but you
should be able to see it by running wgrib on it (add -V for verbose
output).
The matching time window for Point-Stat is configurable in the
"obs_window"
section of the config file. The beginning and ending times of the
window
are defined as a number of seconds relative to the forecast valid
time.
For example, there are 3600 seconds in an hour, so the following would
match any observations +/- 1 hour around the forecast valid time:
obs_window = {
beg = -3600;
end = 3600;
}
Please try setting your time window accordingly and rerun.
Thanks,
John
On Thu, Feb 4, 2016 at 10:10 AM, Samuele Salis via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966 >
>
> Ok John,
>
> I uploaded my files.
>
> Running a level 4 for verbosity I get
>
> DEBUG 2: Processing ACPCP/L0 versus ACPCP/L0, for observation type
ADPSFC,
> over region FULL, for interpolation method NEAREST(1), using 0
pairs.
> DEBUG 3: Number of matched pairs = 0
> DEBUG 3: Observations processed = 16
> DEBUG 3: Rejected: SID exclusion = 0
> DEBUG 3: Rejected: GRIB code = 0
> DEBUG 3: Rejected: valid time = 16
> DEBUG 3: Rejected: bad obs value = 0
> DEBUG 3: Rejected: off the grid = 0
> DEBUG 3: Rejected: level mismatch = 0
> DEBUG 3: Rejected: quality marker = 0
> DEBUG 3: Rejected: message type = 0
> DEBUG 3: Rejected: masking region = 0
> DEBUG 3: Rejected: bad fcst value = 0
> DEBUG 3: Rejected: duplicates = 0
>
> I see all valid times are rejected, I guess something is wrong
there, but
> at the moment I cannot solve the problem.
>
> Thanks,
> Samuele
>
> On 4 February 2016 at 17:36, John Halley Gotway via RT
<met_help at ucar.edu>
> wrote:
>
> > Samuele,
> >
> > When you run Point-Stat, try increasing the verbosity level to 3
or 4 (-v
> > 4) to get more diagnostic information about why counts were or
were not
> > used.
> >
> > At this point, I think it'd be easiest to just have you send me
some
> sample
> > data. I'll run it here and provide you some suggestions.
> >
> > Please post ascii_obs.txt, wrfout.grb, and your Point-Stat
configuration
> > files to our anonymous ftp site following these directions:
> > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >
> > Please let me know when you've posted these files, and I'll go
grab them.
> >
> > Thanks,
> > John
> >
> >
> > On Thu, Feb 4, 2016 at 2:58 AM, Samuele Salis via RT
<met_help at ucar.edu>
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966 >
> > >
> > > Dear John,
> > >
> > > first of all thank you for your kind reply.
> > >
> > > The main error I found is that the APCP field is not included in
my
> > > wrfout.grb file, but I saw the ACPCP parameter, so I changed my
> > > ascii_obs.txt file like this
> > >
> > > ADPSFC AA000A001 20160201_030000 9.458 39.961 813 063 3 813 NA
2.6
> > > ADPSFC AA000A001 20160201_060000 9.458 39.961 813 063 6 813 NA
2.6
> > > ADPSFC AA000A001 20160201_090000 9.458 39.961 813 063 9 813 NA
2.9
> > > ADPSFC AA000A001 20160201_120000 9.458 39.961 813 063 12 813 NA
3.6
> > > ADPSFC AA000A001 20160201_150000 9.458 39.961 813 063 15 813 NA
7.1
> > > ADPSFC AA000A001 20160201_180000 9.458 39.961 813 063 18 813 NA
12.0
> > > ADPSFC AA000A001 20160201_210000 9.458 39.961 813 063 21 813 NA
12.5
> > > ADPSFC AA000A001 20160202_000000 9.458 39.961 813 063 24 813 NA
14.2
> > > ADPSFC AA000A002 20160201_030000 8.832 40.941 6 063 3 6 NA
4.2
> > > ADPSFC AA000A002 20160201_060000 8.832 40.941 6 063 6 6 NA
5.6
> > > ADPSFC AA000A002 20160201_090000 8.832 40.941 6 063 9 6 NA
12.9
> > > ADPSFC AA000A002 20160201_120000 8.832 40.941 6 063 12 6 NA
13.1
> > > ADPSFC AA000A002 20160201_150000 8.832 40.941 6 063 15 6 NA
13.1
> > > ADPSFC AA000A002 20160201_180000 8.832 40.941 6 063 18 6 NA
14.0
> > > ADPSFC AA000A002 20160201_210000 8.832 40.941 6 063 21 6 NA
15.5
> > > ADPSFC AA000A002 20160202_000000 8.832 40.941 6 063 24 6 NA
16.1
> > >
> > > and I generated again the ascii_obs.nc using ascii2nc.
> > >
> > > I modified the PointStatConfig_default file in this way
> > >
> > > fcst = {
> > > wind_thresh = [ NA ];
> > > message_type = [ "ADPSFC" ];
> > > sid_exc = [];
> > >
> > > field = [
> > > {
> > > name = "ACPCP";
> > > level = [ "L000" ];
> > > cat_thresh = [ >1.0 ];
> > > }
> > > ];
> > >
> > > }
> > > obs = fcst;
> > >
> > > and with the command "./point_stat wrfout.grb ascii_obs.nc
> > > PointStatConfig_modified -outdir ./ -v 2"
> > > I obtained apparently no error messages
> > >
> > > DEBUG 1: User Config File: PointStatConfig_modified
> > > GSL_RNG_TYPE=mt19937
> > > GSL_RNG_SEED=18446744072318923606
> > > DEBUG 1: Forecast File: wrfout.grb
> > > DEBUG 1: Observation File: ascii_obs.nc
> > > DEBUG 2:
> > > DEBUG 2:
> > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > DEBUG 2:
> > > DEBUG 2: Reading data for ACPCP/L0.
> > > DEBUG 2: For ACPCP/L0 found 25 forecast levels and 0 climatology
> levels.
> > > DEBUG 2:
> > > DEBUG 2:
> > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > DEBUG 2:
> > > DEBUG 2: Searching 16 observations from 16 messages.
> > > DEBUG 2:
> > > DEBUG 2:
> > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > DEBUG 2:
> > > DEBUG 2: Processing ACPCP/L0 versus ACPCP/L0, for observation
type
> > ADPSFC,
> > > over region FULL, for interpolation method NEAREST(1), using 0
pairs.
> > > DEBUG 2:
> > >
> > > The problem is that the produced statistical point_stat files
are all
> > empty
> > > except the headers.
> > >
> > > I am worry the problem is in my wrfout.grb file, so I will try
with
> > another
> > > one.
> > >
> > >
> > > Thanks,
> > > Samuele
> > >
> > >
> > >
> > > On 3 February 2016 at 18:48, John Halley Gotway via RT <
> > met_help at ucar.edu>
> > > wrote:
> > >
> > > > Samuele,
> > > >
> > > > I see you are having trouble running the Point-Stat tool. The
> relevant
> > > > error here is this:
> > > >
> > > > ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
> > abbreviation
> > > > 'APCP_03' for table version 2
> > > >
> > > > You are passing Point-Stat a forecast GRIB file (wrfout.grb)
but are
> > > > specifying in the config file as if it were a gridded NetCDF
file.
> > Let's
> > > > back up and use wgrib to inventory your GRIB file and see what
> > > > accumulations intervals are present in there:
> > > >
> > > > wgrib wrfout.grb | grep APCP
> > > >
> > > > What do you see? Do you actually have 3-hourly accumulated
precip?
> If
> > > > not, what accumulation interval do you have?
> > > >
> > > > Also, please edit your point observations by specifying a
message
> type
> > in
> > > > the first column. I'd suggest replacing the NA string with
ADPSFC,
> > which
> > > > is PREPBUFR message type convention for surface observations.
> > > >
> > > > Assuming you actually have 3-hourly APCP in your GRIB file,
use this
> in
> > > > your Point-Stat config file:
> > > >
> > > > fcst = {
> > > > sid_exc = [];
> > > > message_type = [ "ADPSFC" ];
> > > >
> > > > field = [
> > > > {
> > > > name = "APCP";
> > > > level = [ "A03" ];
> > > > cat_thresh = [ >1.0 ];
> > > > }
> > > > ];
> > > > }
> > > > obs = fcst;
> > > >
> > > >
> > > > Thanks,
> > > > John Halley Gotway
> > > >
> > > > On Wed, Feb 3, 2016 at 3:43 AM, Samuele Salis via RT <
> > met_help at ucar.edu>
> > > > wrote:
> > > >
> > > > >
> > > > > Wed Feb 03 03:43:55 2016: Request 74966 was acted upon.
> > > > > Transaction: Ticket created by samuelesalis at gmail.com
> > > > > Queue: met_help
> > > > > Subject: Troubles with point_stat
> > > > > Owner: Nobody
> > > > > Requestors: samuelesalis at gmail.com
> > > > > Status: new
> > > > > Ticket <URL:
> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966
> > > >
> > > > >
> > > > >
> > > > > Dear MET users,
> > > > >
> > > > > I have a problem with, I think, the point_stat tool.
> > > > >
> > > > > Here a brief history of what I did, just to check if I made
some
> > former
> > > > > mistake.
> > > > >
> > > > > I made a 24 h simulation with ARW initialised with GFS and
starting
> > > > > 2016020100.
> > > > > I setted the io_form_history=5 to get a grib1 wrfout file
(called
> > > > > wrfout.grb)
> > > > >
> > > > > I created an ascii file containing the accumulated
precipitation in
> > > three
> > > > > hours for
> > > > > the same dates filling the eleven fields described in table
4.4 of
> > > > > Met_user_guide_v5.1
> > > > > just for two stations.
> > > > > Here is the file (ascii_obs.txt)
> > > > > NA AA000A001 20160201_030000 9.458 39.961 813 061 3 813 NA
2.6
> > > > > NA AA000A001 20160201_060000 9.458 39.961 813 061 3 813 NA
2.6
> > > > > NA AA000A001 20160201_090000 9.458 39.961 813 061 3 813 NA
2.9
> > > > > NA AA000A001 20160201_120000 9.458 39.961 813 061 3 813 NA
3.6
> > > > > NA AA000A001 20160201_150000 9.458 39.961 813 061 3 813 NA
7.1
> > > > > NA AA000A001 20160201_180000 9.458 39.961 813 061 3 813 NA
12.0
> > > > > NA AA000A001 20160201_210000 9.458 39.961 813 061 3 813 NA
12.5
> > > > > NA AA000A001 20160202_000000 9.458 39.961 813 061 3 813 NA
14.2
> > > > > NA AA000A002 20160201_030000 8.832 40.941 6 061 3 6 NA
4.2
> > > > > NA AA000A002 20160201_060000 8.832 40.941 6 061 3 6 NA
5.6
> > > > > NA AA000A002 20160201_090000 8.832 40.941 6 061 3 6 NA
12.9
> > > > > NA AA000A002 20160201_120000 8.832 40.941 6 061 3 6 NA
13.1
> > > > > NA AA000A002 20160201_150000 8.832 40.941 6 061 3 6 NA
13.1
> > > > > NA AA000A002 20160201_180000 8.832 40.941 6 061 3 6 NA
14.0
> > > > > NA AA000A002 20160201_210000 8.832 40.941 6 061 3 6 NA
15.5
> > > > > NA AA000A002 20160202_000000 8.832 40.941 6 061 3 6 NA
16.1
> > > > >
> > > > > Afterward I created an nc file using the tool
> > > > > ./ascii2nc ascii_obs.txt ascii_obs.nc -v 2
> > > > >
> > > > > The output of "ncdump -v obs_arr ascii_obs.nc" looks fine.
> > > > >
> > > > > I modified PointStatConfig_default like this
> > > > >
> > > > > fcst = {
> > > > > sid_exc = [];
> > > > >
> > > > > field = [
> > > > > {
> > > > > name = "APCP";
> > > > > level = [ "A03" ];
> > > > > cat_thresh = [ >1.0 ];
> > > > > file_type = [ 1 ] ;
> > > > > init_time = [ "20160201_000000" ];
> > > > > valid_time = [ "20160202_000000" ];
> > > > > }
> > > > > ];
> > > > >
> > > > > }
> > > > > obs = fcst;
> > > > >
> > > > > and using "./point_stat wrfout.grb ascii_obs.nc
> > > > PointStatConfig_modified
> > > > > -outdir ./ -v 2"
> > > > > the program just give me the headers of the statistic files
and the
> > > > > following message error
> > > > >
> > > > > DEBUG 1: User Config File: PointStatConfig_modified
> > > > > GSL_RNG_TYPE=mt19937
> > > > > GSL_RNG_SEED=1087746390
> > > > > DEBUG 1: Forecast File: wrfout_gfs.grib1
> > > > > DEBUG 1: Observation File: prova_ascii_grib1.nc
> > > > > DEBUG 2:
> > > > > DEBUG 2:
> > > > >
> > > > >
> > > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > > > DEBUG 2:
> > > > > DEBUG 2: Reading data for APCP/A03.
> > > > > WARNING:
> > > > > WARNING: process_fcst_climo_files() -> no fields matching
APCP/A03
> > > found
> > > > in
> > > > > file: wrfout_gfs.grib1
> > > > > WARNING:
> > > > > .....
> > > > >
> > > > >
> > > > > I also tried to modify the input file like this
> > > > >
> > > > > fcst = {
> > > > > sid_exc = [];
> > > > >
> > > > > field = [
> > > > > {
> > > > > name = "APCP_03";
> > > > > level = [ "(*,*)" ];
> > > > > cat_thresh = [ >1.0 ];
> > > > > file_type = [ 1 ] ;
> > > > > init_time = [ "20160201_000000" ];
> > > > > valid_time = [ "20160202_000000" ];
> > > > > }
> > > > > ];
> > > > >
> > > > > }
> > > > >
> > > > > obs = {
> > > > > message_type = [ "ADPSFC" ];
> > > > > sid_exc = [];
> > > > >
> > > > > field = [
> > > > > {
> > > > > name = "APCP";
> > > > > level = [ "A03" ];
> > > > > cat_thresh = [ >1.0 ];
> > > > > file_type = [ 1 ] ;
> > > > > init_time = [ "20160201_000000" ];
> > > > > valid_time = [ "20160202_000000" ];
> > > > > }
> > > > > ];
> > > > >
> > > > > }
> > > > >
> > > > > but in this case the message error is even worse and I get
> > > > >
> > > > > DEBUG 1: User Config File: PointStatConfig_modified
> > > > > ERROR :
> > > > > ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
> > > abbreviation
> > > > > 'APCP_03' for table version 2
> > > > > ERROR :
> > > > >
> > > > > without any statistic output file.
> > > > >
> > > > >
> > > > > Any advise is really welcome.
> > > > >
> > > > > Thanks,
> > > > > Samuele
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: Troubles with point_stat
From: Samuele Salis
Time: Fri Feb 05 04:13:44 2016
Ok John,
thanks for your helping.
I enlarged the obs_window until all the data have been included.
Now I get an output for every statistical processing. To check the
goodness of my output now is another step and to do that I have
a lot to study.
Thanks,
Samuele
On 4 February 2016 at 18:54, John Halley Gotway via RT
<met_help at ucar.edu>
wrote:
> Samuele,
>
> Great! I'm glad you were able to diagnose that behavior. I see
that the
> timestamps of your observations range from 20160201 to 20160202.
>
> I don't know the actual timestamp of your forecast GRIB data, but
you
> should be able to see it by running wgrib on it (add -V for verbose
> output).
>
> The matching time window for Point-Stat is configurable in the
"obs_window"
> section of the config file. The beginning and ending times of the
window
> are defined as a number of seconds relative to the forecast valid
time.
> For example, there are 3600 seconds in an hour, so the following
would
> match any observations +/- 1 hour around the forecast valid time:
>
> obs_window = {
> beg = -3600;
> end = 3600;
> }
>
> Please try setting your time window accordingly and rerun.
>
> Thanks,
> John
>
> On Thu, Feb 4, 2016 at 10:10 AM, Samuele Salis via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966 >
> >
> > Ok John,
> >
> > I uploaded my files.
> >
> > Running a level 4 for verbosity I get
> >
> > DEBUG 2: Processing ACPCP/L0 versus ACPCP/L0, for observation type
> ADPSFC,
> > over region FULL, for interpolation method NEAREST(1), using 0
pairs.
> > DEBUG 3: Number of matched pairs = 0
> > DEBUG 3: Observations processed = 16
> > DEBUG 3: Rejected: SID exclusion = 0
> > DEBUG 3: Rejected: GRIB code = 0
> > DEBUG 3: Rejected: valid time = 16
> > DEBUG 3: Rejected: bad obs value = 0
> > DEBUG 3: Rejected: off the grid = 0
> > DEBUG 3: Rejected: level mismatch = 0
> > DEBUG 3: Rejected: quality marker = 0
> > DEBUG 3: Rejected: message type = 0
> > DEBUG 3: Rejected: masking region = 0
> > DEBUG 3: Rejected: bad fcst value = 0
> > DEBUG 3: Rejected: duplicates = 0
> >
> > I see all valid times are rejected, I guess something is wrong
there, but
> > at the moment I cannot solve the problem.
> >
> > Thanks,
> > Samuele
> >
> > On 4 February 2016 at 17:36, John Halley Gotway via RT <
> met_help at ucar.edu>
> > wrote:
> >
> > > Samuele,
> > >
> > > When you run Point-Stat, try increasing the verbosity level to 3
or 4
> (-v
> > > 4) to get more diagnostic information about why counts were or
were not
> > > used.
> > >
> > > At this point, I think it'd be easiest to just have you send me
some
> > sample
> > > data. I'll run it here and provide you some suggestions.
> > >
> > > Please post ascii_obs.txt, wrfout.grb, and your Point-Stat
> configuration
> > > files to our anonymous ftp site following these directions:
> > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > >
> > > Please let me know when you've posted these files, and I'll go
grab
> them.
> > >
> > > Thanks,
> > > John
> > >
> > >
> > > On Thu, Feb 4, 2016 at 2:58 AM, Samuele Salis via RT <
> met_help at ucar.edu>
> > > wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966
>
> > > >
> > > > Dear John,
> > > >
> > > > first of all thank you for your kind reply.
> > > >
> > > > The main error I found is that the APCP field is not included
in my
> > > > wrfout.grb file, but I saw the ACPCP parameter, so I changed
my
> > > > ascii_obs.txt file like this
> > > >
> > > > ADPSFC AA000A001 20160201_030000 9.458 39.961 813 063 3 813
NA 2.6
> > > > ADPSFC AA000A001 20160201_060000 9.458 39.961 813 063 6 813
NA 2.6
> > > > ADPSFC AA000A001 20160201_090000 9.458 39.961 813 063 9 813
NA 2.9
> > > > ADPSFC AA000A001 20160201_120000 9.458 39.961 813 063 12 813
NA 3.6
> > > > ADPSFC AA000A001 20160201_150000 9.458 39.961 813 063 15 813
NA 7.1
> > > > ADPSFC AA000A001 20160201_180000 9.458 39.961 813 063 18 813
NA 12.0
> > > > ADPSFC AA000A001 20160201_210000 9.458 39.961 813 063 21 813
NA 12.5
> > > > ADPSFC AA000A001 20160202_000000 9.458 39.961 813 063 24 813
NA 14.2
> > > > ADPSFC AA000A002 20160201_030000 8.832 40.941 6 063 3 6
NA 4.2
> > > > ADPSFC AA000A002 20160201_060000 8.832 40.941 6 063 6 6
NA 5.6
> > > > ADPSFC AA000A002 20160201_090000 8.832 40.941 6 063 9 6
NA 12.9
> > > > ADPSFC AA000A002 20160201_120000 8.832 40.941 6 063 12 6
NA 13.1
> > > > ADPSFC AA000A002 20160201_150000 8.832 40.941 6 063 15 6
NA 13.1
> > > > ADPSFC AA000A002 20160201_180000 8.832 40.941 6 063 18 6
NA 14.0
> > > > ADPSFC AA000A002 20160201_210000 8.832 40.941 6 063 21 6
NA 15.5
> > > > ADPSFC AA000A002 20160202_000000 8.832 40.941 6 063 24 6
NA 16.1
> > > >
> > > > and I generated again the ascii_obs.nc using ascii2nc.
> > > >
> > > > I modified the PointStatConfig_default file in this way
> > > >
> > > > fcst = {
> > > > wind_thresh = [ NA ];
> > > > message_type = [ "ADPSFC" ];
> > > > sid_exc = [];
> > > >
> > > > field = [
> > > > {
> > > > name = "ACPCP";
> > > > level = [ "L000" ];
> > > > cat_thresh = [ >1.0 ];
> > > > }
> > > > ];
> > > >
> > > > }
> > > > obs = fcst;
> > > >
> > > > and with the command "./point_stat wrfout.grb ascii_obs.nc
> > > > PointStatConfig_modified -outdir ./ -v 2"
> > > > I obtained apparently no error messages
> > > >
> > > > DEBUG 1: User Config File: PointStatConfig_modified
> > > > GSL_RNG_TYPE=mt19937
> > > > GSL_RNG_SEED=18446744072318923606
> > > > DEBUG 1: Forecast File: wrfout.grb
> > > > DEBUG 1: Observation File: ascii_obs.nc
> > > > DEBUG 2:
> > > > DEBUG 2:
> > > >
> > > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > > DEBUG 2:
> > > > DEBUG 2: Reading data for ACPCP/L0.
> > > > DEBUG 2: For ACPCP/L0 found 25 forecast levels and 0
climatology
> > levels.
> > > > DEBUG 2:
> > > > DEBUG 2:
> > > >
> > > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > > DEBUG 2:
> > > > DEBUG 2: Searching 16 observations from 16 messages.
> > > > DEBUG 2:
> > > > DEBUG 2:
> > > >
> > > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > > DEBUG 2:
> > > > DEBUG 2: Processing ACPCP/L0 versus ACPCP/L0, for observation
type
> > > ADPSFC,
> > > > over region FULL, for interpolation method NEAREST(1), using 0
pairs.
> > > > DEBUG 2:
> > > >
> > > > The problem is that the produced statistical point_stat files
are all
> > > empty
> > > > except the headers.
> > > >
> > > > I am worry the problem is in my wrfout.grb file, so I will try
with
> > > another
> > > > one.
> > > >
> > > >
> > > > Thanks,
> > > > Samuele
> > > >
> > > >
> > > >
> > > > On 3 February 2016 at 18:48, John Halley Gotway via RT <
> > > met_help at ucar.edu>
> > > > wrote:
> > > >
> > > > > Samuele,
> > > > >
> > > > > I see you are having trouble running the Point-Stat tool.
The
> > relevant
> > > > > error here is this:
> > > > >
> > > > > ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1 field
> > > abbreviation
> > > > > 'APCP_03' for table version 2
> > > > >
> > > > > You are passing Point-Stat a forecast GRIB file (wrfout.grb)
but
> are
> > > > > specifying in the config file as if it were a gridded NetCDF
file.
> > > Let's
> > > > > back up and use wgrib to inventory your GRIB file and see
what
> > > > > accumulations intervals are present in there:
> > > > >
> > > > > wgrib wrfout.grb | grep APCP
> > > > >
> > > > > What do you see? Do you actually have 3-hourly accumulated
precip?
> > If
> > > > > not, what accumulation interval do you have?
> > > > >
> > > > > Also, please edit your point observations by specifying a
message
> > type
> > > in
> > > > > the first column. I'd suggest replacing the NA string with
ADPSFC,
> > > which
> > > > > is PREPBUFR message type convention for surface
observations.
> > > > >
> > > > > Assuming you actually have 3-hourly APCP in your GRIB file,
use
> this
> > in
> > > > > your Point-Stat config file:
> > > > >
> > > > > fcst = {
> > > > > sid_exc = [];
> > > > > message_type = [ "ADPSFC" ];
> > > > >
> > > > > field = [
> > > > > {
> > > > > name = "APCP";
> > > > > level = [ "A03" ];
> > > > > cat_thresh = [ >1.0 ];
> > > > > }
> > > > > ];
> > > > > }
> > > > > obs = fcst;
> > > > >
> > > > >
> > > > > Thanks,
> > > > > John Halley Gotway
> > > > >
> > > > > On Wed, Feb 3, 2016 at 3:43 AM, Samuele Salis via RT <
> > > met_help at ucar.edu>
> > > > > wrote:
> > > > >
> > > > > >
> > > > > > Wed Feb 03 03:43:55 2016: Request 74966 was acted upon.
> > > > > > Transaction: Ticket created by samuelesalis at gmail.com
> > > > > > Queue: met_help
> > > > > > Subject: Troubles with point_stat
> > > > > > Owner: Nobody
> > > > > > Requestors: samuelesalis at gmail.com
> > > > > > Status: new
> > > > > > Ticket <URL:
> > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966
> > > > >
> > > > > >
> > > > > >
> > > > > > Dear MET users,
> > > > > >
> > > > > > I have a problem with, I think, the point_stat tool.
> > > > > >
> > > > > > Here a brief history of what I did, just to check if I
made some
> > > former
> > > > > > mistake.
> > > > > >
> > > > > > I made a 24 h simulation with ARW initialised with GFS and
> starting
> > > > > > 2016020100.
> > > > > > I setted the io_form_history=5 to get a grib1 wrfout file
> (called
> > > > > > wrfout.grb)
> > > > > >
> > > > > > I created an ascii file containing the accumulated
precipitation
> in
> > > > three
> > > > > > hours for
> > > > > > the same dates filling the eleven fields described in
table 4.4
> of
> > > > > > Met_user_guide_v5.1
> > > > > > just for two stations.
> > > > > > Here is the file (ascii_obs.txt)
> > > > > > NA AA000A001 20160201_030000 9.458 39.961 813 061 3 813 NA
2.6
> > > > > > NA AA000A001 20160201_060000 9.458 39.961 813 061 3 813 NA
2.6
> > > > > > NA AA000A001 20160201_090000 9.458 39.961 813 061 3 813 NA
2.9
> > > > > > NA AA000A001 20160201_120000 9.458 39.961 813 061 3 813 NA
3.6
> > > > > > NA AA000A001 20160201_150000 9.458 39.961 813 061 3 813 NA
7.1
> > > > > > NA AA000A001 20160201_180000 9.458 39.961 813 061 3 813 NA
12.0
> > > > > > NA AA000A001 20160201_210000 9.458 39.961 813 061 3 813 NA
12.5
> > > > > > NA AA000A001 20160202_000000 9.458 39.961 813 061 3 813 NA
14.2
> > > > > > NA AA000A002 20160201_030000 8.832 40.941 6 061 3 6 NA
4.2
> > > > > > NA AA000A002 20160201_060000 8.832 40.941 6 061 3 6 NA
5.6
> > > > > > NA AA000A002 20160201_090000 8.832 40.941 6 061 3 6 NA
12.9
> > > > > > NA AA000A002 20160201_120000 8.832 40.941 6 061 3 6 NA
13.1
> > > > > > NA AA000A002 20160201_150000 8.832 40.941 6 061 3 6 NA
13.1
> > > > > > NA AA000A002 20160201_180000 8.832 40.941 6 061 3 6 NA
14.0
> > > > > > NA AA000A002 20160201_210000 8.832 40.941 6 061 3 6 NA
15.5
> > > > > > NA AA000A002 20160202_000000 8.832 40.941 6 061 3 6 NA
16.1
> > > > > >
> > > > > > Afterward I created an nc file using the tool
> > > > > > ./ascii2nc ascii_obs.txt ascii_obs.nc -v 2
> > > > > >
> > > > > > The output of "ncdump -v obs_arr ascii_obs.nc" looks fine.
> > > > > >
> > > > > > I modified PointStatConfig_default like this
> > > > > >
> > > > > > fcst = {
> > > > > > sid_exc = [];
> > > > > >
> > > > > > field = [
> > > > > > {
> > > > > > name = "APCP";
> > > > > > level = [ "A03" ];
> > > > > > cat_thresh = [ >1.0 ];
> > > > > > file_type = [ 1 ] ;
> > > > > > init_time = [ "20160201_000000" ];
> > > > > > valid_time = [ "20160202_000000" ];
> > > > > > }
> > > > > > ];
> > > > > >
> > > > > > }
> > > > > > obs = fcst;
> > > > > >
> > > > > > and using "./point_stat wrfout.grb ascii_obs.nc
> > > > > PointStatConfig_modified
> > > > > > -outdir ./ -v 2"
> > > > > > the program just give me the headers of the statistic
files and
> the
> > > > > > following message error
> > > > > >
> > > > > > DEBUG 1: User Config File: PointStatConfig_modified
> > > > > > GSL_RNG_TYPE=mt19937
> > > > > > GSL_RNG_SEED=1087746390
> > > > > > DEBUG 1: Forecast File: wrfout_gfs.grib1
> > > > > > DEBUG 1: Observation File: prova_ascii_grib1.nc
> > > > > > DEBUG 2:
> > > > > > DEBUG 2:
> > > > > >
> > > > > >
> > > > >
> > > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > > > > DEBUG 2:
> > > > > > DEBUG 2: Reading data for APCP/A03.
> > > > > > WARNING:
> > > > > > WARNING: process_fcst_climo_files() -> no fields matching
> APCP/A03
> > > > found
> > > > > in
> > > > > > file: wrfout_gfs.grib1
> > > > > > WARNING:
> > > > > > .....
> > > > > >
> > > > > >
> > > > > > I also tried to modify the input file like this
> > > > > >
> > > > > > fcst = {
> > > > > > sid_exc = [];
> > > > > >
> > > > > > field = [
> > > > > > {
> > > > > > name = "APCP_03";
> > > > > > level = [ "(*,*)" ];
> > > > > > cat_thresh = [ >1.0 ];
> > > > > > file_type = [ 1 ] ;
> > > > > > init_time = [ "20160201_000000" ];
> > > > > > valid_time = [ "20160202_000000" ];
> > > > > > }
> > > > > > ];
> > > > > >
> > > > > > }
> > > > > >
> > > > > > obs = {
> > > > > > message_type = [ "ADPSFC" ];
> > > > > > sid_exc = [];
> > > > > >
> > > > > > field = [
> > > > > > {
> > > > > > name = "APCP";
> > > > > > level = [ "A03" ];
> > > > > > cat_thresh = [ >1.0 ];
> > > > > > file_type = [ 1 ] ;
> > > > > > init_time = [ "20160201_000000" ];
> > > > > > valid_time = [ "20160202_000000" ];
> > > > > > }
> > > > > > ];
> > > > > >
> > > > > > }
> > > > > >
> > > > > > but in this case the message error is even worse and I get
> > > > > >
> > > > > > DEBUG 1: User Config File: PointStatConfig_modified
> > > > > > ERROR :
> > > > > > ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1
field
> > > > abbreviation
> > > > > > 'APCP_03' for table version 2
> > > > > > ERROR :
> > > > > >
> > > > > > without any statistic output file.
> > > > > >
> > > > > >
> > > > > > Any advise is really welcome.
> > > > > >
> > > > > > Thanks,
> > > > > > Samuele
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: Troubles with point_stat
From: John Halley Gotway
Time: Fri Feb 05 07:59:06 2016
Samuele,
Great, I'm glad you were able to figure it out. You should run point-
stat
once for each model output time. Typically, users have more than 2
observation locations at which they're verifying. And you'd run
configure
it to only use observations close to the valid time of the forecast.
Just let us know if more questions or issues arise.
Thanks
John
On Friday, February 5, 2016, Samuele Salis via RT <met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966 >
>
> Ok John,
>
> thanks for your helping.
>
> I enlarged the obs_window until all the data have been included.
>
> Now I get an output for every statistical processing. To check the
> goodness of my output now is another step and to do that I have
> a lot to study.
>
> Thanks,
> Samuele
>
> On 4 February 2016 at 18:54, John Halley Gotway via RT
<met_help at ucar.edu
> <javascript:;>>
> wrote:
>
> > Samuele,
> >
> > Great! I'm glad you were able to diagnose that behavior. I see
that the
> > timestamps of your observations range from 20160201 to 20160202.
> >
> > I don't know the actual timestamp of your forecast GRIB data, but
you
> > should be able to see it by running wgrib on it (add -V for
verbose
> > output).
> >
> > The matching time window for Point-Stat is configurable in the
> "obs_window"
> > section of the config file. The beginning and ending times of the
window
> > are defined as a number of seconds relative to the forecast valid
time.
> > For example, there are 3600 seconds in an hour, so the following
would
> > match any observations +/- 1 hour around the forecast valid time:
> >
> > obs_window = {
> > beg = -3600;
> > end = 3600;
> > }
> >
> > Please try setting your time window accordingly and rerun.
> >
> > Thanks,
> > John
> >
> > On Thu, Feb 4, 2016 at 10:10 AM, Samuele Salis via RT
<met_help at ucar.edu
> <javascript:;>>
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966 >
> > >
> > > Ok John,
> > >
> > > I uploaded my files.
> > >
> > > Running a level 4 for verbosity I get
> > >
> > > DEBUG 2: Processing ACPCP/L0 versus ACPCP/L0, for observation
type
> > ADPSFC,
> > > over region FULL, for interpolation method NEAREST(1), using 0
pairs.
> > > DEBUG 3: Number of matched pairs = 0
> > > DEBUG 3: Observations processed = 16
> > > DEBUG 3: Rejected: SID exclusion = 0
> > > DEBUG 3: Rejected: GRIB code = 0
> > > DEBUG 3: Rejected: valid time = 16
> > > DEBUG 3: Rejected: bad obs value = 0
> > > DEBUG 3: Rejected: off the grid = 0
> > > DEBUG 3: Rejected: level mismatch = 0
> > > DEBUG 3: Rejected: quality marker = 0
> > > DEBUG 3: Rejected: message type = 0
> > > DEBUG 3: Rejected: masking region = 0
> > > DEBUG 3: Rejected: bad fcst value = 0
> > > DEBUG 3: Rejected: duplicates = 0
> > >
> > > I see all valid times are rejected, I guess something is wrong
there,
> but
> > > at the moment I cannot solve the problem.
> > >
> > > Thanks,
> > > Samuele
> > >
> > > On 4 February 2016 at 17:36, John Halley Gotway via RT <
> > met_help at ucar.edu <javascript:;>>
> > > wrote:
> > >
> > > > Samuele,
> > > >
> > > > When you run Point-Stat, try increasing the verbosity level to
3 or 4
> > (-v
> > > > 4) to get more diagnostic information about why counts were or
were
> not
> > > > used.
> > > >
> > > > At this point, I think it'd be easiest to just have you send
me some
> > > sample
> > > > data. I'll run it here and provide you some suggestions.
> > > >
> > > > Please post ascii_obs.txt, wrfout.grb, and your Point-Stat
> > configuration
> > > > files to our anonymous ftp site following these directions:
> > > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > > >
> > > > Please let me know when you've posted these files, and I'll go
grab
> > them.
> > > >
> > > > Thanks,
> > > > John
> > > >
> > > >
> > > > On Thu, Feb 4, 2016 at 2:58 AM, Samuele Salis via RT <
> > met_help at ucar.edu <javascript:;>>
> > > > wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966 >
> > > > >
> > > > > Dear John,
> > > > >
> > > > > first of all thank you for your kind reply.
> > > > >
> > > > > The main error I found is that the APCP field is not
included in my
> > > > > wrfout.grb file, but I saw the ACPCP parameter, so I changed
my
> > > > > ascii_obs.txt file like this
> > > > >
> > > > > ADPSFC AA000A001 20160201_030000 9.458 39.961 813 063 3 813
NA
> 2.6
> > > > > ADPSFC AA000A001 20160201_060000 9.458 39.961 813 063 6 813
NA
> 2.6
> > > > > ADPSFC AA000A001 20160201_090000 9.458 39.961 813 063 9 813
NA
> 2.9
> > > > > ADPSFC AA000A001 20160201_120000 9.458 39.961 813 063 12 813
NA
> 3.6
> > > > > ADPSFC AA000A001 20160201_150000 9.458 39.961 813 063 15 813
NA
> 7.1
> > > > > ADPSFC AA000A001 20160201_180000 9.458 39.961 813 063 18 813
NA
> 12.0
> > > > > ADPSFC AA000A001 20160201_210000 9.458 39.961 813 063 21 813
NA
> 12.5
> > > > > ADPSFC AA000A001 20160202_000000 9.458 39.961 813 063 24 813
NA
> 14.2
> > > > > ADPSFC AA000A002 20160201_030000 8.832 40.941 6 063 3 6
NA
> 4.2
> > > > > ADPSFC AA000A002 20160201_060000 8.832 40.941 6 063 6 6
NA
> 5.6
> > > > > ADPSFC AA000A002 20160201_090000 8.832 40.941 6 063 9 6
NA
> 12.9
> > > > > ADPSFC AA000A002 20160201_120000 8.832 40.941 6 063 12 6
NA
> 13.1
> > > > > ADPSFC AA000A002 20160201_150000 8.832 40.941 6 063 15 6
NA
> 13.1
> > > > > ADPSFC AA000A002 20160201_180000 8.832 40.941 6 063 18 6
NA
> 14.0
> > > > > ADPSFC AA000A002 20160201_210000 8.832 40.941 6 063 21 6
NA
> 15.5
> > > > > ADPSFC AA000A002 20160202_000000 8.832 40.941 6 063 24 6
NA
> 16.1
> > > > >
> > > > > and I generated again the ascii_obs.nc using ascii2nc.
> > > > >
> > > > > I modified the PointStatConfig_default file in this way
> > > > >
> > > > > fcst = {
> > > > > wind_thresh = [ NA ];
> > > > > message_type = [ "ADPSFC" ];
> > > > > sid_exc = [];
> > > > >
> > > > > field = [
> > > > > {
> > > > > name = "ACPCP";
> > > > > level = [ "L000" ];
> > > > > cat_thresh = [ >1.0 ];
> > > > > }
> > > > > ];
> > > > >
> > > > > }
> > > > > obs = fcst;
> > > > >
> > > > > and with the command "./point_stat wrfout.grb ascii_obs.nc
> > > > > PointStatConfig_modified -outdir ./ -v 2"
> > > > > I obtained apparently no error messages
> > > > >
> > > > > DEBUG 1: User Config File: PointStatConfig_modified
> > > > > GSL_RNG_TYPE=mt19937
> > > > > GSL_RNG_SEED=18446744072318923606
> > > > > DEBUG 1: Forecast File: wrfout.grb
> > > > > DEBUG 1: Observation File: ascii_obs.nc
> > > > > DEBUG 2:
> > > > > DEBUG 2:
> > > > >
> > > > >
> > > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > > > DEBUG 2:
> > > > > DEBUG 2: Reading data for ACPCP/L0.
> > > > > DEBUG 2: For ACPCP/L0 found 25 forecast levels and 0
climatology
> > > levels.
> > > > > DEBUG 2:
> > > > > DEBUG 2:
> > > > >
> > > > >
> > > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > > > DEBUG 2:
> > > > > DEBUG 2: Searching 16 observations from 16 messages.
> > > > > DEBUG 2:
> > > > > DEBUG 2:
> > > > >
> > > > >
> > > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > > > DEBUG 2:
> > > > > DEBUG 2: Processing ACPCP/L0 versus ACPCP/L0, for
observation type
> > > > ADPSFC,
> > > > > over region FULL, for interpolation method NEAREST(1), using
0
> pairs.
> > > > > DEBUG 2:
> > > > >
> > > > > The problem is that the produced statistical point_stat
files are
> all
> > > > empty
> > > > > except the headers.
> > > > >
> > > > > I am worry the problem is in my wrfout.grb file, so I will
try with
> > > > another
> > > > > one.
> > > > >
> > > > >
> > > > > Thanks,
> > > > > Samuele
> > > > >
> > > > >
> > > > >
> > > > > On 3 February 2016 at 18:48, John Halley Gotway via RT <
> > > > met_help at ucar.edu <javascript:;>>
> > > > > wrote:
> > > > >
> > > > > > Samuele,
> > > > > >
> > > > > > I see you are having trouble running the Point-Stat tool.
The
> > > relevant
> > > > > > error here is this:
> > > > > >
> > > > > > ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1
field
> > > > abbreviation
> > > > > > 'APCP_03' for table version 2
> > > > > >
> > > > > > You are passing Point-Stat a forecast GRIB file
(wrfout.grb) but
> > are
> > > > > > specifying in the config file as if it were a gridded
NetCDF
> file.
> > > > Let's
> > > > > > back up and use wgrib to inventory your GRIB file and see
what
> > > > > > accumulations intervals are present in there:
> > > > > >
> > > > > > wgrib wrfout.grb | grep APCP
> > > > > >
> > > > > > What do you see? Do you actually have 3-hourly
accumulated
> precip?
> > > If
> > > > > > not, what accumulation interval do you have?
> > > > > >
> > > > > > Also, please edit your point observations by specifying a
message
> > > type
> > > > in
> > > > > > the first column. I'd suggest replacing the NA string
with
> ADPSFC,
> > > > which
> > > > > > is PREPBUFR message type convention for surface
observations.
> > > > > >
> > > > > > Assuming you actually have 3-hourly APCP in your GRIB
file, use
> > this
> > > in
> > > > > > your Point-Stat config file:
> > > > > >
> > > > > > fcst = {
> > > > > > sid_exc = [];
> > > > > > message_type = [ "ADPSFC" ];
> > > > > >
> > > > > > field = [
> > > > > > {
> > > > > > name = "APCP";
> > > > > > level = [ "A03" ];
> > > > > > cat_thresh = [ >1.0 ];
> > > > > > }
> > > > > > ];
> > > > > > }
> > > > > > obs = fcst;
> > > > > >
> > > > > >
> > > > > > Thanks,
> > > > > > John Halley Gotway
> > > > > >
> > > > > > On Wed, Feb 3, 2016 at 3:43 AM, Samuele Salis via RT <
> > > > met_help at ucar.edu <javascript:;>>
> > > > > > wrote:
> > > > > >
> > > > > > >
> > > > > > > Wed Feb 03 03:43:55 2016: Request 74966 was acted upon.
> > > > > > > Transaction: Ticket created by samuelesalis at gmail.com
> <javascript:;>
> > > > > > > Queue: met_help
> > > > > > > Subject: Troubles with point_stat
> > > > > > > Owner: Nobody
> > > > > > > Requestors: samuelesalis at gmail.com <javascript:;>
> > > > > > > Status: new
> > > > > > > Ticket <URL:
> > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966
> > > > > >
> > > > > > >
> > > > > > >
> > > > > > > Dear MET users,
> > > > > > >
> > > > > > > I have a problem with, I think, the point_stat tool.
> > > > > > >
> > > > > > > Here a brief history of what I did, just to check if I
made
> some
> > > > former
> > > > > > > mistake.
> > > > > > >
> > > > > > > I made a 24 h simulation with ARW initialised with GFS
and
> > starting
> > > > > > > 2016020100.
> > > > > > > I setted the io_form_history=5 to get a grib1 wrfout
file
> > (called
> > > > > > > wrfout.grb)
> > > > > > >
> > > > > > > I created an ascii file containing the accumulated
> precipitation
> > in
> > > > > three
> > > > > > > hours for
> > > > > > > the same dates filling the eleven fields described in
table 4.4
> > of
> > > > > > > Met_user_guide_v5.1
> > > > > > > just for two stations.
> > > > > > > Here is the file (ascii_obs.txt)
> > > > > > > NA AA000A001 20160201_030000 9.458 39.961 813 061 3 813
NA 2.6
> > > > > > > NA AA000A001 20160201_060000 9.458 39.961 813 061 3 813
NA 2.6
> > > > > > > NA AA000A001 20160201_090000 9.458 39.961 813 061 3 813
NA 2.9
> > > > > > > NA AA000A001 20160201_120000 9.458 39.961 813 061 3 813
NA 3.6
> > > > > > > NA AA000A001 20160201_150000 9.458 39.961 813 061 3 813
NA 7.1
> > > > > > > NA AA000A001 20160201_180000 9.458 39.961 813 061 3 813
NA 12.0
> > > > > > > NA AA000A001 20160201_210000 9.458 39.961 813 061 3 813
NA 12.5
> > > > > > > NA AA000A001 20160202_000000 9.458 39.961 813 061 3 813
NA 14.2
> > > > > > > NA AA000A002 20160201_030000 8.832 40.941 6 061 3 6
NA 4.2
> > > > > > > NA AA000A002 20160201_060000 8.832 40.941 6 061 3 6
NA 5.6
> > > > > > > NA AA000A002 20160201_090000 8.832 40.941 6 061 3 6
NA 12.9
> > > > > > > NA AA000A002 20160201_120000 8.832 40.941 6 061 3 6
NA 13.1
> > > > > > > NA AA000A002 20160201_150000 8.832 40.941 6 061 3 6
NA 13.1
> > > > > > > NA AA000A002 20160201_180000 8.832 40.941 6 061 3 6
NA 14.0
> > > > > > > NA AA000A002 20160201_210000 8.832 40.941 6 061 3 6
NA 15.5
> > > > > > > NA AA000A002 20160202_000000 8.832 40.941 6 061 3 6
NA 16.1
> > > > > > >
> > > > > > > Afterward I created an nc file using the tool
> > > > > > > ./ascii2nc ascii_obs.txt ascii_obs.nc -v 2
> > > > > > >
> > > > > > > The output of "ncdump -v obs_arr ascii_obs.nc" looks
fine.
> > > > > > >
> > > > > > > I modified PointStatConfig_default like this
> > > > > > >
> > > > > > > fcst = {
> > > > > > > sid_exc = [];
> > > > > > >
> > > > > > > field = [
> > > > > > > {
> > > > > > > name = "APCP";
> > > > > > > level = [ "A03" ];
> > > > > > > cat_thresh = [ >1.0 ];
> > > > > > > file_type = [ 1 ] ;
> > > > > > > init_time = [ "20160201_000000" ];
> > > > > > > valid_time = [ "20160202_000000" ];
> > > > > > > }
> > > > > > > ];
> > > > > > >
> > > > > > > }
> > > > > > > obs = fcst;
> > > > > > >
> > > > > > > and using "./point_stat wrfout.grb ascii_obs.nc
> > > > > > PointStatConfig_modified
> > > > > > > -outdir ./ -v 2"
> > > > > > > the program just give me the headers of the statistic
files and
> > the
> > > > > > > following message error
> > > > > > >
> > > > > > > DEBUG 1: User Config File: PointStatConfig_modified
> > > > > > > GSL_RNG_TYPE=mt19937
> > > > > > > GSL_RNG_SEED=1087746390
> > > > > > > DEBUG 1: Forecast File: wrfout_gfs.grib1
> > > > > > > DEBUG 1: Observation File: prova_ascii_grib1.nc
> > > > > > > DEBUG 2:
> > > > > > > DEBUG 2:
> > > > > > >
> > > > > > >
> > > > > >
> > > > >
> > > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > > > > > DEBUG 2:
> > > > > > > DEBUG 2: Reading data for APCP/A03.
> > > > > > > WARNING:
> > > > > > > WARNING: process_fcst_climo_files() -> no fields
matching
> > APCP/A03
> > > > > found
> > > > > > in
> > > > > > > file: wrfout_gfs.grib1
> > > > > > > WARNING:
> > > > > > > .....
> > > > > > >
> > > > > > >
> > > > > > > I also tried to modify the input file like this
> > > > > > >
> > > > > > > fcst = {
> > > > > > > sid_exc = [];
> > > > > > >
> > > > > > > field = [
> > > > > > > {
> > > > > > > name = "APCP_03";
> > > > > > > level = [ "(*,*)" ];
> > > > > > > cat_thresh = [ >1.0 ];
> > > > > > > file_type = [ 1 ] ;
> > > > > > > init_time = [ "20160201_000000" ];
> > > > > > > valid_time = [ "20160202_000000" ];
> > > > > > > }
> > > > > > > ];
> > > > > > >
> > > > > > > }
> > > > > > >
> > > > > > > obs = {
> > > > > > > message_type = [ "ADPSFC" ];
> > > > > > > sid_exc = [];
> > > > > > >
> > > > > > > field = [
> > > > > > > {
> > > > > > > name = "APCP";
> > > > > > > level = [ "A03" ];
> > > > > > > cat_thresh = [ >1.0 ];
> > > > > > > file_type = [ 1 ] ;
> > > > > > > init_time = [ "20160201_000000" ];
> > > > > > > valid_time = [ "20160202_000000" ];
> > > > > > > }
> > > > > > > ];
> > > > > > >
> > > > > > > }
> > > > > > >
> > > > > > > but in this case the message error is even worse and I
get
> > > > > > >
> > > > > > > DEBUG 1: User Config File: PointStatConfig_modified
> > > > > > > ERROR :
> > > > > > > ERROR : VarInfoGrib::set_dict() - unrecognized GRIB1
field
> > > > > abbreviation
> > > > > > > 'APCP_03' for table version 2
> > > > > > > ERROR :
> > > > > > >
> > > > > > > without any statistic output file.
> > > > > > >
> > > > > > >
> > > > > > > Any advise is really welcome.
> > > > > > >
> > > > > > > Thanks,
> > > > > > > Samuele
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: Troubles with point_stat
From: Samuele Salis
Time: Mon Feb 08 13:10:36 2016
Yes, I think the problem is substantially solved. I just have to
prepare
gribs with only one forecast time. This week I am busy with other
tasks in
my job, but I am confident I will have not problems.
Thanks,
Samuele
Il 08/feb/2016 17:40 "John Halley Gotway via RT" <met_help at ucar.edu>
ha
scritto:
> According to our records, your request has been resolved. If you
have any
> further questions or concerns, please respond to this message.
>
------------------------------------------------
Subject: Troubles with point_stat
From: John Halley Gotway
Time: Mon Feb 08 13:37:51 2016
Samuele,
Are you saying that you already have GRIB files containing multiple
output
times?
Or are you saying that you just need to automate the running of UPP
for
each output time?
MET actually can be configured to handle GRIB files which contain
multiple
output times, but the configuration file is a bit messier. Or I could
suggest some wgrib commands you could use to split multiple times from
a
GRIB file into separate files.
Thanks,
John
On Mon, Feb 8, 2016 at 1:10 PM, Samuele Salis via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966 >
>
> Yes, I think the problem is substantially solved. I just have to
prepare
> gribs with only one forecast time. This week I am busy with other
tasks in
> my job, but I am confident I will have not problems.
> Thanks,
> Samuele
> Il 08/feb/2016 17:40 "John Halley Gotway via RT" <met_help at ucar.edu>
ha
> scritto:
>
> > According to our records, your request has been resolved. If you
have any
> > further questions or concerns, please respond to this message.
> >
>
>
------------------------------------------------
Subject: Troubles with point_stat
From: Samuele Salis
Time: Tue Feb 09 02:37:53 2016
Ok John,
at the moment I am running WRF with the option to write directly a
grib1
output with multiple output times.
According to your advice I separated single output times using
wgrib to avoid the difficulties in editing of MET's configuration
time.
I already used some wgrib commands like this one to get (e.g.) ACPC in
the
first hour
wgrib -s wrfout.grib1 | grep ACPCP |grep "0-1hr acc"| wgrib -i -grib
wrfout.grib1 -o wrfout.grib1_ACPCP_01
Now I have to prepare a new file with my observations, create the nc
file
with ascii2nc,
prepare a new PointStatConfig file, run again point_stat and finally
check
my output.
I am still not fond in statistical analysis applied to meteorological
forecasts, so at the
moment I am just interested on the feasibility of the statistical
analysis
using MET and the
data I can access. If I will be sure about this first step, I will
start
studying the theory
to decide which statistical methods are more suitable for the
verification
I want to do
for my meteorological simulations.
Thanks,
Samuele
On 8 February 2016 at 21:37, John Halley Gotway via RT
<met_help at ucar.edu>
wrote:
> Samuele,
>
> Are you saying that you already have GRIB files containing multiple
output
> times?
>
> Or are you saying that you just need to automate the running of UPP
for
> each output time?
>
> MET actually can be configured to handle GRIB files which contain
multiple
> output times, but the configuration file is a bit messier. Or I
could
> suggest some wgrib commands you could use to split multiple times
from a
> GRIB file into separate files.
>
> Thanks,
> John
>
> On Mon, Feb 8, 2016 at 1:10 PM, Samuele Salis via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966 >
> >
> > Yes, I think the problem is substantially solved. I just have to
prepare
> > gribs with only one forecast time. This week I am busy with other
tasks
> in
> > my job, but I am confident I will have not problems.
> > Thanks,
> > Samuele
> > Il 08/feb/2016 17:40 "John Halley Gotway via RT"
<met_help at ucar.edu> ha
> > scritto:
> >
> > > According to our records, your request has been resolved. If you
have
> any
> > > further questions or concerns, please respond to this message.
> > >
> >
> >
>
>
------------------------------------------------
Subject: Troubles with point_stat
From: Samuele Salis
Time: Tue Feb 09 06:09:52 2016
Ok,
I processed the new data creating the new observation file until I
processed it
through point_stat.
The output looks at a first glance good (at the moment I just checked
the
contingency table in file .ctc ).
Later on I will control other outputs, but I think I finally overcame
this
first
step using MET.
Thanks again,
Samuele
On 9 February 2016 at 10:37, Samuele Salis <samuelesalis at gmail.com>
wrote:
> Ok John,
>
> at the moment I am running WRF with the option to write directly a
grib1
> output with multiple output times.
>
> According to your advice I separated single output times using
> wgrib to avoid the difficulties in editing of MET's configuration
time.
>
> I already used some wgrib commands like this one to get (e.g.) ACPC
in the
> first hour
> wgrib -s wrfout.grib1 | grep ACPCP |grep "0-1hr acc"| wgrib -i -grib
> wrfout.grib1 -o wrfout.grib1_ACPCP_01
>
> Now I have to prepare a new file with my observations, create the nc
file
> with ascii2nc,
> prepare a new PointStatConfig file, run again point_stat and finally
check
> my output.
>
> I am still not fond in statistical analysis applied to
meteorological
> forecasts, so at the
> moment I am just interested on the feasibility of the statistical
analysis
> using MET and the
> data I can access. If I will be sure about this first step, I will
start
> studying the theory
> to decide which statistical methods are more suitable for the
verification
> I want to do
> for my meteorological simulations.
>
> Thanks,
> Samuele
>
>
>
> On 8 February 2016 at 21:37, John Halley Gotway via RT
<met_help at ucar.edu>
> wrote:
>
>> Samuele,
>>
>> Are you saying that you already have GRIB files containing multiple
output
>> times?
>>
>> Or are you saying that you just need to automate the running of UPP
for
>> each output time?
>>
>> MET actually can be configured to handle GRIB files which contain
multiple
>> output times, but the configuration file is a bit messier. Or I
could
>> suggest some wgrib commands you could use to split multiple times
from a
>> GRIB file into separate files.
>>
>> Thanks,
>> John
>>
>> On Mon, Feb 8, 2016 at 1:10 PM, Samuele Salis via RT
<met_help at ucar.edu>
>> wrote:
>>
>> >
>> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74966 >
>> >
>> > Yes, I think the problem is substantially solved. I just have to
prepare
>> > gribs with only one forecast time. This week I am busy with other
tasks
>> in
>> > my job, but I am confident I will have not problems.
>> > Thanks,
>> > Samuele
>> > Il 08/feb/2016 17:40 "John Halley Gotway via RT"
<met_help at ucar.edu> ha
>> > scritto:
>> >
>> > > According to our records, your request has been resolved. If
you have
>> any
>> > > further questions or concerns, please respond to this message.
>> > >
>> >
>> >
>>
>>
>
------------------------------------------------
More information about the Met_help
mailing list