子程序调用的问题

ji460600985

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

module data_fixed  !save common data using module 
     integer,parameter :: row=32              
     integer,parameter ::n=row**3
     integer,parameter :: fileid=10             !   unit=fileid  
     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
     character(len=800) ::filenameP="C:\Users\Administrator\Desktop\blend 413K\F0.6 DPD Dynamics (2)\F0.6.00020000.P.Dpd_den"
     character(len=800) ::filenameT="C:\Users\Administrator\Desktop\blend 413K\F0.6 DPD Dynamics (2)\F0.6.00020000.T.Dpd_den"
 end module

 Program Test             !read file from external file
     use data_fixed       ! call the data module
     implicit none

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

     call readfile(filenameF,array_F)   !call the subroutine
     call readfile(filenameP,array_P)
     call readfile(filenameT,array_T)
     write(*,*) array_F(1:10),array_P(1:10),array_T(1:10)

 end program

 subroutine readfile(a,b)  
     use data_fixed        !call the data module
     implicit none

     character(len=80) :: tempstr
     character(len=800) ::a
     real :: b(n)

     integer :: counter=0

     open( unit=fileid , file=a ,access='sequential',status='old',recl=6)
     do while (counter < n )   ! pass the first line data
         counter=counter+1  
         if ( counter ==1 ) then
             read(fileid,"(A80)") tempstr
         end if      
         read(fileid,*) b(counter)
     end do
     close(fileid)
 end subroutine

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

li913

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

open( unit=fileid , file=trim(a) ,status='old')
 read(fileid,*) 
 do i=1,n
 read(fileid,*) b(i)    
 end do
 close(fileid)

ji460600985

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

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

fcode

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

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