grep

grep -a 'WTP' scc.out* |awk -F " " '{print$2" 1 "$3}'>>WTP.dat
grep -a 'STY' scc.out* |awk -F " " '{print$2" 2 "$3}'>>STY.dat
grep -a 'SGS' scc.out* |awk -F " " '{print$2" 3 "$3}'>>SGS.dat
grep -a 'CHN1' scc.out* |awk -F " " '{print$2" 4 "$3}'>>CHN1.dat
grep -a 'YULB' scc.out* |awk -F " " '{print$2" 5 "$3}'>>YULB.dat
grep -a 'ALS' scc.out* |awk -F " " '{print$2" 6 "$3}'>>ALS.dat
grep -a 'TWL' scc.out* |awk -F " " '{print$2" 7 "$3}'>>TWL.dat
grep -a 'TW.TPUB.BHE.sac.f2-8.en.s50.sh1' scc.out* |awk -F " " '{print$2" 8 "$3}'>>TPUB.dat

Response

%sac
sac> r TPUB.BHE.sac
sac> rmean
sac> rtrend
sac> transfer from polezero s s13.pz freqlimits f1 f2 f3 f4
sac> w TPUB.BHE.remove.sac

#2008.03.12.BHE.TPUB.pz(積分到位移)

ZEROS 10
 0.000000E+00     0.000000E+00
 0.000000E+00     0.000000E+00
-5.907000E+03     3.411000E+03
-5.907000E+03    -3.411000E+03
-6.839000E+02     1.755000E+02
-6.839000E+02    -1.755000E+02
-5.551000E+02     0.000000E+00
-2.946000E+02     0.000000E+00
-1.075000E+01     0.000000E+00
POLES 13
-6.909000E+03     9.208000E+03
-6.909000E+03    -9.208000E+03
-6.227000E+03     0.000000E+00
-4.936000E+03     4.713000E+03
-4.936000E+03    -4.713000E+03
-1.391000E+03     0.000000E+00
-5.568000E+02     6.005000E+01
-5.568000E+02    -6.005000E+01
-9.844000E+01     4.428000E+02
-9.844000E+01    -4.428000E+02
-1.095000E+01     0.000000E+00
-3.700000E-02     3.700000E-02
-3.700000E-02    -3.700000E-02
-2.551000E+02     0.000000E+00
CONSTANT  5.906669E+26

Jrdseed manual

Note: 


rdseed 解seed波形到sac格式
rdseed -d -o 1 -f *.seed
解出pole zero & Response
rdseed -p -R -f  *.seed


以下為手冊

INSTALLATION AND USAGE

You can download Jrdseed here

For those not familiar with the 'rdseed' interface, please refer to the rdseed online manual at:
http://www.iris.edu/manuals/rdseed.htm

The rdseed manual will tell you how to use Jrdseed, aside from the feature differences mentioned above. The other critical difference
is how you start Jrdseed. Since it uses Java, you have to use the Java launcher on the .jar file itself:
java -jar JrdseedVer0.10.1.jar (Java 1.5 or later is required.)

How do I run Jrdseed?
Jrdseed must be run from the command line.

In Windows:
Go to Start-->Run...
Type 'cmd' in the prompt. You should get a terminal window.

For other operating systems, use a Terminal window.

enter: cd 'your_Jrdseed_directory' (substitute with your directory for Jrdseed)

enter:
java -Xmx512m -jar JrdseedVer0.10.1.jar  (for a 512 MB memory allocation)
or:
java -Xmx1024m -jar JrdseedVer0.10.1.jar   (for a 1 GB memory allocation)


A sample session with Jrdseed looks like:

