X 方向の運動方程式の粘性拡散項を評価する. 音波に関連しない項として扱うために長い時間ステップで評価
use dc_trace, only: BeginSub, EndSub use gridset, only: DimXMin, DimXMax, DimZMin, DimZMax, DelXZ use differentiate, only: fs_dx_ss, ss_dx_fs, fs_dz_ff, ff_dx_sf, & & fs_dz_ff, ff_dz_fs, fs_dx_ss 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) :: fs_TurbVelX_n(DimXMin:DimXMax, DimZMin:DimZMax) !粘性拡散項