常用模型 June

#waxian velocity
from rsf.proj import *
import math
par = {
    'nz':400,  'dz':2.0,   'oz':0,
    'nx':400,  'dx':10.0,  'ox':0
    }
   
def reflector(formula,par):
    return '''
    math n1=%d o1=%g d1=%g n2=1 output="%s"
    ''' % (par['nx']/2,par['ox'],par['dx'],formula)

Flow('left',None,reflector('202',par))
Flow('meda',None,reflector('402',par))
Flow('rigt',None,reflector('202',par))

Flow('spk','left meda rigt',
     '''
     cat ${SOURCES[1:3]} axis=1 |
     unif3 n1=%(nz)d d1=%(dz)g v00=1,2
     ''' % par)
    
End()

#waxian velocity
from rsf.proj import *
import math
par = {
    'nz':400,  'dz':2.0,   'oz':0,
    'nx':400,  'dx':10.0,  'ox':0
    }
   
def reflector(formula,par):
    return '''
    math n1=%d o1=%g d1=%g n2=1 output="%s"
    ''' % (par['nx']/2,par['ox'],par['dx'],formula)

Flow('left',None,reflector('201',par))
Flow('meda',None,reflector('x1*0.1+201',par))
Flow('rigt',None,reflector('400',par))

Flow('spk','left meda rigt',
     '''
     cat ${SOURCES[1:3]} axis=1 |
     unif3 n1=%(nz)d d1=%(dz)g v00=1,2
     ''' % par)
    
End()