HE-VI 法を用いたエクスナー関数の計算.
離散化する際, 上下境界条件として鉛直速度が零を与えている.
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
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) !無次元圧力
real(8), intent(out) :: ss_Exner_a(DimXMin:DimXMax, DimZMin:DimZMax) !無次元圧力 (1 ステップ後)