sconstruct

国外免费地震资料处理软件包

from rsfproj import *

Flow('data1',None,'spike n1=30 k1=15 label1=" " ')
Flow('data2','data1','causint')
Flow('data3','data1','causint adj=y')

Flow('data4',None,'spike n1=30 nsp=2 k1=10,15 mag=1,-1 label1=" " ')
Flow('data5','data4','causint')
Flow('data6','data4','causint adj=y')
Flow('data7','data5','causint')

data = []
for i in range(1,8):
    data.append('data%d' % i)

def plot(label):
    return '''cat axis=2 ${SOURCES[1:%d]} |
    dots labels="%s"
    yreverse=1 overlap=.8 constsep=1 connect=0
    ''' % (len(data),label)

Result('causint',data,plot("in1:C in1:C\` in1:in2:C in2:C\` in2:CC in2"))
Result('quiz',   data,plot("in1:C in1:C\` in1:in2:C in2:C\` in2:      "))

###########################################################################
End()