dvm_mainprogram.f90

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

program DVMmain
use varible

call setup()
call spini()
call para_sigma()
!call circle()
!write(*,"(10f8.4)") (sigma(i),i=1,nall)
open(unit=10,file='Gx.xls')
open(unit=20,file='Gy.xls')
!open(unit=30,file='pvx.xls')
!open(unit=40,file='pvy.xls')
!open(unit=50,file='px.xls')
!open(unit=60,file='py.xls')
write(10,"(1f12.6)") ((gx(i,j),i=1,nlimit),j=1,knp)
!write(*,*)"================="
write(20,"(1f12.6)") ((gy(i,j),i=1,nlimit),j=1,knp)
open(unit=30,file='bv1.txt')
open(unit=40,file='bv2.txt')
open(unit=50,file='px.txt')
open(unit=60,file='py.txt')

!write(30,"(1f8.4)")(bv(i,1),i=1,nall)
!write(40,"(1f8.4)")(bv(i,2),i=1,nall)
write(30,"(1f12.6)") ((pvx(i,j),i=1,nlimit),j=1,1)
write(40,"(1f12.6)") ((pvy(i,j),i=1,nlimit),j=1,1)
write(50,"(1f12.6)") ((px(i,j),i=1,nlimit),j=1,pnum)
write(60,"(1f12.6)") ((py(i,j),i=1,nlimit),j=1,pnum)
!write(*,*)"================="
!write(*,*) ((ndp(i,j),i=1,ns+1),j=1,knp)
!write(10,*)"============================================="
!write(*,*)((ndp(k,l),k=1,ns),L=1,knp)
!write(*,*)"done"
end program