最近看孙文科教授的相关文章。要获得位错产生的地表形变，要完成以下任务：

1、计算表面自由边界边值问题－与计算固体潮汐一致的微分方程组和边界条件

2、计算表面负荷边界边值问题－与计算负荷潮汐一致的微分方程组和边界条件

3、计算表面剪切边界边值问题－这个，，，这个？？？，反正就是边界条件不一样而方程一样。（通过上述三种边界的解的组合可以获得任意边界的解）

4、利用okubo的等效原理（俺的翻译，呵呵reciprocity theorem），内部位错在地表产生的变形等效于上述三个结果在位错深度处的线性组合。其中负荷边界和剪切边界由位错参数确定。

5、计算4种独立位错类型产生的变形的计算（任何位错都是这四种的线性组合）

6、利用坐标转换将北极处的位错（方位角为0）转到实际位错发生的地点（利用经纬度，位错方位角），计算地表变形(变形包括位移倾斜重力应变等)

UdS=2.49e+13m^3  ,dip angle=60degree, rake=30degree,depth=32km.

                                               （近场和远场地震同震垂直位移）