cgnslib.f90

国外大名顶顶的“台风”并行计算流体力学CFD软件的早期版本的源代码

! Fortran 90 version of cgnslib.h : MODULE

module CGNSLIB

  integer, parameter :: cgnslen = 32
!* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - *\
!*      modes for cgns file                                          *
!* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
  integer, parameter :: MODE_READ   = 0
  integer, parameter :: MODE_WRITE  = 1
  integer, parameter :: MODE_CLOSED = 2
  integer, parameter :: MODE_MODIFY = 3

!* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - *\
!*      some error code                                               *
!* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
  integer, parameter :: ALL_OK         = 0
  integer, parameter :: ERROR          = 1
  integer, parameter :: NODE_NOT_FOUND = 2
  integer, parameter :: INCORRECT_PATH = 3

!* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - *\
!*      Dimensional Units                                                *
!* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
  integer, parameter :: Null = 0
  integer, parameter :: UserDefined = 1
  integer, parameter :: Kilogram  = 2
  integer, parameter :: Gram      = 3
  integer, parameter :: Slug      = 4
  integer, parameter :: PoundMass = 5

  character(len=cgnslen), parameter :: MassUnitsName(0:5) = &
    (/'Null                            ','UserDefined                     ',&
      'Kilogram                        ','Gram                            ',&
      'Slug                            ','PoundMass                       '/)

  integer, parameter :: Meter      = 2
  integer, parameter :: Centimeter = 3 
  integer, parameter :: Millimeter = 4
  integer, parameter :: Foot       = 5
  integer, parameter :: Inch       = 6
  character(len=cgnslen), parameter :: LengthUnitsName(0:6) = &
    (/'Null                            ','UserDefined                     ',&
      'Meter                           ','Centimeter                      ',&
      'Millimeter                      ','Foot                            ',&
      'Inch                            ' /)

  integer, parameter :: Second = 2
  character(len=cgnslen), parameter :: TimeUnitsName(0:2) = &
    (/'Null                            ','UserDefined                     ',&
      'Second                          ' /)

  integer, parameter :: Kelvin     = 2
  integer, parameter :: Celcius    = 3
  integer, parameter :: Rankine    = 4
  integer, parameter :: Fahrenheit = 5
  character(len=cgnslen), parameter :: TemperatureUnitsName(0:5) = &
    (/'Null                            ','UserDefined                     ',&
      'Kelvin                          ','Celcius                         ',&
      'Rankine                         ','Fahrenheit                      '/)

  integer, parameter :: Degree = 2
  integer, parameter :: Radian = 3
  character(len=cgnslen), parameter :: AngleUnitsName(0:3) = &
    (/'Null                            ','UserDefined                     ',&
      'Degree                          ','Radian                          '/)

!* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - *\
!*      Data Class                                                       *
!* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
  character(len=cgnslen), parameter :: DataClassName(0:6) = &
    (/'Null                            ','UserDefined                     ',&
      'Dimensional                     ','NormalizedByDimensional         ',&
      'NormalizedByUnknownDimensional  ','NondimensionalParameter         ',&
      'DimensionlessConstant           ' /)
  integer, parameter :: Dimensional                    = 2
  integer, parameter :: NormalizedByDimensional        = 3
  integer, parameter :: NormalizedByUnknownDimensional = 4
  integer, parameter :: NondimensionalParameter        = 5
  integer, parameter :: DimensionlessConstant          = 6


!* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - *\
!*      Grid Location                                                    *
!* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */

  character(len=cgnslen), parameter :: GridLocationName(0:8) = &
    (/'Null                            ','UserDefined                     ',&
      'Vertex                          ','CellCenter                      ',&
      'FaceCenter                      ','IFaceCenter                     ',&
      'JFaceCenter                     ','KFaceCenter                     ',&
      'EdgeCenter                      '/)
  integer, parameter :: Vertex      = 2
  integer, parameter :: CellCenter  = 3
  integer, parameter :: FaceCenter  = 4
  integer, parameter :: IFaceCenter = 5
  integer, parameter :: JFaceCenter = 6
  integer, parameter :: KFaceCenter = 7
  integer, parameter :: EdgeCenter  = 8

!* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - *\
!*      Grid Connectivity Types                                          *
!* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */

  character(len=cgnslen), parameter :: GridConnectivityTypeName(0:4) = &
    (/'Null                            ','UserDefined                     ',&
      'Overset                         ','Abutting                        ',&
      'Abutting1to1                    '/)
  integer, parameter :: Overset      = 2
  integer, parameter :: Abutting     = 3
  integer, parameter :: Abutting1to1 = 4

!* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - *\
!*      Point Set Types                                                  *
!* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */

  character(len=cgnslen), parameter :: PointSetTypeName(0:8) = &
    (/'Null                            ','UserDefined                     ',&
      'PointList                       ','PointListDonor                  ',&
      'PointRange                      ','PointRangeDonor                 ',&
      'ElementRange                    ','ElementList                     ',&
      'CellListDonor                   ' /)
  integer, parameter :: PointList       = 2
  integer, parameter :: PointListDonor  = 3
  integer, parameter :: PointRange      = 4
  integer, parameter :: PointRangeDonor = 5
  integer, parameter :: ElementRange    = 6
  integer, parameter :: ElementList     = 7
  integer, parameter :: CellListDonor   = 8

!* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - *\
!*      Governing Equations and Physical Models Types                    *
!* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */

  integer, parameter :: FullPotential             = 2
  integer, parameter :: Euler                     = 3
  integer, parameter :: NSLaminar                 = 4
  integer, parameter :: NSTurbulent               = 5
  integer, parameter :: NSLaminarIncompressible   = 6
  integer, parameter :: NSTurbulentIncompressible = 7
  character(len=cgnslen), parameter :: GoverningEquationsTypeName(0:7) = &
    (/'Null                            ','UserDefined                     ',&
      'FullPotential                   ','Euler                           ',&
      'NSLaminar                       ','NSTurbulent                     ',&
      'NSLaminarIncompressible         ','NSTurbulentIncompressible       '/)

! Gas Model
  integer, parameter :: Ideal                       = 2
  integer, parameter :: VanderWaals                 = 3
! Viscosity model
  integer, parameter :: Constant                    = 4
! Viscosity and thermal conductivity model
  integer, parameter :: PowerLaw                    = 5
  integer, parameter :: SutherlandLaw               = 6
! Thermal conductivity model only
  integer, parameter :: ConstantPrandtl             = 7
! Turbulence closure
  integer, parameter :: EddyViscosity               = 8
  integer, parameter :: ReynoldsStress              = 9
  integer, parameter :: ReynoldsStressAlgebraic     = 10
! Turbulence model
  integer, parameter :: Algebraic_BaldwinLomax      = 11
  integer, parameter :: Algebraic_CebeciSmith       = 12
  integer, parameter :: HalfEquation_JohnsonKing    = 13