Class | NumDiffusion |
In: |
util/numdiffusion.f90
|
数値拡散項の計算モジュール 現在は 2 次精度中心差分を利用
Subroutine : |
NumDiffusion モジュールの初期化ルーチン
subroutine NumDiffusion_init() ! ! NumDiffusion モジュールの初期化ルーチン ! !暗黙の型宣言禁止 implicit none ! CReSS マニュアルでは, 2 次精度中心差分の場合, Alpha < 1/8 くらいが適当と述べている. ! deepconv では NuH として 500 以上の値が入っていたそうだ, NuH = Alpha * ( DelX ** 2.0d0 ) / DelTimeLong NuV = Alpha * ( DelZ ** 2.0d0 ) / DelTimeLong !確認 call MessageNotify( "M", "NumDiffusion_init", "NuH = %f", d=(/NuH/) ) call MessageNotify( "M", "NumDiffusion_init", "NuV = %f", d=(/NuV/) ) end subroutine NumDiffusion_init
Function : | |||
pz_NumDiffVelX(DimXMin:DimXMax, DimZMin:DimZMax) : | real(8)
| ||
pz_VarX(DimXMin:DimXMax, DimZMin:DimZMax) : | real(8), intent(in)
|
z 方向に半格子ずれた点での数値拡散項を評価
function pz_NumDiffVelX(pz_VarX) ! ! z 方向に半格子ずれた点での数値拡散項を評価 ! !暗黙の型宣言禁止 implicit none !変数定義 real(8), intent(in) :: pz_VarX(DimXMin:DimXMax, DimZMin:DimZMax) !物理量 real(8) :: pz_NumDiffVelX(DimXMin:DimXMax, DimZMin:DimZMax) !数値拡散 ! pz_NumDiffVelX = 0.0d0 !初期化 pz_NumDiffVelX = NuH * ( pz_dx_xz( xz_dx_pz( pz_VarX ) ) ) + NuV * ( pz_dz_pr( pr_dz_pz( pz_VarX ) ) ) end function pz_NumDiffVelX
Function : | |||
xr_NumDiffVelZ(DimXMin:DimXMax, DimZMin:DimZMax) : | real(8)
| ||
xr_VarZ(DimXMin:DimXMax, DimZMin:DimZMax) : | real(8), intent(in)
|
x 方向に半格子ずれた点での数値拡散項を評価
function xr_NumDiffVelZ(xr_VarZ) ! ! x 方向に半格子ずれた点での数値拡散項を評価 ! !暗黙の型宣言禁止 implicit none !変数定義 real(8), intent(in) :: xr_VarZ(DimXMin:DimXMax, DimZMin:DimZMax) !物理量 real(8) :: xr_NumDiffVelZ(DimXMin:DimXMax, DimZMin:DimZMax) !数値拡散 ! xr_NumDiffVelZ = 0.0d0 xr_NumDiffVelZ = NuH * ( xr_dx_pr( pr_dx_xr( xr_VarZ ) ) ) + NuV * ( xr_dz_xz( xz_dz_xr( xr_VarZ ) ) ) end function xr_NumDiffVelZ
Function : | |||
xz_NumDiffKm(DimXMin:DimXMax, DimZMin:DimZMax) : | real(8)
| ||
xz_Scalar(DimXMin:DimXMax, DimZMin:DimZMax) : | real(8), intent(in)
|
x, z 方向に半格子ずれた点での数値拡散項を評価
function xz_NumDiffKm(xz_Scalar) ! ! x, z 方向に半格子ずれた点での数値拡散項を評価 ! !暗黙の型宣言禁止 implicit none !変数定義 real(8), intent(in) :: xz_Scalar(DimXMin:DimXMax, DimZMin:DimZMax) !スカラー量 real(8) :: xz_NumDiffKm(DimXMin:DimXMax, DimZMin:DimZMax) !水平方向の数値拡散 xz_NumDiffKm = NuH * (xz_dx_pz(pz_dx_xz( xz_Scalar ))) + NuV * (xz_dz_xr(xr_dz_xz( xz_Scalar ))) end function xz_NumDiffKm
Function : | |||
xz_NumDiffScalar(DimXMin:DimXMax, DimZMin:DimZMax) : | real(8)
| ||
xz_Scalar(DimXMin:DimXMax, DimZMin:DimZMax) : | real(8), intent(in)
|
x, z 方向に半格子ずれた点での数値拡散項を評価
function xz_NumDiffScalar(xz_Scalar) ! ! x, z 方向に半格子ずれた点での数値拡散項を評価 ! !暗黙の型宣言禁止 implicit none !変数定義 real(8), intent(in) :: xz_Scalar(DimXMin:DimXMax, DimZMin:DimZMax) !スカラー量 real(8) :: xz_NumDiffScalar(DimXMin:DimXMax, DimZMin:DimZMax) !水平方向の数値拡散 ! xz_NumDiffScalar = 0.0d0 !初期化 xz_NumDiffScalar = NuH * (xz_dx_pz(pz_dx_xz( xz_Scalar ))) + NuV * (xz_dz_xr(xr_dz_xz( xz_Scalar ))) call StorePotTempDiff( xz_NumDiffScalar ) end function xz_NumDiffScalar
Function : | |||
xza_NumDiffScalar(DimXMin:DimXMax, DimZMin:DimZMax, SpcNum) : | real(8)
| ||
xza_Scalar(DimXMin:DimXMax, DimZMin:DimZMax, SpcNum) : | real(8), intent(in)
|
x, z 方向に半格子ずれた点での数値拡散項を評価
function xza_NumDiffScalar(xza_Scalar) ! ! x, z 方向に半格子ずれた点での数値拡散項を評価 ! !暗黙の型宣言禁止 implicit none !変数定義 real(8), intent(in) :: xza_Scalar(DimXMin:DimXMax, DimZMin:DimZMax, SpcNum) !スカラー量 real(8) :: xza_NumDiffScalar(DimXMin:DimXMax, DimZMin:DimZMax, SpcNum) !水平方向の数値拡散 integer :: s do s = 1, SpcNum xza_NumDiffScalar(:,:,s) = NuH * (xz_dx_pz(pz_dx_xz( xza_Scalar(:,:,s) ))) + NuV * (xz_dz_xr(xr_dz_xz( xza_Scalar(:,:,s) ))) end do call StoreMixRtDiff( xza_NumDiffScalar ) end function xza_NumDiffScalar