Subroutine NumDiffusionVelZ

• Developer: SUGIYAMA Ko-ichiro (sugiyama(at)gfd-dennou.org)
• Version: $Id: numdiffusionvelz.f90,v 1.3 2005/02/02 09:53:48 odakker Exp $
• Tag Name: $Name: $
• Change History:

Overview

数値粘性項を評価する.

Error Handling

Known Bugs

Note

粘性係数は arareprm 中で, \nu = \frac{\alpha \DeltaX}{\Delta t} で決める. \alpha < 1/8 である. CReSS マニュアル参照.

Future Plans

Dependency

use dc_trace, only: BeginSub, EndSub
use gridset,  only: DimXMin, DimXMax, DimZMin, DimZMax
use differentiate,  only: sf_dx_ff, ff_dx_sf, sf_dz_ss, ss_dz_sf
use arareprm, only: Cm, NuH, NuV

Input

real(8), intent(in) :: sf_VelZ_n(DimXMin:DimXMax, DimZMin:DimZMax) !速度 w

Output

real(8), intent(out) :: sf_NumDiffVelZ_n(DimXMin:DimXMax, DimZMin:DimZMax) !s数値粘性項