unit 7的文件(efdc.inp)打开失败.

    Program NHM
    INTEGER*2 IC,JC,LINE,Depth_type,SurElv_Type,Max_Iter,
     &          PreCon_type,AdvecType,Step_Iterval
    REAL*4    OX,OY,DX,DY,IniDepth,IniSurElv,n,StartTime,EndTime,
     &          DT,RSQM
C
    OPEN(UNIT=11,FILE='EFDC.INP')
    DO LINE=1,10
        READ(11,*)
    END DO
    READ(11,*)OX,OY,IC,JC,DX,DY
    PRINT*,OX,OY,IC,JC,DX,DY
C
    DO LINE=1,13
        READ(11,*)
    END DO
    READ(11,*)Depth_type,SurElv_Type,IniDepth,IniSurElv,n
    PRINT*,Depth_type,SurElv_Type,IniDepth,IniSurElv,n
C
    DO LINE=1,7
        READ(11,*)
    END DO
    READ(11,*)StartTime,EndTime,DT
    PRINT*,StartTime,EndTime,DT
C
    DO LINE=1,8
        READ(11,*)
    END DO
    READ(11,*)Max_Iter,RSQM,PreCon_type
    PRINT*,Max_Iter,RSQM,PreCon_type
C
    DO LINE=1,8
        READ(11,*)
    END DO
    READ(11,*)AdvecType
    PRINT*,AdvecType
C
    DO LINE=1,5
        READ(11,*)
    END DO
    READ(11,*)Step_Iterval
    PRINT*,Step_Iterval
C

    STOP
    END
请问,上面的代码在读取文件时,怎么提示end-of-file during read,这个文件的内容如下所示(包括第一行的---线):
---------------------------------------------------------------
* Card 01 ************* Grid Generation
*         OX: Original Point's X Co-ordinate
*         OY: Original Point's Y Co-ordinate
*         IC: Number of cells in X-direction
*         JC: Number of cells in Y-direction
*         DX: Length of the cells in X-direction. Unit: m
*         DY: Length of teh cells in Y-direction. Unit: m
*
* Card 01   OX  OY   IC   JC   DX   DY
10.0   10.0  100  100   0.1   0.1
---------------------------------------------------------------
* Card 02 ************* Read bathymetry & surface elevation 
*         Depth_type : 1, Set the initial depth for all
*                      2, Read spatial bathymetry data from bathymetry.inp file
*
*         SurElv_Type: 1, Set the initial water surface elevation for all
*                      2, Read spatial water surface elevation from  IniSurElv.inp file
*
*         IniDepth   : value, when Depth_type =1. Unit:m 
*         IniSurElv  : value, when SurElv_Type=1. Unit:m
*         n          : value, Roughness height in the Manning Equations
*
* Card 02    Depth_type    SurElv_Type     IniDepth        IniSurElv        n
1  1  10  0  0
---------------------------------------------------------------
* Card 03 ************* Run time 
*         StartTime  : The beginning time,s
*         EndTime    : The end time,s
*         DT         : Time step,s
*
* Card 03    StartTime        EndTime        DT
0  10  0.05
---------------------------------------------------------------
* Card 04 ************* Computational Options 
*         Max_Iter   : Maximum number of iterations in Bi-CGSTAB
*         RSQM       : Maximum sum of residual squared in Bi-CGSTAB
*         PreCon_type: 1, Using the ILU(0) as the Preconditioner of Bi-CGSTAB
*                      2, No using the ILU(0) as the Preconditioner of Bi-CGSTAB
*
* Card 04    Max_Iter       RSQM         PreCon_type
200  0.000000001  1 
---------------------------------------------------------------
* Card 05 ************* Convective/Advective terms Options 
*        AdvecType   : 1, explicit centeral scheme
*                      2, explicit semi-momentum scheme
*                      3, explicit upwind scheme
*                      4, Euler-Lagrange method
*
* Card 05     AdvecType 
1  
---------------------------------------------------------------
* Card 06 ************* Output Options  
*        Step_Iterval: the frequency output the results
*         
* Card 06     Step_Iterval 
200