!=================================================================
! 一次元灰色大気放射対流平衡計算プログラム; 
! 黒体輻射の計算
!
!  1999/05/27  小高正嗣
!  1999/10/18  小高正嗣
!=================================================================

      subroutine CLBBRAD
     I   ( TEMP0,
     O     BBRAD,
     D     NZ   , NZB )     

!-----------------------------------------------------------------      
      implicit real*8 ( a-h,o-z )
!-----------------------------------------------------------------      

      dimension
     \   TEMP0(-2:NZB), BBRAD(-2:NZB)

      parameter( NU = 667.0D0 )

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

      do iz = 0, NZB
         TEMPM       = ( TEMP0( iz ) + TEMP0( iz - 1 ) ) * 0.5D0
         
         BBRAD( iz ) = ( 1.1911D-8 * NU**3 ) /
     \        ( EXP( 1.439D0 * NU / TEMPM ) - 1.0D0 )

      end do
      
      BBRAD(-1) = BBRAD(0)
      BBRAD(-2) = BBRAD(0)
      
      return

      end
