新浪博客

全球岩石圈模型Lith1.0

2014-08-07 11:38阅读：

http://blog.sina.cn/dpool/blog/u/3887657512

基本信息
LITHO 1.0 v1是一个2014年公布的全球岩石圈模型，横向为1°的完全嵌合节点，垂向上分为一系列地球物理层，如水，冰，沉积层，结晶地壳，岩石圈盖层，软流圈。该模型是以crust1.0模型为初始模型，拟合全球5-40mHz面波数据，将crust1.0模型扩展到上地幔岩石圈和下部的软流圈。

模型分层
每个节点模型可以分为11层，分别为
1. 水层；
2. 冰层；
3. 上沉积层；
4. 中沉积层；
5. 下沉积层；
6.
上地壳；
7. 中地壳；
8. 下地壳；
9. 岩石圈地幔
10.软流圈
11.ak135模型
官方下载
http://igppweb.ucsd.edu/~gabi/litho1.0.html
内容
内含说明文档README一份；
内含模型文件40962个；
内含Cpp代码1个，用于读取模型，编译文件makefile一个；
如何编译
1解压后在LITHO1.0中新建一个bin文件夹
2进入progs文件夹
3把Makefile中相关变量修改为自己系统的文件夹
BIN = /home/gecan/FedoraShare/opt/LITHO1.0/bin
CPP = g++
FLAGS = -DMODELLOC=\'/home/gecan/FedoraShare/opt/LITHO1.0/litho_model\'
4把access_litho.cc源文件中
#include string 修改为#include string.h (尖括号打不出来)
避免调用strcpy出现strcpy' was not declared in this scope这样的错误。
5运行 make access_litho进行编译
6将新程序加入系统路径中
[gecan@localhost bin]$ vi ~/.bashrc
在文件末尾加入
#litho1.0
export PATH=$PATH:/home/gecan/FedoraShare/opt/LITHO1.0/bin
保存并激活
[gecan@localhost bin]$ source ~/.bashrc
如何使用
$access_litho -h
access_litho -p lat lon [ -d depth] [-l level] [-e] [-h]
-h help
-p lat lon (runs in profile mode)
-d depth (runs in point mode)
-l level

$ access_litho -p 12. 34.
60781. 3300.00 7874.65 4279.70 0.00 70.00 7874.65 4279.70 1.00000 ASTHENO-TOP
60781. 3300.00 8035.35 4578.55 0.00 200.00 8035.35 4578.55 1.00000 LID-BOTTOM
43521. 3300.00 8035.35 4578.55 0.00 200.00 8035.35 4578.55 1.00000 LID-TOP
43521. 3034.12 7180.76 4096.07 0.00 600.00 7180.76 4096.07 1.00000 CRUST3-BOTTOM
29011. 3034.12 7180.76 4096.07 0.00 600.00 7180.76 4096.07 1.00000 CRUST3-TOP
29011. 2852.08 6573.94 3782.54 0.00 600.00 6573.94 3782.54 1.00000 CRUST2-BOTTOM
14490. 2852.08 6573.94 3782.54 0.00 600.00 6573.94 3782.54 1.00000 CRUST2-TOP
14490. 2750.94 6270.52 3580.27 0.00 600.00 6270.52 3580.27 1.00000 CRUST1-BOTTOM
-471. 2750.94 6270.52 3580.27 0.00 600.00 6270.52 3580.27 1.00000 CRUST1-TOP
-471. 2110.00 2500.00 1070.00 0.00 600.00 2500.00 1070.00 1.00000 SEDS1-BOTTOM
-479. 2110.00 2500.00 1070.00 0.00 600.00 2500.00 1070.00 1.00000 SEDS1-TOP

The output is:
depth(m) density(kg/m3) Vp(m/s) Vs(m/s) Qkappa Qmu Vp2(m/s) Vs2(m/s) eta layername
参考文献
Pasyanos, M.E., T.G. Masters, G. Laske, and Z. Ma (2014). LITHO1.0: An updated crust and lithospheric model of the Earth, J. Geophys. Res., DOI: 10.1002/2013JB010626.
http://onlinelibrary.wiley.com/doi/10.1002/2013JB010626/full

我的更多文章

下载客户端阅读体验更佳

APP专享