2.5 次の乱流クロージャーを利用し, 次の時刻の渦粘性係数と渦拡散係数を求めるサブルーチン.
use dc_trace, only: BeginSub, EndSub use gridset, only: DimXMin, DimXMax, DimZMin, DimZMax, DelXZ use bcset, only: ss_BC use physset, only: Grav use basicset, only: ss_PotTempBasicZ use arareset, only: Cm use average, only: ss_avr_fs, ss_avr_ff, ss_avr_sf use differentiate_center4, only: fs_dx_ss, ss_dx_fs, ff_dx_sf, & & ss_dz_sf, sf_dz_ss, ff_dz_fs
real(8), intent(in) :: DelTime real(8), intent(in) :: ss_Km_in(DimXMin:DimXMax, DimZMin:DimZMax) real(8), intent(in) :: fs_VelX_adv(DimXMin:DimXMax, DimZMin:DimZMax) real(8), intent(in) :: sf_VelZ_adv(DimXMin:DimXMax, DimZMin:DimZMax) real(8), intent(in) :: ss_Km_adv(DimXMin:DimXMax, DimZMin:DimZMax) real(8), intent(in) :: fs_VelX_dif(DimXMin:DimXMax, DimZMin:DimZMax) real(8), intent(in) :: sf_VelZ_dif(DimXMin:DimXMax, DimZMin:DimZMax) real(8), intent(in) :: ss_PotTemp_dif(DimXMin:DimXMax, DimZMin:DimZMax) real(8), intent(in) :: ss_Km_dif(DimXMin:DimXMax, DimZMin:DimZMax)
real(8), intent(out) :: ss_Km_out(DimXMin:DimXMax, DimZMin:DimZMax) real(8), intent(out) :: ss_Kh_out(DimXMin:DimXMax, DimZMin:DimZMax)