def geometry(ss,rr,qq,par):
# source coordinates
Plot(ss,fdmod.ssplot('',par))
Plot('win'+ss,ss,winss('',par))
# receiver coordinates
Plot(rr,fdmod.rrplot('plotfat=12',par))
# image coordinates
fdmod.makebox(qq,par['zminq'],par['zmaxq'],par['xminq'],par['xmaxq'],par)
Plot(qq,fdmod.bbplot('plotcol=1',par))
def geometry(ss, rr, qq, par):
# source coordinates
Plot(ss, fdmod.ssplot('', par))
Plot('win' + ss, ss, winss('', par))
# receiver coordinates
Plot(rr, fdmod.rrplot('plotfat=12', par))
# image coordinates
fdmod.makebox(qq, par['zminq'], par['zmaxq'], par['xminq'], par['xmaxq'],
par)
Plot(qq, fdmod.bbplot('plotcol=1', par))
def model(rr,par):
Flow('zero-2d',None,
'''
spike nsp=1 mag=0.0
n1=%(nz)d o1=%(oz)g d1=%(dz)g
n2=%(nx)d o2=%(ox)g d2=%(dx)g |
put label1=%(lz)s label2=%(lx)s unit1=%(uz)s unit2=%(ux)s
''' % par)
Flow('vz-2d','zero-2d',
'''
spike nsp=5
nsp=5 k1=10,40,70,100,130 l1=39,69,99,129,%(nz)d mag=0.2,0.4,0.6,0.8,1.0
n1=%(nz)d o1=%(oz)g d1=%(dz)g
n2=%(nx)d o2=%(ox)g d2=%(dx)g |
put label1=%(lz)s label2=%(lx)s unit1=%(uz)s unit2=%(ux)s |
add add=%(vp)f
''' % par)
Flow('fault-2d','zero-2d',
'''
spike nsp=1 k1=40 mag=1.0 l1=%(nz)d k2=60 l2=%(nx)d p2=1
n1=%(nz)d o1=%(oz)g d1=%(dz)g
n2=%(nx)d o2=%(ox)g d2=%(dx)g |
put label1=%(lz)s label2=%(lx)s unit1=%(uz)s unit2=%(ux)s
''' % par)
Flow('const-2d','zero-2d',
'''
spike nsp=1 mag=1.0 k1=40 l1=%(nz)d k2=1 l2=59
n1=%(nz)d o1=%(oz)g d1=%(dz)g
n2=%(nx)d o2=%(ox)g d2=%(dx)g |
put label1=%(lz)s label2=%(lx)s unit1=%(uz)s unit2=%(ux)s
''' % par)
Flow('vp-2d','vz-2d','window')
Flow('ro-2d','zero-2d','math output="%(ro)g"' %par)
fdmod.makebox('box',0.2,0.4,0.4,0.9,par)
Plot('box',fdmod.bbplot('',par))
Plot('vp-2d',fdmod.cgrey('allpos=y pclip=100 bias=1.5 ',par))
Plot('ro-2d',fdmod.cgrey('bias=2. allpos=y',par))
Result('vp-2d','vp-2d ss-2d rr-2d box','Overlay')
Result('ro-2d','ro-2d ss-2d','Overlay')