Subroutine Exner_HeVi

Overview

HE-VI 法を用いたエクスナー関数の計算.

Error Handling

Known Bugs

Note

離散化する際, 上下境界条件として鉛直速度が零を与えている.

Future Plans

Dependency

use gridset,  only: DimXMin, DimXMax, DimZMin, DimZMax, &
     &              RegXMax, RegXMin, RegZMin, RegZMax, DelZ
use timeset,  only: DelTimeShort
use basicset, only: ss_VelSoundBasicZ, ss_CpBasicZ, &
     &              ss_DensBasicZ, ss_PotTempBasicZ
use arareprm, only: alpha, beta
use linlib,   only: LinSolv
use average,  only: sf_avr_ss
use differentiate, only: sf_dz_ss, ss_dz_sf, ss_dx_fs

Input

real(8), intent(in)      :: ss_DivVel_n(DimXMin:DimXMax, DimZMin:DimZMax)  !速度の収束
real(8), intent(in)      :: fs_VelX_a(DimXMin:DimXMax, DimZMin:DimZMax)    !速度 u
real(8), intent(in)      :: sf_VelZ_n(DimXMin:DimXMax, DimZMin:DimZMax)    !速度 w
real(8), intent(in)      :: sf_Fz_n(DimXMin:DimXMax, DimZMin:DimZMax)      !Z 方向の外力項
real(8), intent(in)      :: ss_Exner_n(DimXMin:DimXMax, DimZMin:DimZMax)   !無次元圧力

Output

real(8), intent(out)     :: ss_Exner_a(DimXMin:DimXMax, DimZMin:DimZMax)   !無次元圧力 (1 ステップ後)