Class | CFLCheck_3d |
In: |
util/cflcheck_3d.f90
|
CFL 条件のチェックをするためのパッケージ型モジュール
* 音波に対して CFL 条件をチェック * 入力された速度に対して CFL 条件をチェック
Subroutine : | |
pyz_VelX(DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax) : | real(DP), intent(in) |
水平速度に対して CFL 条件をチェック.
subroutine CFLCheckTimeLongVelX( pyz_VelX ) ! !水平速度に対して CFL 条件をチェック. ! !暗黙の型宣言禁止 implicit none !内部変数 real(DP), intent(in) :: pyz_VelX (DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax) real(DP) :: CFL !CFL 条件を求める CFL = 2.0d0*DelTimeLong/ max(abs(maxval(xyz_dx/pyz_VelX)), abs(minval(xyz_dx/pyz_VelX))) !メッセージ出力 call MessageNotify( "M", "CFLCheckTimeLongVelX", "Courant number of VelX for DelTimeLong = %f", d=(/CFL/) ) end subroutine CFLCheckTimeLongVelX
Subroutine : | |
xqz_VelY(DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax) : | real(DP), intent(in) |
水平速度に対して CFL 条件をチェック.
subroutine CFLCheckTimeLongVelY( xqz_VelY ) ! !水平速度に対して CFL 条件をチェック. ! !暗黙の型宣言禁止 implicit none !内部変数 real(DP), intent(in) :: xqz_VelY (DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax) real(DP) :: CFL !CFL 条件を求める CFL = 2.0d0*DelTimeLong/ max(abs(maxval(xyz_dy/xqz_VelY)), abs(minval(xyz_dy/xqz_VelY))) !メッセージ出力 call MessageNotify( "M", "CFLCheckTimeLongVelY", "Courant number of VelY for DelTimeLong = %f", d=(/CFL/) ) end subroutine CFLCheckTimeLongVelY
Subroutine : | |
xyr_VelZ(DimXMin:DimXMax,DimYMin:DImYMax,DimZMin:DimZMax) : | real(DP), intent(in) |
水平速度に対して CFL 条件をチェック.
subroutine CFLCheckTimeLongVelZ( xyr_VelZ ) ! !水平速度に対して CFL 条件をチェック. ! !暗黙の型宣言禁止 implicit none !内部変数 real(DP), intent(in) :: xyr_VelZ (DimXMin:DimXMax,DimYMin:DImYMax,DimZMin:DimZMax) real(DP) :: CFL !CFL 条件を求める CFL = 2.0d0*DelTimeLong/ max(abs(maxval(xyz_dz/xyr_VelZ)), abs(minval(xyz_dz/xyr_VelZ))) !メッセージ出力 call MessageNotify( "M", "CFLCheckTimeLongVelZ", "Courant number of VelZ for DelTimeLong = %f", d=(/CFL/) ) end subroutine CFLCheckTimeLongVelZ
Subroutine : | |||
xyz_VelSound(DimXMin:DimXMax,DimYMin:DImYMax, DimZMin:DimZMax) : | real(DP), intent(in)
|
音波に対して CFL 条件をチェック
subroutine CFLCheckTimeShort( xyz_VelSound ) ! !音波に対して CFL 条件をチェック ! !暗黙の型宣言禁止 implicit none !変数定義 real(DP), intent(in) :: xyz_VelSound (DimXMin:DimXMax,DimYMin:DImYMax, DimZMin:DimZMax) !音速 real(DP) :: CFL !クーラン数 real(DP) :: DxyzMin !最小格子間隔 !音速と CFL 条件を求める DxyzMin = min(min(minval(xyz_dx), minval(xyz_dy)), minval(xyz_dz)) CFL = DelTimeShort * maxval(xyz_VelSound) / DxyzMin !メッセージ call MessageNotify( "M", "CFLCheckTimeShort", "Sound Wave Velocity = %f", d=(/maxval(xyz_VelSound)/) ) call MessageNotify( "M", "CFLCheckTimeShort", "min(DelX, DelY, DelZ) = %f", d=(/DxyzMin/) ) call MessageNotify( "M", "CFLCheckTimeShort", "DelTimeShort = %f", d=(/DelTimeSHort/) ) !警告メッセージ if ( CFL >= 1.0) then call MessageNotify( "E", "CFLCheckTimeShort", "CFL Condition is broken, DelTimeShort * VelSound > min(DelX, DelZ)") else call MessageNotify( "M", "CFLCheckTimeShort", "Courant number for DelTimeSort = %f", d=(/CFL/) ) end if end subroutine CFLCheckTimeShort