鉛直速度自由表面境界条件の変更について

自由表面における境界条件を検討し、GMS1.8.0でGMSソルバーの実装を変更しました。

 

変更前

uz1(i,j,0)=uz1(i,j,1)+(vp*vp-2.*vs*vs)/vp/vp*(

-dtdx1/dtdz1*(ux1(i ,j ,1 )-ux1(i-1,j ,1 ))

-dtdy1/dtdz1*(uy1(i ,j ,1 )-uy1(i ,j-1,1 )) )

 

変更後

uz1(i,j,0)=uz1(i,j,1)+(vp*vp-2.*vs*vs)/vp/vp*(

+dtdx1/dtdz1*(ux1(i ,j ,1 )-ux1(i-1,j ,1 ))

+dtdy1/dtdz1*(uy1(i ,j ,1 )-uy1(i ,j-1,1 )) )