!c Description: !c ▽・u の計算 !c !c Current Code Owner: !c sugiyama@gfd-dennou.org !c !c Copyright (C) SUGIYAMA Ko-ichiro, 2004, All rights reserved subroutine VelDiv(fs_VelX, sf_VelZ, ss_VelDiv) !--- モジュールの読み込み use gridset, only: DimXMin, DimXMax, DimZMin, DimZMax, ss_Bc use if_diff, only: ss_dx_fs, ss_dz_sf use if_boundary !--- 暗黙の型宣言禁止 implicit none !--- 入出力変数 real(8), intent(in) :: fs_Velx(DimXMin:DimXMax, DimZMin:DimZMax) real(8), intent(in) :: sf_Velz(DimXMin:DimXMax, DimZMin:DimZMax) real(8), intent(out) :: ss_VelDiv(DimXMin:DimXMax, DimZMin:DimZMax) !--- 速度 u ss_VelDiv = ss_dx_fs( fs_VelX ) + ss_dz_sf( sf_VelZ ) !--- 境界条件 call boundary(ss_Bc, ss_VelDiv) end subroutine VelDiv