!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