星期四, 七月 22, 2010

fortran 数据格式转换的又一次教训

(E_M_max/E_M_P)**(k/N_E)这个实际上是错的!(其中k和N_E都是整型)
原因是(k/N_E)这个也会是整型,结果如果k 哎,又一次教训啊,还是ftnchek帮忙查出来的。
应该改为如下:
(E_M_max/E_M_P)**(k*1.D0/N_E)

没有评论: