同一段包含su与mada命令的代码两种实现方式(不同部分粗体标注了一下)
方法1.
from rsf.proj import *
os.system("suplane>plane.su")
os.system("segyhdrs<plane.su |segywrite tape=plane.segy")
Flow('plane tfile hfile bfile','plane.segy','segyread endian=1 tfile=${TARGETS[1]} hfile=${TARGETS[2]} bfile=${TARGETS[3]}')
End()
方法2.
from rsf.proj import *
import rsf.suproj
suproject = rsf.suproj.SUProject()
suproject.Flow('plane.su',None,'suplane')
suproject.Flow('plane.segy','plane.su','segyhdrs |segywrite tape=$TARGET')
Flow('plane tfile hfile bfile','plane.segy','segyread endian=1 tfile=${TARGETS[1]} hfile=${TARGETS[2]} bfile=${TARGETS[3]}')
End()