オイラー法を用いた X 方向の速度の計算. deepconv/arare では時間積分として He-VI 法を利用しているので, X 方向の速度は陽解法で解く.
use dc_trace, only: BeginSub, EndSub use gridset, only: DimXMin, DimXMax, DimZMin, DimZMax use basicset, only: ss_CpBasicZ, ss_PotTempBasicZ use bcset, only: fs_BC use arareset, only: alpha use average, only: fs_avr_ss use differentiate_center2, only: fs_dx_ss, ss_dx_fs, ss_dz_sf
real(8), intent(in) :: DelTime real(8), intent(in) :: fs_VelX_in(DimXMin:DimXMax, DimZMin:DimZMax) real(8), intent(in) :: ss_Exner(DimXMin:DimXMax, DimZMin:DimZMax) real(8), intent(in) :: fs_VelX(DimXMin:DimXMax, DimZMin:DimZMax) real(8), intent(in) :: sf_VelZ(DimXMin:DimXMax, DimZMin:DimZMax) real(8), intent(in) :: fs_AdvDiffX(DimXMin:DimXMax, DimZMin:DimZMax)
real(8), intent(out) :: fs_VelX_out(DimXMin:DimXMax, DimZMin:DimZMax)