gen psrayp_like_file for SRF

[JouSF-1409]

在将psrayp迁移到SRF以及其他震相上，有几个东西不是很理解：

1. taup_time， 自己测试了下转换波的taup_time。taup -ph P40s,P60s -deg 60 -h 10这样的结果，实际上只是得到了35km处间断面Ps波的信息。是否需要根据layers的划分单独生成一个model来做（）。
2. 从taup获取的layer数组与深度关系。 psrayp.py line 21用的是np.where所以 fake_layers应该是浅于11km的层，我们不单独计算他们的rayp。但是在line 48里直接 在数组切片里layer -11，有点不太理解。

第一个问题很好解决。 这个我们可以使用DepModel和layers做插值后，生成tvel格式的速度文件给taup算。 但是还是得明确与深度的明确关系，以及怎么生成格式正确的rayplib。