速度の発散 ▽・u の計算
use dc_trace, only: BeginSub, EndSub use gridset, only: DimXMin, DimXMax, DimZMin, DimZMax use differentiate, only: ss_dx_fs, ss_dz_sf
real(8), intent(in) :: fs_Velx(DimXMin:DimXMax, DimZMin:DimZMax) !速度 u real(8), intent(in) :: sf_Velz(DimXMin:DimXMax, DimZMin:DimZMax) !速度 w
real(8), intent(out) :: ss_DivVel(DimXMin:DimXMax, DimZMin:DimZMax) !速度の収束