Beispiel #1
0
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))
Beispiel #2
0
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))
Beispiel #3
0
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')