sconstruct

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

from rsfproj import *

Fetch('WGstack.H','book/iee')
Flow('dgulf','WGstack.H','dd form=native | window n1=1000 f2=1020 n2=200')

Fetch('wz.27.H','wz')
Flow('yc','wz.27.H','dd form=native | window n1=400 n2=64| tpow tpow=2')
Flow('yc2','yc','window f2=21')
Flow('yc27','yc yc2','window n2=20 | pad n2=21 | cat axis=2 ${SOURCES[1]}')

Flow('sigmoid',None,
     '''
     sigmoid d1=.004 n1=200 d2=.008 n2=200 |
     smooth rect1=3 diff1=1 |
     smooth rect1=3
     ''')

Flow('conflict',None,'conflict -')

for example in ('sigmoid','conflict','yc27','dgulf'):
    lom = 'lom-'+example
    lag = 'lag-'+example
    mpa = 'mpa-'+example
    Flow([lom,lag],example,
         'lopef a=7,2 w=20,%d lag=${TARGETS[1]}' %
         [8,4][int(example=='conflict')])
    if example=='conflict':
        rect1=1
    elif example=='dgulf':
        rect1=4
    else:
        rect1=3
    Flow(mpa,[example,lom],'decon filt=${SOURCES[1]} rect1=%d' % rect1)
    if example=='dgulf':
        pclip=96
    elif example=='yc27':
        pclip=99.5
    else:
        pclip=100
    Plot(example,'grey pclip=%g title=INPUT' % pclip)
    Plot(mpa,'grey pclip=%g title=LOMOPLAN' % pclip)
    Result(example,[example,mpa],'SideBySideAniso')
    
End()