!================================================================= ! 一次元 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