# ---------------------------------------------------------------------------
# pyCGNS - Python package for CFD General Notation System -
# See license.txt file in the root directory of this Python module source
# ---------------------------------------------------------------------------
#
from __future__ import unicode_literals
import CGNS.PAT.cgnslib as C
import CGNS.PAT.cgnskeywords as K
import CGNS.PAT.cgnserrors as E
import numpy as N
#
#
data = C.newArbitraryGridMotion(None, '{ArbitraryGridMotion}')
C.newRind(data, N.array([0, 0, 0, 0, 1, 1]))
C.newGridLocation(data)
C.newDataArray(data, K.GridVelocityX_s)
C.newDataArray(data, K.GridVelocityY_s)
C.newDataArray(data, K.GridVelocityZ_s)
C.newDataClass(data)
C.newDimensionalUnits(data)
C.newUserDefinedData(data, '{UserDefinedData}')
C.newDescriptor(data, '{Descriptor}')
#
status = '11.3'
comment = 'Full SIDS with all optionals'
pattern = [data, status, comment]
#
m = CGL.newArbitraryGridMotion(z, '{ArbitraryGridMotion}')
d = CGL.newDataClass(m, CGK.NondimensionalParameter_s)
x = CGL.newZoneGridConnectivity(z)
x = CGL.newGridConnectivity(x, '{GridConnectivity}', z[0])
CGL.newPointRange(x,
value=NPY.array([[1, 1], [1, 1], [1, 1]],
dtype=NPY.int32,
order='Fortran'))
p = CGL.newGridConnectivityProperty(x)
m = CGL.newPeriodic(p)
d = CGL.newDataClass(m, CGK.NondimensionalParameter_s)
w = CGL.newConvergenceHistory(b)
d = CGL.newDataClass(w, CGK.NondimensionalParameter_s)
i = CGL.newIntegralData(b, '{IntegralData}')
d = CGL.newDataClass(i, CGK.NondimensionalParameter_s)
i = CGL.newUserDefinedData(b, '{UserDefinedData}')
d = CGL.newDataClass(i, CGK.NondimensionalParameter_s)
i = CGL.newGravity(b)
d = CGL.newDataClass(i, CGK.NondimensionalParameter_s)
f = CGL.newFlowEquationSet(b)
d = CGL.newDataClass(f, CGK.NondimensionalParameter_s)
m = CGL.newGasModel(f)
d = CGL.newDataClass(m, CGK.NondimensionalParameter_s)
m = CGL.newViscosityModel(f)
d = CGL.newDataClass(m, CGK.NondimensionalParameter_s)
m = CGL.newThermalConductivityModel(f)
d = CGL.newDataClass(m, CGK.NondimensionalParameter_s)
m = CGL.newThermalRelaxationModel(f)
d = CGL.newDataClass(m, CGK.NondimensionalParameter_s)
m = CGL.newChemicalKineticsModel(f)
d = CGL.newDataClass(m, CGK.NondimensionalParameter_s)
# pyCGNS - Python package for CFD General Notation System -
# See license.txt file in the root directory of this Python module source
# -------------------------------------------------------------------------
#
from __future__ import unicode_literals
import CGNS.PAT.cgnslib as CGL
import CGNS.PAT.cgnsutils as CGU
import CGNS.PAT.cgnskeywords as CGK
import numpy as NPY
TESTS = []
# -------------------------------------------------------------------------
tag = 'base level family #1'
diag = True
T = CGL.newCGNSTree()
b = CGL.newBase(T, '{Base#001}', 3, 3)
CGL.newFamily(b, '{Family#001}')
TESTS.append((tag, T, diag))
# -------------------------------------------------------------------------
tag = 'family elsewhere than base'
diag = False
T = CGL.newCGNSTree()
b = CGL.newBase(T, '{Base#001}', 3, 3)
u = CGL.newUserDefinedData(b, '{UD#001}')
CGL.newFamily(u, '{Family#001}')
TESTS.append((tag, T, diag))
# -------------------------------------------------------------------------
NPY.array([[5, 4, 0], [7, 6, 0], [9, 8, 0]], order='F'))
g = CGL.newGridCoordinates(z, 'GridCoordinates')
g = CGL.newGridCoordinates(z, '{Grid#002}')
g = CGL.newGridCoordinates(z, '{Grid#003}')
g = CGL.newGridCoordinates(z, '{Grid#004}')
TESTS.append((tag, T, diag))
# -------------------------------------------------------------------------
tag = 'grid empty #2'
diag = False
T = CGL.newCGNSTree()
b = CGL.newBase(T, '{Base}', 3, 3)
z = CGL.newZone(b, '{Zone}',
NPY.array([[5, 4, 0], [7, 6, 0], [9, 8, 0]], order='F'))
g = CGL.newGridCoordinates(z, 'GridCoordinates')
r = CGL.newUserDefinedData(g, '{UserDefinedData}')
TESTS.append((tag, T, diag))
# -------------------------------------------------------------------------
tag = 'grid bad dims #1'
diag = False
T = CGL.newCGNSTree()
b = CGL.newBase(T, '{Base}', 3, 3)
z = CGL.newZone(b, '{Zone}',
NPY.array([[5, 4, 0], [7, 6, 0], [9, 8, 0]], order='F'))
g = CGL.newGridCoordinates(z, 'GridCoordinates')
d = CGL.newDataArray(g, CGK.CoordinateX_s,
NPY.ones((4, 7, 9), dtype='float64', order='F'))
d = CGL.newDataArray(g, CGK.CoordinateY_s,
NPY.ones((5, 7, 9), dtype='float64', order='F'))
d = CGL.newDataArray(g, CGK.CoordinateZ_s,
