原因判明

MPI部分が原因ではなく、それ以外のところが原因だった。
今自分が解いている問題に依存した事であることが判明。
具体的には、円柱座標系で解いているのだが、軸（z軸）上の
値を周りの点の平均として計算している。それにより
軸上でも軸方向成分以外のベクトル成分が0でない値を
持てるようにしている。その部分が悪さをしている模様
（今日はそこまでしか分からなかったが、大きな進歩）。
明日にでもゆっくり考えよう。

codeのコアな部分（並列化部分を含む）のバグじゃなかった
のには一安心。しかし逆にそういう細かい所のバグで、
ずーっと悩んでいたのは悔しい。Fortran77→Fortran90の
書き替えを決心した不可解な振る舞いもこれに関係ある
のか？まあどちらにしろFortran90に少し詳しくなったし、
よしとしよう。