[Met_help] [rt.rap.ucar.edu #84858] History for How to manage APCP with different accumulation period in MET
John Halley Gotway via RT
met_help at ucar.edu
Tue Jul 9 12:03:46 MDT 2019
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
I am currently using UPP with hourly output of WRF (24 frames in a file) in
order to be able to validate my WRF simulation with MET tool. So I get
hourly grib1 WRFPRS files, each one containig a field with APCP at the
level A01.
I would like to validate this field with SYNOP data. In these data, I have
only APCP data at the level A03 (cumul over 3hrs).
I use point_stat and try to test the hourly WRF field:
{
name = "APCP";
level = [ "A01","A03","A06" ] ;
},
I use the same field definition for fcst and obs.
In the log, I see that some forecast data are found for APCP/A01 or
APCP/A03.
But I do not get data for APCP/A01 for every hour. It seems strange because
I work with hourly WRF output. Is this normal ? How are APCP with no
precipitation (0 mm) processed in point_stat ?
And how do I get forecast data for APCP/A03 while I'm running point_stat
for each time file separately ???
Could you explain me how point_stat is processing APCP ? There is very
little information in the doc.
Thanks a lot for your help
PS: I send you in copy the pointstat_config file and the log after
point_stat
*Blandine L'Hévéder*
Météo - Océan -Climat
blandine.lheveder at gmail.com
tel. 06 62 53 32 37
https://fr.linkedin.com/in/blandine-l-heveder-493171a9
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: How to manage APCP with different accumulation period in MET
From: John Halley Gotway
Time: Mon Apr 23 14:37:13 2018
Hello Blandine,
I see that you have questions about using Point-Stat to compute
statistics
for accumulated precipitation. Thanks for sending your Point-Stat
config
file and log from Point-Stat. My first suggestion would be increasing
the
logging level when you run Point-Stat. Try passing in "-v 3" on the
command line to increase the logging level to 3. That'll include log
messages that look something like this:
DEBUG 3: Number of matched pairs = 132
DEBUG 3: Observations processed = 89893
DEBUG 3: Rejected: SID exclusion = 0
DEBUG 3: Rejected: obs type = 79788
DEBUG 3: Rejected: valid time = 0
DEBUG 3: Rejected: bad obs value = 0
DEBUG 3: Rejected: off the grid = 5
DEBUG 3: Rejected: level mismatch = 9219
DEBUG 3: Rejected: quality marker = 0
DEBUG 3: Rejected: message type = 344
DEBUG 3: Rejected: masking region = 405
DEBUG 3: Rejected: bad fcst value = 0
DEBUG 3: Rejected: duplicates = 0
For each verification task, you'll see these counts describing the
number
of matched pairs that were used to compute statistics... as well as
counts
for the different reasons why observations were *not* used for this
task.
Since precipitation is accumulated in time, it is more complicated to
verify. Ultimately, you're limited to verifying only at the
accumulation
intervals for which you have observations. So the real question is,
what
point observation datasets are you using and accumulation intervals
are
present in them?
If you'd like, you could post a sample observation file to our
anonymous
ftp site along with the PB2NC (presumably) config file (or perhaps
they're
ascii?) you're using to process it:
https://dtcenter.org/met/users/support/met_help.php#ftp
Thanks,
John
On Mon, Apr 23, 2018 at 4:53 AM, Blandine Lheveder via RT
<met_help at ucar.edu
> wrote:
>
> Mon Apr 23 04:53:34 2018: Request 84858 was acted upon.
> Transaction: Ticket created by blandine.lheveder at gmail.com
> Queue: met_help
> Subject: How to manage APCP with different accumulation period
in MET
> Owner: Nobody
> Requestors: blandine.lheveder at gmail.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84858 >
>
>
> I am currently using UPP with hourly output of WRF (24 frames in a
file) in
> order to be able to validate my WRF simulation with MET tool. So I
get
> hourly grib1 WRFPRS files, each one containig a field with APCP at
the
> level A01.
>
> I would like to validate this field with SYNOP data. In these data,
I have
> only APCP data at the level A03 (cumul over 3hrs).
>
> I use point_stat and try to test the hourly WRF field:
> {
> name = "APCP";
> level = [ "A01","A03","A06" ] ;
> },
> I use the same field definition for fcst and obs.
>
> In the log, I see that some forecast data are found for APCP/A01 or
> APCP/A03.
> But I do not get data for APCP/A01 for every hour. It seems strange
because
> I work with hourly WRF output. Is this normal ? How are APCP with no
> precipitation (0 mm) processed in point_stat ?
> And how do I get forecast data for APCP/A03 while I'm running
point_stat
> for each time file separately ???
>
> Could you explain me how point_stat is processing APCP ? There is
very
> little information in the doc.
>
> Thanks a lot for your help
>
> PS: I send you in copy the pointstat_config file and the log after
> point_stat
>
> *Blandine L'Hévéder*
> Météo - Océan -Climat
>
> blandine.lheveder at gmail.com
> tel. 06 62 53 32 37
> https://fr.linkedin.com/in/blandine-l-heveder-493171a9
>
>
------------------------------------------------
Subject: How to manage APCP with different accumulation period in MET
From: Blandine Lheveder
Time: Tue Apr 24 04:36:10 2018
Dear John,
Thank you a lot for your help.
I have increased the login level to 3 and found the problem.
Point_stat could not read the valid_time in the forecast file (output
of
pcp_combine) because I had put a wrong argument " file_type =
NETCDF_NCCF;"
for the fcst.
Now it is OK :-)
I have another question:
I use SYNOP observations to validate with MET my WRF simulation. I
would
like to validate the 10m wind. In the SYNOP I have only the wind
direction
(0-360°, GRIB code 031) and wind speed (m/s, GRIB code 032). In the
wrf
output, after the UPP post-processing, I do not get these fields but
only
the UGRD and VGRD at Z10.
I have seen that in stat_analysis it is possible to compute statistics
on
wind direction, but only if pair of variables UGRD and VGRD are
processed
in the Point_stat step.
Do you have an idea how to proceed in this case ?
Thanks a lot
Best regards
Blandine
*Blandine L'Hévéder*
Météo - Océan -Climat
blandine.lheveder at gmail.com
tel. 06 62 53 32 37
https://fr.linkedin.com/in/blandine-l-heveder-493171a9
2018-04-23 22:37 GMT+02:00 John Halley Gotway via RT
<met_help at ucar.edu>:
> Hello Blandine,
>
> I see that you have questions about using Point-Stat to compute
statistics
> for accumulated precipitation. Thanks for sending your Point-Stat
config
> file and log from Point-Stat. My first suggestion would be
increasing the
> logging level when you run Point-Stat. Try passing in "-v 3" on the
> command line to increase the logging level to 3. That'll include
log
> messages that look something like this:
>
> DEBUG 3: Number of matched pairs = 132
> DEBUG 3: Observations processed = 89893
> DEBUG 3: Rejected: SID exclusion = 0
> DEBUG 3: Rejected: obs type = 79788
> DEBUG 3: Rejected: valid time = 0
> DEBUG 3: Rejected: bad obs value = 0
> DEBUG 3: Rejected: off the grid = 5
> DEBUG 3: Rejected: level mismatch = 9219
> DEBUG 3: Rejected: quality marker = 0
> DEBUG 3: Rejected: message type = 344
> DEBUG 3: Rejected: masking region = 405
> DEBUG 3: Rejected: bad fcst value = 0
> DEBUG 3: Rejected: duplicates = 0
>
> For each verification task, you'll see these counts describing the
number
> of matched pairs that were used to compute statistics... as well as
counts
> for the different reasons why observations were *not* used for this
task.
>
> Since precipitation is accumulated in time, it is more complicated
to
> verify. Ultimately, you're limited to verifying only at the
accumulation
> intervals for which you have observations. So the real question is,
what
> point observation datasets are you using and accumulation intervals
are
> present in them?
>
> If you'd like, you could post a sample observation file to our
anonymous
> ftp site along with the PB2NC (presumably) config file (or perhaps
they're
> ascii?) you're using to process it:
> https://dtcenter.org/met/users/support/met_help.php#ftp
>
> Thanks,
> John
>
> On Mon, Apr 23, 2018 at 4:53 AM, Blandine Lheveder via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > Mon Apr 23 04:53:34 2018: Request 84858 was acted upon.
> > Transaction: Ticket created by blandine.lheveder at gmail.com
> > Queue: met_help
> > Subject: How to manage APCP with different accumulation
period in
> MET
> > Owner: Nobody
> > Requestors: blandine.lheveder at gmail.com
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84858 >
> >
> >
> > I am currently using UPP with hourly output of WRF (24 frames in a
file)
> in
> > order to be able to validate my WRF simulation with MET tool. So I
get
> > hourly grib1 WRFPRS files, each one containig a field with APCP at
the
> > level A01.
> >
> > I would like to validate this field with SYNOP data. In these
data, I
> have
> > only APCP data at the level A03 (cumul over 3hrs).
> >
> > I use point_stat and try to test the hourly WRF field:
> > {
> > name = "APCP";
> > level = [ "A01","A03","A06" ] ;
> > },
> > I use the same field definition for fcst and obs.
> >
> > In the log, I see that some forecast data are found for APCP/A01
or
> > APCP/A03.
> > But I do not get data for APCP/A01 for every hour. It seems
strange
> because
> > I work with hourly WRF output. Is this normal ? How are APCP with
no
> > precipitation (0 mm) processed in point_stat ?
> > And how do I get forecast data for APCP/A03 while I'm running
point_stat
> > for each time file separately ???
> >
> > Could you explain me how point_stat is processing APCP ? There is
very
> > little information in the doc.
> >
> > Thanks a lot for your help
> >
> > PS: I send you in copy the pointstat_config file and the log after
> > point_stat
> >
> > *Blandine L'Hévéder*
> > Météo - Océan -Climat
> >
> > blandine.lheveder at gmail.com
> > tel. 06 62 53 32 37
> > https://fr.linkedin.com/in/blandine-l-heveder-493171a9
> >
> >
>
>
------------------------------------------------
Subject: How to manage APCP with different accumulation period in MET
From: John Halley Gotway
Time: Tue Apr 24 09:55:08 2018
Blandine,
OK, so your SYNOP observations include wind speed and direction rather
than
U and V. And you're GRIB output from UPP contain U and V. For wind
speed,
everything is fine. Just request that "WIND" be verified... that's
the
GRIB code abbreviation for wind speed. When you request WIND from a
GRIB
file, MET first searches for a matching record. If not found, it then
finds the corresponding U and V and automatically derives wind speed
for
you.
Unfortunately though, I do not have a good solution for wind
direction. In
the latest release, met-7.0 we added a new VCNT line type which
includes
vector wind statistics and correspond to the vector partial sums
included
in the VL1L2 line type. But those are computed by verifying U and V,
not
wind speed and direction.
Point-Stat and Grid-Stat do not compute statistics on wind direction
directly. And that's because it's a "circular" quantity. 1 degree
and 359
degrees are only 2 degrees different, but MET doesn't include special
logic
for handling that.
I'm not familiar with the SYNOP observations, but one potential option
would be using wind speed and direction to derive the corresponding U
and
V. Sorry I don't have a better solution for you.
Thanks,
John
On Tue, Apr 24, 2018 at 4:36 AM, Blandine Lheveder via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84858 >
>
> Dear John,
>
> Thank you a lot for your help.
> I have increased the login level to 3 and found the problem.
> Point_stat could not read the valid_time in the forecast file
(output of
> pcp_combine) because I had put a wrong argument " file_type =
NETCDF_NCCF;"
> for the fcst.
>
> Now it is OK :-)
>
> I have another question:
> I use SYNOP observations to validate with MET my WRF simulation. I
would
> like to validate the 10m wind. In the SYNOP I have only the wind
direction
> (0-360°, GRIB code 031) and wind speed (m/s, GRIB code 032). In the
wrf
> output, after the UPP post-processing, I do not get these fields but
only
> the UGRD and VGRD at Z10.
> I have seen that in stat_analysis it is possible to compute
statistics on
> wind direction, but only if pair of variables UGRD and VGRD are
processed
> in the Point_stat step.
>
> Do you have an idea how to proceed in this case ?
>
> Thanks a lot
>
> Best regards
>
> Blandine
>
> *Blandine L'Hévéder*
> Météo - Océan -Climat
>
> blandine.lheveder at gmail.com
> tel. 06 62 53 32 37
> https://fr.linkedin.com/in/blandine-l-heveder-493171a9
>
> 2018-04-23 22:37 GMT+02:00 John Halley Gotway via RT
<met_help at ucar.edu>:
>
> > Hello Blandine,
> >
> > I see that you have questions about using Point-Stat to compute
> statistics
> > for accumulated precipitation. Thanks for sending your Point-Stat
config
> > file and log from Point-Stat. My first suggestion would be
increasing
> the
> > logging level when you run Point-Stat. Try passing in "-v 3" on
the
> > command line to increase the logging level to 3. That'll include
log
> > messages that look something like this:
> >
> > DEBUG 3: Number of matched pairs = 132
> > DEBUG 3: Observations processed = 89893
> > DEBUG 3: Rejected: SID exclusion = 0
> > DEBUG 3: Rejected: obs type = 79788
> > DEBUG 3: Rejected: valid time = 0
> > DEBUG 3: Rejected: bad obs value = 0
> > DEBUG 3: Rejected: off the grid = 5
> > DEBUG 3: Rejected: level mismatch = 9219
> > DEBUG 3: Rejected: quality marker = 0
> > DEBUG 3: Rejected: message type = 344
> > DEBUG 3: Rejected: masking region = 405
> > DEBUG 3: Rejected: bad fcst value = 0
> > DEBUG 3: Rejected: duplicates = 0
> >
> > For each verification task, you'll see these counts describing the
number
> > of matched pairs that were used to compute statistics... as well
as
> counts
> > for the different reasons why observations were *not* used for
this task.
> >
> > Since precipitation is accumulated in time, it is more complicated
to
> > verify. Ultimately, you're limited to verifying only at the
accumulation
> > intervals for which you have observations. So the real question
is, what
> > point observation datasets are you using and accumulation
intervals are
> > present in them?
> >
> > If you'd like, you could post a sample observation file to our
anonymous
> > ftp site along with the PB2NC (presumably) config file (or perhaps
> they're
> > ascii?) you're using to process it:
> > https://dtcenter.org/met/users/support/met_help.php#ftp
> >
> > Thanks,
> > John
> >
> > On Mon, Apr 23, 2018 at 4:53 AM, Blandine Lheveder via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > >
> > > Mon Apr 23 04:53:34 2018: Request 84858 was acted upon.
> > > Transaction: Ticket created by blandine.lheveder at gmail.com
> > > Queue: met_help
> > > Subject: How to manage APCP with different accumulation
period in
> > MET
> > > Owner: Nobody
> > > Requestors: blandine.lheveder at gmail.com
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84858
> >
> > >
> > >
> > > I am currently using UPP with hourly output of WRF (24 frames in
a
> file)
> > in
> > > order to be able to validate my WRF simulation with MET tool. So
I get
> > > hourly grib1 WRFPRS files, each one containig a field with APCP
at the
> > > level A01.
> > >
> > > I would like to validate this field with SYNOP data. In these
data, I
> > have
> > > only APCP data at the level A03 (cumul over 3hrs).
> > >
> > > I use point_stat and try to test the hourly WRF field:
> > > {
> > > name = "APCP";
> > > level = [ "A01","A03","A06" ] ;
> > > },
> > > I use the same field definition for fcst and obs.
> > >
> > > In the log, I see that some forecast data are found for APCP/A01
or
> > > APCP/A03.
> > > But I do not get data for APCP/A01 for every hour. It seems
strange
> > because
> > > I work with hourly WRF output. Is this normal ? How are APCP
with no
> > > precipitation (0 mm) processed in point_stat ?
> > > And how do I get forecast data for APCP/A03 while I'm running
> point_stat
> > > for each time file separately ???
> > >
> > > Could you explain me how point_stat is processing APCP ? There
is very
> > > little information in the doc.
> > >
> > > Thanks a lot for your help
> > >
> > > PS: I send you in copy the pointstat_config file and the log
after
> > > point_stat
> > >
> > > *Blandine L'Hévéder*
> > > Météo - Océan -Climat
> > >
> > > blandine.lheveder at gmail.com
> > > tel. 06 62 53 32 37
> > > https://fr.linkedin.com/in/blandine-l-heveder-493171a9
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: How to manage APCP with different accumulation period in MET
From: Blandine Lheveder
Time: Thu Apr 26 06:55:14 2018
Hi John,
Thanks for your help. Now I can validate wind speed.
I have received that Met-help history exchange. Could it work in my
case or
do I need to have pairs of U and V for forecast AND Obs ?
> Hi Jeff.
>
> Neither point_stat nor grid_stat evaluate wind direction directly
because
> it is a circular variable, however, here is a way that you could get
this
> type of information using stat_analysis, which may be useful to you:
>
> stat_analysis -lookin point_stat_360000L_20070331_120000V_mpr.met-
7.0.txt
> -job aggregate_stat -line_type MPR -out_line_type WDIR -by
> OBS_SID,OBS_VALID_BEG,OBTYPE,OBS_HGT,VX_MASK,INTERP_MTHD,INTERP_PNTS
> -out wdir_job.out
>
> In the output, wdir_job.out, you'll find the wind direction forecast
value
> in "FBAR" and the wind direction observation value in "OBAR".
>
> I have attached the input file I used (one for met-7.0 and one for
met-6.0
> depending on what version you are using) and the output file for you
to see
> the results.
>
> I hope this will help.
>
> Thanks,
> Julie
*Blandine L'Hévéder*
Météo - Océan -Climat
blandine.lheveder at gmail.com
tel. 06 62 53 32 37
https://fr.linkedin.com/in/blandine-l-heveder-493171a9
2018-04-24 17:55 GMT+02:00 John Halley Gotway via RT
<met_help at ucar.edu>:
> Blandine,
>
> OK, so your SYNOP observations include wind speed and direction
rather than
> U and V. And you're GRIB output from UPP contain U and V. For wind
speed,
> everything is fine. Just request that "WIND" be verified... that's
the
> GRIB code abbreviation for wind speed. When you request WIND from a
GRIB
> file, MET first searches for a matching record. If not found, it
then
> finds the corresponding U and V and automatically derives wind speed
for
> you.
>
> Unfortunately though, I do not have a good solution for wind
direction. In
> the latest release, met-7.0 we added a new VCNT line type which
includes
> vector wind statistics and correspond to the vector partial sums
included
> in the VL1L2 line type. But those are computed by verifying U and
V, not
> wind speed and direction.
>
> Point-Stat and Grid-Stat do not compute statistics on wind direction
> directly. And that's because it's a "circular" quantity. 1 degree
and 359
> degrees are only 2 degrees different, but MET doesn't include
special logic
> for handling that.
>
> I'm not familiar with the SYNOP observations, but one potential
option
> would be using wind speed and direction to derive the corresponding
U and
> V. Sorry I don't have a better solution for you.
>
> Thanks,
> John
>
>
> On Tue, Apr 24, 2018 at 4:36 AM, Blandine Lheveder via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84858 >
> >
> > Dear John,
> >
> > Thank you a lot for your help.
> > I have increased the login level to 3 and found the problem.
> > Point_stat could not read the valid_time in the forecast file
(output of
> > pcp_combine) because I had put a wrong argument " file_type =
> NETCDF_NCCF;"
> > for the fcst.
> >
> > Now it is OK :-)
> >
> > I have another question:
> > I use SYNOP observations to validate with MET my WRF simulation. I
would
> > like to validate the 10m wind. In the SYNOP I have only the wind
> direction
> > (0-360°, GRIB code 031) and wind speed (m/s, GRIB code 032). In
the wrf
> > output, after the UPP post-processing, I do not get these fields
but only
> > the UGRD and VGRD at Z10.
> > I have seen that in stat_analysis it is possible to compute
statistics on
> > wind direction, but only if pair of variables UGRD and VGRD are
processed
> > in the Point_stat step.
> >
> > Do you have an idea how to proceed in this case ?
> >
> > Thanks a lot
> >
> > Best regards
> >
> > Blandine
> >
> > *Blandine L'Hévéder*
> > Météo - Océan -Climat
> >
> > blandine.lheveder at gmail.com
> > tel. 06 62 53 32 37
> > https://fr.linkedin.com/in/blandine-l-heveder-493171a9
> >
> > 2018-04-23 22:37 GMT+02:00 John Halley Gotway via RT
<met_help at ucar.edu
> >:
> >
> > > Hello Blandine,
> > >
> > > I see that you have questions about using Point-Stat to compute
> > statistics
> > > for accumulated precipitation. Thanks for sending your Point-
Stat
> config
> > > file and log from Point-Stat. My first suggestion would be
increasing
> > the
> > > logging level when you run Point-Stat. Try passing in "-v 3" on
the
> > > command line to increase the logging level to 3. That'll
include log
> > > messages that look something like this:
> > >
> > > DEBUG 3: Number of matched pairs = 132
> > > DEBUG 3: Observations processed = 89893
> > > DEBUG 3: Rejected: SID exclusion = 0
> > > DEBUG 3: Rejected: obs type = 79788
> > > DEBUG 3: Rejected: valid time = 0
> > > DEBUG 3: Rejected: bad obs value = 0
> > > DEBUG 3: Rejected: off the grid = 5
> > > DEBUG 3: Rejected: level mismatch = 9219
> > > DEBUG 3: Rejected: quality marker = 0
> > > DEBUG 3: Rejected: message type = 344
> > > DEBUG 3: Rejected: masking region = 405
> > > DEBUG 3: Rejected: bad fcst value = 0
> > > DEBUG 3: Rejected: duplicates = 0
> > >
> > > For each verification task, you'll see these counts describing
the
> number
> > > of matched pairs that were used to compute statistics... as well
as
> > counts
> > > for the different reasons why observations were *not* used for
this
> task.
> > >
> > > Since precipitation is accumulated in time, it is more
complicated to
> > > verify. Ultimately, you're limited to verifying only at the
> accumulation
> > > intervals for which you have observations. So the real question
is,
> what
> > > point observation datasets are you using and accumulation
intervals are
> > > present in them?
> > >
> > > If you'd like, you could post a sample observation file to our
> anonymous
> > > ftp site along with the PB2NC (presumably) config file (or
perhaps
> > they're
> > > ascii?) you're using to process it:
> > > https://dtcenter.org/met/users/support/met_help.php#ftp
> > >
> > > Thanks,
> > > John
> > >
> > > On Mon, Apr 23, 2018 at 4:53 AM, Blandine Lheveder via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > >
> > > > Mon Apr 23 04:53:34 2018: Request 84858 was acted upon.
> > > > Transaction: Ticket created by blandine.lheveder at gmail.com
> > > > Queue: met_help
> > > > Subject: How to manage APCP with different accumulation
period
> in
> > > MET
> > > > Owner: Nobody
> > > > Requestors: blandine.lheveder at gmail.com
> > > > Status: new
> > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> Ticket/Display.html?id=84858
> > >
> > > >
> > > >
> > > > I am currently using UPP with hourly output of WRF (24 frames
in a
> > file)
> > > in
> > > > order to be able to validate my WRF simulation with MET tool.
So I
> get
> > > > hourly grib1 WRFPRS files, each one containig a field with
APCP at
> the
> > > > level A01.
> > > >
> > > > I would like to validate this field with SYNOP data. In these
data, I
> > > have
> > > > only APCP data at the level A03 (cumul over 3hrs).
> > > >
> > > > I use point_stat and try to test the hourly WRF field:
> > > > {
> > > > name = "APCP";
> > > > level = [ "A01","A03","A06" ] ;
> > > > },
> > > > I use the same field definition for fcst and obs.
> > > >
> > > > In the log, I see that some forecast data are found for
APCP/A01 or
> > > > APCP/A03.
> > > > But I do not get data for APCP/A01 for every hour. It seems
strange
> > > because
> > > > I work with hourly WRF output. Is this normal ? How are APCP
with no
> > > > precipitation (0 mm) processed in point_stat ?
> > > > And how do I get forecast data for APCP/A03 while I'm running
> > point_stat
> > > > for each time file separately ???
> > > >
> > > > Could you explain me how point_stat is processing APCP ? There
is
> very
> > > > little information in the doc.
> > > >
> > > > Thanks a lot for your help
> > > >
> > > > PS: I send you in copy the pointstat_config file and the log
after
> > > > point_stat
> > > >
> > > > *Blandine L'Hévéder*
> > > > Météo - Océan -Climat
> > > >
> > > > blandine.lheveder at gmail.com
> > > > tel. 06 62 53 32 37
> > > > https://fr.linkedin.com/in/blandine-l-heveder-493171a9
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: How to manage APCP with different accumulation period in MET
From: John Halley Gotway
Time: Thu Apr 26 09:12:15 2018
Blandine,
Unfortunately, that logic for wind direction evaluation in STAT-
Analysis is
done using the U/V matched pairs. So I don't that that functionality
will
be useful to you without U/V pairs.
Are you able to send me a sample file of the SYNOP observations that
you're
using? I'm curious to see them.
Thanks,
John
On Thu, Apr 26, 2018 at 6:55 AM, Blandine Lheveder via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84858 >
>
> Hi John,
>
> Thanks for your help. Now I can validate wind speed.
> I have received that Met-help history exchange. Could it work in my
case or
> do I need to have pairs of U and V for forecast AND Obs ?
>
> > Hi Jeff.
> >
> > Neither point_stat nor grid_stat evaluate wind direction directly
> because
> > it is a circular variable, however, here is a way that you could
get
> this
> > type of information using stat_analysis, which may be useful to
you:
> >
> > stat_analysis -lookin point_stat_360000L_20070331_120000V_mpr.met-
> 7.0.txt
> > -job aggregate_stat -line_type MPR -out_line_type WDIR -by
> >
OBS_SID,OBS_VALID_BEG,OBTYPE,OBS_HGT,VX_MASK,INTERP_MTHD,INTERP_PNTS
> > -out wdir_job.out
> >
> > In the output, wdir_job.out, you'll find the wind direction
forecast
> value
> > in "FBAR" and the wind direction observation value in "OBAR".
> >
> > I have attached the input file I used (one for met-7.0 and one for
> met-6.0
> > depending on what version you are using) and the output file for
you
> to see
> > the results.
> >
> > I hope this will help.
> >
> > Thanks,
> > Julie
>
> *Blandine L'Hévéder*
> Météo - Océan -Climat
>
> blandine.lheveder at gmail.com
> tel. 06 62 53 32 37
> https://fr.linkedin.com/in/blandine-l-heveder-493171a9
>
> 2018-04-24 17:55 GMT+02:00 John Halley Gotway via RT
<met_help at ucar.edu>:
>
> > Blandine,
> >
> > OK, so your SYNOP observations include wind speed and direction
rather
> than
> > U and V. And you're GRIB output from UPP contain U and V. For
wind
> speed,
> > everything is fine. Just request that "WIND" be verified...
that's the
> > GRIB code abbreviation for wind speed. When you request WIND from
a GRIB
> > file, MET first searches for a matching record. If not found, it
then
> > finds the corresponding U and V and automatically derives wind
speed for
> > you.
> >
> > Unfortunately though, I do not have a good solution for wind
direction.
> In
> > the latest release, met-7.0 we added a new VCNT line type which
includes
> > vector wind statistics and correspond to the vector partial sums
included
> > in the VL1L2 line type. But those are computed by verifying U and
V, not
> > wind speed and direction.
> >
> > Point-Stat and Grid-Stat do not compute statistics on wind
direction
> > directly. And that's because it's a "circular" quantity. 1
degree and
> 359
> > degrees are only 2 degrees different, but MET doesn't include
special
> logic
> > for handling that.
> >
> > I'm not familiar with the SYNOP observations, but one potential
option
> > would be using wind speed and direction to derive the
corresponding U and
> > V. Sorry I don't have a better solution for you.
> >
> > Thanks,
> > John
> >
> >
> > On Tue, Apr 24, 2018 at 4:36 AM, Blandine Lheveder via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=84858 >
> > >
> > > Dear John,
> > >
> > > Thank you a lot for your help.
> > > I have increased the login level to 3 and found the problem.
> > > Point_stat could not read the valid_time in the forecast file
(output
> of
> > > pcp_combine) because I had put a wrong argument " file_type =
> > NETCDF_NCCF;"
> > > for the fcst.
> > >
> > > Now it is OK :-)
> > >
> > > I have another question:
> > > I use SYNOP observations to validate with MET my WRF simulation.
I
> would
> > > like to validate the 10m wind. In the SYNOP I have only the wind
> > direction
> > > (0-360°, GRIB code 031) and wind speed (m/s, GRIB code 032). In
the wrf
> > > output, after the UPP post-processing, I do not get these fields
but
> only
> > > the UGRD and VGRD at Z10.
> > > I have seen that in stat_analysis it is possible to compute
statistics
> on
> > > wind direction, but only if pair of variables UGRD and VGRD are
> processed
> > > in the Point_stat step.
> > >
> > > Do you have an idea how to proceed in this case ?
> > >
> > > Thanks a lot
> > >
> > > Best regards
> > >
> > > Blandine
> > >
> > > *Blandine L'Hévéder*
> > > Météo - Océan -Climat
> > >
> > > blandine.lheveder at gmail.com
> > > tel. 06 62 53 32 37
> > > https://fr.linkedin.com/in/blandine-l-heveder-493171a9
> > >
> > > 2018-04-23 22:37 GMT+02:00 John Halley Gotway via RT <
> met_help at ucar.edu
> > >:
> > >
> > > > Hello Blandine,
> > > >
> > > > I see that you have questions about using Point-Stat to
compute
> > > statistics
> > > > for accumulated precipitation. Thanks for sending your Point-
Stat
> > config
> > > > file and log from Point-Stat. My first suggestion would be
> increasing
> > > the
> > > > logging level when you run Point-Stat. Try passing in "-v 3"
on the
> > > > command line to increase the logging level to 3. That'll
include log
> > > > messages that look something like this:
> > > >
> > > > DEBUG 3: Number of matched pairs = 132
> > > > DEBUG 3: Observations processed = 89893
> > > > DEBUG 3: Rejected: SID exclusion = 0
> > > > DEBUG 3: Rejected: obs type = 79788
> > > > DEBUG 3: Rejected: valid time = 0
> > > > DEBUG 3: Rejected: bad obs value = 0
> > > > DEBUG 3: Rejected: off the grid = 5
> > > > DEBUG 3: Rejected: level mismatch = 9219
> > > > DEBUG 3: Rejected: quality marker = 0
> > > > DEBUG 3: Rejected: message type = 344
> > > > DEBUG 3: Rejected: masking region = 405
> > > > DEBUG 3: Rejected: bad fcst value = 0
> > > > DEBUG 3: Rejected: duplicates = 0
> > > >
> > > > For each verification task, you'll see these counts describing
the
> > number
> > > > of matched pairs that were used to compute statistics... as
well as
> > > counts
> > > > for the different reasons why observations were *not* used for
this
> > task.
> > > >
> > > > Since precipitation is accumulated in time, it is more
complicated to
> > > > verify. Ultimately, you're limited to verifying only at the
> > accumulation
> > > > intervals for which you have observations. So the real
question is,
> > what
> > > > point observation datasets are you using and accumulation
intervals
> are
> > > > present in them?
> > > >
> > > > If you'd like, you could post a sample observation file to our
> > anonymous
> > > > ftp site along with the PB2NC (presumably) config file (or
perhaps
> > > they're
> > > > ascii?) you're using to process it:
> > > > https://dtcenter.org/met/users/support/met_help.php#ftp
> > > >
> > > > Thanks,
> > > > John
> > > >
> > > > On Mon, Apr 23, 2018 at 4:53 AM, Blandine Lheveder via RT <
> > > > met_help at ucar.edu
> > > > > wrote:
> > > >
> > > > >
> > > > > Mon Apr 23 04:53:34 2018: Request 84858 was acted upon.
> > > > > Transaction: Ticket created by blandine.lheveder at gmail.com
> > > > > Queue: met_help
> > > > > Subject: How to manage APCP with different accumulation
period
> > in
> > > > MET
> > > > > Owner: Nobody
> > > > > Requestors: blandine.lheveder at gmail.com
> > > > > Status: new
> > > > > Ticket <URL: https://rt.rap.ucar.edu/rt/
> > Ticket/Display.html?id=84858
> > > >
> > > > >
> > > > >
> > > > > I am currently using UPP with hourly output of WRF (24
frames in a
> > > file)
> > > > in
> > > > > order to be able to validate my WRF simulation with MET
tool. So I
> > get
> > > > > hourly grib1 WRFPRS files, each one containig a field with
APCP at
> > the
> > > > > level A01.
> > > > >
> > > > > I would like to validate this field with SYNOP data. In
these
> data, I
> > > > have
> > > > > only APCP data at the level A03 (cumul over 3hrs).
> > > > >
> > > > > I use point_stat and try to test the hourly WRF field:
> > > > > {
> > > > > name = "APCP";
> > > > > level = [ "A01","A03","A06" ] ;
> > > > > },
> > > > > I use the same field definition for fcst and obs.
> > > > >
> > > > > In the log, I see that some forecast data are found for
APCP/A01 or
> > > > > APCP/A03.
> > > > > But I do not get data for APCP/A01 for every hour. It seems
strange
> > > > because
> > > > > I work with hourly WRF output. Is this normal ? How are APCP
with
> no
> > > > > precipitation (0 mm) processed in point_stat ?
> > > > > And how do I get forecast data for APCP/A03 while I'm
running
> > > point_stat
> > > > > for each time file separately ???
> > > > >
> > > > > Could you explain me how point_stat is processing APCP ?
There is
> > very
> > > > > little information in the doc.
> > > > >
> > > > > Thanks a lot for your help
> > > > >
> > > > > PS: I send you in copy the pointstat_config file and the log
after
> > > > > point_stat
> > > > >
> > > > > *Blandine L'Hévéder*
> > > > > Météo - Océan -Climat
> > > > >
> > > > > blandine.lheveder at gmail.com
> > > > > tel. 06 62 53 32 37
> > > > > https://fr.linkedin.com/in/blandine-l-heveder-493171a9
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
More information about the Met_help
mailing list