1*59599516SKenneth E. Jansen subroutine tnanq (u, n, arrname) 2*59599516SKenneth E. Jansen 3*59599516SKenneth E. Jansen include "common.h" 4*59599516SKenneth E. Jansen 5*59599516SKenneth E. Jansen dimension u(nshg,n),rnan(2) 6*59599516SKenneth E. Jansen character*8 arrname 7*59599516SKenneth E. Jansen 8*59599516SKenneth E. Jansen nnanq = 0 9*59599516SKenneth E. Jansen nlarge = 0 10*59599516SKenneth E. Jansen DO j = 1,n 11*59599516SKenneth E. Jansen DO i = 1,nshg 12*59599516SKenneth E. Jansen if (abs(u(i,j)).gt.1.0e10) nlarge=nlarge+1 13*59599516SKenneth E. Jansen if (u(i,j) .ne. u(i,j)) then 14*59599516SKenneth E. Jansen nnanq = nnanq + 1 15*59599516SKenneth E. Jansen u(i,j)=9.876543e21 16*59599516SKenneth E. Jansen endif 17*59599516SKenneth E. Jansen ENDDO 18*59599516SKenneth E. Jansen ENDDO 19*59599516SKenneth E. Jansen rnan(1)=nnanq 20*59599516SKenneth E. Jansen rnan(2)=nlarge 21*59599516SKenneth E. Jansen call sumgatN(rnan,2,summed,1) 22*59599516SKenneth E. Jansen if (summed.ge.1) then 23*59599516SKenneth E. Jansen call write_restart(myrank,9876543,nshg,n,u,u) 24*59599516SKenneth E. Jansen call error('tnanq ',arrname,nnanq) 25*59599516SKenneth E. Jansen endif 26*59599516SKenneth E. Jansen return 27*59599516SKenneth E. Jansen end 28