1*59599516SKenneth E. Jansen subroutine tnanqe (u, n, arrname) 2*59599516SKenneth E. Jansen 3*59599516SKenneth E. Jansen include "common.h" 4*59599516SKenneth E. Jansen 5*59599516SKenneth E. Jansen dimension u(npro,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,npro 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 write(*,*) myrank, i,j 15*59599516SKenneth E. Jansen nnanq = nnanq + 1 16*59599516SKenneth E. Jansen u(i,j)=9.876543e21 17*59599516SKenneth E. Jansen endif 18*59599516SKenneth E. Jansen ENDDO 19*59599516SKenneth E. Jansen ENDDO 20*59599516SKenneth E. Jansen rnan(1)=nnanq 21*59599516SKenneth E. Jansen rnan(2)=nlarge 22*59599516SKenneth E. Jansen call sumgatN(rnan,2,summed,1) 23*59599516SKenneth E. Jansen if (summed.ge.1) then 24*59599516SKenneth E. Jansen do i=1,npro 25*59599516SKenneth E. Jansen write(8+myrank,245) (u(i,j), j=1,n) 26*59599516SKenneth E. Jansen enddo 27*59599516SKenneth E. Jansen call error('tnanqe ',arrname,nnanq) 28*59599516SKenneth E. Jansen endif 29*59599516SKenneth E. Jansen245 format(10(e14.7,2x)) 30*59599516SKenneth E. Jansen return 31*59599516SKenneth E. Jansen end 32