[Met_help] [rt.rap.ucar.edu #77843] History for MET Tools help

[John Halley Gotway via RT]

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

Hello Met Help,

 

I am trying to compare Stage iv precipitation with HWRF ACPC dataset. I ran
pcp_combine to get the accumulated precip for a 24 hr period and ran
grid_stat to compare HWRF output with Stage 4 precipitation dataset.

 

I have a bunch of statistics file and trying to plot the data with the R
script (on the DTC website) - "plot_mpr.R". This is the command that I used
and I get an error using this. 

 

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

[/glade/scratch/subram/Precip/grid_stat]% Rscript plot_mpr.R
grid_stat_120000L_20150616_120000V.stat

[1] "Reading: grid_stat_120000L_20150616_120000V.stat"

Error in read.table(tmp_file, header = FALSE) :

  no lines available in input

[1] "Found: 0 MPR lines"

Error in `colnames<-`(`*tmp*`, value = c("VERSION", "MODEL", "FCST_LEAD",  :

  attempt to set colnames on object with less than two dimensions

Execution halted

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

 

All my data is in Yellowstone (/glade/scratch/subram/Precip/grid_stat).
MET_TOOL executables are in /glade/scratch/subram/MET/met-5.2/bin

 

Can you please help?

 

Thank you.

 

Regards

Subashini

 

--

Subashini Subramanian

Graduate Student

Dept. of Earth, Atmospheric and Planetary Sciences

Purdue University

Ph: (765)337-1027

 

 

 



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

Subject: MET Tools help
From: John Halley Gotway
Time: Fri Sep 09 09:47:50 2016

Subashini,

I see that you have been able to run the pcp_combine and grid_stat
tools to
generate statistics verifying precipitation forecasts against StageIV
analysis.

And you are now to the point where you need to visualize those output
stats.  You tried running the "plot_mpr.R" script from the MET
website.
That script is designed to plot a very specific line type which is
produced
by the Point-Stat tool.  MPR stands for matched pair data and contains
individual paired forecast/observation values.  That script will not
be
useful to you since that's not the type of data you have.  Presumably,
you
have continuous (CNT line type) and categorical (CTS line type)
statistics
produced by grid_stat.

Now, how should you visualize the output statistics created by MET?
We do
not support a comprehensive suite of scripts for this and instead
leave it
up to the user to visualize their data.  You're welcome to use
whatever
plotting/scripting tools you find most convenient... other have used
Python, NCL, IDL, Matlab, and R.

But...

- The plot_cnt.R script (
http://www.dtcenter.org/met/users/downloads/analysis_scripts.php)
provides
an example of reading the continuous statistics (CNT) line type and
making
some plots.   However, the header it assumes is from METv3.0, so that
may
need some updating.  You're welcome to use that script and adapt it to
do
whatever you need.

- Within the DTC itself, we have developed a nice user interface for
interactively making plots with the output of MET.  It is called
METViewer:
   http://www.dtcenter.org/met/metviewer/metviewer1.jsp

I know you formerly were a DTC visitor... perhaps you could contact
whomever you worked with and ask them about loading your MET output in
METViewer.  Although using METViewer is a pretty big learning curve
and may
be overkill for what you're trying to do.

Thanks,
John Halley Gotway



On Fri, Sep 9, 2016 at 7:05 AM, Subashini Subramanian via RT <
met_help at ucar.edu> wrote:

>
> Fri Sep 09 07:05:20 2016: Request 77843 was acted upon.
> Transaction: Ticket created by subram13 at purdue.edu
>        Queue: met_help
>      Subject: MET Tools help
>        Owner: Nobody
>   Requestors: subram13 at purdue.edu
>       Status: new
>  Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=77843 >
>
>
> Hello Met Help,
>
>
>
> I am trying to compare Stage iv precipitation with HWRF ACPC
dataset. I ran
> pcp_combine to get the accumulated precip for a 24 hr period and ran
> grid_stat to compare HWRF output with Stage 4 precipitation dataset.
>
>
>
> I have a bunch of statistics file and trying to plot the data with
the R
> script (on the DTC website) - "plot_mpr.R". This is the command that
I used
> and I get an error using this.
>
>
>
> ------------------------------------------------------------
> ----------------
> ---------
>
> [/glade/scratch/subram/Precip/grid_stat]% Rscript plot_mpr.R
> grid_stat_120000L_20150616_120000V.stat
>
> [1] "Reading: grid_stat_120000L_20150616_120000V.stat"
>
> Error in read.table(tmp_file, header = FALSE) :
>
>   no lines available in input
>
> [1] "Found: 0 MPR lines"
>
> Error in `colnames<-`(`*tmp*`, value = c("VERSION", "MODEL",
"FCST_LEAD",
> :
>
>   attempt to set colnames on object with less than two dimensions
>
> Execution halted
>
> ------------------------------------------------------------
> ----------------
> ----------------
>
>
>
> All my data is in Yellowstone
(/glade/scratch/subram/Precip/grid_stat).
> MET_TOOL executables are in /glade/scratch/subram/MET/met-5.2/bin
>
>
>
> Can you please help?
>
>
>
> Thank you.
>
>
>
> Regards
>
> Subashini
>
>
>
> --
>
> Subashini Subramanian
>
> Graduate Student
>
> Dept. of Earth, Atmospheric and Planetary Sciences
>
> Purdue University
>
> Ph: (765)337-1027
>
>
>
>
>
>
>
>
>

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