GSAC

#!/bin/sh

infile=sac2/2012.083.09.25.52.0195.II.WRAB.00.BHZ.M.SAC psfile=fig/wrab_2012

#infile=sac2/2013.160.14.22.47.0195.II.WRAB.00.BHZ.M.SAC #psfile=fig/wrab_2013

gsac << EOF bg plt xlim 40 320 r $infile pctl xlen 4.0 ylen 2.0 grid on p xlim off taper fft plotsp cut 60 120 r $infile p taper fft plotsp cut 140 160 r $infile p taper fft plotsp cut 190 240 r $infile p taper fft plotsp quit EOF

for file in 001 002 003 004 005 006 007 008 do echo P$file.PLT $psfile.$file.ps plotnps -F5 -S0.5 -A4 -R -W3 -H30 -KB < P$file.PLT > $psfile.$file.ps done

convert ps file to jpeg

for file in ls fig/*ps do ps2raster -A -P $file done