作者PTT0207445 (我是乡民)
看板Fortran
标题Re: [问题] 档案读取问题
时间Thu Jan 17 09:10:25 2019
program main
implicit none
integer,
parameter :: X
=5, N
=40, N2
=N
/2
real,
dimension(N) :: ex
real,
dimension(N2) :: ans
integer :: i
open(
30,
status='old',
file='INP.TXT')
open(
31,
file='OPT.TXT')
do i
=1,N
read(
30,
*) ex(i)
enddo
ans
=ex(
1:N2)
*ex(N2
+1:N)
do i
=1,N2
write(
31,
*) ans(i)
enddo
close(
30)
close(
31)
end program main
※ 引述《btoker (阿瑞)》之铭言：
: 有一个变数X=5 会产生下面的文字档
: 文字档共有40列 如下
: 10.66
: 8.43
: 72.47
: 50.29
: 29.91
: 10.66
: 8.43
: 72.47
: 50.29
: 29.91
: 99.99
: 8.43
: 72.47
: 50.29
: 29.91
: 99.99
: 8.43
: 72.47
: 50.29
: 29.91
: 1.
: 2.
: 3.
: 4.
: 5.
: 1.
: 2.
: 3.
: 4.
: 5.
: 6.
: 7.
: 8.
: 9.
: 10.
: 6.
: 7.
: 8.
: 9.
: 10.
: 可以明显的发现这个数列是5个5个一组的
: 我想把第1列的10.66与第21列的1.0相乘
: 第2列的 8.43与第22列的2.0相乘 输出到另一个档案
: 初步错误的程式码如下
: OPEN(UNIT=30,STATUS="OLD",FILE="INP.TXT")
: OPEN(UNIT=31,FILE="OPT.TXT")
: DO I=1,X*4,1
: READ(30,*)EX
: DO J=1,X*4,1
: READ(30,*)AX
: END DO
: ANS=EX*AX
: WRITE(31,*)ANS
: END DO
: CLOSE(30)
: CLOSE(31)
: 请问该如何修正
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 86.156.211.166
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Fortran/M.1547687427.A.D9C.html