[Fortran] 纯文本查看 复制代码
subroutine sub2()
implicit none
integer,parameter:: i0=20
integer,parameter:: j0=20
integer,parameter:: k0=20
integer i
integer j
integer k
real X
real Y
real T
real::a(i0,j0,k0)
real::b(i0,j0,k0)
real::c(i0,j0,k0)
real::d(i0,j0,k0)
real::ans(i0,j0,k0)
call sub1(X,Y,T,a,b,c,d,ans)
open(unit=11,file='data2.txt')
do i=1,i0
do j=1,j0
do k=1,k0
if (ans(i,j,k)>-0.01.and.ans(i,j,k)<0.01) then
write(11,"('(',f5.2,',',f5.2,','f5.2')=',(2X,F9.6),(2X,F9.6),(2X,F9.6),(2X,F9.6),(2X,F9.6))") X,Y,T,a(i,j,k),b(i,j,k),c(i,j,k),d(i,j,k),ans(i,j,k)
end if
end do
end do
end do
return
end