myserver% java -jar JrdseedVer0.10.1.jar

        << Java SEED Reader, Release 0.10.1 237 >>
        Input File (/dev/nrst0) or 'Quit' to Exit: waveforms_99.seed
        Output File (stdout) :
        Volume # [(1)-N] :
        Options [acCsSpRtde] : d
        Station List (ALL) :
        Channel List (ALL) : BH?,LH?
        Network List (ALL) :
        Loc Ids (ALL ["--" for spaces]) :
        Output Format [(1=SAC), 2=AH, 3=CSS, 4=mini seed, 5=seed, 6=sac ascii,
        7=SEGY] : 4
        Select Data Type [(E=Everything), D=Data of Undetermined State, R=Raw
        waveform Data, Q=QC'd data] : E
        Start Time(s) YYYY,DDD,HH:MM:SS.FFFF : 1999,027,06:24:59.0000
        End Time(s) YYYY,DDD,HH:MM:SS.FFFF :
        Extract Responses [Y/(N)] : N
        Input File (/dev/nrst0) or 'Quit' to Exit: q


Usage:  Jrdseed -f inputfile -{a | d [list] | l | s | t}
 where inputfile = name of input file or device;
  -a   retrieve abbreviation dictionaries
  -c   retrieve volume table of contents
  -C   retrieve the comments where arg is either STN or CHN
  -d   read data from tape
  -e   extract event/station data
  -f   input filename
  -g   alternate response file(s)
  -i   ignore location codes
  -k   strip records from output which have a zero sample count. (miniseed and seed only)
  -l   list contents of each block in volume
  -o   specify output format (1=SAC, 2=AH, 3=CSS, 4=MINISEED, 5=SEED [default=1])
  -q   specify the output directory
  -Q   specify data type (E=Everything, D=Data of undetermined state, R=Raw waveform data, Q=QC'd data
  -R   print response data
  -s   retrieve all station header information
  -S   retrieve station summary information
  -t   show start/stop times, start blocks of events
  -u   usage (this list)

foreach

#!/bin/csh

foreach fi (`ls *.HHZ`)
# fi為變數名稱，括號裡面加上想要變成變數的東西，例如這邊就是把整   個資料夾裡面帶有HHZ的檔案都變成變數
set se = `grep Sensitivity: $fi | awk '{print $3}'`
#設定se為一個新的變數，且變數內容是從grep 對應到fi這個變數的檔案中，尋找"Sensitivity"這個字串所在的行數，並print其中的第三欄而來。
echo $se >> test.txt
#叫出上面這個叫做se的變數，然後寫進去一個叫做text.txt的檔案中
echo "123"
#叫出這個自己打的字串~記得字串必須前後用" "括起來，如果是變數就不用括
end
#結束~開始跑scrip!


飯粒 -- 一次處理一大堆的sac檔(一整天的連續資料如果讓sac一次讀進去，記憶體會不足，所以必須批次處理，這時候foreach就是好幫手囉!)
 #!/bin/csh

foreach sacp (`ls *.sac`)
echo "sac << !" >>sacproce.sh
echo "r "$sacp  >>sacproce.sh
echo "interpolate delta 0.05;wh" >> sacproce.sh
echo "q" >>sacproce.sh
echo "!"  >>sacproce.sh
echo "sac << !" >>sacproce.sh
echo "r "$sacp >>sacproce.sh
echo "bp p 2 n 4 c 2 8;w append .f2-8" >>sacproce.sh
echo "r "$sacp".f2-8" >>sacproce.sh
echo "envelope;w append .en" >>sacproce.sh
echo "q" >>sacproce.sh
echo "!" >>sacproce.sh
end

飯粒--用foreach 處理從BATS下載的sac檔，更名以及ch O
 #!/bin/csh
foreach batsp (`ls *SAC`)
set YER=`echo $batsp | awk -F "." '{print$1}'`
set DAY=`echo $batsp | awk -F "." '{print$2}'`
echo "sac >>!" >> batsprocess.sh
echo "r "$batsp >> batsprocess.sh
echo "ch O GMT "$YER $DAY" 00 00 00 000">>batsprocess.sh
echo "ch allt ( sub 0.0 &1,O ) IZTYPE IO">>batsprocess.sh
echo "wh" >>batsprocess.sh
echo "q">>batsprocess.sh
echo "!" >>batsprocess.sh

echo $batsp |awk -F "." '{print "mv "$0" "$7"."$8"."$10".sac"}'>>batsprocess.sh
end