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

代码如下：set dyn=on large i$p2am
8f  
ini xvel=0 yvel=0 zvel=0 Cse`MP  
ini xdis 0 ydis 0 zdis 0 }aVzr}!  
free x y z range y -.1 .1 ,`Z4fz:  
def setup u=qaz7E  
omega = 4.0 * pi G u_\ySV/y  
alfa = 2.2 4eB oR%2o  
beta = 3.75e-1 <`"  
gamma = 8.0 rQlQ^W$=?  
end M;b3- i  
def wave Dq<DW2It>  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) g 8uq6U  
end rv(N0p/  
setup t_xK?``  
apply xacc=1.0 hist=wave  range  y -.1 .1 q2J
|koT  
apply ff />}zB![(K  
~]RfOpq^w  
hist gp xvel 140,0,0 boCi*]  
hist gp xvel 200,98,0 wu eDedz\  
hist gp yvel 200,98,0 DnJ `]r  
hist gp xdis 122,41,0 [;7zg@Sa  
hist gp xdis 122,41,40 j;b>~_ U%  
hist gp ydis 122,41,0 ,SNrcwv  
hist gp ydis 122,41,40 3M+rFB}tS  
hist dytime [e ztu9  
hist wave 0{%@"Fb0O  
set dyn damp local 0.15 ;(approx 5% damping) LxGD=b  
solve age 4 !4D?X\~"%  
save dam3.sav 3u[8;1}7Q  
?wl
RHVZ  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 ]|PTZ1?j  

声明下y是竖直方向