fk算法绘制层状介质理论地震图

经过了一个学期理论地震学推公式写程序的折磨,作为小白,从对Linux系统一窍不通到多文件编译运行,能正演出了第一张理论地震图还是非常开心的!

朱露培老师的fk3.2程序包下载:http://www.eas.slu.edu/People/LZhu/home.html

首先在shell中将程序解压,编译,得到的主要文件有:
fk.f 计算层状介质中爆炸源/单力/力偶的格林函数(.grn文件)
syn.c 用已知的震源机制和fk.f得到的格林函数合成理论地震图(.z文件)
trav.f 计算到时
hk 一维波速模型参数设置(hk模型)
fk.pl Perl语言写成的程序包,可以直接运行
注意:由于有makefile,编译只需make一下,(对于没有makefile的工程大家可以参考我的另一篇博客:
如何编译含有多个文件的fortran程序包(Linux多文件编译)https://blog.csdn.net/weixin_44498266/article/details/106772137

绘制理论地震图的过程:
第一步,获得格林函数.grn文件
方法1:使用perl包,用法只需要运行fk.pl即可查询
方法2:运行fk文件,输入参数参考sample_input文件

第二步,用grn文件合成理论地震图.z文件
方法:运行syn.c文件,用法只需运行syn.c即可查询

第三步,将.z文件绘制理论地震图
方法1:sac绘图
方法2:gmt里的pssac绘制

最后要感谢李老师和张老师一学期对我们的帮助,入学以来以来收获最大的课程,没有之一。po一张自己的时距曲线成图,供学弟学妹参考。
学弟学妹加油冲鸭

Logo

更多推荐