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 rnan(1)=nnanq 20 rnan(2)=nlarge 21 call sumgatN(rnan,2,summed,1) 22 if (summed.ge.1) then 23 call write_restart(myrank,9876543,nshg,n,u,u) 24 call error('tnanq ',arrname,nnanq) 25 endif 26 return 27 end 28