!====================================================================
!  拡散方程式解法ルーチン(陰解法)
!  初期値計算
!
!  履歴  2000/01/07  小高正嗣
!====================================================================
      subroutine INIT
     I   ( ZGRID,
     O     TEMP, CDTURB,
     D     NZ, NZB )
!--------------------------------------------------------------------
      implicit real*8 ( a-h,o-z )
!--------------------------------------------------------------------

      dimension TEMP(-2:NZB), CDTURB(-2:NZB), ZGRID(-2:NZB)

!--------------------------------------------------------------------

      do iz = 0, NZ/2-2
         CDTURB(iz) = 1.0D0
      end do

      do iz = NZ/2-1, NZ
         CDTURB(iz) = 0.0D0
      end do

      do iz = 0, NZ/2-1
         TEMP(iz)  = 2.0D0 * ZGRID(iz)
      end do

      do iz = NZ/2, NZ-1
         TEMP(iz)  = 2.0D0 * ( 1.0D0 - ZGRID(iz) )
      end do

      TEMP(-2)  =0.0D0
      TEMP(-1)  =0.0D0
      TEMP(NZ)  =0.0D0
      TEMP(NZ+1)=0.0D0
      TEMP(NZ+2)=0.0D0

      call BOUND( CDTURB, NZ, NZB )

!--------------------------------------------------------------------
      return
      end
