!=====================================================================
!  鉛直 1 次元地表面熱収支計算プログラム; 大気上端の日射の計算
!
!  履歴  2000/09/20  小高正嗣; 作成開始
!        2000/09/21  小高正嗣
!        2000/09/25  小高正嗣
!=====================================================================
      subroutine CLFXSOL
     I   ( S00  , ALS   , YLAT , TIME, TIMES,
     O     S0   , COSZ  ,
     C     EXCTL, OBLIC , OMGA0, DAY , PI   )
!---------------------------------------------------------------------
      implicit real*8 (a-h,o-z)
!---------------------------------------------------------------------

      TIMEL = ( TIME + TIMES ) * 2.0 * PI / DAY

      COSZ = MAX ( 0.0D0, 
     \          ( sin( YLAT ) * sin( OBLIC ) * sin( ALS ) + 
     \            cos( TIMEL - PI ) * cos( YLAT ) *
     \            sqrt( 1.0D0 - (sin( OBLIC )*sin( ALS ))**2.0D0 ) ) )

      S0 = S00 * ( ( 1.0D0 + EXCTL * cos( ALS + OMGA0 ) )
     \     / ( 1.0D0 - EXCTL**2.0D0 ) )**2.0D0 

      return

      end
