An approach is described for creating random complex screens to be used in computer simulations of arbitrary Schell-model beams with a prescribed far-field intensity distribution. Simulation examples including beam profiles with reflection symmetry and rotational symmetry, flat-top, and pyramidal shapes are presented to verify the proposed approach. A more general scenario with a nonsymmetric far-field beam shape is illustrated to demonstrate the evolution in the free-space propagation from the source plane to the far zone.