[Met_help] [rt.rap.ucar.edu #94140] History for Threshold NA in Stat Anal Output
John Halley Gotway via RT
met_help at ucar.edu
Wed Feb 12 11:15:18 MST 2020
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
John, when I run the stat_anal job below, the out_thresh field is set to NA
instead of the threshold I have selected. How can I get this field to
display my threshold?
/h/MET/bin/stat_analysis -lookin
/h/data/global/WXQC/data/met/ptstat/stats/compare_pairs/WS/point_stat_galwem
10_WS_SFC_1320000L_20191228_060000V.stat -out
/h/data/global/WXQC/data/met/summary/WS_gt15_cts_galwem10_conus_comp -model
galwem10 -fcst_var WIND -vx_mask FULL -mask_poly
/h/data/global/WXQC/data/met/poly/CONUS.poly -tmp_dir
/h/data/global/WXQC/data/met/tmp -job aggregate_stat -line_type MPR
-out_line_type CTS -out_thresh gt15 -by OUT_THRESH -fcst_valid_beg
20191217_000000 -fcst_valid_end 20191231_180000 -fcst_lev Z10 -by FCST_LEAD
-config
/h/data/global/WXQC/data/met/config/met_config/STATAnalysisConfig_bootstrap_
aggstat -v 4.0
Thanks
Bob
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Threshold NA in Stat Anal Output
From: John Halley Gotway
Time: Tue Feb 11 13:30:29 2020
Bob,
Short answer... "-by OUT_THRESH" doesn't work and please start using
the
"-out_stat" option.
Longer answer...
Thanks for sending you sample data file and listing the command you're
running. I used it to run the following command:
stat_analysis \
-lookin point_stat_galwem10_WS_SFC_1320000L_20191228_060000V.stat \
-model galwem10 -fcst_var WIND -fcst_lev Z10 -vx_mask FULL \
-fcst_valid_beg 20191217_000000 -fcst_valid_end 20191231_180000 \
-job aggregate_stat -line_type MPR -out_line_type CTS -out_thresh gt15
\
-mask_poly CONUS.poly -by FCST_LEAD -out_stat agg_MPR_to_CTS.stat
Let me make a few comments about this:
(1) In the command you sent, "-by OUT_THRESH" has not effect. The
arguments for the "-by" option are the column names from the input
.stat
files being read. There is no column named "OUT_THRESH". The "-by"
option
is used to subset the input data into cases... and the job specified
is run
for each of those cases. So "-by FCST_LEAD -by FCST_LEV" would work
well
(assuming you had different lead times and vertical levels in the
.stat
files you're passing to STAT-Analysis).
(2) Your specific question is why doesn't ">15" show up in the output?
The
original behavior of stat-analysis was to print the job output to the
screen, and it does still do that. In one of the releases along the
way,
we added the "-out_stat" option which I'm using in the job listed
above.
That writes a full .stat file, including all the header columns. The
">15"
threshold does indeed appear in the FCST_THRESH and OBS_THRESH header
columns in that file.
See attached.
Hope that helps.
Thanks,
John
On Tue, Feb 11, 2020 at 12:44 PM Randy Bullock via RT
<met_help at ucar.edu>
wrote:
>
> Tue Feb 11 12:43:38 2020: Request 94140 was acted upon.
> Transaction: Given to johnhg (John Halley Gotway) by bullock
> Queue: met_help
> Subject: Threshold NA in Stat Anal Output
> Owner: johnhg
> Requestors: robert.craig.2 at us.af.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94140 >
>
>
> This transaction appears to have no content
>
------------------------------------------------
Subject: RE: [Non-DoD Source] Re: [rt.rap.ucar.edu #94140] Threshold NA in Stat Anal Output
From: robert.craig.2 at us.af.mil
Time: Tue Feb 11 14:17:55 2020
Okay, that worked for me.
Thanks
-----Original Message-----
From: John Halley Gotway via RT <met_help at ucar.edu>
Sent: Tuesday, February 11, 2020 2:30 PM
To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #94140] Threshold NA in
Stat Anal Output
Bob,
Short answer... "-by OUT_THRESH" doesn't work and please start using
the "-out_stat" option.
Longer answer...
Thanks for sending you sample data file and listing the command you're
running. I used it to run the following command:
stat_analysis \
-lookin point_stat_galwem10_WS_SFC_1320000L_20191228_060000V.stat \
-model galwem10 -fcst_var WIND -fcst_lev Z10 -vx_mask FULL \
-fcst_valid_beg 20191217_000000 -fcst_valid_end 20191231_180000 \ -job
aggregate_stat -line_type MPR -out_line_type CTS -out_thresh gt15 \
-mask_poly CONUS.poly -by FCST_LEAD -out_stat agg_MPR_to_CTS.stat
Let me make a few comments about this:
(1) In the command you sent, "-by OUT_THRESH" has not effect. The
arguments for the "-by" option are the column names from the input
.stat files being read. There is no column named "OUT_THRESH". The
"-by" option is used to subset the input data into cases... and the
job specified is run for each of those cases. So "-by FCST_LEAD -by
FCST_LEV" would work well (assuming you had different lead times and
vertical levels in the .stat files you're passing to STAT-Analysis).
(2) Your specific question is why doesn't ">15" show up in the output?
The original behavior of stat-analysis was to print the job output to
the screen, and it does still do that. In one of the releases along
the way, we added the "-out_stat" option which I'm using in the job
listed above.
That writes a full .stat file, including all the header columns. The
">15"
threshold does indeed appear in the FCST_THRESH and OBS_THRESH header
columns in that file.
See attached.
Hope that helps.
Thanks,
John
On Tue, Feb 11, 2020 at 12:44 PM Randy Bullock via RT
<met_help at ucar.edu>
wrote:
>
> Tue Feb 11 12:43:38 2020: Request 94140 was acted upon.
> Transaction: Given to johnhg (John Halley Gotway) by bullock
> Queue: met_help
> Subject: Threshold NA in Stat Anal Output
> Owner: johnhg
> Requestors: robert.craig.2 at us.af.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94140
> >
>
>
> This transaction appears to have no content
>
------------------------------------------------
More information about the Met_help
mailing list