C----------------------------------------------- C MONITOR.f C An INCLUDE code C Monitor transient results. C Currently at increment = INCR+1 C----------------------------------------------- C ------------------------------------ C Calculate difference between current C PHI and the steady-state value. C ------------------------------------ E = PHI(4) - (3*XORD(4) -5*YORD(4)) if(abs(E).le.0.0001) MNI = 0 C ---------------------------------- C Save values of PHI(4) for plotting C ---------------------------------- write(75,1001) PHI(4),E 1001 format(2e14.5)