sconstruct

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

from rsfproj import *

Fetch('midpts.hh','midpts')

Flow('cdp','midpts.hh','window n3=1 j1=2 | dd form=native')
Flow('mutter','cdp','mutter half=n v0=1.4') 

Flow('vscan','mutter','vscan half=n v0=1.3 nv=81 dv=0.02 semblance=y')
Flow('vrms','vscan','pick rect1=50')

cases = ['cdp','mutter']
nmos = []

clip=4141100

for case in cases:
    nmo = 'nmo'+case
    nmos.append(nmo)
    
    Plot(case,'window n1=300 | grey wanttitle=n')
    Flow(nmo,[case,'vrms'],'nmo half=n velocity=${SOURCES[1]} str=0')
    Plot(nmo,'grey wanttitle=n clip=%g' % clip)
        
Result('mutter',cases,'SideBySideAniso',vppen='txscale=1.5')
Result('nmogath',nmos,'SideBySideAniso',vppen='txscale=1.5')

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