[交流] 【求助】FORTRAN编程添加原子 已有1人参与

【求助】在MS导出的文件中，按照规则添加原子 上个帖子我说了，可以把晶体结构保存为PDB格式，比如说，我要在一个沸石结构上表面添加氨基，那么像这样的程序为什么运行不出来？     OPEN(10,FILE="MFI.PDB",POSITION="ASIS",ACTION="READWRITE",STATUS="OLD"         OPEN(60,FILE="MFINEW.PDB",STATUS="REPLACE"         READ(10,"(A80)"ATOMNUM,XXX,YYY,ZZZ C     %计算新加入的第一个原子与原来骨架距离%              DO  XXX=1,7776,1          DO  YYY=1,7776,1           DO  ZZZ=1,7776,1            CALL RAND(XNEW)            CALL RAND(YNEW)            CALL RAND(ZNEW) 20           IF SQRT(((XNEW-XXX)**2+(YNEW-YYY)**2+(ZNEW-ZZZ)**2).LT.0.05)             THEN                 WRITE(20,“(4A7)"ATOMNUM,XXNEW,YYNEW,ZZNEW                 ELSE                           END IF           END DO          END DO         END DO C      %计算新加入的第二个原子与原来骨架距离%         DO  XXX=1,7777,1          DO  YYY=1,7777,1           DO  ZZZ=1,7777,1            CALL RAND(XXNEW)            CALL RAND(YYNEW)            CALL RAND(ZZNEW) 30           IF SQRT(((XXNEW-XXX)**2+(YYNEW-YYY)**2+(ZZNEW-ZZZ)**2).LT.0.05)             THEN                 WRITE(30,“(4A7)"ATOMNUM,XXNEW,YYNEW,ZZNEW                 ELSE                                    END IF           END DO          END DO         END DO C     %计算新加入的第三个原子与原来骨架距离%         DO  XXX=1,7778,1          DO  YYY=1,7778,1           DO  ZZZ=1,7778,1            CALL RAND(XXXNEW)            CALL RAND(YYYNEW)            CALL RAND(ZZZNEW) 40          IF SQRT(((XXXNEW-XXX)**2+(YYYNEW-YYY)**2+(ZZZNEW-ZZZ)**2).LT.0.0           5 )THEN           WRITE(40,“(4A7)"ATOMNUM,XNEW,YNEW,ZNEW           ELSE           END IF           END DO          END DO         END DO       WRITE(60,“(10A7)"ATOMNUM,XNEW,YNEW,ZNEW,XXNEW,YYNEW,ZZNEW,XXXNE         W,YYNEW,ZZNEW,XXXNEW,YYYNEW,ZZZNEW       END                SUBROUTINE RAND(M)         DOUBLE PRECISION S,U,V,R         S=65536.0         U=2053.0         V=13849.0         M=R/S         R=R-M*S         R=U*R+V         M=R/S         R=R-M*S         NRND1=R/S         RETURN         END 我的思想是读入文件，然后再骨架中，引用随机数，随即的插入原子，只有当心加入的原子与原来原子的距离小于0.05的时候，才认为成键了，不知道这样的做法对不对？在FORTRAN编译器中，读入文件，对文件操作需要把文件和程序放在一个文件夹里面。

回复此楼