On p. 2381, in the definition of the collisionality parameter $\nu$, 4 lines before Eq. 2, the factor of 2.1 should be replaced by 2.5.
On p. 2381, in the definition of R/L_n^* in the 3rd line after Eq. 2, R/L_n^* is defined as max(6,R/Ln). It should be defined as min(6,R/Ln).
Both of these are typos in the paper only. These corrections bring the paper into agreement with the subroutine that has been used to evaluate and test the model. The IFS-PPPL parameterization in this paper is fairly complicated, and there are many opportunities for typos if you try to code this up from scratch. Contact us and we will send you a copy of the fortran subroutine used to evaluate the formulas.
Clarifications/addition:
In calculating D_ML = D_mixing = max(gamma/k_perp^2) as described on the bottom of p. 2382, for practicality the range of k_theta*rho scanned was typically from approximately 0.1 to 0.8.
Though it is only used for intermediate calculations and not for the final form of chi, we define here the "W_NL" function used on page 2382. This definition is taken from two emails (the second one corrects a factor of 6.6 to be 13) from: Date: Tue, 13 May 1997 00:25:29 From: Bill Dorland bdorland@peaches.ph.utexas.edu To: Multiple recipients of list cyclone@oerhp01.er.doe.gov *************************************************************** The precise relationship between chi and D_ML in the IFS-PPPL model is the following: chi_i = W * D_ML where D_ML == max(gamma/k_perp**2) where the maximum is taken over all k_theta, and k_perp**2 is defined in our PoP paper, Vol. 2, page 2382, 1995, and where W == 13 * (1 + s_hat**2.5)/(1 + s_hat) / G(R/LT - R/LTcrit) where G(x) = 1 for x < 1, and sqrt(x) for x > 1. Here, s_hat is the usual magnetic shear parameter. The strong increase of W with s_hat is not representative of the dependence of chi on s_hat, since D ~ 1/(1+s_hat**2.5), according to our calculations. The final expression for chi_i that constitutes the IFS/PPPL is constructed from this W and a parameterization of the D_ML from Kotschenreuther's GK codes. Thus, the differences between the observed nonlinear chi from the GF ITG simulations and the mixing length estimate was parameterized by three things: a) a factor of 13 b) weaker magnetic shear dependence c) weaker dependence on R/LT far from marginal stability.