def movie(wfl, vel, par):
if (not par.has_key('vbias')): par['vbias'] = 1.5
fdmod.wom(wfl + 'm', wfl, vel, par['vbias'], par)
Plot(wfl + 'm', fdmod.wgrey('pclip=99', par))
for i in range(0, par['nt'] / par['jsnap'], 1):
fdmod.wframe(wfl + '-' + str(i), wfl + 'm', i, 'pclip=99', par)
Result(wfl + '-' + str(i), 'Overlay')
def movie(wfl,vel,par):
if(not par.has_key('vbias')): par['vbias']=1.5
fdmod.wom(wfl+'m',wfl,vel,par['vbias'],par)
Plot( wfl+'m',fdmod.wgrey('pclip=99',par))
for i in range(0,par['nt']/par['jsnap'],1):
fdmod.wframe(wfl+'-'+str(i),wfl+'m',i,'pclip=99',par)
Result(wfl+'-'+str(i),'Overlay')
def modeling(wav, vel, den, ss, rr, qq, dat, wfl, par):
fdmod.awefd(dat + '-tmp', wfl, wav, vel, den, ss, rr, 'free=y', par)
Flow(dat, dat + '-tmp',
'window f2=%(kt)d | pad end2=%(kt)d | put o2=%(ot)g' % par)
# wavefield movie
fdmod.wom(wfl + 'om', wfl, vel, 4.0, par)
Plot(wfl + 'om', fdmod.wgrey('pclip=99', par), view=1)
for f in range(0, par['nt'] / par['jsnap'], 2):
tag = "-%02d" % f
fdmod.wframe(wfl + 'om-' + tag, wfl + 'om', f, 'pclip=99', par)
Result(wfl + tag, [wfl + 'om-' + tag, ss, rr, qq], 'Overlay')
def modeling(wav,vel,den,ss,rr,qq,dat,wfl,par):
fdmod.awefd(dat+'-tmp',
wfl,
wav,
vel,
den,
ss,rr,'free=y',par)
Flow(dat,
dat+'-tmp',
'window f2=%(kt)d | pad end2=%(kt)d | put o2=%(ot)g' %par)
# wavefield movie
fdmod.wom(wfl+'om',wfl,vel,4.0,par)
Plot( wfl+'om',
fdmod.wgrey('pclip=99',par),view=1)
for f in range(0,par['nt']/par['jsnap'],2):
tag = "-%02d"%f
fdmod.wframe( wfl+'om-'+tag,wfl+'om',f,'pclip=99',par)
Result(wfl+tag,[wfl+'om-'+tag,ss,rr,qq],'Overlay')
