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')