c**************************************************
c when diffusive matrix has negative overflow
	subroutine bd_stop(n)
         implicit double precision(a-h, o-z)
      parameter(maxa=500,maxn=maxa*3,bolz=8.31d-3/4.18,one3=1.d0/3.d0)
      integer nout, nn,ntime,n
      double precision gc(maxn),xold(maxa,3),gcnew(maxn),fcnew,
     +   y(maxn),rran(maxn),fold(maxa,3),vold(maxa,3),uold(maxa,3)
     +   , winv(maxn),ftor(maxa),phiold(maxa)
      common/diffu_tensor/diffu(maxn,maxn),w(maxn,maxn),varr(maxn)
     
	write(999,*) 'diffu'
            do i=1,n
	     do j=1,n
	     write(999,*) diffu(i,j),i,j
	     enddo
	    enddo
	write(999,*) 'w(i,j)'
	     do i=1,n
                do j=1,n
                write(999,*) w(i,j), i,j
        	enddo
	     enddo

	return
	end