[热点讨论]free-field边界问题 [复制链接]

代码如下：set dyn=on large tF;aB*  
ini xvel=0 yvel=0 zvel=0 NR3IeTd  
ini xdis 0 ydis 0 zdis 0 }\
irr9,  
free x y z range y -.1 .1 8(+X0}  
def setup 5<S1,u5  
omega = 4.0 * pi Psv-y  
alfa = 2.2 D2hvf^g'*  
beta = 3.75e-1 ]:f1r8<3p  
gamma = 8.0 AfY(+w6!K  
end .3Ap+V8?  
def wave :@p`E}1r{  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) kBT cND|  
end 'IgtBd|K>  
setup c.XLEjV|  
apply xacc=1.0 hist=wave  range  y -.1 .1 a@X'oV`(2b  
apply ff @e s
lF  
h(zi$V  
hist gp xvel 140,0,0 ^8\pJg_0  
hist gp xvel 200,98,0 1"e=Zqn$)  
hist gp yvel 200,98,0 G(4k#jB  
hist gp xdis 122,41,0 InDR\=o  
hist gp xdis 122,41,40 SHaZ-d  
hist gp ydis 122,41,0 N7e^XUG   
hist gp ydis 122,41,40 vuK 5DG4  
hist dytime LD5'4,%-  
hist wave S
Y{J  
set dyn damp local 0.15 ;(approx 5% damping) <.AIVp  
solve age 4 PK~okz4b  
save dam3.sav Zdak))7  
EYQ!ELuF  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 O8lOr(|l  

声明下y是竖直方向