Question

fortran有限元编程的若干问题

105 FORMAT(2X,'NN,NE,ND.NFIX,E,ANU,T,GM,NTYPE',1/4I4,E10.4,2F7.3,E10.4,I3) READ(5,*) (LOC(I,1),LOC(I,2),LOC(I,3),I=1,NE) READ(5,*) (CX(J),CY(J),J=1.NN) READ(5,*) (IFIX(K),K=1,NFIX)C- - - - - - - - - - - - - - - - - - - - - - DO 10 I=1,ND10 F(I)=0.0 F(2)=-1.0C- - - - - - - - - - - - - - - - - - - - - - CALL CST(LOC,CX,CY,IFIX,F,GK,STRES,BAK) WRITE(6,110)110 FORMAT(/4X,'NODE',5X,'X-DISP',8X,'Y-DISP') WRITE(6,115) (I,F(2*I-1),F(2*I),I=1,NN)115 FORMAT(2X,I5,2E15.6)C- - - - - - - - - - - - - - - - - - - - - - WRITE(6,120)120 FORMAT(/4X,'ELEMENT',4X,'X-STR',8X,'Y-STR',8X,'XY-STR') WRITE(6,125) (I,(STRES(I,J),J=1,3),I=1,NE)125 FORMAT(2X,I4,3E15.6) STOP END
我调试这段fortran代码的时候，发现很多的错误，我是第一次使用fortran，所以不知道的地方很多。
我贴上来的图

是我运行的错误的地方，请不吝赐教





Answer 1

看起来你的代码应该是固定格式的代码。

而你使用了 source1.f90 这个文件名，大多数编译器会把他视为自由格式。

追问

其中C- - - - - - -这种肉眼可以认为是功能分割，但是代码可以这样直接敲进去吗？

追答

C---------- 表示注释行。这是固定格式的注释方法，如果你使用 .f90 扩展名，表示是自由格式。

自由格式的注释行是不一样的，这样基本的常识如果你还欠缺，建议你买一本教材系统学习。