Z 方向の運動方程式の粘性拡散項を評価する. 音波に関連しない項として扱うため長い時間ステップで評価.
use dc_trace, only: BeginSub, EndSub use gridset, only: DimXMin, DimXMax, DimZMin, DimZMax, DelXZ use differentiate, only: sf_dx_ff, ff_dz_fs, ff_dx_sf, & & sf_dz_ss, ss_dz_sf use average, only: ff_avr_ss use arareprm, only: Cm
real(8), intent(in) :: fs_VelX_n(DimXMin:DimXMax, DimZMin:DimZMax) !速度 u real(8), intent(in) :: sf_VelZ_n(DimXMin:DimXMax, DimZMin:DimZMax) !速度 w real(8), intent(in) :: ss_Km_n(DimXMin:DimXMax, DimZMin:DimZMax) !渦粘性係数
real(8), intent(out) :: sf_TurbVelZ_n(DimXMin:DimXMax, DimZMin:DimZMax) !粘性拡散項