!=================================================================
!  地表面温度計算プログラム
!
!  1999/07/21  小高正嗣
!  1999/10/26  小高正嗣
!=================================================================

      subroutine CLTSFC
     I   ( FHEAT , FXRADD, FXSOL,
     O     TSFC  ,
     C     SIGMA ,
     D     NX    , NXB   , 
     D     NZ    , NZB            )

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

      dimension
     \   FHEAT ( 2 )

      dimension
     \   FXRADD( -2:NZB )

      dimension
     \   TSFC ( -2:NXB )

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

      do ix = -2, NX+2
         FXRADU_SFC = FXSOL + FXRADD( 0 ) - FHEAT( 1 ) 
         TSFC( ix ) = ( FXRADU_SFC / SIGMA )**0.25D0
      end do

      return

      end

