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

      subroutine SETCSG
     I   ( DAY   ,
     O     ALBEDO, DIFG, CPGRND, DENSG, TURBG, SKDPS )

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

      ALBEDO = 0.25D0           ! 地表面アルベド
      TI     = 272.0D0          ! 熱慣性(Wm^-2sec^0.5K^-1)
      
      CPGRND = 588.0D0          ! 比熱(Jkg^-1K^-1)
      DENSG  = 1650.0D0         ! 土壌密度(kgm^-3)

      DIFG = TI**2.0 / CPGRND / DENSG
      TURBG = DIFG / CPGRND / DENSG
      
      SKDPS = sqrt( TURBG * DAY )

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

      return
      end
      
      
