接下来简单介绍一下程序：

1）        main.f90 是主程序，主要是调用其他程序

2）        主要的计算程序在thermodynamics.f90, tool.f90 and solver.f90

3）        Mpi.f90是communication between processes

4）        Ghostnode.f90是边界条件

5）        Input.f90主要是从input file里读取数据，output.f90就是输出计算结果，inti.f90主要是allocate数组

6）        Global.f90是定义的全局变量

现在我怀疑是不是MPI communication的时候出现问题，因为NaN的值出现在第一个processor的边界上，而且这个程序在自己的电脑（windows）下串行运行是OK的但是当我在超级计算机上必行运算时，程序就不行，然后我用DDT debug，发现在这里出现问题，如下图，然后我进一步check，发现变量xp有四个值是NaN, xp被计算用subroutineconvert_composition(ix,iy)。

2，By the way, 如果要运行这个程序的话需要:

Intel fotransompiler, intelmkl and open MPI.

3.source code 在src  file里， input file 里是需要读入的数据，结果会写在results file里。Makefile是用来compile code的。run.sh是PBS file.

程序在这里：http://pan.baidu.com/share/link?shareid=433689&uk=3305225013

十分感谢！

备注：问题也在这儿  http://emuch.net/bbs/viewthread.php?tid=5846918