!=================================================================
!  一次元 CO2 大気放射対流平衡計算プログラム; 光路長の計算
!
!  1999/05/25  小高正嗣
!  1999/06/04  小高正嗣
!  2000/05/03  小高正嗣; 射出率スキームへ(圧力スケーリングを適用)
!=================================================================

      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 ( SFACT = 0.75D0 )

!-----------------------------------------------------------------
      OPL(0) = PSFC**( 1.0D0 + SFACT ) / 
     \     ( 1.0D0 + SFACT ) / (1000D+2)**SFACT / GRAV / 1.96D0 

      do iz = 1, NZ
         PRES = ( PRES0( iz ) + PRES0( iz - 1 ) ) * 0.5D0 
         OPL( iz ) = PRES**( 1.0D0 + SFACT ) / (1000D+2)**SFACT / 
     \        ( 1.0D0 + SFACT ) / GRAV / 1.96D0 
      end do

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

!      OPL(NZ+1) = OPL(NZ)
!      OPL(NZ+2) = OPL(NZ)

      return
      
      end
