自由表面における境界条件を検討し、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 )) )