常用模型

最新推荐文章于 2020-08-29 22:38:30 发布

weixin_30276935

最新推荐文章于 2020-08-29 22:38:30 发布

阅读量115

点赞数

原文链接：http://www.cnblogs.com/zebra/archive/2011/04/11/2012038.html

版权

# waxian velocity
from rsf.proj import *
import math
par = {
' nz ' : 400 , ' dz ' : 2.0 , ' oz ' :0,
' nx ' : 400 , ' dx ' : 10.0 , ' ox ' :0
}

def reflector(formula,par):
return '''
math n1=%d o1=%g d1=%g n2=1 output="%s"
''' % (par[ ' nx ' ] / 2 ,par[ ' ox ' ],par[ ' dx ' ],formula)

Flow( ' left ' ,None,reflector( ' 202 ' ,par))
Flow( ' meda ' ,None,reflector( ' 402 ' ,par))
Flow( ' rigt ' ,None,reflector( ' 202 ' ,par))

Flow( ' spk ' , ' left meda rigt ' ,
'''
cat ${SOURCES[1:3]} axis=1 |
unif3 n1=%(nz)d d1=%(dz)g v00=1,2
''' % par)

End()

# waxian velocity
from rsf.proj import *
import math
par = {
' nz ' : 400 , ' dz ' : 2.0 , ' oz ' :0,
' nx ' : 400 , ' dx ' : 10.0 , ' ox ' :0
}

def reflector(formula,par):
return '''
math n1=%d o1=%g d1=%g n2=1 output="%s"
''' % (par[ ' nx ' ] / 2 ,par[ ' ox ' ],par[ ' dx ' ],formula)

Flow( ' left ' ,None,reflector( ' 201 ' ,par))
Flow( ' meda ' ,None,reflector( ' x1*0.1+201 ' ,par))
Flow( ' rigt ' ,None,reflector( ' 400 ' ,par))

Flow( ' spk ' , ' left meda rigt ' ,
'''
cat ${SOURCES[1:3]} axis=1 |
unif3 n1=%(nz)d d1=%(dz)g v00=1,2
''' % par)

End()