Subroutine EddyViscosity

Overview

2.5 次の乱流クロージャーを利用し, 次の時刻の渦粘性係数と渦拡散係数を求めるサブルーチン.

Error Handling

Known Bugs

Note

Future Plans

Dependency

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

Input

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)

Output

real(8), intent(out) :: ss_Km_out(DimXMin:DimXMax, DimZMin:DimZMax)
real(8), intent(out) :: ss_Kh_out(DimXMin:DimXMax, DimZMin:DimZMax)