A computer package has been developed to predict the pattern of SH-wave amplitudes in the vicinity of an earthquake. The phase-front parabolic approximation, which can be applied where there are significant lateral variations in the elastic constants within the Earth, has been adapted so that it also allows for surface topography. Examples are presented from a study of shaking in the Hutt Valley and the surrounding region as a result of movement on the Wellington Fault.