gCAP

1. 下载

   # gCAP1.0包
   1 $ wget http://www.eas.slu.edu/People/LZhu/downloads/gcap1.0.tar
   # 所需NR函数包	
   2 $ wget https://blog.seisman.info/downloads/gcap_utils.tar.gz
   # 编译好的Makefile
   3 $ wget https://blog.seisman.info/downloads/Makefile.gCAP
   

2. 解压安装

   1 $ tar -xvf gcap1.0.tar
   2 $ tar -xvf gcap_utils.tar.gz --strip-components 1 -C /home/wangjq/software/gcap
   3 $ mv Makefile.gCAP /home/wangjq/software/gcap/Makefile
   4 $ mv gcap /home/wangjq/opt/gcap
   5 $ cd ~/opt/gcap
   

3. 修改Makefile

   1 $ vi Makefile
   # 修改SAC路径
   	SACHOME = /home/wangjq/opt/sac
   

4. 编译

   1 $ make
   

5. 添加环境变量

   1 $ sudo vi ~/.bashrc
   	export PATH=/home/wangjq/opt/gcap:${PATH}
   2 $ source ~/.bashrc
   

6. 修改cap.pl/cap_plt.pl

   1 $ vi cap.pl
   	# 修改如下语句
   	require "/home/wangjq/opt/gcap/cap_plt.pl"
   	$green = "/home/wangjq/data/models/Glib"
   2 $ vi cap_plt.pl
   	# 在 cap_plt.pl 的第 6 行之后加上如下语句
   	system "gmtset MEASURE_UNIT inch";
   	system "gmtset PAGE_ORIENTATION portrait";
   

7. 执行

   1 $ cap.pl
   ===== CAP seismic source tensor inversion using seismic waveforms ====
       Ref: Zhu and Helmberger, 1996, BSSA 86, 1645-1641.
            Zhu and Ben-Zion, 2013, GJI, submitted.
   

8. 用 fk 构建格林函数库

   1 # 切换到格林函数库所在目录
   2 $ cd ~/data/models/Glib
   3 # 为 hk 模型新建目录
   4 $ mkdir hk
   5 $ cd hk
   6 # 将 fk 自带的 hk 模型文件复制到当前目录下
   7 $ cp ~/opt/fk/hk ~/data/models/Glib/hk/hk
   8 # 用 fk 计算格林函数库
   9 $ fk.pl -Mhk/15/k -N512/0.2 -S2 05 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85
   90 95 100 105 110 115 120 125 130 135 140 145 150 155 160 165 170 175 180 185
   190 195 200 205 210 215 220 225 230 235 240 245 250 255 260 265 270 275 280
   285 290 295 300 305 310 315 320 325 330 335 340 345 350 355 360 365 370 375
   380 385 390 395 400 405 410 415
   10 $ fk.pl -Mhk/15/k -N512/0.2 -S0 05 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85
   90 95 100 105 110 115 120 125 130 135 140 145 150 155 160 165 170 175 180 185
   190 195 200 205 210 215 220 225 230 235 240 245 250 255 260 265 270 275 280
   285 290 295 300 305 310 315 320 325 330 335 340 345 350 355 360 365 370 375
   380 385 390 395 400 405 410 415
   # -M 选项中 15 代表震源深度是 15 km，你可以从 1 km 一直算到 40 km 甚至更深，此处只算了 15 km 以供示例使用。
   

9. 运行 gcap 示例

   1 # 切换回 gcap 目录
   2 $ cd ~/opt/gcap
   3 # 为了防止误操作，先将自带的示例备份一下
   4 $ cp -r 20080418093700 20080418093700.bak
   5 # 运行 gcap 示例
   6 $ cap.pl -H0.2 -P0.3 -S2/5/0 -T35/70 -F -D2/1/0.5 -C0.05/0.3/0.02/0.1 -W1 -X10 -Mhk_15/5.0 20080418093700