159599516SKenneth E. Jansen subroutine tnanq (u, n, arrname) 259599516SKenneth E. Jansen 359599516SKenneth E. Jansen include "common.h" 459599516SKenneth E. Jansen 559599516SKenneth E. Jansen dimension u(nshg,n),rnan(2) 69798c0ddSCameron Smith character(len=*) arrname 759599516SKenneth E. Jansen 859599516SKenneth E. Jansen nnanq = 0 959599516SKenneth E. Jansen nlarge = 0 1059599516SKenneth E. Jansen DO j = 1,n 1159599516SKenneth E. Jansen DO i = 1,nshg 1259599516SKenneth E. Jansen if (abs(u(i,j)).gt.1.0e10) nlarge=nlarge+1 1359599516SKenneth E. Jansen if (u(i,j) .ne. u(i,j)) then 1459599516SKenneth E. Jansen nnanq = nnanq + 1 1559599516SKenneth E. Jansen u(i,j)=9.876543e21 1659599516SKenneth E. Jansen endif 1759599516SKenneth E. Jansen ENDDO 1859599516SKenneth E. Jansen ENDDO 19*513954efSKenneth E. Jansen 2059599516SKenneth E. Jansen rnan(1)=nnanq 2159599516SKenneth E. Jansen rnan(2)=nlarge 2259599516SKenneth E. Jansen call sumgatN(rnan,2,summed,1) 23*513954efSKenneth E. Jansen 2459599516SKenneth E. Jansen if (summed.ge.1) then 25*513954efSKenneth E. Jansen close(1001) !Hack to close the varts files and flush the buffers. 2659599516SKenneth E. Jansen call write_restart(myrank,9876543,nshg,n,u,u) 2759599516SKenneth E. Jansen call error('tnanq ',arrname,nnanq) 2859599516SKenneth E. Jansen endif 2959599516SKenneth E. Jansen return 30*513954efSKenneth E. Jansen 3159599516SKenneth E. Jansen end 32