Class NumDiffusion
In: util/numdiffusion.f90

    Copyright (C) GFD Dennou Club, 2005. All rights reserved.

Module NumDiffusion

  * 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:

Overview

数値拡散項の計算モジュール

Error Handling

Bugs

Note

Future Plans

Methods

Included Modules

gridset timeset differentiate_center4

Public Instance methods

[Source]

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

[Validate]