2025-01-01实际波形和理论波形对比图画法-GMT版

（a）按-En画在一起的；

（b）(c) 按-En分开画的；

（d）按方位角画在一起的；

重点就是把波形名字、其他信息放在一个文件里面，调用文件画图

脚本如下：注意区分实际波形和理论波形的名称
   #!bin/bash

   # cp ./pick_obstimeT/*.t .

   # cp ./syn/*.t .

   # printf “r **.t
rtr;rmean;taper
bp c 0.02 0.08 n 4 p 2
w over
quit
” | sac

   ls ??.???.00.t | awk {split($1,a,”.”);print 160,NR,$1,a[2],a[2]”_101.8200000_32.2500000.t”} > stalist.txt

   saclst dist f ??.???.00.t | sort -gk2 | awk {split($1,a,”.”);printNR,$2,$1,a[2],a[2]”_101.8200000_32.2500000.t”} > stalist_az.txt

   ##################

   gmt begin EKall jpg

   cat stalist.txt | awk {print $3} |gmt sac -BWS -Bx+l”time/s” -By+l”aligned by number” -JX6c/20c -R-10/180/0/26 -Baf -En1 -M0.1i -W1p,black

   cat stalist.txt | awk {print $5} |gmt sac -JX6c/20c -R-10/180/0/26 -Baf -En0.8 -M0.1i -W1p,red

   cat stalist.txt | awk {print $1,$2-0.2,$4} | gmt text -JX6c/20c -R-10/180/0/26 -Baf -F+f10p

   echo 0 25.5 (a) | gmt text -JX6c/20c -R-10/180/0/26 -Baf -F+f10p

   ##################

   cat stalist.txt | awk {print $3} |gmt sac -BWSen -Bx+l”obs time/s” -By+l”aligned by number” -JX6c/20c -R-10/180/0/26 -En1 -M0.1i -W1p,black -X4i

   cat stalist.txt | awk {print $1,$2-0.2,$4} | gmt text -JX6c/20c -R-10/180/0/26 -Baf -F+f10p

   echo 0 25.5 (b) | gmt text -JX6c/20c -R-10/180/0/26 -Baf -F+f10p

   ###################

   cat stalist.txt | awk {print $5} |gmt sac -BWSen -Bx+l”osyn time/s” -By+l”aligned by number” -JX6c/20c -R-10/180/0/26 -En1 -M0.1i -W1p,red -X3.2i

   cat stalist.txt | awk {print $1,$2-0.2,$4} | gmt text -JX6c/20c -R-10/180/0/26 -Baf -F+f10p

   echo 0 25.5 (c) | gmt text -JX6c/20c -R-10/180/0/26 -Baf -F+f10p

   ##################

   gmt basemap -JX6c/20c -R-10/180/0/26 -B+n -X4i

   cat stalist_az.txt | awk {print $3,-10,$2} |gmt sac -JX6c/20c -R-10/180/0/370 -BWS -Bx+l”time/s” -By+l”azimuth/degree” -M0.1i -W1p,black

   cat stalist_az.txt | awk {print $5,-10,$2-2} |gmt sac -JX6c/20c -R-10/180/0/370  -M0.1i -W1p,red

   cat stalist_az.txt | awk {print 160,$2-0.2,$4} | gmt text -JX6c/20c -R-10/180/0/370 -F+f10p

   echo 0 365 (d) | gmt text -JX6c/20c -R-10/180/0/26 -Baf -F+f10p