General joint conditional simulations using a fast Fourier transform method

A procedure for generating joint statistically homogeneous random fields is examined. The method is based on the spectral representation theorem. It handles large fields easily and is both rapid and flexible. Algorithm development and examples are presented. The procedure is adapted further to include the possibility of generating fields that are jointly conditioned on data from two related fields.