chaunceyyou 发表于 2018-5-18 09:50
非常谢谢你，我刚接触fortran没多久，书上的程序用的版本比较老，再继续时一定注重规范 ...

客气 祝学习愉快

吉大渣渣-固体 发表于 2018-5-17 23:04
坏在您的主程序中CALL PUTIN(NJ,N,NE,NM,NPJ,NPF,JN,X,Y,JE,JEAI,EAI,JPJ,PJ,JPF,PF)这个过程中子例行程 ...

非常谢谢你，我刚接触fortran没多久，书上的程序用的版本比较老，再继续时一定注重规范

chaunceyyou 发表于 2018-5-17 22:16
2
EXAMPLE----(5-2)
5,7,7,2,3,0

坏在您的主程序中CALL PUTIN(NJ,N,NE,NM,NPJ,NPF,JN,X,Y,JE,JEAI,EAI,JPJ,PJ,JPF,PF)这个过程中子例行程序PUTIN中I,I1，IE,IP,J，KK这些变量都没有初始化,以Intel visual Fortran编译器为例，这些变量不初始化直接拿来算编译器默认给的值是-858993460.有这些不确定性的隐患存在，导致后面EAI的值不稳定是正常的。

衷心给您的敬告：
1.尽管您的程序以Release模式可以得过且过运行出结果并正常退出，但是不初始化就拿来算这个隐患希望您心里有数。
2.不写implicit none,不重视初始化的程序都是在耍流氓。

吉大渣渣-固体 发表于 2018-5-17 21:54
请提供文件"PFSAP.IN".在没有此文件的情况下，直觉建议您：
1.主程序和子例行程序加上 implicit none;
2.详 ...

2
EXAMPLE----(5-2)
5,7,7,2,3,0
0,1,0,0.0,0.0,2,3,0,150.0,35.0
4,5,0,300.0,70.0,6,7,0,150.0,160.0
0,0,0,300.0,320.0
1,2,1,2,3,1,1,4,1,2,4,2,3,4,2
4,5,2,3,5,2
3.0E6,144.0,0.0
3.0E6,180.0,0.0
1,1610.0,3,3220.0,5,1610.0
0

请提供文件"PFSAP.IN".在没有此文件的情况下，直觉建议您：
1.主程序和子例行程序加上 implicit none;
2.详细检查所有变量、数组的声明和初始化