A new methodology for computing synthetic seismograms, complete of the main direct, refracted, converted phases and surface waves, in three \u2013 dimensional anelastic lateral heterogeneous media is applied to the modeling of the local records of the 2003 Bam Mw 6.6 Earthquake, South-Eastern Iran. The method is based on the combination of the Modal Summation (MS) technique with the Asymptotic Ray Theory (ART). The three \u2013 dimensional models are determined by a set of vertically heterogeneous sections (1D structures) that are juxtaposed on a regular grid. The distribution of these sections in the grid is done in such a way to satisfy the condition of applicability of the WKBJ approximation, i.e. the lateral variation of all the elastic parameters has to be small with respect to the prevailing wavelength. The procedure, described very briefly, has been applied for the area of destructive 26 December 2003 Bam earthquake Mw= 6.6