The Gaussian distribution is a stable distribution, and so, a Gaussian random field with any desired energy spectrum may be generated quite easily using Fourier-space filtering. A log-normal simulation may then be generated by exponentiating this Gaussian random field. However, the difficulty arises because a log-normal field generated in this manner will not have the same energy spectrum as the original Gaussian field. Ideally, it would be possible to analytically determine the energy spectrum for the Gaussian field, such that exponentiation of this field results in a log-normal field with the desired energy spectrum. This approach is implemented in the method of Evans et al. (1999), which is valid for the case beta<1.

We extend the range of validity to 1 < beta < 3 by using a simple iterative method to numerically determine the appropriate energy spectrum for the Gaussian field. The iteration produces a sequence of Gaussian random fields, where the difference between the desired energy spectrum (for the log-normal field) and the energy spectrum of the exponentiation of a field in the sequence is used to generate the energy spectrum for the next Gaussian field in the sequence (i.e. the next iterate). The main assumption of the method is that the exponentiation of the Gaussian field does not substantially alter the form of the energy spectrum. This is a valid assumption when the mean and the variance of the Gaussian field are not large, which is the case of interest for cloud statistics. We present results from this technique and discuss its performance, including the variance range for which convergence occurs.

Supplementary URL: http://www.eos.ubc.ca/research/clouds/radjp416