- character没有指明长度
program main
implicit none
integer::it,i,j,tmax=20,ny=100,nx=100
real::u(100,100),v(100,100),theta(100,100),A=-1/3,B=1/3
character::filename
do it=1,tmax
do j=2,ny-1
do i=2,nx-1
u(i,j)=Au(i,j)+B((u(i-1,j)+u(i+1,j))+u(i,j-1)u(i,j+1))
enddo
enddo
do j=ny-1,2,-1
do i=nx-1,2,-1
v(i,j)=Av(i,j)+B*(v(i-1,j)+v(i+1,j)+v(i,j-1)+v(i,j+1))
enddo
enddo
theta=(v+u)/2
write(filename,"(a,i2,a)") ‘theta’,it,’.txt’
open(unit=11,file=filename)
do i=1,nx
do j=1,ny
write(11,"(f8.4)") theta(i,j)
end do
enddo
close(11)
enddo
end
f:end of record的n种原因
最新推荐文章于 2021-12-30 12:56:49 发布