c*******************************************************
c init_diffu(n)
c initialiate w(i,j)
      subroutine init_diffu(n)
      implicit double precision(a-h, o-z)
      parameter(maxa=500,maxn=maxa*3)
	integer n3
       common/diffu_tensor/diffu(maxn,maxn),w(maxn,maxn),var(maxn)
      
      n3=n*3
      do i=1,n3
         do j=1,n3
            w(i,j)=0.d0
         enddo
      enddo

      return 
      end