[Met_help] [rt.rap.ucar.edu #98772] History for How to use stat_analysis

[John Halley Gotway via RT]

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

Hello,

I have a series of .stat files from a point_stat run containing CNT and MPR  line_types for airtmp and geopht for taus ranging from 0 to 48h at 12 h increments.  These lines are valid at a single dtg (i.e., a stat file for 2019072500 will contain 7/23 0Z tau48, 7/23 12Z tau36, 7/24 0Z tau24, etc.).  I am trying to get the ME and RMSE for each of these taus over all stat files (so all taus from 2019072500 to 2019080300) and am currently executing:

stat_analysis -lookin ./ -job aggregate_stat -line_type MPR -out_line_type CNT -fcst_var airtmp -obs_var t -by airtmp,500,24 -out_stat 24h_out.stat

This is being done for 500 mb airtemp, tau 24.  The resulting file, 24h_out.stat, shows my FCST_LEAD being 48, even though I specified 24h.  Any tips on how to complete this?

Thanks,
Justin


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

Subject: How to use stat_analysis
From: John Halley Gotway
Time: Fri Feb 19 16:28:18 2021

Justin,

I think you’re really close. I notice a small problem in your job
command:
 -by airtmp,500,24

The -by option should specify the NAME of the columns which define the
case
information and NOT specific values within those columns.

Please try:
  -by FCST_VAR,FCST_LEV,FCST_LEAD

Instead

Thanks
John HG

On Fri, Feb 19, 2021 at 4:08 PM Tsu, Mr. Justin via RT
<met_help at ucar.edu>
wrote:

>
> Fri Feb 19 16:08:16 2021: Request 98772 was acted upon.
> Transaction: Ticket created by justin.tsu at nrlmry.navy.mil
>        Queue: met_help
>      Subject: How to use stat_analysis
>        Owner: Nobody
>   Requestors: justin.tsu at nrlmry.navy.mil
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=98772 >
>
>
> Hello,
>
> I have a series of .stat files from a point_stat run containing CNT
and
> MPR  line_types for airtmp and geopht for taus ranging from 0 to 48h
at 12
> h increments.  These lines are valid at a single dtg (i.e., a stat
file for
> 2019072500 will contain 7/23 0Z tau48, 7/23 12Z tau36, 7/24 0Z
tau24,
> etc.).  I am trying to get the ME and RMSE for each of these taus
over all
> stat files (so all taus from 2019072500 to 2019080300) and am
currently
> executing:
>
> stat_analysis -lookin ./ -job aggregate_stat -line_type MPR
-out_line_type
> CNT -fcst_var airtmp -obs_var t -by airtmp,500,24 -out_stat
24h_out.stat
>
> This is being done for 500 mb airtemp, tau 24.  The resulting file,
> 24h_out.stat, shows my FCST_LEAD being 48, even though I specified
24h.
> Any tips on how to complete this?
>
> Thanks,
> Justin
>
>

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