







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())
小家伙已经完全接受去上幼儿园啦。
先准备多段数据文件,用> Z作为数据头,后面跟随具体数值,如:
> -Z0.500000
如果使用GMT6.0绘图,代码可以写为:
gmt plot prov_poly.dat -Ccolor.cpt -G+z
上面一行是用来填写的,如果不带-G,则默认描边。
有网页中提出的codesign命令并没有解决这个问题
直接在终端下进入/Applications/MATLAB2023a/Contents/MacOS/
./MATLAB
一定要终端进入,鼠标点击不行。
主要是units出错,经查