subroutine tnanq (u, n, arrname)

      include "common.h"

      dimension   u(nshg,n),rnan(2)
      character(len=*) arrname

      nnanq = 0
      nlarge = 0
      DO j = 1,n
	DO i = 1,nshg
	  if (abs(u(i,j)).gt.1.0e10)  nlarge=nlarge+1
	  if (u(i,j) .ne. u(i,j)) then
	     nnanq = nnanq + 1
	     u(i,j)=9.876543e21
	  endif 
	ENDDO
      ENDDO
  	rnan(1)=nnanq
	rnan(2)=nlarge
	call sumgatN(rnan,2,summed,1)
      if (summed.ge.1) then
	 call write_restart(myrank,9876543,nshg,n,u,u)
         call error('tnanq   ',arrname,nnanq)
      endif
      return
      end