!c Description:
!c   Z 方向の運動方程式において長い時間ステップで評価する項の計算
!c
!c Current Code Owner:
!c   sugiyama@gfd-dennou.org
!c
!c Copyright (C) SUGIYAMA Ko-ichiro, 2004, All rights reserved

subroutine DiffuseZ(fs_VelX_B, sf_VelZ_B, sf_DiffuseZ_B)

  !--- モジュールの読み込み
  use gridset,  only: DimXMin, DimXMax, DimZMin, DimZMax, sf_Bc
  use if_avr,   only: sf_avr_fs, sf_avr_ss
  use if_diff,  only: sf_dx_sf, sf_dz_sf
  use if_boundary

  !--- 暗黙の型宣言禁止
  implicit none

  !--- 入出力変数
  real(8), intent(in)     :: fs_VelX_B(DimXMin:DimXMax, DimZMin:DimZMax)
  real(8), intent(in)     :: sf_VelZ_B(DimXMin:DimXMax, DimZMin:DimZMax)
  real(8), intent(out)    :: sf_DiffuseZ_B(DimXMin:DimXMax, DimZMin:DimZMax)

  !--- 拡散項 (現在は 0.0 とする)
  sf_DiffuseZ_B = 0.0d0

  !--- 境界条件
  call boundary(sf_Bc, sf_DiffuseZ_B)

end subroutine DiffuseZ