Question

fortran输出txt文件

比如我想用用fortran输出一个txt文件，里面的内容写着
压力 温度 密度
500 30 1000

求程序代码和解释。

Answer 1

给一个代码做个输出（你可以看到代码里只有计算，没有输出）。
在/claster/book1/user/data 这个路径中输出一个TXT文件。
输出源代码中所计算出的几个数据。

原代码如下：

c this just plots
integer sz,N,i,j
parameter(sz=500)
real tm(sz),tp(sz),jmn,jmx
real dms(sz),dkin(sz),dpot(sz)
real dcom(sz,3),dvel(sz,3),dang(sz,3)

c read in the disk data
jmn=100000000000.
jmx=0.
N=0
open(unit2,file='ints.star',status='old')
do i=1,sz
read(2,*,end=10)tm(i),dms(i),dkin(i),dpot(i),
dcom(i,1),dcom(i,2),dcom(i,3),
dvel(i,1),dvel(i,2),veel(i,3),
dang(i,1),dang(i,2),dang(i,3)
do j=1,3
jmn=min(jmn,dang(i,j))
jmx=max(jmx,dang(i,j))
end do
N=N+1
end do
10 close(2)

c Jz
call PGPAGE
call PGSLW(3)
call PGSVP(0.2,0.8,0.2,0.8)
call PGWINDOW(tm(1)-0.1,tm(N)+0.1,jmn,jmx*1.05)
call PGBOX('BCNST'0.,0,'BCNST'0.,0)
CALL PGLAB('t','J\\dz\\u')

do i=1,N
tp(i)=dang(i,3)
end do
call PGLINE(N,tm,tp)

c Jx
call PGPAGE
call PGSLW(3)
call (0.2,0.8,0.2,0.8)
call PGWINDOW(tm(1)-0.1.tm (N)+0.1,jmn,jmx*3.)
call PGBOX('BCNST',0.,0,'BCNST',0.,0)
call PGLAB（'t','J\\dx\\u,''）

do i=1,N
tp(i)=dang(i,1)
end PGLINE(N,tm,tp)

c Jy
call PGPAGE
call PGSLW(3)
call (0.2,0.8,0.2,0.8)
call PGWINDOW(tm(1)-0.1.tm (N)+0.1,jmn,jmx*1.05.)
call PGBOX('BCNST',0.,0,'BCNST',0.,0)
call PGLAB('t','J\\dx\\u','')
do i=1,N
tp(i)=dang(i,2)
end do
call PGLINE(N,tm,tp)

c Jtot
call PGPAGE
call PGSLW(3)
call (0.2,0.8,0.2,0.8)
call PGWINDOW(tm(1)-0.1.tm (N)+0.1,0,,jmx*4.)
call PGBOX('BCNST',0.,0,'BCNST',0.,0)
call PGLAB('t','J\\dx\\u','')
do i=1,N
tp(i)=sqrt(dang(i,1)**2+dang(i,2)**2+dang(i,3)**2)
end do
call PGLINE(N,tm,tp)
return
end

c end of subroutine PlotJz

Answer 2

character(len = 20) :: filename = "homework.txt" !设置文件名
integer, parameter :: fileid = 10 !设置write输出位置
。。。。。。
open(unit = fileid, file = filename)
。。。。。。
write (fileid, "('压力 温度 密度')")
。。。。。。
write (fileid, "(3I4)")500，30，1000
。。。。。。
close(fileid)

Answer 3

open（12，file="result.txt",status='new'）
write(12,'(a)') "压力 温度 密度"
write(12,'(3i10)')500,30,1000
close(12)