[Met_help] [rt.rap.ucar.edu #86716] History for FW: MPR File Format For HIRA

[John Halley Gotway via RT]

----------------------------------------------------------------
  Initial Request
----------------------------------------------------------------

John, I think this might not have made it to you yesterday so I am resending.

Thanks Bob

-----Original Message-----
From: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN 
Sent: Tuesday, August 21, 2018 3:35 PM
To: 'met_help at ucar.edu' <met_help at ucar.edu>
Subject: MPR File Format For HIRA

John, since MET can't handle rotated grids, we are generating the MPR files
with neighborhoods outside of MET and then having stat_analysis read them to
generate pstd stats (brier scrores and brier skill scores).    Attached is a
MPR file.  When I run stat_analysis on this file I get a brier score but the
briercl and bss are NA.  We are using persistence as the climo reference and
I have been putting a 1 or 0 in the climo mean field indicating if the
persistence forecast had a precipitation value greater than or equal to 1
hundreth of an inch ( the threshold we are using).  Since the climo columns
in the pstd file are all NA, obviously I am not presenting the climo value
to MET properly.   What should the MPR file contain in the climo mean,
stdev, or cdf column to get MET to us it properly?  The brier scores
calculate out close to what they should be so the problem seems to be just
with the climo representation.

Thanks
Bob Craig


----------------------------------------------------------------
  Complete Ticket History
----------------------------------------------------------------

Subject: FW: MPR File Format For HIRA
From: John Halley Gotway
Time: Wed Aug 22 16:26:46 2018

Bob,

FYI, we are working on adding support for rotated lat/lon grids for
the
next release of met, met-8.0.  And we recently made progress on our
Fortify
scans.  We're now able to run them in an automated way every night.
So
we're now ready to get going addressing the Fortify issues.  But we
won't
be able to finish before the met-8.0 code freeze next week.  All of
this is
to say that met-8.1 will include all the Fortify fixes... and support
for
rotated lat/lon grids.

Thanks for sending your sample file.  That helps with the debugging.

OK, so good news and bad news.

The bad news is that this is a bug.  In the file named
"parse_stat_line.cc"
in the function "parse_mpr_line()", we are looking for a column of
data
named "CLIMO" to get the climatology mean value:
   m_data.climo    = atof(l.get_item("CLIMO"));

In met-6.0 and earlier versions, that was the correct column name.  In
met-6.1 and future versions, we switched to writing:
   CLIMO_MEAN CLIMO_STDEV CLIMO_CDF

Now the good news.  With the existing code, if you change the column
name
in your .stat file from "CLIMO_MEAN" to "CLIMO" then it all works.

I need to think some more about how the code should handle this
difference
in column names.

Thanks,
John

On Wed, Aug 22, 2018 at 9:37 AM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:

>
> Wed Aug 22 09:37:21 2018: Request 86716 was acted upon.
> Transaction: Ticket created by robert.craig.2 at us.af.mil
>        Queue: met_help
>      Subject: FW: MPR File Format For HIRA
>        Owner: Nobody
>   Requestors: robert.craig.2 at us.af.mil
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=86716 >
>
>
> John, I think this might not have made it to you yesterday so I am
> resending.
>
> Thanks Bob
>
> -----Original Message-----
> From: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
> Sent: Tuesday, August 21, 2018 3:35 PM
> To: 'met_help at ucar.edu' <met_help at ucar.edu>
> Subject: MPR File Format For HIRA
>
> John, since MET can't handle rotated grids, we are generating the
MPR files
> with neighborhoods outside of MET and then having stat_analysis read
them
> to
> generate pstd stats (brier scrores and brier skill scores).
Attached is
> a
> MPR file.  When I run stat_analysis on this file I get a brier score
but
> the
> briercl and bss are NA.  We are using persistence as the climo
reference
> and
> I have been putting a 1 or 0 in the climo mean field indicating if
the
> persistence forecast had a precipitation value greater than or equal
to 1
> hundreth of an inch ( the threshold we are using).  Since the climo
columns
> in the pstd file are all NA, obviously I am not presenting the climo
value
> to MET properly.   What should the MPR file contain in the climo
mean,
> stdev, or cdf column to get MET to us it properly?  The brier scores
> calculate out close to what they should be so the problem seems to
be just
> with the climo representation.
>
> Thanks
> Bob Craig
>
>

------------------------------------------------