Subroutine EddyViscosity

• Developer: SUGIYAMA Ko-ichiro (sugiyama(at)gfd-dennou.org)
• Version: $Id: eddyviscosity.f90,v 1.2 2005/01/28 12:41:21 sugiyama Exp $
• Tag Name: $Name: $
• Change History:

Overview

乱流エネルギーを定常とした場合の渦拡散係数, 渦粘性係数を求める. この乱流パラメタリゼーションは, Mellor and Yamada (1974) の Level 1 Closure に対応するが, Level1 Closure に存在する \bar{\theta^2} は無視されている.

Error Handling

Known Bugs

Note

乱流エネルギーが負になったら, 零を入れる.

Future Plans

Dependency

use timeset,  only: DelTimeLong
use gridset,  only: DimXMin, DimXMax, DimZMin, DimZMax, DelX
use average,  only: ss_avr_ff, ss_avr_sf
use arareprm, only: Cm, Grav
use differentiate, only: ss_dx_fs, sf_dz_ss, ff_dx_sf, ff_dz_fs, ss_dz_sf
use basicset, only: ss_PotTempBasicZ
use bcset,    only: ss_BC

Input

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_PotTemp_n(DimXMin:DimXMax, DimZMin:DimZMax) !温位

Output

real(8), intent(out) :: &
     & ss_Km_n(DimXMin:DimXMax, DimZMin:DimZMax) !渦粘性係数
real(8), intent(out) :: &
     & ss_Kh_n(DimXMin:DimXMax, DimZMin:DimZMax) !渦拡散係数