Class | NumDiffusion |
In: |
util/numdiffusion.f90
|
Copyright (C) GFD Dennou Club, 2005. All rights reserved.
* Developer: SUGIYAMA Ko-ichiro * Version: $Id: numdiffusion.f90,v 1.1.1.1 2006/04/25 03:43:58 deepconv Exp $ * Tag Name: $Name: $ * Change History:
数値拡散項の計算モジュール
subroutine NumDiffusion_init () ! ! NumDiffusion モジュールの初期化ルーチン ! !暗黙の型宣言禁止 implicit none ! 2 次精度中心差分の場合 ! CReSS マニュアルでは, Alpha < 1/8 くらいが適当と述べている. NuH = Alpha * ( DelX ** 2.0d0 ) / DelTimeLong NuV = Alpha * ( DelZ ** 2.0d0 ) / DelTimeLong ! 4 次精度中心差分の場合は以下を利用. ! CReSS マニュアルでは, Alpha = 0.001 くらいが適当と述べている. ! NuH = Alpha * ( DelX ** 4.0d0 ) / DelTimeLong ! NuV = Alpha * ( DelZ ** 4.0d0 ) / DelTimeLong !確認 write(*,*) "NumDiffusion_init, NuH: ", NuH write(*,*) "NumDiffusion_init, NuV: ", NuV end subroutine