Using teleseismic data, the spatial and temporal slip distribution of the 2003 Puerto Plata earthquake in the Dominican Republic, with the recent method, was determined. We obtained the slip distribution over the fault plane, with the following source parameters: (strike, dip and rake) = (110.0º, 8.0º and 76.6 º); seismic moment of M0 = 1.446 1019 Nm (Mw6.7), source duration = 14 sec. These results obtained from the inversion analysis were used as parameters for simulating strong ground motion in the Puerto Plata area and other cities by empirical attenuation formulas for Peak Ground Acceleration (PGA) and Peak Ground Velocity (PGV) with the shortest distance to the coseismic rupture and its moment magnitude. The fault area was estimated to be 30 km x 27 km = 810km2 from the aftershock distribution obtained by a local seismic network. These PGA and PGV values are converted to the modified Mercalli intensity by using the empirical formula. The obtained values are quite consistent with an isoseismal map with the registered intensity that was felt during the Puerto Plata earthquake.