back.txt

圆柱绕流问题的涡方法数值模拟（fortran源码）

!--------------------壁面镜相-------------------

!   do k=1,kk
!      do L=1,knp
!	     if((xl2-0.02)<=dx(k,L).and.dx(k,L)<=xl2.and.(-yl2+0.02)<=dy(k,L).and.dy(k,L)<=(yl2-0.02))then
!		   dx(k,L)=xl2-dx(k,L)+xl2
!		 elseif(-xl2<=dx(k,L).and.dx(k,L)<=(-xl2+0.02).and.(-yl2+0.02)<=dy(k,L).and.dy(k,L)<=(yl2-0.02))then
!		   dx(k,L)=-xl2+(-xl2-dx(k,L))
!		 elseif((yl2-0.02)<dy(k,L).and.dy(k,L)<yl2.and.(-xl2+0.02)<=dx(k,L).and.dx(k,L)<(xl2-0.02))then
!		   dy(k,L)=yl2-dy(k,L)+yl2
!		 elseif(-yl2<dy(k,L).and.dy(k,L)<(-yl2+0.02).and.(-xl2+0.02)<=dx(k,L).and.dx(k,L)<(xl2-0.02))then
!		   dy(k,L)=-yl2+(-yl2-dy(k,L))
!		 elseif((xl2-0.02)<dx(k,L).and.dx(k,L)<xl2.and.(yl2-0.02)<dy(k,L).and.dy(k,L)<yl2)then
!		   dy(k,L)=yl2-dy(k,L)+yl2
!		   dx(k,L)=xl2-dx(k,L)+xl2
!		 elseif((xl2-0.02)<dx(k,L).and.dx(k,L)<xl2.and.-yl2<dy(k,L).and.dy(k,L)<(-yl2+0.02))then
!		   dy(k,L)=-yl2+(-yl2-dy(k,L))
!		   dx(k,L)=xl2-dx(k,L)+xl2
!		 elseif(-xl2<dx(k,L).and.dx(k,L)<(-xl2+0.02).and.(yl2-0.02)<dy(k,L).and.dy(k,L)<yl2)then
!		   dy(k,L)=yl2-dy(k,L)+yl2
!		   dx(k,L)=-xl2+(-xl2-dx(k,L))
!		 elseif(-xl2<dx(k,L).and.dx(k,L)<(-xl2+0.02).and.-yl2<dy(k,L).and.dy(k,L)<(-yl2+0.02))then
!		   dx(k,L)=-xl2+(-xl2-dx(k,L))
!		   dy(k,L)=-yl2+(-yl2-dy(k,L))
!		 else
!		 continue
!		 end if
!	  end do
 ! end do

!--------------------受力分析-------------------

 nstep=nstep+1

!end do
 open(unit=10,file='dx.txt')
 open(unit=20,file='dy.txt')
 write(10,"(1f8.4)")((dx(i,j),i=1,nlimit),j=1,knp)
 write(20,"(1f8.4)")((dy(i,j),i=1,nlimit),j=1,knp)
!write(*,"(8f8.4)")(btao(i),i=1,nall)
 close(10)
 close(20)
  !write(*,*)((ndp(i,j),i=1,nlimit),j=1,knp)