sconstruct

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

from rsfproj import *

Flow('time',None,
     'math n1=200 d1=1 o1=1 output="(x1-1)*0.4" ')
Flow('head','time',"window n1=41 | math output='exp(0.25*input)+11'")
Flow('alias','head',"math output='sin(input)*0.5' | put head=$SOURCE")

doty = '''
dots dots=0 connect=0 gaineach=0 constsep=1 strings=0 clip=1.2 label1=Sample
'''

Flow('data','alias','bin1 nx=200 dx=0.4 x0=0')
Plot('data',doty + ' strings=1 d1=1 o1=1')

doty = doty + ' connect=1 d1=1 o1=1'

for size in [3,5]:
    sine = 'subsine%d' % size
    Flow(sine,'alias','levint nx=200 dx=0.4 x0=0 na=%d' % size)
    Plot(sine,doty)
    Result(sine,['data',sine],'Overlay')

###########################################################################

End()