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一张自己的时距曲线成图,供学弟学妹参考。

fk算法绘制层状介质理论地震图相关推荐

  1. OpenGL——使用Bresenham算法绘制圆

    Bresenham算法是计算机图形学中为了"显示器(屏幕或打印机)系由像素构成"的这个特性而设计出来的算法,使得在求直线各点的过程中全部以整数来运算,因而大幅度提升计算速度.--摘 ...

  2. 如何用算法绘制一张上海外滩夜景图

    突发奇想:一个数学界的未解之谜 无意间想到一个有趣的数学问题,也是数学界困扰人类很久的一个谜题: 克拉茨猜想 (图片来源 [数学狂]如何形象地展示克拉茨猜想 -- Collatz Conjecture ...

  3. TIT 计算机图形学 实验三 使用重心坐标算法绘制颜色渐变的正六面体

    TIT 计算机图形学 实验三 使用重心坐标算法绘制颜色渐变的正六面体 前言 参考视频计算机图形学全套算法讲解和C++编码实现(共23讲配套源码),计算机图形学案例视频讲解以及主页相关算法.孔老师是我的 ...

  4. 基于中点算法和Bresenham算法绘制椭圆

     方法一 一.设计思路 该实验利用中点算法和Bresenham算法寻找最逼近椭圆的理想像素集. 对于,判断下一点取还是, 先取中点,判断该中点在圆内还是圆外, 若在圆外,说明椭圆上的点离下方的像素点更 ...

  5. 在Unity中使用四叉树算法绘制地形

    四叉树算法在游戏中获得了广泛的应用,前几年3D引擎实现的地形绘制大部分都是用四叉树生成的,因为移动端在硬件方面的限制,我们的地形使用的是美术自己制作的地形,对于程序来说省去了不少工作量,但是作为程序开 ...

  6. Wu反走样算法绘制圆(C++/MFC实现)

    Wu反走样圆 原理:参考Bresenham算法,在主位移过程中计算出离理想圆最近的两个点,赋予不同的亮度值,绘制像素点即可! MFC 中CXXXView类中添加函数: //Wu算法画反走样圆 void ...

  7. OpenGL实现Hermite算法绘制三次曲线

    首先是推导:节省功夫我就直接贴照片了. 程序加了鼠标的监听器,可以移动控制点和型值点. 注意:图片中矩阵第二行第二列应该为3,当时笔误 程序效果: 代码如下: #include<gl/glut. ...

  8. 利用python进行tf-idf算法绘制词云图_利用python实现通过TF-IDF和BM25提取文章关键词...

    使用TF-IDF和BM25提取文章关键词 评估方法: 人工从文章中提取1-5个关键词,和机器提取的关键词做比较 召回 = 机器提词∩人工提词 / 人工提词 准确 = 机器提词∩人工提词 / 机器提词 ...

  9. Wu反走样算法绘制直线段

    Wu反走样算法 原理:在我看来,Wu反走样算法是在Bresenham算法基础上改进了一番,它给最靠近理想直线/曲线的两个点以不同的亮度值,以达到模糊锯齿的效果.因为人眼看到的是线附近亮度的平均值. M ...

最新文章

  1. Linux用命令修改dpi,Ubuntu17.10通过dpi更改系统字体大小比例的方法
  2. 计算机中职高考,中职计算机高考中的应用
  3. 突发!贝壳董事长左晖因疾病意外去世
  4. chrome inspect 离线调试-工具包 怎么使用
  5. 网站社区类产品管理经验
  6. linux dns配置bind9,DNS服务(bind9)配置过程
  7. python批量上传 服务器_Python Tornado批量上传图片并显示功能
  8. 选择列表和可多选的选择列表
  9. 与殿堂级大师隔空对话——一张纸是怎么限制了我们的想象力
  10. react 日历组件_anujs1.5.1支持React.Suspense与lazy
  11. BT3入门之中文语言支持
  12. MongoDb学习(四)--Repository
  13. win和linux生成的密文不一样,AES加密解密Windows下跟linux下结果不同的解决方案
  14. 解决sese9 安装时多个屏幕
  15. 3dmax2017卸载/安装失败/如何彻底卸载清除干净3dmax2017注册表和文件的方法
  16. exsi rh2288hv5 驱动_华为2288H V5阵列卡驱动下载|
  17. js 实现浏览器下载视频
  18. 跨越opengl和d3d的鸿沟(四):完结篇,平台和未来
  19. 嗅探(窃听网络上流经的数据包)
  20. 怎样把两个excel表格合并

热门文章

  1. mini2440 led驱动
  2. 无Internet,安全的解决方法
  3. linux下挂载新的磁盘
  4. Swiper + 图片懒加载
  5. RabbitMQ 安装及常用操作命令
  6. SOC的定义、适用性和组成
  7. 最后完美解决pip没法用的问题
  8. 浏览器中的data类型的Url格式,data:image/png,data:image/jpeg!
  9. spring入门例子
  10. List元素自动排序