kk=0
do ii=1,60
do jj=2,61
!射线向右 判断射线与边界有几条交点
if ((pointx(ii)>xm .and. pointxx(jj)>xm) .or. (pointx(ii)>xm .and. pointxx(jj)<xm) .or. (pointx(ii)<xm .and. pointxx(jj)>xm)) then
if (pointy(ii)>=ym .and. pointyy(jj)<ym) .or. (pointy(ii)<=ym .and. pointyy(jj)>ym) then
if (MOD(kk+1,2.)/=0.) then
kk=kk+1
ux(i,j,k)=0.
uy(i,j,k)=0.
uz(i,j,k)=0.
ep1(i,j,k)=1.
end if
end if
end if
stop
end do
end do
return