|
在我这里,并没有出现你截图中的 end-of-file 错误。而是浮点数错误。
究其原因,在于
womiga0 = abs(acos(-tan(delta)*tan(fai)))
此处的 tan(delta)*tan(fai) 无法保证在 [ -1 , 1 ] 区间内。当 delta 大于 -1.5049 时,tan(delta)*tan(fai) 的结果为 1.0288
至此,无法计算 acos
你可以在此句前面加上 write(*,*) delta,tan(delta),tan(delta)*tan(fai) 以便观察。
最后,对于这种运行时错误,建议你学会自行debug,可参考此文:http://debug.w.fcode.cn
|
|