关于子程序读取txt文件

freshman · 发表于 2019-2-19 19:36:07

本帖最后由 freshman 于 2019-2-19 21:32 编辑

我想用子程序来读取一个文件中的数据并且储存到数组中，然后另外一个子程序会使用该数组中的元素进行计算。
现在代码跑不通跪求大神帮忙看看

[Fortran] 纯文本查看 复制代码

1MODULE inputpower
2    IMPLICIT NONE
3    REAL k(50,152)
4END MODULE

[Fortran] 纯文本查看 复制代码

01SUBROUTINE getpow(xx,yy)
02    USE inputpower
03    IMPLICIT NONE
04    INTEGER i,j
05    OPEN(unit=99,file='test1.txt',status='old',form = 'formatted')
06        DO i=1,xx
07             READ(99,*) (k(j,i) ,j=1,yy)
08             !print*, (k(j,i) ,j=1,yy)
09        END DO
10    CLOSE(99)
11END SUBROUTINE getpow(xx,yy)

[Fortran] 纯文本查看 复制代码

01PROGRAM inpow  
02    USE inputpower
03    IMPLICIT NONE
04    INTEGER na,nr
05    na = 7
06    nr = 7
07    !read the input data
08    CALL getpow(na,nr) 
09     
10    END PROGRAM inpow

text文件内容如下
1 1 1 1 1 1 1
2 1 1 1 1 1 1
3 1 1 1 1 1 1
4 1 1 1 1 1 1
5 1 1 1 1 1 1
6 1 1 1 1 1 1
7 1 1 1 1 1 1

问题出在了ps：自学没几天，请大佬推荐比较好的fortran90书籍

freshman · 发表于 2019-2-19 19:45:43

可能代码特别不规范主要是参照了一些fortran77的代码

freshman · 发表于 2019-2-19 21:57:08

SUBROUTINE getpow(xx,yy)
USE inputpower
IMPLICIT NONE
INTEGER i,j
OPEN(unit=99,file='test1.txt',status='old',form = 'formatted')
      DO i=1,xx
         READ(99,*) (k(j,i) ,j=1,yy)
         !print*, (k(j,i) ,j=1,yy)
      END DO
CLOSE(99)
END SUBROUTINE getpow(xx,yy)
用排除法得到是红色的两行代码有问题

fcode · 发表于 2019-2-19 22:02:10

[Fortran] 纯文本查看 复制代码

01MODULE inputpower
02    IMPLICIT NONE
03    REAL k(50,152)
04END MODULE
05 
06SUBROUTINE getpow(xx,yy)
07    USE inputpower
08    IMPLICIT NONE
09    INTEGER i,j,xx,yy!//XX,YY需要定义
10    OPEN(unit=99,file='test1.txt',status='old',form = 'formatted')
11        DO i=1,xx
12             READ(99,*) (k(j,i) ,j=1,yy)
13             !print*, (k(j,i) ,j=1,yy)
14        END DO
15    CLOSE(99)
16END SUBROUTINE getpow!//这里不写括号
17 
18PROGRAM inpow  
19    USE inputpower
20    IMPLICIT NONE
21    INTEGER na,nr
22    na = 7
23    nr = 7
24    !read the input data
25    CALL getpow(na,nr) 
26     
27    END PROGRAM inpow

如果还有问题，请给出错误提示描述。

freshman · 发表于 2019-2-19 22:27:15

fcode 发表于 2019-2-19 22:02
[mw_shl_code=fortran,true]MODULE inputpower
IMPLICIT NONE
REAL k(50,152)

解决了
非常感谢！
谢谢您！！！！！！

[文件读写] 关于子程序读取txt文件

美女勋章

热心勋章

星光勋章

新人勋章

贡献勋章

管理勋章

帅哥勋章

爱心勋章

规矩勋章

元老勋章

水王勋章