Description Usage Arguments Note References
View source: R/willoughby2006eqs.R
This function uses the NewtonRaphson method to solve equation 3 (the dualexponential profile version) for ξ in Willoughby et al. (2006). This value of ξ can then be used to determine R1 for that storm observation.
1  solve_for_xi(xi0 = 0.5, eq3_right, eps = 0.001, itmax = 100)

xi0 
A numeric value giving the starting guess for ξ 
eq3_right 
A numerical value with the righthand side of Willoughby
et al. (2006), Eqn. 3, the dualexponential version. This value is
calculated at each storm observation point using the

eps 
The convergence threshold for determining if the algorithm has converged. 
itmax 
The maximum number of iterations to try before deciding that the algorithm did not converge. 
If this algorithm does not converge, the function returns a missing value for ξ.
Jones O, Maillardet R, and Robinson A. 2009. Introduction to Scientific Programming and Simulation Using R. Boca Raton, FL: Chapman & Hall/CRC Press.
Press WH, Teukolsky SA, Vetterling WT, and Flannery BP. 2002. Numerical Recipes in C++: The Art of Scientific Computing. 2nd ed. Cambridge, UK: Cambridge University Press.
Willoughby HE, Darling RWR, and Rahn ME. 2006. Parametric representation of the primary hurricane vortex. Part II: A new family of sectionally continuous profiles. Monthly Weather Review 134(4):11021120.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.