!=================================================================
! 一次元灰色大気放射対流平衡計算プログラム; 加熱率の計算
! 
!
!  1999/05/31  小高正嗣
!=================================================================

      subroutine CLQRAD
     I   ( FXRADN, FDZ, FDZM, DENS0,
     O     QRAD  ,
     C     CPDRY , DX ,
     D     NZ    , NZB              )

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

      dimension
     \   DENS0(-2:NZB) , FDZ(-2:NZB), FDZM(-2:NZB),
     \   FXRADN(-2:NZB),QRAD(-2:NZB)

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

      do iz = 0, NZ
         QRAD( iz ) = - ( FXRADN( iz + 1 ) - FXRADN ( iz ) ) /
     \        ( CPDRY * DENS0( iz ) * FDZ( iz ) * DX ) * 86400.0D0
      end do

      QRAD( -1 ) = QRAD ( 0 )
      QRAD( -2 ) = QRAD ( 0 )
      QRAD( NZ+1 ) = QRAD ( NZ )
      QRAD( NZ+2 ) = QRAD ( NZ )

      end 
