xumi1993 / seispy

Python module of seismology and receiver functions

Home Page:https://seispy.xumijian.me

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

station naming with network or not

ShihChi90 opened this issue · comments

Describe the bug
The naming of station name seem not consistant in different modules
for example in rfcorrect.py the staname is {net}.{station}
but in ccpprofile.py the staname is just station
since they both read station from same list (in ccp.cfg), this would cause some inconvinience

For example i have a station name TR01 in net TW, and my station list is also naming it TR01.
When saving data in rfdepth, it would be named as TW.TR01. But when doing ccpstacking, it would say that the station TR01 is not in rfdepth structure. Then I would have to rename my station in list to TW.TR01, it would work as expected.

Expected behavior
using same config (station list) for rfdepth and ccpprofile should not cause this error

Error raised
2024-05-01 10:01:10,591 [CCP] WARNING: TR01 does not in RFdepth structure
image

Platform (please complete the following information):

  • Linux
  • seispy 1.3.6+ dev

Additional context
I think when saving data in rfdepth, just use the station name without the network.

The CCP module will read full station name from stalist or stackfile, so you need to set up the station name to net.sta as you have mentioned. I recommend to use commands rf2depth and ccp_profile to avoid this issue, instead of invoking modules in a python script.

Please refer to here for the standard procedure for CCP stacking.

I understand your concern. A better way is separately saving net and sta in the RFdepth file, but modifying the structure will influence the use by other users.