1、Introducing ekional source model to account for some of the problems mentioned.
2、Rupture which is initially allowed at any point on the infinite fault plane is firsr restricted to this zone. Secondly, we allow rupture only happen within a certain radius R away from the previously defined point source location labeled center point. The radius is the only variabe in our model controlling the size of the rupture surface.
Kiwi has mentioned the disadvantages of traditional dscretization includes these errors:
1、 Each point in the sub-fault incorrectly radiates the same Green's function waveforms; (它会影响最低频率的滤波)
2、The temporal shifts of the signal radiated by different parts of the subfaults are neglected;（难道不是体现为破裂传播速度的不同吗？）
3、Suppression of local inhomogeneities of the moment tensor density inside the sub-fault.
Give a simplerule about the value of rupture velocity; And the size of sub-faullt is made dependent on the maximum frequency and the maximum slowness.
这里的Greenfunction内插方法使用的是bilinear interpolation, module is 'interpolation'; The API to read and write green's funtion database and the Green's function catche is implemented in 'gfdb';
Conduction of filtering and tapering: