[Met_help] [rt.rap.ucar.edu #76825] History for Range Check Error

Julie Prestopnik via RT met_help at ucar.edu
Tue Jul 19 10:06:36 MDT 2016 

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

John, thanks to you help,  I am progressing with HIRA and have moved on to the next piece -> temperature.  For temperature, instead of generating a Brier Skill Score, I am generating a CRPSS.  
So in this case I am writing the neighborhood data as the orank line type (each neighborhood data point is considered an ensemble member).  I sort the "members" and get the rank of the observation.  I don't know what PIT is so I am putting NA  for it.  The members are included in the orank line per the format in the user guide.  The persistence value is used for climo.

When running through MET, I get the infamous "range check error".   The command line I used is below:

/home/qcteam/met-5.1/bin/stat_analysis -lookin /h/data/global/WXQC/data/met/mdlob_pairs/TT -out /h/data/global/WXQC/data/met/summary/GALWEM_TMP_24hr_9_RHIST_0Z -job aggregate_stat -line_type ORANK -out_line_type RHIST -fcst_lead 240000 -v 6 -dump_row /h/data/global/WXQC/data/met/filter_job.stat

I placed the orank files on the ftp server.   Below is the error I get:


DEBUG 3: Processing STAT file "/h/data/global/WXQC/data/met/mdlob_pairs/TT/point_stat_9_TT_galwem_240000L_20160604_000000V.stat" ... 6 of 12

DEBUG 3: Processing STAT file "/h/data/global/WXQC/data/met/mdlob_pairs/TT/point_stat_9_TT_galwem_480000L_20160603_000000V.stat" ... 7 of 12

DEBUG 3: Processing STAT file "/h/data/global/WXQC/data/met/mdlob_pairs/TT/point_stat_9_TT_galwem_240000L_20160603_000000V.stat" ... 8 of 12

DEBUG 3: Processing STAT file "/h/data/global/WXQC/data/met/mdlob_pairs/TT/point_stat_9_TT_galwem_720000L_20160606_000000V.stat" ... 9 of 12

DEBUG 3: Processing STAT file "/h/data/global/WXQC/data/met/mdlob_pairs/TT/point_stat_9_TT_galwem_480000L_20160605_000000V.stat" ... 10 of 12

DEBUG 3: Processing STAT file "/h/data/global/WXQC/data/met/mdlob_pairs/TT/point_stat_9_TT_galwem_240000L_20160602_000000V.stat" ... 11 of 12

DEBUG 3: Processing STAT file "/h/data/global/WXQC/data/met/mdlob_pairs/TT/point_stat_9_TT_galwem_480000L_20160604_000000V.stat" ... 12 of 12

DEBUG 2: STAT Lines read     = 130649
DEBUG 2: STAT Lines retained = 32548
DEBUG 4: 
DEBUG 4: Initializing Job 1 to default job: "-job aggregate_stat -fcst_lead 240000 -line_type ORANK -dump_row /h/data/global/WXQC/data/met/filter_job.stat -out_line_type RHIST "
DEBUG 4: 
DEBUG 4: Amending Job 1 with options: "-job aggregate_stat -line_type ORANK -out_line_type RHIST -fcst_lead 240000 -dump_row /h/data/global/WXQC/data/met/filter_job.stat"
DEBUG 4: 
DEBUG 4: Amending Job 1 with command line options: "-job aggregate_stat -line_type ORANK -out_line_type RHIST -fcst_lead 240000 -dump_row /h/data/global/WXQC/data/met/filter_job.stat"
DEBUG 2: 
DEBUG 2: Processing Job 1: -job aggregate_stat -fcst_lead 240000 -line_type ORANK -dump_row /h/data/global/WXQC/data/met/filter_job.stat -out_line_type RHIST 
DEBUG 1: Creating dump row output file "/h/data/global/WXQC/data/met/filter_job.stat"
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=18446744071690283351
ERROR  : 
ERROR  : NumArray::operator[](int) const -> range check error
ERROR  :

What causes this error?

Thanks
Bob



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

Subject: Range Check Error
From: Julie Prestopnik
Time: Tue Jun 21 15:43:32 2016

Hi Bob.  Thank you for putting your files on the ftp server.

The RANK values in your input files range from 0 - 8, however the RANK
column is expected to be 1 base rank.  MET is expecting to see values
of
the range 1 - (n+1),  In this case, with a 0 base rank, MET is using 0
- 1
(or -1) to index into the array, which causes the error you are
getting.
If you change your RANK values to be 1 base rank, this problem should
go
away.

I hope that helps.

Julie

On Tue, Jun 21, 2016 at 3:03 PM, robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:

>
> Tue Jun 21 15:03:06 2016: Request 76825 was acted upon.
> Transaction: Ticket created by robert.craig.2 at us.af.mil
>        Queue: met_help
>      Subject: Range Check Error
>        Owner: Nobody
>   Requestors: robert.craig.2 at us.af.mil
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76825 >
>
>
> John, thanks to you help,  I am progressing with HIRA and have moved
on to
> the next piece -> temperature.  For temperature, instead of
generating a
> Brier Skill Score, I am generating a CRPSS.
> So in this case I am writing the neighborhood data as the orank line
type
> (each neighborhood data point is considered an ensemble member).  I
sort
> the "members" and get the rank of the observation.  I don't know
what PIT
> is so I am putting NA  for it.  The members are included in the
orank line
> per the format in the user guide.  The persistence value is used for
climo.
>
> When running through MET, I get the infamous "range check error".
The
> command line I used is below:
>
> /home/qcteam/met-5.1/bin/stat_analysis -lookin
> /h/data/global/WXQC/data/met/mdlob_pairs/TT -out
> /h/data/global/WXQC/data/met/summary/GALWEM_TMP_24hr_9_RHIST_0Z -job
> aggregate_stat -line_type ORANK -out_line_type RHIST -fcst_lead
240000 -v 6
> -dump_row /h/data/global/WXQC/data/met/filter_job.stat
>
> I placed the orank files on the ftp server.   Below is the error I
get:
>
>
> DEBUG 3: Processing STAT file
>
"/h/data/global/WXQC/data/met/mdlob_pairs/TT/point_stat_9_TT_galwem_240000L_20160604_000000V.stat"
> ... 6 of 12
>
> DEBUG 3: Processing STAT file
>
"/h/data/global/WXQC/data/met/mdlob_pairs/TT/point_stat_9_TT_galwem_480000L_20160603_000000V.stat"
> ... 7 of 12
>
> DEBUG 3: Processing STAT file
>
"/h/data/global/WXQC/data/met/mdlob_pairs/TT/point_stat_9_TT_galwem_240000L_20160603_000000V.stat"
> ... 8 of 12
>
> DEBUG 3: Processing STAT file
>
"/h/data/global/WXQC/data/met/mdlob_pairs/TT/point_stat_9_TT_galwem_720000L_20160606_000000V.stat"
> ... 9 of 12
>
> DEBUG 3: Processing STAT file
>
"/h/data/global/WXQC/data/met/mdlob_pairs/TT/point_stat_9_TT_galwem_480000L_20160605_000000V.stat"
> ... 10 of 12
>
> DEBUG 3: Processing STAT file
>
"/h/data/global/WXQC/data/met/mdlob_pairs/TT/point_stat_9_TT_galwem_240000L_20160602_000000V.stat"
> ... 11 of 12
>
> DEBUG 3: Processing STAT file
>
"/h/data/global/WXQC/data/met/mdlob_pairs/TT/point_stat_9_TT_galwem_480000L_20160604_000000V.stat"
> ... 12 of 12
>
> DEBUG 2: STAT Lines read     = 130649
> DEBUG 2: STAT Lines retained = 32548
> DEBUG 4:
> DEBUG 4: Initializing Job 1 to default job: "-job aggregate_stat
> -fcst_lead 240000 -line_type ORANK -dump_row
> /h/data/global/WXQC/data/met/filter_job.stat -out_line_type RHIST "
> DEBUG 4:
> DEBUG 4: Amending Job 1 with options: "-job aggregate_stat
-line_type
> ORANK -out_line_type RHIST -fcst_lead 240000 -dump_row
> /h/data/global/WXQC/data/met/filter_job.stat"
> DEBUG 4:
> DEBUG 4: Amending Job 1 with command line options: "-job
aggregate_stat
> -line_type ORANK -out_line_type RHIST -fcst_lead 240000 -dump_row
> /h/data/global/WXQC/data/met/filter_job.stat"
> DEBUG 2:
> DEBUG 2: Processing Job 1: -job aggregate_stat -fcst_lead 240000
> -line_type ORANK -dump_row
/h/data/global/WXQC/data/met/filter_job.stat
> -out_line_type RHIST
> DEBUG 1: Creating dump row output file
> "/h/data/global/WXQC/data/met/filter_job.stat"
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=18446744071690283351
> ERROR  :
> ERROR  : NumArray::operator[](int) const -> range check error
> ERROR  :
>
> What causes this error?
>
> Thanks
> Bob
>
>
>


--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu

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

More information about the Met_help mailing list