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