[Met_help] [rt.rap.ucar.edu #84002] History for point_stat can deal with less 46 stations?
John Halley Gotway via RT
met_help at ucar.edu
Tue Jul 9 12:03:54 MDT 2019
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Dear Sir
When I used point_stat to verify WRF forcast 6hr cumulative precipitation
in a small domain(20kmX30km),which has 46 obs stations, point_stat find 0
pairs.
WRF spatial distance is 3km,while the distance of obs stations is less 2km,
Does obs stations distance less WRF spatial distance can lead to the
following trouble?
Searching 46 observations from 46 PrepBufr messages.
----------------------------------------------------------------------------
----
Processing APCP_06(*,*) versus APCP_06/A6, for observation type ADPSFC, over
region FULL, for interpolation method DW_MEAN(9), using 0 pairs.
----------------------------------------------------------------------------
----
Output file:
/data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800/point_stat_APCP_6H
R_AWS_300000L_20170809_060000V.stat
Output file:
/data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800/point_stat_APCP_6H
R_AWS_300000L_20170809_060000V_fho.txt
Output file:
/data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800/point_stat_APCP_6H
R_AWS_300000L_20170809_060000V_ctc.txt
Output file:
/data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800/point_stat_APCP_6H
R_AWS_300000L_20170809_060000V_cts.txt
Output file:
/data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800/point_stat_APCP_6H
R_AWS_300000L_20170809_060000V_mpr.txt
Forecast File:
/data6t/MET-results/fcst/wrfrain/RMAPS/wrfrain_d02_f2017080800_v2017080912_a
cc6
Climatology File: none
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: point_stat can deal with less 46 stations?
From: John Halley Gotway
Time: Fri Feb 09 09:22:36 2018
Hello Hgao,
I see that you're getting 0 matched pairs from Point-Stat. There are
many
possible reasons why Point-Stat found 0 matched pairs. Please rerun
Point-Stat but turn up the verbosity level to 3 using the "-v 3"
command
line option. That will dump out a list of counts for reasons *why*
the
observations were not used. And that should point you in the right
direction.
Thanks,
John Halley Gotway
On Fri, Feb 9, 2018 at 12:09 AM, hgao via RT <met_help at ucar.edu>
wrote:
>
> Fri Feb 09 00:09:08 2018: Request 84002 was acted upon.
> Transaction: Ticket created by hgao at ium.cn
> Queue: met_help
> Subject: point_stat can deal with less 46 stations?
> Owner: Nobody
> Requestors: hgao at ium.cn
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84002 >
>
>
> Dear Sir
>
> When I used point_stat to verify WRF forcast 6hr cumulative
> precipitation
> in a small domain(20kmX30km),which has 46 obs stations, point_stat
find 0
> pairs.
> WRF spatial distance is 3km,while the distance of obs stations is
less 2km,
> Does obs stations distance less WRF spatial distance can lead to the
> following trouble?
>
>
>
> Searching 46 observations from 46 PrepBufr messages.
>
> ------------------------------------------------------------
> ----------------
> ----
>
> Processing APCP_06(*,*) versus APCP_06/A6, for observation type
ADPSFC,
> over
> region FULL, for interpolation method DW_MEAN(9), using 0 pairs.
>
> ------------------------------------------------------------
> ----------------
> ----
>
> Output file:
> /data6t/MET-results/scores/point_stat/RMAPS-ST/
> 2017080800/point_stat_APCP_6H
> R_AWS_300000L_20170809_060000V.stat
> Output file:
> /data6t/MET-results/scores/point_stat/RMAPS-ST/
> 2017080800/point_stat_APCP_6H
> R_AWS_300000L_20170809_060000V_fho.txt
> Output file:
> /data6t/MET-results/scores/point_stat/RMAPS-ST/
> 2017080800/point_stat_APCP_6H
> R_AWS_300000L_20170809_060000V_ctc.txt
> Output file:
> /data6t/MET-results/scores/point_stat/RMAPS-ST/
> 2017080800/point_stat_APCP_6H
> R_AWS_300000L_20170809_060000V_cts.txt
> Output file:
> /data6t/MET-results/scores/point_stat/RMAPS-ST/
> 2017080800/point_stat_APCP_6H
> R_AWS_300000L_20170809_060000V_mpr.txt
> Forecast File:
> /data6t/MET-results/fcst/wrfrain/RMAPS/wrfrain_d02_
> f2017080800_v2017080912_a
> cc6
> Climatology File: none
>
>
------------------------------------------------
Subject: point_stat can deal with less 46 stations?
From: John Halley Gotway
Time: Fri Feb 09 09:24:04 2018
I also see that you submitted another ticket showing the configuration
file
you used for MET version 3.0.1. There have been many, many updates
since
that version. The current version is 6.1.
Before doing anything else, I'd recommend upgrading to the current
version
of the software.
Thanks,
John
On Fri, Feb 9, 2018 at 9:22 AM, John Halley Gotway <johnhg at ucar.edu>
wrote:
> Hello Hgao,
>
> I see that you're getting 0 matched pairs from Point-Stat. There
are many
> possible reasons why Point-Stat found 0 matched pairs. Please rerun
> Point-Stat but turn up the verbosity level to 3 using the "-v 3"
command
> line option. That will dump out a list of counts for reasons *why*
the
> observations were not used. And that should point you in the right
> direction.
>
> Thanks,
> John Halley Gotway
>
> On Fri, Feb 9, 2018 at 12:09 AM, hgao via RT <met_help at ucar.edu>
wrote:
>
>>
>> Fri Feb 09 00:09:08 2018: Request 84002 was acted upon.
>> Transaction: Ticket created by hgao at ium.cn
>> Queue: met_help
>> Subject: point_stat can deal with less 46 stations?
>> Owner: Nobody
>> Requestors: hgao at ium.cn
>> Status: new
>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84002 >
>>
>>
>> Dear Sir
>>
>> When I used point_stat to verify WRF forcast 6hr cumulative
>> precipitation
>> in a small domain(20kmX30km),which has 46 obs stations, point_stat
find 0
>> pairs.
>> WRF spatial distance is 3km,while the distance of obs stations is
less
>> 2km,
>> Does obs stations distance less WRF spatial distance can lead to
the
>> following trouble?
>>
>>
>>
>> Searching 46 observations from 46 PrepBufr messages.
>>
>> ------------------------------------------------------------
>> ----------------
>> ----
>>
>> Processing APCP_06(*,*) versus APCP_06/A6, for observation type
ADPSFC,
>> over
>> region FULL, for interpolation method DW_MEAN(9), using 0 pairs.
>>
>> ------------------------------------------------------------
>> ----------------
>> ----
>>
>> Output file:
>> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
>> <(201)%20708-0800>/point_stat_APCP_6H
>> R_AWS_300000L_20170809_060000V.stat
>> Output file:
>> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
>> <(201)%20708-0800>/point_stat_APCP_6H
>> R_AWS_300000L_20170809_060000V_fho.txt
>> Output file:
>> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
>> <(201)%20708-0800>/point_stat_APCP_6H
>> R_AWS_300000L_20170809_060000V_ctc.txt
>> Output file:
>> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
>> <(201)%20708-0800>/point_stat_APCP_6H
>> R_AWS_300000L_20170809_060000V_cts.txt
>> Output file:
>> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
>> <(201)%20708-0800>/point_stat_APCP_6H
>> R_AWS_300000L_20170809_060000V_mpr.txt
>> Forecast File:
>> /data6t/MET-results/fcst/wrfrain/RMAPS/wrfrain_d02_f20170808
>> 00_v2017080912_a
>> cc6
>> Climatology File: none
>>
>>
>
------------------------------------------------
Subject: point_stat can deal with less 46 stations?
From: hgao
Time: Sat Feb 10 05:02:12 2018
Dear Dr. John,
Thank you very much for your rapid advice , I tried with -
v 3,got the following log :
Reading records for APCP_06(*,*).
For APCP_06(*,*) found 1 forecast levels and 0 climatology
levels.--------------------------------------------------------------------------------
Searching 46 observations from 46 PrepBufr
messages.--------------------------------------------------------------------------------
Processing APCP_06(*,*) versus APCP_06/A6, for observation type
ADPSFC, over region FULL, for interpolation method DW_MEAN(9), using 0
pairs.
Number of matched pairs = 0
Observations processed = 46
Rejected: GRIB code = 0
Rejected: valid time = 46
Rejected: bad obs value = 0
Rejected: off the grid = 0
Rejected: level mismatch = 0
Rejected: message type = 0
Rejected: masking region = 0
Rejected: bad fcst value =
0--------------------------------------------------------------------------------
On the other hand, I 've tried point_stat of version 6.0,but
point_stat can't accept Netcdf file of WRF forcast . I 'll try v6.1
next .
------------------------------------------------------------------
发件人：John Halley Gotway via RT <met_help at ucar.edu>发送时间：2018年2月10日(星期六)
00:24收件人：高华 <hgao at ium.cn>主 题：Re: [rt.rap.ucar.edu #84002] point_stat
can deal with less 46
stations? I also see that you submitted another ticket showing the configuration file
you used for MET version 3.0.1. There have been many, many updates since
that version. The current version is 6.1.
Before doing anything else, I'd recommend upgrading to the current version
of the software.
Thanks,
John
On Fri, Feb 9, 2018 at 9:22 AM, John Halley Gotway <johnhg at ucar.edu> wrote:
> Hello Hgao,
>
> I see that you're getting 0 matched pairs from Point-
Stat. There are many
> possible reasons why Point-Stat found 0 matched pairs. Please rerun
> Point-Stat but turn up the verbosity level to 3 using the "-
v 3" command
> line option. That will dump out a list of counts for reasons *why* the
> observations were not used. And that should point you in the right
> direction.
>
> Thanks,
> John Halley Gotway
>
> On Fri, Feb 9, 2018 at 12:09 AM, hgao via RT <met_help at ucar.edu> wrote:
>
>>
>> Fri Feb 09 00:09:08 2018: Request 84002 was acted upon.
>> Transaction: Ticket created by hgao at ium.cn
>> Queue: met_help
>> Subject: point_stat can deal with less 46 stations?
>> Owner: Nobody
>> Requestors: hgao at ium.cn
>> Status: new
>> Ticket <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84002 >
>>
>>
>> Dear Sir
>>
>> When I used point_stat to verify WRF forcast 6hr cumulative
>> precipitation
>> in a small domain(20kmX30km),which has 46 obs stations, point_stat find 0
>> pairs.
>> WRF spatial distance is 3km,while the distance of obs stations is less
>> 2km,
>> Does obs stations distance less WRF spatial distance can lead to the
>> following trouble?
>>
>>
>>
>> Searching 46 observations from 46 PrepBufr messages.
>>
>> ------------------------------------------------------------
>> ----------------
>> ----
>>
>> Processing APCP_06(*,*) versus APCP_06/A6, for observation type ADPSFC,
>> over
>> region FULL, for interpolation method DW_MEAN(9), using 0 pairs.
>>
>> ------------------------------------------------------------
>> ----------------
>> ----
>>
>> Output file:
>> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
>> <(201)%20708-0800>/point_stat_APCP_6H
>> R_AWS_300000L_20170809_060000V.stat
>> Output file:
>> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
>> <(201)%20708-0800>/point_stat_APCP_6H
>> R_AWS_300000L_20170809_060000V_fho.txt
>> Output file:
>> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
>> <(201)%20708-0800>/point_stat_APCP_6H
>> R_AWS_300000L_20170809_060000V_ctc.txt
>> Output file:
>> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
>> <(201)%20708-0800>/point_stat_APCP_6H
>> R_AWS_300000L_20170809_060000V_cts.txt
>> Output file:
>> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
>> <(201)%20708-0800>/point_stat_APCP_6H
>> R_AWS_300000L_20170809_060000V_mpr.txt
>> Forecast File:
>> /data6t/MET-results/fcst/wrfrain/RMAPS/wrfrain_d02_f20170808
>> 00_v2017080912_a
>> cc6
>> Climatology File: none
>>
>>
>
------------------------------------------------
Subject: point_stat can deal with less 46 stations?
From: hgao
Time: Sat Feb 10 20:33:25 2018
Dr. John
I run point_stat of v6.1_bugfix,it can't read Netcdf WRF file
by UPP:[iumnwp at Ec226 bin]$ ./point_stat
./wrfrain_d02_f2017080800_v2017080812_acc6.nc
./20170808120000.06_new.nc ../scripts/config/PointStatConfig -outdir
./ -v 2
DEBUG 1: Default Config File: /home/iumnwp/met-
6.1_bugfix/share/met/config/PointStatConfig_default
DEBUG 1: User Config File: ../scripts/config/PointStatConfig
ERROR :
ERROR : VarInfoGrib::add_grib_code() -> unrecognized GRIB1 field
abbreviation 'APCP_06' for table version 2
ERROR :
------------------------------------------------------------------
发件人：John Halley Gotway via RT <met_help at ucar.edu>发送时间：2018年2月10日(星期六)
00:24收件人：高华 <hgao at ium.cn>主 题：Re: [rt.rap.ucar.edu #84002] point_stat
can deal with less 46
stations? I also see that you submitted another ticket showing the configuration file
you used for MET version 3.0.1. There have been many, many updates since
that version. The current version is 6.1.
Before doing anything else, I'd recommend upgrading to the current version
of the software.
Thanks,
John
On Fri, Feb 9, 2018 at 9:22 AM, John Halley Gotway <johnhg at ucar.edu> wrote:
> Hello Hgao,
>
> I see that you're getting 0 matched pairs from Point-
Stat. There are many
> possible reasons why Point-Stat found 0 matched pairs. Please rerun
> Point-Stat but turn up the verbosity level to 3 using the "-
v 3" command
> line option. That will dump out a list of counts for reasons *why* the
> observations were not used. And that should point you in the right
> direction.
>
> Thanks,
> John Halley Gotway
>
> On Fri, Feb 9, 2018 at 12:09 AM, hgao via RT <met_help at ucar.edu> wrote:
>
>>
>> Fri Feb 09 00:09:08 2018: Request 84002 was acted upon.
>> Transaction: Ticket created by hgao at ium.cn
>> Queue: met_help
>> Subject: point_stat can deal with less 46 stations?
>> Owner: Nobody
>> Requestors: hgao at ium.cn
>> Status: new
>> Ticket <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84002 >
>>
>>
>> Dear Sir
>>
>> When I used point_stat to verify WRF forcast 6hr cumulative
>> precipitation
>> in a small domain(20kmX30km),which has 46 obs stations, point_stat find 0
>> pairs.
>> WRF spatial distance is 3km,while the distance of obs stations is less
>> 2km,
>> Does obs stations distance less WRF spatial distance can lead to the
>> following trouble?
>>
>>
>>
>> Searching 46 observations from 46 PrepBufr messages.
>>
>> ------------------------------------------------------------
>> ----------------
>> ----
>>
>> Processing APCP_06(*,*) versus APCP_06/A6, for observation type ADPSFC,
>> over
>> region FULL, for interpolation method DW_MEAN(9), using 0 pairs.
>>
>> ------------------------------------------------------------
>> ----------------
>> ----
>>
>> Output file:
>> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
>> <(201)%20708-0800>/point_stat_APCP_6H
>> R_AWS_300000L_20170809_060000V.stat
>> Output file:
>> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
>> <(201)%20708-0800>/point_stat_APCP_6H
>> R_AWS_300000L_20170809_060000V_fho.txt
>> Output file:
>> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
>> <(201)%20708-0800>/point_stat_APCP_6H
>> R_AWS_300000L_20170809_060000V_ctc.txt
>> Output file:
>> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
>> <(201)%20708-0800>/point_stat_APCP_6H
>> R_AWS_300000L_20170809_060000V_cts.txt
>> Output file:
>> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
>> <(201)%20708-0800>/point_stat_APCP_6H
>> R_AWS_300000L_20170809_060000V_mpr.txt
>> Forecast File:
>> /data6t/MET-results/fcst/wrfrain/RMAPS/wrfrain_d02_f20170808
>> 00_v2017080912_a
>> cc6
>> Climatology File: none
>>
>>
>
------------------------------------------------
Subject: point_stat can deal with less 46 stations?
From: John Halley Gotway
Time: Wed Feb 14 09:33:43 2018
OK, I want to circle back to your email from 4 days ago. You sent me
this
output from Point-Stat:
Number of matched pairs = 0
Observations processed = 46
Rejected: GRIB code = 0
Rejected: valid time = 46
Rejected: bad obs value = 0
Rejected: off the grid = 0
Rejected: level mismatch = 0
Rejected: message type = 0
Rejected: masking region = 0
Rejected: bad fcst value = 0
This tells you the reason why Point-Stat is not finding any matched
pairs.
All 46 observations were skipped due to the valid time. None of them
were
close enough in time to the valid time of the forecast for them to be
used.
This matching time window is defined in the Point-Stat configuration
file
using the "obs_window" setting:
obs_window = {
beg = -5400;
end = 5400;
}
The "beg" and "end" settings specify the boundaries of the matching
time
window as the number of seconds relative to the forecast valid time.
The +/- 5400 seconds listed above means that the observations have to
be
within 1.5 hours of the forecast valid time to be used. You can make
that
matching time window very, very large and then see if any matches show
up.
For example, the following would match any observations within 1 year
of
the forecast!
obs_window = {
beg = -60*60*24*365;
end = 60*60*24*365;
}
Thanks,
John
On Sat, Feb 10, 2018 at 8:33 PM, hgao via RT <met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84002 >
>
> Dr. John
> I run point_stat of v6.1_bugfix,it can't read Netcdf WRF file
> by UPP:[iumnwp at Ec226 bin]$ ./point_stat ./wrfrain_d02_f2017080800_
> v2017080812_acc6.nc ./20170808120000.06_new.nc
../scripts/config/PointStatConfig
> -outdir ./ -v 2
> DEBUG 1: Default Config File: /home/iumnwp/met-6.1_bugfix/
> share/met/config/PointStatConfig_default
> DEBUG 1: User Config File: ../scripts/config/PointStatConfig
> ERROR :
> ERROR : VarInfoGrib::add_grib_code() -> unrecognized GRIB1 field
> abbreviation 'APCP_06' for table version 2
> ERROR :
------------------------------------------------------------------发件人：John
> Halley Gotway via RT <met_help at ucar.edu>发送时间：2018年2月10日(星期六)
00:24收件人：高华 <
> hgao at ium.cn>主 题：Re: [rt.rap.ucar.edu #84002] point_stat can deal
with
> less 46 stations? I also see that you submitted another ticket
> showing the configuration file
> you used for MET version 3.0.1. There have been many, many updates
since
> that version. The current version is 6.1.
>
> Before doing anything else, I'd recommend upgrading to the current
version
> of the software.
>
> Thanks,
> John
>
> On Fri, Feb 9, 2018 at 9:22 AM, John Halley Gotway <johnhg at ucar.edu
> > wrote:
>
> > Hello Hgao,
> >
> > I see that you're getting 0 matched pairs from Point-Stat.
> There are many
> > possible reasons why Point-Stat found 0 matched pairs. Please
rerun
> > Point-Stat but turn up the verbosity level to 3 using the "-v 3"
command
> > line option. That will dump out a list of counts for reasons
*why* the
> > observations were not used. And that should point you in the
right
> > direction.
> >
> > Thanks,
> > John Halley Gotway
> >
> > On Fri, Feb 9, 2018 at 12:09 AM, hgao via RT <met_help at ucar.edu>
wrote:
> >
> >>
> >> Fri Feb 09 00:09:08 2018: Request 84002 was acted upon.
> >> Transaction: Ticket created by hgao at ium.cn
> >> Queue: met_help
> >> Subject: point_stat can deal with less 46 stations?
> >> Owner: Nobody
> >> Requestors: hgao at ium.cn
> >> Status: new
> >> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84002 >
> >>
> >>
> >> Dear Sir
> >>
> >> When I used point_stat to verify WRF forcast 6hr cumulative
> >> precipitation
> >> in a small domain(20kmX30km),which has 46 obs
> stations, point_stat find 0
> >> pairs.
> >> WRF spatial distance is 3km,while the distance of obs stations is
less
> >> 2km,
> >> Does obs stations distance less WRF spatial distance can lead to
the
> >> following trouble?
> >>
> >>
> >>
> >> Searching 46 observations from 46 PrepBufr messages.
> >>
> >> ------------------------------------------------------------
> >> ----------------
> >> ----
> >>
> >> Processing APCP_06(*,*) versus APCP_06/A6, for observation type
ADPSFC,
> >> over
> >> region FULL, for interpolation method DW_MEAN(9), using 0 pairs.
> >>
> >> ------------------------------------------------------------
> >> ----------------
> >> ----
> >>
> >> Output file:
> >> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
> >> <(201)%20708-0800>/point_stat_APCP_6H
> >> R_AWS_300000L_20170809_060000V.stat
> >> Output file:
> >> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
> >> <(201)%20708-0800>/point_stat_APCP_6H
> >> R_AWS_300000L_20170809_060000V_fho.txt
> >> Output file:
> >> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
> >> <(201)%20708-0800>/point_stat_APCP_6H
> >> R_AWS_300000L_20170809_060000V_ctc.txt
> >> Output file:
> >> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
> >> <(201)%20708-0800>/point_stat_APCP_6H
> >> R_AWS_300000L_20170809_060000V_cts.txt
> >> Output file:
> >> /data6t/MET-results/scores/point_stat/RMAPS-ST/2017080800
> >> <(201)%20708-0800>/point_stat_APCP_6H
> >> R_AWS_300000L_20170809_060000V_mpr.txt
> >> Forecast File:
> >> /data6t/MET-results/fcst/wrfrain/RMAPS/wrfrain_d02_f20170808
> >> 00_v2017080912_a
> >> cc6
> >> Climatology File: none
> >>
> >>
> >
>
>
------------------------------------------------
More information about the Met_help
mailing list