[Met_help] [rt.rap.ucar.edu #95065] History for PJC output

John Halley Gotway via RT met_help at ucar.edu
Mon Apr 27 16:27:13 MDT 2020


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

Folks - Where are the OY_TP_i, ON_TP_i, and LIKELIHOOD_i values for probability of 1.0 included? On my PJC output files, I only see values for THRESH_n equal to 0.9. Thanks.

Efren A. Serra (Contractor)
Physicist

DeVine Consulting, Inc.
Naval Research Laboratory
Marine Meteorology Division
7 Grace Hopper Ave., STOP 2
Monterey, CA 93943
Code 7542
Mobile: 408-425-5027



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

Subject: PJC output
From: John Halley Gotway
Time: Mon Apr 27 15:48:43 2020

Efren,

The pairs with probabilities of 1.0 *should* be included in the last
probability bin.  When you set "cat_thresh = [ ==0.1 ];", MET
automatically
expands that into 10 probability bins:
bin1: 0.0 <= p < 0.1
bin2: 0.1 <= p < 0.2
bin3: 0.2 <= p < 0.3
bin4: 0.3 <= p < 0.4
bin5: 0.4 <= p < 0.5
bin6: 0.5 <= p < 0.6
bin7: 0.6 <= p < 0.7
bin8: 0.7 <= p < 0.8
bin9: 0.8 <= p < 0.9
bin10: 0.9 <= p <= 1.0

This should be occurring in the function named "value_to_row()" in
this
file:
https://github.com/NCAR/MET/blob/master_v9.0/met/src/libcode/vx_statistics/contable_nx2.cc

Specifically, on line 234, we check to see if the probability value
equals
the last threshold (e.g. 1.0), and then assign it to the last bin:
if ( is_eq(t, Thresholds[Nrows]) ) return ( Nrows - 1 );

But it sounds like that's not the behavior you're seeing. If you send
me
some sample data to illustrate with a detailed description of what's
not
working as expected, I could take a look.

FYI, if needed you can post data to our anonymous ftp site following
these
instructions:
https://dtcenter.org/community-code/model-evaluation-tools-met/met-
help-desk#ftp

Thanks,
John Halley Gotway

On Mon, Apr 27, 2020 at 1:22 PM efren.serra.ctr at nrlmry.navy.mil via RT
<
met_help at ucar.edu> wrote:

>
> Mon Apr 27 13:22:52 2020: Request 95065 was acted upon.
> Transaction: Ticket created by efren.serra.ctr at nrlmry.navy.mil
>        Queue: met_help
>      Subject: PJC output
>        Owner: Nobody
>   Requestors: efren.serra.ctr at nrlmry.navy.mil
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=95065 >
>
>
> Folks - Where are the OY_TP_i, ON_TP_i, and LIKELIHOOD_i values for
> probability of 1.0 included? On my PJC output files, I only see
values for
> THRESH_n equal to 0.9. Thanks.
>
> Efren A. Serra (Contractor)
> Physicist
>
> DeVine Consulting, Inc.
> Naval Research Laboratory
> Marine Meteorology Division
> 7 Grace Hopper Ave., STOP 2
> Monterey, CA 93943
> Code 7542
> Mobile: 408-425-5027
>
>
>

------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #95065] PJC output
From: efren.serra.ctr at nrlmry.navy.mil
Time: Mon Apr 27 16:11:18 2020

Buck - Does this answer your question?

-----Original Message-----
From: John Halley Gotway via RT <met_help at ucar.edu>
Sent: Monday, April 27, 2020 2:49 PM
To: Serra, Mr. Efren, Contractor, Code 7531
<efren.serra.ctr at nrlmry.navy.mil>
Subject: Re: [rt.rap.ucar.edu #95065] PJC output

Efren,

The pairs with probabilities of 1.0 *should* be included in the last
probability bin.  When you set "cat_thresh = [ ==0.1 ];", MET
automatically expands that into 10 probability bins:
bin1: 0.0 <= p < 0.1
bin2: 0.1 <= p < 0.2
bin3: 0.2 <= p < 0.3
bin4: 0.3 <= p < 0.4
bin5: 0.4 <= p < 0.5
bin6: 0.5 <= p < 0.6
bin7: 0.6 <= p < 0.7
bin8: 0.7 <= p < 0.8
bin9: 0.8 <= p < 0.9
bin10: 0.9 <= p <= 1.0

This should be occurring in the function named "value_to_row()" in
this
file:
https://github.com/NCAR/MET/blob/master_v9.0/met/src/libcode/vx_statistics/contable_nx2.cc

Specifically, on line 234, we check to see if the probability value
equals the last threshold (e.g. 1.0), and then assign it to the last
bin:
if ( is_eq(t, Thresholds[Nrows]) ) return ( Nrows - 1 );

But it sounds like that's not the behavior you're seeing. If you send
me some sample data to illustrate with a detailed description of
what's not working as expected, I could take a look.

FYI, if needed you can post data to our anonymous ftp site following
these
instructions:
https://dtcenter.org/community-code/model-evaluation-tools-met/met-
help-desk#ftp

Thanks,
John Halley Gotway

On Mon, Apr 27, 2020 at 1:22 PM efren.serra.ctr at nrlmry.navy.mil via RT
< met_help at ucar.edu> wrote:

>
> Mon Apr 27 13:22:52 2020: Request 95065 was acted upon.
> Transaction: Ticket created by efren.serra.ctr at nrlmry.navy.mil
>        Queue: met_help
>      Subject: PJC output
>        Owner: Nobody
>   Requestors: efren.serra.ctr at nrlmry.navy.mil
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=95065
> >
>
>
> Folks - Where are the OY_TP_i, ON_TP_i, and LIKELIHOOD_i values for
> probability of 1.0 included? On my PJC output files, I only see
values
> for THRESH_n equal to 0.9. Thanks.
>
> Efren A. Serra (Contractor)
> Physicist
>
> DeVine Consulting, Inc.
> Naval Research Laboratory
> Marine Meteorology Division
> 7 Grace Hopper Ave., STOP 2
> Monterey, CA 93943
> Code 7542
> Mobile: 408-425-5027
>
>
>



------------------------------------------------
Subject: PJC output
From: Sampson, Mr. Buck
Time: Mon Apr 27 16:15:08 2020

It certainly does.  It was just a question that I figured Efren could
verify
by testing a field with a lot of 100% in it, and we have a lot of
those.
Thanks for the assist, John.  Very helpful.


Buck Sampson
NRL Monterey
7 Grace Hopper Ave, Monterey CA 93943
Ph: 831 656-4714
Cell: 831 359-8114
https://scholar.google.com/citations?hl=en&pli=1&user=nhNwujgAAAAJ



-----Original Message-----
From: Serra, Mr. Efren, Contractor, Code 7531
<efren.serra.ctr at nrlmry.navy.mil>
Sent: Monday, April 27, 2020 3:11 PM
To: met_help at ucar.edu; Sampson, Mr. Buck
<Buck.Sampson at nrlmry.navy.mil>
Subject: RE: [rt.rap.ucar.edu #95065] PJC output

Buck - Does this answer your question?

-----Original Message-----
From: John Halley Gotway via RT <met_help at ucar.edu>
Sent: Monday, April 27, 2020 2:49 PM
To: Serra, Mr. Efren, Contractor, Code 7531
<efren.serra.ctr at nrlmry.navy.mil>
Subject: Re: [rt.rap.ucar.edu #95065] PJC output

Efren,

The pairs with probabilities of 1.0 *should* be included in the last
probability bin.  When you set "cat_thresh = [ ==0.1 ];", MET
automatically
expands that into 10 probability bins:
bin1: 0.0 <= p < 0.1
bin2: 0.1 <= p < 0.2
bin3: 0.2 <= p < 0.3
bin4: 0.3 <= p < 0.4
bin5: 0.4 <= p < 0.5
bin6: 0.5 <= p < 0.6
bin7: 0.6 <= p < 0.7
bin8: 0.7 <= p < 0.8
bin9: 0.8 <= p < 0.9
bin10: 0.9 <= p <= 1.0

This should be occurring in the function named "value_to_row()" in
this
file:
https://github.com/NCAR/MET/blob/master_v9.0/met/src/libcode/vx_statistics/contable_nx2.cc

Specifically, on line 234, we check to see if the probability value
equals the
last threshold (e.g. 1.0), and then assign it to the last bin:
if ( is_eq(t, Thresholds[Nrows]) ) return ( Nrows - 1 );

But it sounds like that's not the behavior you're seeing. If you send
me some
sample data to illustrate with a detailed description of what's not
working as
expected, I could take a look.

FYI, if needed you can post data to our anonymous ftp site following
these
instructions:
https://dtcenter.org/community-code/model-evaluation-tools-met/met-
help-desk#ftp

Thanks,
John Halley Gotway

On Mon, Apr 27, 2020 at 1:22 PM efren.serra.ctr at nrlmry.navy.mil via RT
<
met_help at ucar.edu> wrote:

>
> Mon Apr 27 13:22:52 2020: Request 95065 was acted upon.
> Transaction: Ticket created by efren.serra.ctr at nrlmry.navy.mil
>        Queue: met_help
>      Subject: PJC output
>        Owner: Nobody
>   Requestors: efren.serra.ctr at nrlmry.navy.mil
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=95065
> >
>
>
> Folks - Where are the OY_TP_i, ON_TP_i, and LIKELIHOOD_i values for
> probability of 1.0 included? On my PJC output files, I only see
values
> for THRESH_n equal to 0.9. Thanks.
>
> Efren A. Serra (Contractor)
> Physicist
>
> DeVine Consulting, Inc.
> Naval Research Laboratory
> Marine Meteorology Division
> 7 Grace Hopper Ave., STOP 2
> Monterey, CA 93943
> Code 7542
> Mobile: 408-425-5027
>
>
>


------------------------------------------------
Subject: PJC output
From: efren.serra.ctr at nrlmry.navy.mil
Time: Mon Apr 27 16:17:44 2020

I'm on to the testing of that case; thanks for reply John.

-----Original Message-----
From: Sampson, Mr. Buck <Buck.Sampson at nrlmry.navy.mil>
Sent: Monday, April 27, 2020 3:15 PM
To: Serra, Mr. Efren, Contractor, Code 7531
<efren.serra.ctr at nrlmry.navy.mil>;
met_help at ucar.edu
Subject: RE: [rt.rap.ucar.edu #95065] PJC output

It certainly does.  It was just a question that I figured Efren could
verify
by testing a field with a lot of 100% in it, and we have a lot of
those.
Thanks for the assist, John.  Very helpful.


Buck Sampson
NRL Monterey
7 Grace Hopper Ave, Monterey CA 93943
Ph: 831 656-4714
Cell: 831 359-8114
https://scholar.google.com/citations?hl=en&pli=1&user=nhNwujgAAAAJ



-----Original Message-----
From: Serra, Mr. Efren, Contractor, Code 7531
<efren.serra.ctr at nrlmry.navy.mil>
Sent: Monday, April 27, 2020 3:11 PM
To: met_help at ucar.edu; Sampson, Mr. Buck
<Buck.Sampson at nrlmry.navy.mil>
Subject: RE: [rt.rap.ucar.edu #95065] PJC output

Buck - Does this answer your question?

-----Original Message-----
From: John Halley Gotway via RT <met_help at ucar.edu>
Sent: Monday, April 27, 2020 2:49 PM
To: Serra, Mr. Efren, Contractor, Code 7531
<efren.serra.ctr at nrlmry.navy.mil>
Subject: Re: [rt.rap.ucar.edu #95065] PJC output

Efren,

The pairs with probabilities of 1.0 *should* be included in the last
probability bin.  When you set "cat_thresh = [ ==0.1 ];", MET
automatically
expands that into 10 probability bins:
bin1: 0.0 <= p < 0.1
bin2: 0.1 <= p < 0.2
bin3: 0.2 <= p < 0.3
bin4: 0.3 <= p < 0.4
bin5: 0.4 <= p < 0.5
bin6: 0.5 <= p < 0.6
bin7: 0.6 <= p < 0.7
bin8: 0.7 <= p < 0.8
bin9: 0.8 <= p < 0.9
bin10: 0.9 <= p <= 1.0

This should be occurring in the function named "value_to_row()" in
this
file:
https://github.com/NCAR/MET/blob/master_v9.0/met/src/libcode/vx_statistics/contable_nx2.cc

Specifically, on line 234, we check to see if the probability value
equals the
last threshold (e.g. 1.0), and then assign it to the last bin:
if ( is_eq(t, Thresholds[Nrows]) ) return ( Nrows - 1 );

But it sounds like that's not the behavior you're seeing. If you send
me some
sample data to illustrate with a detailed description of what's not
working as
expected, I could take a look.

FYI, if needed you can post data to our anonymous ftp site following
these
instructions:
https://dtcenter.org/community-code/model-evaluation-tools-met/met-
help-desk#ftp

Thanks,
John Halley Gotway

On Mon, Apr 27, 2020 at 1:22 PM efren.serra.ctr at nrlmry.navy.mil via RT
<
met_help at ucar.edu> wrote:

>
> Mon Apr 27 13:22:52 2020: Request 95065 was acted upon.
> Transaction: Ticket created by efren.serra.ctr at nrlmry.navy.mil
>        Queue: met_help
>      Subject: PJC output
>        Owner: Nobody
>   Requestors: efren.serra.ctr at nrlmry.navy.mil
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=95065
> >
>
>
> Folks - Where are the OY_TP_i, ON_TP_i, and LIKELIHOOD_i values for
> probability of 1.0 included? On my PJC output files, I only see
values
> for THRESH_n equal to 0.9. Thanks.
>
> Efren A. Serra (Contractor)
> Physicist
>
> DeVine Consulting, Inc.
> Naval Research Laboratory
> Marine Meteorology Division
> 7 Grace Hopper Ave., STOP 2
> Monterey, CA 93943
> Code 7542
> Mobile: 408-425-5027
>
>
>


------------------------------------------------
Subject: PJC output
From: John Halley Gotway
Time: Mon Apr 27 16:27:07 2020

Sounds good.  I'll go ahead and resolve this issue.

Just let us know what other questions or issues arise in your use of
MET.

Thanks,
John

On Mon, Apr 27, 2020 at 4:18 PM efren.serra.ctr at nrlmry.navy.mil via RT
<
met_help at ucar.edu> wrote:

>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=95065 >
>
> I'm on to the testing of that case; thanks for reply John.
>
> -----Original Message-----
> From: Sampson, Mr. Buck <Buck.Sampson at nrlmry.navy.mil>
> Sent: Monday, April 27, 2020 3:15 PM
> To: Serra, Mr. Efren, Contractor, Code 7531 <
> efren.serra.ctr at nrlmry.navy.mil>;
> met_help at ucar.edu
> Subject: RE: [rt.rap.ucar.edu #95065] PJC output
>
> It certainly does.  It was just a question that I figured Efren
could
> verify
> by testing a field with a lot of 100% in it, and we have a lot of
those.
> Thanks for the assist, John.  Very helpful.
>
>
> Buck Sampson
> NRL Monterey
> 7 Grace Hopper Ave, Monterey CA 93943
> Ph: 831 656-4714
> Cell: 831 359-8114
> https://scholar.google.com/citations?hl=en&pli=1&user=nhNwujgAAAAJ
>
>
>
> -----Original Message-----
> From: Serra, Mr. Efren, Contractor, Code 7531
> <efren.serra.ctr at nrlmry.navy.mil>
> Sent: Monday, April 27, 2020 3:11 PM
> To: met_help at ucar.edu; Sampson, Mr. Buck
<Buck.Sampson at nrlmry.navy.mil>
> Subject: RE: [rt.rap.ucar.edu #95065] PJC output
>
> Buck - Does this answer your question?
>
> -----Original Message-----
> From: John Halley Gotway via RT <met_help at ucar.edu>
> Sent: Monday, April 27, 2020 2:49 PM
> To: Serra, Mr. Efren, Contractor, Code 7531 <
> efren.serra.ctr at nrlmry.navy.mil>
> Subject: Re: [rt.rap.ucar.edu #95065] PJC output
>
> Efren,
>
> The pairs with probabilities of 1.0 *should* be included in the last
> probability bin.  When you set "cat_thresh = [ ==0.1 ];", MET
automatically
> expands that into 10 probability bins:
> bin1: 0.0 <= p < 0.1
> bin2: 0.1 <= p < 0.2
> bin3: 0.2 <= p < 0.3
> bin4: 0.3 <= p < 0.4
> bin5: 0.4 <= p < 0.5
> bin6: 0.5 <= p < 0.6
> bin7: 0.6 <= p < 0.7
> bin8: 0.7 <= p < 0.8
> bin9: 0.8 <= p < 0.9
> bin10: 0.9 <= p <= 1.0
>
> This should be occurring in the function named "value_to_row()" in
this
> file:
>
>
https://github.com/NCAR/MET/blob/master_v9.0/met/src/libcode/vx_statistics/contable_nx2.cc
>
> Specifically, on line 234, we check to see if the probability value
equals
> the
> last threshold (e.g. 1.0), and then assign it to the last bin:
> if ( is_eq(t, Thresholds[Nrows]) ) return ( Nrows - 1 );
>
> But it sounds like that's not the behavior you're seeing. If you
send me
> some
> sample data to illustrate with a detailed description of what's not
> working as
> expected, I could take a look.
>
> FYI, if needed you can post data to our anonymous ftp site following
these
> instructions:
>
> https://dtcenter.org/community-code/model-evaluation-tools-met/met-
help-desk#ftp
>
> Thanks,
> John Halley Gotway
>
> On Mon, Apr 27, 2020 at 1:22 PM efren.serra.ctr at nrlmry.navy.mil via
RT <
> met_help at ucar.edu> wrote:
>
> >
> > Mon Apr 27 13:22:52 2020: Request 95065 was acted upon.
> > Transaction: Ticket created by efren.serra.ctr at nrlmry.navy.mil
> >        Queue: met_help
> >      Subject: PJC output
> >        Owner: Nobody
> >   Requestors: efren.serra.ctr at nrlmry.navy.mil
> >       Status: new
> >  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=95065
> > >
> >
> >
> > Folks - Where are the OY_TP_i, ON_TP_i, and LIKELIHOOD_i values
for
> > probability of 1.0 included? On my PJC output files, I only see
values
> > for THRESH_n equal to 0.9. Thanks.
> >
> > Efren A. Serra (Contractor)
> > Physicist
> >
> > DeVine Consulting, Inc.
> > Naval Research Laboratory
> > Marine Meteorology Division
> > 7 Grace Hopper Ave., STOP 2
> > Monterey, CA 93943
> > Code 7542
> > Mobile: 408-425-5027
> >
> >
> >
>
>
>

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


More information about the Met_help mailing list