set terminal png size 600,600 set output "fig_TRS.png" set size ratio 1.0 set title "Tripartite response spectrum" set xlabel "Period (sec)" set ylabel "Velocity (kine)" set logscale x set logscale y set xrange [0.02:5] set yrange [1:500] set xtics nomirror (0.02,0.05,0.1,0.2,0.5,1,2,5) set ytics nomirror (1,2,5,10,20,50,100,200,500) # set x2label "Acceleration (gal)" set logscale x2 set x2range [0.02:5] set x2tics nomirror ("" 0.02,"50000" 0.0628,"20000" 0.157,"10000" 0.314,"5000" 0.628,"2000" 1.57,"1000" 3.14,"" 5) rotate by 45 # set y2label "Displacement (cm)" set logscale y2 set y2range [1:500] set y2tics nomirror ("" 1,"1" 1.256,"2" 2.512,"5" 6.28,"10" 12.56,"20" 25.12,"50" 62.8,"100" 125.6,"200" 251.2,"" 500) rotate by -45 # set grid xtics ytics mxtics mytics # a0=1 ;f0(x)=a0*x/(2.0*pi) a1=2 ;f1(x)=a1*x/(2.0*pi) a2=5 ;f2(x)=a2*x/(2.0*pi) a3=10 ;f3(x)=a3*x/(2.0*pi) a4=20 ;f4(x)=a4*x/(2.0*pi) a5=50 ;f5(x)=a5*x/(2.0*pi) a6=100 ;f6(x)=a6*x/(2.0*pi) a7=200 ;f7(x)=a7*x/(2.0*pi) a8=500 ;f8(x)=a8*x/(2.0*pi) a9=1000 ;f9(x)=a9*x/(2.0*pi) a10=2000 ;f10(x)=a10*x/(2.0*pi) a11=5000 ;f11(x)=a11*x/(2.0*pi) a12=10000;f12(x)=a12*x/(2.0*pi) a13=20000;f13(x)=a13*x/(2.0*pi) a14=50000;f14(x)=a14*x/(2.0*pi) # d0=0.01 ;g0(x)=2.0*pi*d0/x; d1=0.02 ;g1(x)=2.0*pi*d1/x; d2=0.05 ;g2(x)=2.0*pi*d2/x; d3=0.1 ;g3(x)=2.0*pi*d3/x; d4=0.2 ;g4(x)=2.0*pi*d4/x; d5=0.5 ;g5(x)=2.0*pi*d5/x; d6=1 ;g6(x)=2.0*pi*d6/x; d7=2 ;g7(x)=2.0*pi*d7/x; d8=5 ;g8(x)=2.0*pi*d8/x; d9=10 ;g9(x)=2.0*pi*d9/x; d10=20 ;g10(x)=2.0*pi*d10/x; d11=50 ;g11(x)=2.0*pi*d11/x; d12=100 ;g12(x)=2.0*pi*d12/x; d13=200 ;g13(x)=2.0*pi*d13/x; plot "inp_test2.txt" using 1:3 with lines linewidth 2 linetype -1 notitle, \ f1(x) with lines linewidth 1 linetype 3 notitle, \ f2(x) with lines linewidth 1 linetype 3 notitle, \ f3(x) with lines linewidth 1 linetype 3 notitle, \ f4(x) with lines linewidth 1 linetype 3 notitle, \ f5(x) with lines linewidth 1 linetype 3 notitle, \ f6(x) with lines linewidth 1 linetype 3 notitle, \ f7(x) with lines linewidth 1 linetype 3 notitle, \ f8(x) with lines linewidth 1 linetype 3 notitle, \ f9(x) with lines linewidth 1 linetype 3 notitle, \ f10(x) with lines linewidth 1 linetype 3 notitle, \ f11(x) with lines linewidth 1 linetype 3 notitle, \ f12(x) with lines linewidth 1 linetype 3 notitle, \ f13(x) with lines linewidth 1 linetype 3 notitle, \ f14(x) with lines linewidth 1 linetype 3 notitle, \ g0(x) with lines linewidth 1 linetype 1 notitle, \ g1(x) with lines linewidth 1 linetype 1 notitle, \ g2(x) with lines linewidth 1 linetype 1 notitle, \ g3(x) with lines linewidth 1 linetype 1 notitle, \ g4(x) with lines linewidth 1 linetype 1 notitle, \ g5(x) with lines linewidth 1 linetype 1 notitle, \ g6(x) with lines linewidth 1 linetype 1 notitle, \ g7(x) with lines linewidth 1 linetype 1 notitle, \ g8(x) with lines linewidth 1 linetype 1 notitle, \ g9(x) with lines linewidth 1 linetype 1 notitle, \ g10(x) with lines linewidth 1 linetype 1 notitle, \ g11(x) with lines linewidth 1 linetype 1 notitle, \ g12(x) with lines linewidth 1 linetype 1 notitle, \ g13(x) with lines linewidth 1 linetype 1 notitle