!=================================================================
!  一次元 CO2 大気放射対流平衡計算プログラム; 光路長の計算
!
!  1999/05/25  小高正嗣
!  1999/06/04  小高正嗣
!  2000/05/03  小高正嗣; 射出率スキームへ(圧力スケーリングを適用)
!  2000/05/09  小高正嗣; 単位を g/cm^2 に変更, 拡散因子を考慮
!=================================================================

      subroutine SETOPL
     I   ( PRES0, PSFC  ,
     O     OPL  , 
     C     GRAV , 
     D     NZ   , NZB     )

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

      dimension 
     /   PRES0(-2:NZB), OPL(-2:NZB)

      parameter ( diffc = 1.67D0 )

!-----------------------------------------------------------------
      OPL(0) = PSFC / GRAV * diffc * 1.0D-1

      do iz = 1, NZ
         OPL( iz ) = ( PRES0( iz ) + PRES0( iz - 1 ) ) * 0.5D0 
     \        / GRAV * diffc * 1.0D-1
      end do

      OPL(-2)   = OPL(0)
      OPL(-1)   = OPL(0)
      OPL(NZ+1) = 0.0D0
      OPL(NZ+2) = 0.0D0

      return
      
      end
