\HeaderA{cyclstat}{Cyclone statistics}{cyclstat}
\keyword{manip}{cyclstat}
\begin{Description}\relax
Produce storm statistics/analysis. This function was used to produce the figures in the paper Benestad, Sorteberg and Chen, 'The use of a Calculus-based Cyclone Identification method for generating storm statistics' submitted to Climate Dynamics 2004.
\end{Description}
\begin{Usage}
\begin{verbatim}
cyclstat(fname=NULL,psl0=1000,topo="data/etopo60.Rdata",cyclone=TRUE,
         x.rng=c(5,35),y.rng=c(55,72),cmp=FALSE,mon=c(12,1,2))
\end{verbatim}
\end{Usage}
\begin{Arguments}
\begin{ldescription}
\item[\code{fname}] Filename of cyclone stats data.
\item[\code{psl0}] Treshold value
\item[\code{topo}] Name of topography file.
\item[\code{cyclone}] Treshold: TRUE -> less than, FALSE -> greater than
\item[\code{x.rng}] Interval of longitudes of interest (deg E)
\item[\code{y.rng}] Interval of latitudes of interest (deg N)
\item[\code{cmp}] Flag set for comparing two storm data sets: one given a fname and the default in cyclones ('data(storms)')
\item[\code{mon}] Months to plot a histogram for.
\end{ldescription}
\end{Arguments}
\begin{Value}
\end{Value}
\begin{Author}\relax
R.E. Benestad
\end{Author}
\begin{Examples}
\begin{ExampleCode}
## Not run: 
library(cyclones)
load("data/cyclonesERA40.Rdata")
results.era40 <- results

europe <- cyclstat(results.era40,cmp=TRUE,x.rng=c(-10,20),y.rng=c(40,55))
## End(Not run)
\end{ExampleCode}
\end{Examples}

