!=====================================================================
!  鉛直 1 次元地表面熱収支計算プログラム; 物理定数の設定
!
!  履歴  2000/09/17  小高正嗣; 作成開始
!
!=====================================================================

      subroutine SETGGD
     I   ( SKDPS, GGRID,
     O     GDZ  , GDZM,
     D     NZG  , NZGB        )

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

      dimension
     \     GGRID(-1:NZGB),      ! 地面鉛直座標(無次元)
     \     GDZ  (0:NZGB),       ! 鉛直格子間隔
     \     GDZM (0:NZG)         ! 鉛直格子間隔(半整数格子上)

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

      do iz = 0, NZGB
         GDZ(iz) = ( GGRID(iz) - GGRID(iz-1) ) * SKDPS
      end do

      do iz = 0, NZG
         GDZM(iz) = 0.5 * ( GDZ(iz+1) + GDZ(iz) )
      end do

      return
      end