Hi,
Finally, I concluded that the comment ( http://software.intel.com/sites/default/files/article/164389/fp-consistency-... ) which Ian suggested, is right, that is, "the random results is caused by variations in the starting address and alignment of the global stack".
So, in order to avoid this problems, I rewrote from "apply_dissipation.F77" to "apply_dissipation.F90" using Fortran 90 style. (I'm attaching the modified files.) Then, we can avoid random results at least my test environments.
If you don't have any objection, could you commit these change to repository?
Kentaro TAKAMI