from rsf.proj import *
Flow('vel7',None,
'''
spike n1=513 d1=5 k1=98 n2=513 d2=5 k2=98 |
math output="550+0.00015*(x1-800)*(x1-800)+0.0001*(x2-500)*(x2-500)"
''')
Plot('vel7',
'''
grey transp=n pclip=100 screenratio=1 color=i
label1="X" label2="Z" scalebar=y poly=y bias=1200
unit1=m unit2=m title=a barlabel=Velocity barunit=m/s
''')
Result('vel7',
'''
grey transp=n pclip=100 screenratio=1 color=j
label1="X" label2="Z" scalebar=y poly=y mean=y
unit1=m unit2=m title=Velocity barlabel=Velocity barunit=m/s
''')
nt=1250
dt=0.002
nv=2
Flow('source7',None,
'''
spike n1=%d d1=%g k1=100 |
ricker1 frequency=20
'''%(nt,dt))
Result('source7','graph title="Source Wavelet" ')
Flow('wavp7',['source7','vel7'],
'''
pspifcos2bs vel=${SOURCES[1]} dt=%g nt=%d nv=%d isx=256 isz=256 nb=0 c=0.01
'''%(dt,nt,nv) )
Plot('wavp7',
'''
window j3=20 |
grey poly=y label2="Z" label1="X" title="pspif2" unit1="m" unit2="m"
yreverse=y transp=n gainpanel=each pclip=100
''',view=1)
Flow('wavf7',['source7','vel7'],
'''
ffd2dehf vel=${SOURCES[1]} dt=%g nt=%d isx=256 isz=256 nbb=0 nbl=0 nbt=0 nbr=0 ax=1.0 az=1.0 factor=0.9
'''%(dt,nt) )
Plot('wavf7',
'''
window j3=20 |
grey poly=y label2="Z" label1="X" title="ffd2"
yreverse=y transp=n gainpanel=each pclip=100 unit1="m" unit2="m"
''',view=0)
Flow('wavfd7',['source7','vel7'],
'''
fd2bs vel=${SOURCES[1]} dt=%g nt=%d isx=256 isz=256 nb=0 c=0.01
'''%(dt,nt) )
Plot('wavfd7',
'''
window j3=20 |
grey poly=y label2="Z" label1="X" title="fd2" unit1="m" unit2="m"
yreverse=y transp=n gainpanel=each pclip=100
''',view=0)
Flow('wavsp7',['source7','vel7'],
'''
sp vel=${SOURCES[1]} dt=%g nt=%d isx=256 isz=256 nb=0 c=0.01 ax=1.0 az=1.0 factor=0.7
'''%(dt,nt) )
Plot('wavfd7',
'''
window j3=20 |
grey poly=y label2="Z" label1="X" title="fd2" unit1="m" unit2="m"
yreverse=y transp=n gainpanel=each pclip=100
''',view=0)
Flow('wavfd_7','wavfd7','window n3=1 f3=800')
Flow('wavffd_7','wavf7','window n3=1 f3=800')
Flow('wavpspi_7','wavp7','window n3=1 f3=800')
Flow('wavsp_7','wavsp7','window n3=1 f3=800')
Flow('compare_7',['wavfd_7','wavffd_7','wavpspi_7'],
'''
cat ${SOURCES[1]} ${SOURCES[2]} axis=3
''')
Plot('wavfd_7',
'''
grey poly=y label2="Z" label1="X" title=b screenratio=1
yreverse=y transp=n gainpanel=all pclip=99 unit1=m unit2=m
''' )
Plot('wavffd_7',
'''
grey poly=y label2="Z" label1="X" title=c screenratio=1
yreverse=y transp=n gainpanel=all pclip=99 unit1=m unit2=m
''' )
Plot('wavpspi_7',
'''
grey poly=y label2="Z" label1="X" title=d screenratio=1
yreverse=y transp=n gainpanel=all pclip=99 unit1=m unit2=m
''' )
Plot('wavsp_7',
'''
grey poly=y label2="Z" label1="X" title=d screenratio=1
yreverse=y transp=n gainpanel=all pclip=99 unit1=m unit2=m labelsz=6
''' )
Result('wavfd','wavfd_7',
'''
grey poly=y label2="Depth" label1="Distance" title= screenratio=1
yreverse=y transp=n gainpanel=all pclip=99 unit1=m unit2=m labelsz=6
''' )
Result('wavffd','wavffd_7',
'''
grey poly=y label2="Depth" label1="Distance" title= screenratio=1
yreverse=y transp=n gainpanel=all pclip=99 unit1="m" unit2="m" labelsz=6
''' )
Result('wavpspi','wavpspi_7',
'''
grey poly=y label2="Depth" label1="Distance" title= screenratio=1
yreverse=y transp=n gainpanel=all pclip=99 unit1="m" unit2="m" labelsz=6
''' )
Result('wavsp','wavsp_7',
'''
grey poly=y label2="Depth" label1="Distance" title= screenratio=1
yreverse=y transp=n gainpanel=all pclip=99 unit1="m" unit2="m" labelsz=6
''' )
End() |