程序运行时出现 run-time error F6501错误

木头 · 发表于 2017-2-16 15:10:35

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

001C----------------------------------------------------------------------
002C     形成单元节点的岩土强度指标FA(I),C(I),单元联系数据,输入节点应力(平面问题）
003        IMPLICIT REAL*8 (A-H,O-Z)
004      COMMON FA(10000),C(10000),NUMNP,NENM,NUM
005        DIMENSION NODMJ(8,10000),NENMJ(16,10000),FAM(100),CM(100)
006      DIMENSION STRESS(4),NSTR(10000),NPAR(20,50)
007        CHARACTER*30 INFNAME
008C******************************************************************
009        INFNAME='STRESS-CHANGE-4.OUT'
010C******************************************************************
011        DO I=1,30
012           IF(INFNAME(I:I).EQ.'.') GOTO 10
013        END DO
01410    INFNAME(I+1:I+4)='CON '
015      OPEN(UNIT=11,FILE=INFNAME)
016      INFNAME(I+1:I+4)='NSTR'
017      OPEN(UNIT=13,FILE=INFNAME)
018C   *********************************************************************
019C    单元平均节点应力由有限元输出获得.
020C    文件格式:NODE(I),SGXX,SGYY,SGZZ,SGYZ,节点编号应依此增大.
021 
022      INFNAME(I+1:I+4)='STR '
023      OPEN(23,FILE=INFNAME)
024C   *********************************************************************
025 
026      INFNAME(I+1:I+4)='ENM '
027      OPEN(22,FILE=INFNAME,FORM='UNFORMATTED',ACCESS='DIRECT',RECL=120)
028 
029C    输入单元数据,节点数NUMNP，单元组数NUM,单元数NUMN
030        READ(11,*) NUMNP,NUM,NENM
031        DO I=1,NUM
032           READ(11,*) (NPAR(J,I),J=1,10)
033           READ(11,*) (NPAR(J,I),J=11,20)
034        END DO
035C******************************************************************
036C    直接输入单元组的岩土强度指标,FAM(I),CM(I).
037      FAM(1)=34.0
038        CM(1)=35000.0
039      FAM(2)=34.0
040        CM(2)=35000.0
041      FAM(3)=34.0
042        CM(3)=35000.0
043      FAM(4)=25.0
044        CM(4)=25000.0
045C******************************************************************
046 
047C    形成单元联系单元 NENMJ(8,NUMN),输出到*.NSTR文件中
048      IELN=8
049        DO I=1,NENM
050           READ(22,REC=I) IDEGEN,IELD,ISOCOR,(NODMJ(K,I),K=1,IELN)
051        END DO
052        CLOSE(22)
053 
054        DO I=1,NENM
055           DO J=1,16
056              NENMJ(J,I)=0
057           END DO
058           MM=0
059           DO J=1,NENM
060              IF(J.NE.I) THEN
061                      CALL NENMJD2(NODMJ,I,J,MM,NENMJ)
062              END IF
063           END DO
064           WRITE(13,'(16I6)') (NENMJ(K,I),K=1,16)
065        END DO
066 
067C    形成单元节点应力 STRESSJ(4,NUNMP), 岩土强度指标FA(NUNMP),C(NUNMP),输出到*.NSTR文件中.
068        DO I=1,NUMNP
069           NSTR(I)=0
070           FA(I)=0.0
071           C(I)=0
072        END DO
073 
074        DO I=1,NENM
075           NM=0
076C     第MM组单元的单元数NPAR(2,MM)
077           DO MM=1,NUM
078              IF(I.GT.NM.AND.I.LE.NM+NPAR(2,MM)) GOTO 200
079                     NM=NM+NPAR(2,MM)
080           END DO
081C     求单元各节点的岩土强度指标和.
082200      DO J=1,8
083              IF(NODMJ(J,I).NE.0) THEN
084                FA(NODMJ(J,I))=FA(NODMJ(J,I))+FAM(MM)
085                C(NODMJ(J,I))=C(NODMJ(J,I))+CM(MM)
086                NSTR(NODMJ(J,I))=NSTR(NODMJ(J,I))+1
087              END IF
088           END DO
089        END DO
090        ZZ=ATAN(1.0)/45.0
091        DO I=1,NUMNP
092C     求单元各节点的岩土强度指标平均值.  NSTR(I)为节点I的联系单元数
093           FA(I)=FA(I)/NSTR(I)
094           C(I)=C(I)/NSTR(I)
095 
096           READ(23,*) NODE,(STRESS(J),J=1,4)
097C    形成单元节点应力 STRESSJ(4), 岩土强度指标FA(I),C(I),输出到*.NSTR文件中.
098           WRITE(13,2000) (-STRESS(J),J=1,4),FA(I)*ZZ,C(I)
099        END DO
100        CLOSE(23)
101        CLOSE(13)
102 
1031000  FORMAT(1X,'INPUT FILE NAME :',$)
1041001  FORMAT(A12)
1052000  FORMAT(6E14.6)
1062100  FORMAT(16I5)
107        END
108 
109C----------------------------------------------------------------------
110C   形成单元联系标记NENMJ(20,I), 单元M1和节点号NODMJ(4,M1),联系单元M2和节点号NODMJ(4,M2)
111        SUBROUTINE NENMJD2(NODMJ,M1,M2,MM,NENMJ)
112        IMPLICIT REAL*8 (A-H,O-Z)
113        DIMENSION NENMJ(16,10000),NODMJ(8,10000)
114        N=1
115        DO I=1,4
116           DO J=1,4
117                     IF(NODMJ(I,M1).EQ.NODMJ(J,M2)) GOTO 100
118               END DO
119        END DO
120        N=0
121100   IF(N.GE.1) THEN
122        MM=MM+1       
123          NENMJ(MM,M1)=M2
124           END IF
125        RETURN
126        END

木头 · 发表于 2017-2-16 15:11:57

本人刚刚接触Fortran，导师给了一些编好的程序，让我计算，计算时出现这样的错误，不知道该怎么解决，想请教各位大神一下

fcode · 发表于 2017-2-16 15:15:32

这是读取文件时遇到了文件的结束。例如，代码要求数据有3行数据，而实际输入文件只有2行。在某些时候，甚至输入文件根本不存在或是空白文件。此外，数据文件缺失某些行某些列也是较大的可能。当然，也可能是循环读取文件时没控制住。比如想读100行，结果写错了，读了1000行。有时候，可能是文件路径错误，导致实际读取了另一个空白的文件。或者文件扩展名被隐藏了，比如读取 a.txt，但实际的文件名却是 a.txt.txt