本帖最后由 楚香饭 于 2015-1-20 22:49 编辑
做人呢,一定要长记性。
代码写错不丢人,丢人的是,同样的错误多次犯。
http://bbs.fcode.cn/thread-53-1-1.html
这是你上次的帖子吧??
上次你错在,Common没有对齐。这次你犯了同样的错。
主程序中的 common:
[Fortran] 纯文本查看 复制代码 common a,b,rr,epson
common ua0,uw0
common m1kw,m2w,m1ka,m2a,m1ks,m2s
common Q1,Q2,q12,q21
common Cw,Cvw,Ca,Cva
子程序 SET 中的 Common:
[Fortran] 纯文本查看 复制代码 common a,b,r,epson
common n0,Sl0,Sg0,ua0,uw0,rw,kw,ka,g,Temp,M,Uatm !// 这里多了很多
common m1kw,m2w,m1ka,m2a,m1ks,m2s
common Q1,Q2,q12,q21
这就导致 m2w = 0,在
Cw=(m1kw-m2w)/m2w 时,除数为 0,开始出错。后面各种浮点数错误。NaN 等等等等。
叫我说你什么好?
请一定要理解:COMMON不是根据变量名称进行对应,而是根据变量的顺序进行对应!!这就是为什么不提倡使用 COMMON 的原因。
|