This code reads the database containing median spectra and simulates response spectra by using monte carlo simulation on the basis of the correlation coefficient between logarithmic spectral acceleration at two different periods as proposed by Jayaram and Baker
Journal Paper Baker, J., & Jayaram, N. (2008). Correlation of spectral acceleration values from NGA ground motion models. Earthquake Spectra(24(1)), 299-317