fluidcompute.f90

来自「visual fortran 源码」· F90 代码 · 共 26 行

subroutine fluidcompute
use varible
use non_common
open (unit=101,file='simpleresult.txt')
call grid()
!write(*,*)"xu"
!write(*,"(F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4)")(xu(I),I=1,NI)
!write(*,*)"yv"
!write(*,"(F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4)")(yv(I),I=1,Nj)
call setup1()
!write(*,"(F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4)")(x(I),I=1,NI)
!write(*,*)
!write(*,"(F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4)")(y(I),I=1,Nj)
call start()
!write(101,"(F8.1,F8.1,F8.1,F8.1,F8.1,F8.1,F8.1,F8.1,F8.1,F8.1,F8.1,F8.1,F8.1,F8.1)")((t(I,j),I=1,NI),j=1,nj)
do while(lstop==.FALSE.)
call dense()
call bound()
!write(101,"(F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4,F8.4)")((aim(I,J),I=1,NI),J=1,NJ)
call setup2()

!write(101,*)"next"
end do
!write(*,"(F8.2,F8.2,F8.2,F8.2,F8.2,F8.2,F8.2,F8.2,F8.2,F8.2,F8.2,F8.2)")((v(I,j),I=1,NI),j=1,nj)
write(101,"(F8.2,F8.2,F8.2,F8.2,F8.2,F8.2,F8.2)")((t(I,j),I=1,NI),j=1,nj)
end subroutine fluidcompute