基本信息

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