StatAnalysis to add new column FCST_THRESH and OBS_THRESH when -out_thresh is used

[pratiman-91]

Hi METplus team,

I am using StatAnalaysis for CTS (-out_line_type) with MPR (-line_type) with multiple thresholds (-out_thresh >=1.0 -out_thresh >=10.0, etc). I am able to run it, however, the final output does not create "FCST_THRESH and/or OBS_THRESH" columns. It becomes difficult to know the threshold used in the aggregated outputs.

Please let me know if you need more information.

Thanks!
-Pratiman

[j-opatz]

Hello Pratiman,

And thank you for your question. It seems like you're having difficulties getting the threshold values you are using in Stat-Analysis to appear in the output files.

By default, the Stat-Analysis tool will print out the column information to the screen; that information can also be redirected to the log file using the -out command line option. But as you found, the column header information is not included through these default behaviors.

Our suggestion would be to utilize the -out_stat setting, which will create a full stat output file (rather than a log file that contains the statistical information) for your aggregated runs. The output stat files have all of the same statistical information, but will also include the column header information like FCST_THRESH and OBS_THRESH.

Please let us know if this solves your problem, or if you need additional guidance on this topic.

[pratiman-91]

Thank you very much! -out_stat works.

[pratiman-91]

Hi @j-opatz ,

When using -out_line_type WDIR and -out_stat. The output stat file is empty, which is not expected. I can see the output is being generated but not written.

[JohnHalleyGotway]

FYI, yes, unfortunately, the WDIR output is specific to Stat-Analysis and does NOT correspond to an actual .stat line type generated by a MET statistics tool. So that is why -out_line_type WDIR and -out_stat does NOT produce any actual output.

I'll note however that the WDIR output pre-dates the addition of the vector continuous statistics line type which includes similar information. So I'd recommend trying -out_line_type VCNT -out_stat vcnt_out.stat.