月度归档： 2025 年 2 月

import subprocess

p=subprocess.Popen(['sac'],stdin=subprocess.PIPE)
s=""
s+="echo on \n"
s+="r *.D.SAC\n"
s+="rmean\n"
s+="rtr\n"
s+="taper\n"
s+="bp co 0.5 12 n 4 p 2\n"
s+="decimate 4\n"
s+="w over\n"
s+="q\n"
p.communicate(s.encode())