An efficient recently developed method, the Inverse Lévy Measure (ILM) algorithm, is presented for drawing random samples from gamma, skewed stable and other nonnegative independent-increment random fields, which we call Lévy random fields. The method is useful for computing posterior distributions in nonparametric hierarchical Bayesian statistical analysis. Algorithms are illustrated through prototype implementations in S-plus.