子程序调用的问题

ji460600985

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

查看代码

打印?

01	module data_fixed  !save common data using module

02	integer,parameter :: row=32

03	integer,parameter ::n=row**3

04	integer,parameter :: fileid=10             !   unit=fileid

05	character(len=800) ::filenameF="C:\Users\Administrator\Desktop\blend 413K\F0.6 DPD Dynamics (2)\F0.6.00020000.F.Dpd_den"  !the site of the file that you want to load

06	character(len=800) ::filenameP="C:\Users\Administrator\Desktop\blend 413K\F0.6 DPD Dynamics (2)\F0.6.00020000.P.Dpd_den"

07	character(len=800) ::filenameT="C:\Users\Administrator\Desktop\blend 413K\F0.6 DPD Dynamics (2)\F0.6.00020000.T.Dpd_den"

08	end module

09

10	Program Test             !read file from external file

11	use data_fixed       ! call the data module

12	implicit none

13

14	real :: array_F(n),array_P(n),array_T(n)

15

16	call readfile(filenameF,array_F)   !call the subroutine

17	call readfile(filenameP,array_P)

18	call readfile(filenameT,array_T)

19	write(*,*) array_F(1:10),array_P(1:10),array_T(1:10)

20

21	end program

22

23	subroutine readfile(a,b)

24	use data_fixed        !call the data module

25	implicit none

26

27	character(len=80) :: tempstr

28	character(len=800) ::a

29	real :: b(n)

30

31	integer :: counter=0

32

33	open( unit=fileid , file=a ,access='sequential',status='old',recl=6)

34	do while (counter < n )   ! pass the first line data

35	counter=counter+1

36	if ( counter ==1 ) then

37	read(fileid,"(A80)") tempstr

38

end if

39	read(fileid,*) b(counter)

40

end do

41	close(fileid)

42	end subroutine

为什么array_F调用成功，而array_P和array_T调用出错啊？

fcode

你运行后遇到了什么错误？请告知错误信息。

如果你要跳过第一行，按照沙发代码那样写不是更好吗？

ji460600985

li913 发表于 2014-12-7 20:59
open( unit=fileid , file=trim(a) ,status='old')
read(fileid,*)
do i=1,n

我是从第二行开始读取数据的？所以要跳过第一行！

li913

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

查看代码

打印?

1	open( unit=fileid , file=trim(a) ,status='old')

2	read(fileid,*)

3

do i=1,n

4	read(fileid,*) b(i)

5

end do

6	close(fileid)