Hi,
I compiled Your code under Linux and found some bugs.
In particular,
1. In "[login to view URL]" variable EPS is not defined.
2. In "[login to view URL]" in the loop 'DO TI=2,TI_MAX-3,1' array bounds are exceeded
for PHI(X,TI-2), GAMMA(X,TI-2) /TI-2<1/, and BETA(X,TI+1) /TI+1>NN/.
Would be glad to further help You.
Thanks.
Bogdan Yavorskyy.