Fortran求解线性方程组，怎么样调用子程序

landaifeng

各位大神，请问我有求解线性方程组的计算MAXA数组的子程序，也有数组，我读入了数组a(7,7),想计算MAXA数组，下面的程序调用这个FMAXA的时候报错了，
（1）warning FOR4227: argument LMT (number 3) in reference to procedure FMAXA from procedure MAIN22 incorrect: is not an array（2） warning FOR4227: argument MAXA (number 4) in reference to procedure FMAXA from procedure MAIN22 incorrect: is not an array
请问大家，应该怎么调试？

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

查看代码

打印?

01	program main22

02	dimension a(7,7)

03	open (1,file='E:\1.txt')

04	do 10 i=1,7

05	read(1,*)(a(i,j),j=1,7)

06	10 continue

07	write(*,*)((a(i,j),j=1,7),i=1,7)

08

09

10	call FMAXA(NN1, NE, LMT, MAXA, NWK, NPF, NDF)

11

END

12

13	SUBROUTINE FMAXA(NN1, NE, LMT, MAXA, NWK, NPF, NDF)

14	!This program forms the MDE address matrix MAXA of [K]

15	IMPLICIT REAL*8(A-H,O-Z)

16	IMPLICIT INTEGER*4(I-N)

17	DIMENSION MAXA(NPF), LMT(NDF,NE)

18	MAXA=0; NWK=0

19

MAXA(1)=1

20	DO I=2, NN1

21

IP=I-1

22

IG=IP

23	DO IE=1, NE

24	DO J=1, NDF

25	IF(LMT(J,IE).EQ.IP)THEN

26	DO K=1, NDF

27	IF(LMT(K,IE).GT.0.AND.LMT(K,IE).LE.IG) IG=LMT(K,IE)

28

ENDDO

29

END IF

30

ENDDO

31

ENDDO

32	MAXA(I)=MAXA(I-1)+IP-IG+1

33

ENDDO

34	NWK=MAXA(NN1)-1

35

RETURN

36	write(*,*)MAXA(I)

37

END

楚香饭

好好看看书吧。你这错误太低级了。

你完全没有搞懂什么是变量名，什么是数组，什么是实参虚参。

你读入了一个数组，叫 a，但是传递的实参全部没有定义，也没有与 a 有任何关系。

其他的虚参，你也没有传入与其意义对应的实参。