# coding: utf-8
import numpy as np
import matplotlib.pyplot as pl

data1 = np.loadtxt('dat_spc_to_EW2.csv',delimiter=',', skiprows=3, usecols=(0,1))
data2 = np.loadtxt('dat_spc_to_NS2.csv',delimiter=',', skiprows=3, usecols=(0,1))
data3 = np.loadtxt('dat_spc_to_UD2.csv',delimiter=',', skiprows=3, usecols=(0,1))

fig=pl.figure(num=1,figsize=(6, 6), dpi=100)
pl.subplot(1,1,1,aspect=1.0)
pl.xlabel('Period (sec)',fontsize=12)
pl.ylabel('Response acceleration (gal)',fontsize=12)
pl.xlim(0.01, 10)
pl.ylim(10.0, 10000)
pl.loglog(data1[:,0],data1[:,1],color='red'  ,linewidth=2.0,linestyle='-',label='EW2')
pl.loglog(data2[:,0],data2[:,1],color='blue' ,linewidth=2.0,linestyle='-',label='NS2')
pl.loglog(data3[:,0],data3[:,1],color='green',linewidth=2.0,linestyle='-',label='UD2')
pl.legend(shadow=True, prop={'size':12},loc='upper right')
pl.grid(True,which='both',ls='-', color='0.65')
#pl.show()
pl.savefig("fig_spc.png", dpi=100)