[Met_help] [rt.rap.ucar.edu #88776] History for Range Check Error in Stat Analysis

John Halley Gotway via RT met_help at ucar.edu
Wed Feb 6 10:19:31 MST 2019


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

John, I am running stat analysis on some point stat (not hira) files and
looking to get CTS results for precip type.  Precip type ranges from 0-7 in
my model and observation files.  Then I run the code I get:

 

DEBUG 3: Processing STAT file
"/h/data/global/WXQC/data/met/ptstat/mdlob_pairs/galwemf/PT/point_stat_galwe
mf_all_150000L_20190115_150000V.stat" ... 81 of 84

DEBUG 3: Processing STAT file
"/h/data/global/WXQC/data/met/ptstat/mdlob_pairs/galwemf/PT/point_stat_galwe
mf_all_270000L_20190115_030000V.stat" ... 82 of 84

DEBUG 3: Processing STAT file
"/h/data/global/WXQC/data/met/ptstat/mdlob_pairs/galwemf/PT/point_stat_galwe
mf_all_210000L_20190112_210000V.stat" ... 83 of 84

DEBUG 3: Processing STAT file
"/h/data/global/WXQC/data/met/ptstat/mdlob_pairs/galwemf/PT/point_stat_galwe
mf_all_330000L_20190114_090000V.stat" ... 84 of 84

DEBUG 2: STAT Lines read     = 186613

DEBUG 2: STAT Lines retained = 2298

DEBUG 4: 

DEBUG 4: Initializing Job 1 to default job: "-job aggregate_stat
-fcst_valid_beg 20190112_000000 -fcst_valid_end 20190112_000000 -line_type
MPR -by DESC -by FCST_THRESH -by FCST_LEAD -dump_row
/h/data/global/WXQC/data/met/filter_job -out_line_type CTS -out_alpha
0.05000 -boot_interval 1 -boot_rep_prop 1.00000 -n_boot_rep 1000 -boot_rng
mt19937 -boot_seed '' "

DEBUG 4: 

DEBUG 4: Amending Job 1 with options: "-job aggregate_stat -line_type MPR
-out_line_type CTS -fcst_valid_beg 20190112_000000 -fcst_valid_end
20190112_000000 -by desc -by FCST_THRESH -by FCST_LEAD -dump_row
/h/data/global/WXQC/data/met/filter_job"

DEBUG 4: 

DEBUG 4: Amending Job 1 with command line options: "-job aggregate_stat
-line_type MPR -out_line_type CTS -fcst_valid_beg 20190112_000000
-fcst_valid_end 20190112_000000 -by desc -by FCST_THRESH -by FCST_LEAD
-dump_row /h/data/global/WXQC/data/met/filter_job"

DEBUG 2: 

DEBUG 2: Processing Job 1: -job aggregate_stat -fcst_valid_beg
20190112_000000 -fcst_valid_end 20190112_000000 -line_type MPR -by DESC -by
FCST_THRESH -by FCST_LEAD -dump_row /h/data/global/WXQC/data/met/filter_job
-out_line_type CTS -out_alpha 0.05000 -boot_interval 1 -boot_rep_prop
1.00000 -n_boot_rep 1000 -boot_rng mt19937 -boot_seed '' 

DEBUG 1: Creating dump row output file
"/h/data/global/WXQC/data/met/filter_job"

GSL_RNG_TYPE=mt19937

GSL_RNG_SEED=18446744072848627804

DEBUG 2: Computing output for 1 case(s).

ERROR  : 

ERROR  : ThreshArray::operator[](int) const -> range check error!

ERROR  :

 

I know we have seen this before but don't remember what triggers it.  It
would be a  nice addition to the users manual to explain possible causes for
errors like this.  Attached is the MPR file I am reading.

 

Thanks

Bob

 

 



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

Subject: FW: [Non-DoD Source] [rt.rap.ucar.edu #88776] AutoReply: Range Check Error in Stat Analysis
From: robert.craig.2 at us.af.mil
Time: Tue Feb 05 14:42:00 2019

John, found a work around.

Bob

-----Original Message-----
From: met_help at ucar.edu via RT <met_help at ucar.edu>
Sent: Monday, February 4, 2019 3:56 PM
To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
Subject: [Non-DoD Source] [rt.rap.ucar.edu #88776] AutoReply: Range
Check Error in Stat Analysis

Greetings,

This message has been automatically generated in response to the
creation of a trouble ticket regarding:
"Range Check Error in Stat Analysis",
a summary of which appears below.

There is no need to reply to this message right now.  Your ticket has
been assigned an ID of [rt.rap.ucar.edu #88776].

Please include the string:

         [rt.rap.ucar.edu #88776]

in the subject line of all future correspondence about this issue. To
do so, you may reply to this message.

For more information, please see:

MET Online Tutorial:
   https://www.dtcenter.org/met/users/support/online_tutorial/index.php

MET Users Guide:
   https://www.dtcenter.org/met/users/docs/overview.php

MET FAQs:
   https://www.dtcenter.org/met/users/support/faqs/index.php

MET-Help Email Archive:
   http://mailman.ucar.edu/pipermail/met_help

                        Thank you,
                        met_help at ucar.edu

-------------------------------------------------------------------------
John, I am running stat analysis on some point stat (not hira) files
and looking to get CTS results for precip type.  Precip type ranges
from 0-7 in my model and observation files.  Then I run the code I
get:



DEBUG 3: Processing STAT file
"/h/data/global/WXQC/data/met/ptstat/mdlob_pairs/galwemf/PT/point_stat_galwe
mf_all_150000L_20190115_150000V.stat" ... 81 of 84

DEBUG 3: Processing STAT file
"/h/data/global/WXQC/data/met/ptstat/mdlob_pairs/galwemf/PT/point_stat_galwe
mf_all_270000L_20190115_030000V.stat" ... 82 of 84

DEBUG 3: Processing STAT file
"/h/data/global/WXQC/data/met/ptstat/mdlob_pairs/galwemf/PT/point_stat_galwe
mf_all_210000L_20190112_210000V.stat" ... 83 of 84

DEBUG 3: Processing STAT file
"/h/data/global/WXQC/data/met/ptstat/mdlob_pairs/galwemf/PT/point_stat_galwe
mf_all_330000L_20190114_090000V.stat" ... 84 of 84

DEBUG 2: STAT Lines read     = 186613

DEBUG 2: STAT Lines retained = 2298

DEBUG 4:

DEBUG 4: Initializing Job 1 to default job: "-job aggregate_stat
-fcst_valid_beg 20190112_000000 -fcst_valid_end 20190112_000000
-line_type MPR -by DESC -by FCST_THRESH -by FCST_LEAD -dump_row
/h/data/global/WXQC/data/met/filter_job -out_line_type CTS -out_alpha
0.05000 -boot_interval 1 -boot_rep_prop 1.00000 -n_boot_rep 1000
-boot_rng
mt19937 -boot_seed '' "

DEBUG 4:

DEBUG 4: Amending Job 1 with options: "-job aggregate_stat -line_type
MPR -out_line_type CTS -fcst_valid_beg 20190112_000000 -fcst_valid_end
20190112_000000 -by desc -by FCST_THRESH -by FCST_LEAD -dump_row
/h/data/global/WXQC/data/met/filter_job"

DEBUG 4:

DEBUG 4: Amending Job 1 with command line options: "-job
aggregate_stat -line_type MPR -out_line_type CTS -fcst_valid_beg
20190112_000000 -fcst_valid_end 20190112_000000 -by desc -by
FCST_THRESH -by FCST_LEAD -dump_row
/h/data/global/WXQC/data/met/filter_job"

DEBUG 2:

DEBUG 2: Processing Job 1: -job aggregate_stat -fcst_valid_beg
20190112_000000 -fcst_valid_end 20190112_000000 -line_type MPR -by
DESC -by FCST_THRESH -by FCST_LEAD -dump_row
/h/data/global/WXQC/data/met/filter_job
-out_line_type CTS -out_alpha 0.05000 -boot_interval 1 -boot_rep_prop
1.00000 -n_boot_rep 1000 -boot_rng mt19937 -boot_seed ''

DEBUG 1: Creating dump row output file
"/h/data/global/WXQC/data/met/filter_job"

GSL_RNG_TYPE=mt19937

GSL_RNG_SEED=18446744072848627804

DEBUG 2: Computing output for 1 case(s).

ERROR  :

ERROR  : ThreshArray::operator[](int) const -> range check error!

ERROR  :



I know we have seen this before but don't remember what triggers it.
It would be a  nice addition to the users manual to explain possible
causes for errors like this.  Attached is the MPR file I am reading.



Thanks

Bob








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


More information about the Met_help mailing list