sconstruct

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

from rsfproj import *
import fdmod,gfield

par = {
    'nx':1200, 'ox':0, 'dx':1.0, 'lx':'x','ux':'m',
    'nz':800,  'oz':0, 'dz':1.0, 'lz':'z','uz':'m',
    'ru':30,'rv':6,'aa':1,'ff':45
    }
fdmod.param(par)

# ------------------------------------------------------------

gfield.execute('gg',par)
Result('gg',fdmod.cgrey('color=E',par))

# ------------------------------------------------------------

Result('ggz','gg','stack norm=y axis=1 | graph')
Result('ggx','gg','stack norm=y axis=2 | graph')

# ------------------------------------------------------------

Flow(  'vo','gg','math output="1500+x1"')
Result('vo',fdmod.cgrey('bias=1500 allpos=y',par))

Flow('nu','gg','scale rscale=0.1')
Flow('vv','vo nu',
     '''
     math v=${SOURCES[0]} n=${SOURCES[1]}
     output="v/sqrt(1+n)"
     ''' % par)
Result('vv',fdmod.cgrey('bias=1500 allpos=y',par))

Flow('ivo','vo','math output=1/input^2 | stack')
Flow('ivv','vv','math output=1/input^2 | stack')

Result('ivall','ivo ivv',
       '''
       cat axis=2 space=n ${SOURCES[1]} |
       graph
       ''')

End()