1 subroutine tnanq (u, n, arrname) 2 3 include "common.h" 4 5 dimension u(nshg,n),rnan(2) 6 character(len=*) arrname 7 8 nnanq = 0 9 nlarge = 0 10 DO j = 1,n 11 DO i = 1,nshg 12 if (abs(u(i,j)).gt.1.0e10) nlarge=nlarge+1 13 if (u(i,j) .ne. u(i,j)) then 14 nnanq = nnanq + 1 15 u(i,j)=9.876543e21 16 endif 17 ENDDO 18 ENDDO 19 20 rnan(1)=nnanq 21 rnan(2)=nlarge 22 call sumgatN(rnan,2,summed,1) 23 24 if (summed.ge.1) then 25 close(1001) !Hack to close the varts files and flush the buffers. 26 call write_restart(myrank,9876543,nshg,n,u,u) 27 call error('tnanq ',arrname,nnanq) 28 endif 29 return 30 31 end 32